Acessibilidade / Reportar erro

Variation of rainfall features across the Plain of Reeds under the impacts of climate change1 1 Research developed at Faculty of Social Science Education, Vietnam

Variação das características pluviométricas na Planície de Reeds sob os impactos das mudanças climáticas

ABSTRACT

Globally, there is enough evidence to affirm that climate change is one of the main factors causing the changing trends in rainfall features (CTRFs). The aim of this study, therefore, is to assess the CTRFs across the rice cultivation regions (RCRs) of the Plain of Reeds in Vietnam. Vietnam is a prime example of the impacts of climate change (ICC) on agricultural activities. Daily rainfall data series from 16 observation stations during the period 1984-2015 were appraised for quality and then the CTRFs were detected by applying the Mann_Kendall test and Sen’s slope estimator. The quality of rainfall data at all observation stations is homogeneous with a significance level of 0.05. Over a year, the trends in rainfall features slightly increased at most obsrvation stations, except for Hong Ngu and Moc Hoa. Additionally, the trends in rainfall features slightly decreased in rainny season rainfall (RSR) at 10 out of 16 observation stations, while insignificant to significant increasesing in dry season rainfall (DSR) were detected at all observation stations.

Key words:
agriculture; Mann-Kendall test; Pettitt; rainfall decline; rice

RESUMO

Globalmente, há evidências suficientes para afirmar que a mudança climática é um dos principais fatores que causam as tendências de mudança nas características das chuvas (TMCCs). O objetivo do estudo, portanto, é avaliar as TMCCs nas regiões de cultivo de arroz (RCAs) da Planície de Reeds. Vietnã é um excelente exemplo dos impactos das mudanças climáticas (IMC) nas atividades agrícolas. Para realização deste estudo, as séries diárias de dados de precipitação em 16 estações de observação durante o período de 1984-2015 foram avaliadas quanto à qualidade e, em seguida, foram detectadas as TMCCs aplicando o teste de Mann-Kendall e o estimador de inclinação de Sen. A qualidade dos dados pluviométricos em todas as estações de observação é homogênea com nível de significância de 0,05. Na escala anual, houve uma tendência crescente insignificante na maioria das estações de observação, exceto Hong Ngu e Moc Hoa, e uma tendência insignificante decrescente na precipitação da estação chuvosa foi detectada em 10 das 16 estações de observação, enquanto tendências crescentes insignificantes a significativas foram detectadas em todas estações de observação. Pode-se confirmar que a área de estudo está passando por um estágio crescente da precipitação na estação seca (PES). Os achados do estudo contribuem positivamente para as atividades agropecuárias, ampliando o entendimento do IMC sobre as características pluviométricas.

Palavras-chave:
agricultura; teste de Mann-Kendall; Pettitt; declínio de chuva; arroz

HIGHLIGHTS:

The Plain of Reeds befallen a slight change tends in annual rainfall and rainy season rainfall.

An increasing trend in dry season rainfall is a living proof of the impact of climate change on rainfall factors.

An increase in the dry season rainfall is expected as a positive highlight for the agricultural sector.

Introduction

In recent decades, climate change has strongly affected aspects of life, with its consequences thought to be a major theme of humanity in the 21st century (Bartels et al., 2020Bartels, R. J.; Black, A. W.; Keim, B. D. Trends in precipitation days in the United States. International Journal of Climatology, v.40, p.1038-1048, 2020. https://doi.org/10.1002/joc.6254
https://doi.org/10.1002/joc.6254...
; Dang et al., 2021Dang, T. A.; Nguyen, V. H.; Mai, P. N. Utilizing rainfed supply and irrigation as a climate variability adaptation solution for coastal lowland areas in Vietnam. Agriculture and natural resources, v.55, p.485-495, 2021. https://doi.org/10.34044/j.anres.2021.55.3.19
https://doi.org/10.34044/j.anres.2021.55...
). According to Alhathlou et al. (2021Alhathloul, S. H.; Khan, A. A.; Mishra, A. K. Trend analysis and change point detection of annual and seasonal horizontal visibility trends in Saudi Arabia. Theoretical and Applied Climatology, v.144, p.127-146, 2021. https://doi.org/10.1007/s00704-021-03533-z
https://doi.org/10.1007/s00704-021-03533...
), intensifying human activities have caused serious effects on the global water cycle, and extreme weather events (EWEs) are an inevitable consequence of climate change. The EWEs are, therefore, seen as a major factor contributing to the loss of crops (Lee & Dang, 2019Lee, S. K.; Dang, T. A. Spatio-temporal variations in meteorology drought over the Mekong River Delta of Vietnam in the recent decades. Paddy and Water Environment, v.17, p.35-44, 2019. https://doi.org/10.1007/s10333-018-0681-8
https://doi.org/10.1007/s10333-018-0681-...
; AlSubih et al., 2021AlSubih, M.; Kumar, M.; Mallick, J.; Ramakrishnan, R.; Islam, S.; Singh, C. K. Time series trend analysis of rainfall in last five decades and its quantification in Aseer Region of Saudi Arabia. Arabian Journal of Geosciences, v.14, p.1-15, 2021. https://doi.org/10.1007/s12517-021-06935-5
https://doi.org/10.1007/s12517-021-06935...
). Relevant studies have pointed out that the frequency of EWEs is increasing as a part of the impacts of climate change (ICC) (Dang et al., 2021Dang, T. A.; Nguyen, V. H.; Mai, P. N. Utilizing rainfed supply and irrigation as a climate variability adaptation solution for coastal lowland areas in Vietnam. Agriculture and natural resources, v.55, p.485-495, 2021. https://doi.org/10.34044/j.anres.2021.55.3.19
https://doi.org/10.34044/j.anres.2021.55...
; Dinh & Dang, 2022Dinh, T. K. H.; Dang, T. A. Potential risks of climate variability on rice cultivation regions in the Mekong Delta, Vietnam. Revista Brasileira de Engenharia Agrícola e Ambiental , v.26, p.348-355, 2022. https://doi.org/10.1590/1807-1929/agriambi.v26n5p348-355
https://doi.org/10.1590/1807-1929/agriam...
). Rainfall is considered a meteorological factor that closely relates to many aspects of life and where rainfed is considered as the main source of irrigation water provided (Nikumbh et al., 2019Nikumbh, A. C.; Chakraborty, A.; Bhat, G. S. Recent spatial aggregation tendency of rainfall extremes over India. Scientifc Reports, v.9, p.1-7, 2019. https://doi.org/10.1038/s41598-019-46719-2
https://doi.org/10.1038/s41598-019-46719...
; Balcha et al., 2022Balcha, Y. A.; Malcherek, A.; Alamirew, T. Understanding future climate in the upper Awash Basin (UASB) with selected climate model outputs under CMIP. Climate, v.10, p.1-28, 2022. https://doi.org/10.3390/cli10120185
https://doi.org/10.3390/cli10120185...
).

The changing trends in rainfall features (CTRFs) can influence many aspects of life, especially in the agricultural sector (Bartels et al., 2020Bartels, R. J.; Black, A. W.; Keim, B. D. Trends in precipitation days in the United States. International Journal of Climatology, v.40, p.1038-1048, 2020. https://doi.org/10.1002/joc.6254
https://doi.org/10.1002/joc.6254...
; Dang et al., 2021Dang, T. A.; Nguyen, V. H.; Mai, P. N. Utilizing rainfed supply and irrigation as a climate variability adaptation solution for coastal lowland areas in Vietnam. Agriculture and natural resources, v.55, p.485-495, 2021. https://doi.org/10.34044/j.anres.2021.55.3.19
https://doi.org/10.34044/j.anres.2021.55...
). Accordingly, an increase in the rainy season rainfall (RSR) may introduce potential risks, such as seed rot or rice falling near harvest time, leading to a grain yield decline or even a yield loss. A decreasing trend in the dry season rainfall (DSR) can be strongly linked to meteorological drought and saltwater intrusion (Balcha et al., 2022Balcha, Y. A.; Malcherek, A.; Alamirew, T. Understanding future climate in the upper Awash Basin (UASB) with selected climate model outputs under CMIP. Climate, v.10, p.1-28, 2022. https://doi.org/10.3390/cli10120185
https://doi.org/10.3390/cli10120185...
; Lee & Dang, 2019Lee, S. K.; Dang, T. A. Spatio-temporal variations in meteorology drought over the Mekong River Delta of Vietnam in the recent decades. Paddy and Water Environment, v.17, p.35-44, 2019. https://doi.org/10.1007/s10333-018-0681-8
https://doi.org/10.1007/s10333-018-0681-...
).

Globally, numerous studies related to the CTRFs under the ICC have attracted the attention of the scientific community (Lee & Dang, 2020Lee, S. K.; Dang, T. A. Extreme rainfall trends over the Mekong Delta under the impacts of climate change. International Journal of Climate Change Strategies and Management, v.12, p.639-652, 2020. https://doi.org/10.1108/IJCCSM-04-2020-0032
https://doi.org/10.1108/IJCCSM-04-2020-0...
; Dinh & Dang, 2022Dinh, T. K. H.; Dang, T. A. Potential risks of climate variability on rice cultivation regions in the Mekong Delta, Vietnam. Revista Brasileira de Engenharia Agrícola e Ambiental , v.26, p.348-355, 2022. https://doi.org/10.1590/1807-1929/agriambi.v26n5p348-355
https://doi.org/10.1590/1807-1929/agriam...
). Nowadays, studies on the CTRFs have been commonly deployed at large cultivation regions, such as the conterminous United States by Gershunov et al. (2019Gershunov, A.; Shulgina, T.; Clemesha, R. E. S.; Guirguis, K.; Pierce, D. W.; Dettinger, M. D.; Lavers, D. A.; Cayan, D. R.; Polade, S. D.; Kalansky, J.; Ralph, F. M. Precipitation regime change in Western North America: the role of atmospheric rivers. Scientific Reports, v.9, p.1-11, 2019.) and Risser et al. (2023Risser, M. D.; Collins, W. D.; Wehner, M. F.; O’Brien, T. A.; Paciorek, C. J.; O’Brien, J. P.; Patricola, C. M.; Huang, H.; Ullrich, P. A.; Loring, B. A framework for detection and attribution of regional precipitation change: Application to the United States historical record. Climate Dynamics, v.60, p.705-741, 2023. https://doi.org/10.1007/s00382-022-06321-1
https://doi.org/10.1007/s00382-022-06321...
); in India by Gupta et al. (2021Gupta, N.; Banerjee, A.; Gupta, S. K. Spatio-temporal trend analysis of climatic variables over Jharkhand, India. Earth Systems and Environment, v.5, p.71-86, 2021. https://doi.org/10.1007/s41748-021-00204-x
https://doi.org/10.1007/s41748-021-00204...
) and Praveen et al. (2020Praveen, B.; Talukdar, S.; Shahfahad; Susanta, M.; Jayanta, M.; Pritee, S.; Abu, R. M.; Towfiqul, I.; Atiqur, R. M. Analyzing trend and forecasting of rainfall changes in India using non-parametrical and machine learning approaches. Scientifc Reports, v.10, p.1-21, 2020. https://doi.org/10.1038/s41598-020-67228-7
https://doi.org/10.1038/s41598-020-67228...
); in Pakistan by Iqbal et al. (2019Iqbal, Z.; Shahid, S.; Ahmed, K.; Ismail, T.; Nawaz, N. Spatial distribution of the trends in precipitation and precipitation extremes in the sub-Himalayan region of Pakistan. Theoretical and Applied Climatology , v.137, p.2755-2769, 2019. https://doi.org/10.1007/s00704-019-02773-4
https://doi.org/10.1007/s00704-019-02773...
); in Iran by Ahmadi et al. (2018Ahmadi, F.; Nazeri, T. M.; Mirabbasi, R.; Khalili, K.; Jhajharia, D. Spatiotemporal trend and abrupt change analysis of temperature in Iran. Meteorological Applications, v.25, p.314-321, 2018. https://doi.org/10.1002/met.1694
https://doi.org/10.1002/met.1694...
), and in Saudi Arabia by Alhathloul et al. (2021Alhathloul, S. H.; Khan, A. A.; Mishra, A. K. Trend analysis and change point detection of annual and seasonal horizontal visibility trends in Saudi Arabia. Theoretical and Applied Climatology, v.144, p.127-146, 2021. https://doi.org/10.1007/s00704-021-03533-z
https://doi.org/10.1007/s00704-021-03533...
).

Thus, a study was conducted to assess the CTRFs across the rice cultivation regions (RCRs) of the Plain of Reeds, in Vietnam, applying the Mann_Kendall test and Sen’s slope estimator based on collected rainfall data series from 16 observation stations, during the period 1984-2015.

Material and Methods

The Plain of Reeds is located in the northwest region of the Mekong Delta (10o 04’ 56’’ to 11o 00’ N latitude and 105° 03’ 15” to 106° 09’ 33” E longitude), stretching a part of the border of three provinces of Long An, Tien Giang, and Dong Thap with a total agricultural land area of approximately of 700,000 hectares and an average annual yield of rice up to 560,000 tons (Lee & Dang, 2020Lee, S. K.; Dang, T. A. Extreme rainfall trends over the Mekong Delta under the impacts of climate change. International Journal of Climate Change Strategies and Management, v.12, p.639-652, 2020. https://doi.org/10.1108/IJCCSM-04-2020-0032
https://doi.org/10.1108/IJCCSM-04-2020-0...
; Dang, 2021). The area has an altitude varying from 0.5 to 2.0 m above mean sea level (Figure 1). However, the area is currently facing a lack of irrigation water due to the CTRFs, as a part of the ICC (Lee & Dang, 2019; Dang et al., 2021).

Figure 1
Map of the Plain of Reeds with rainfall observation stations marked red circles

Unlike the northern plains of Vietnam, which are controlled by four distinct seasons of the year according to the climatology features, the Plain of Reeds is dominated by the two major circulations in the Northeast and Southwest monsoons, creating the dry season from October to April and the rainy season from May to November (Dang, 2021Dang, T. A. Grain yield optimization in the Plain of Reeds in the context of climate variability. Revista Brasileira de Engenharia Agrícola e Ambiental, v.25, p.591-596, 2021. https://doi.org/10.1590/1807-1929/agriambi.v25n9p591-596
https://doi.org/10.1590/1807-1929/agriam...
). Annually, the area receives an average rainfall of approximately 1,465 mm, out of which 85% of rainfall comes from rainy season months (Figure 2).

Figure 2
Illustration of monthly rainfall across the study area in the period 1984 - 2015

To conduct this research, daily rainfall data series from 16 observation stations across the Plain of Reeds (Figure 1) were collected from the National Centre for Hydrometeorological Forecasting (NCHMF) between the dates 1984-2015. Accordingly, the quality of the rainfall data series was appraised by applying the homogeneity tests, namely the Standard Normal Homogeneity Test (SNHT) and Pettitt tests, based on XLSTAT software. This was then synthesized to total rainfall at seasonal and annual scales for further trend analysis.

To avoid interruption of data series, the present study applied the SNHT to expertise observed rainfall series across the study area (Patakamuri et al., 2020Patakamuri, S. K.; Muthiah, K.; Sridhar, V. Long-term homogeneity, trend, and change-point analysis of rainfall in the arid district of Ananthapuramu, Andhra Pradesh State, India. Water, v.12, p.1-22, 2020. https://doi.org/10.3390/w12010211
https://doi.org/10.3390/w12010211...
; Praveen et al., 2020Praveen, B.; Talukdar, S.; Shahfahad; Susanta, M.; Jayanta, M.; Pritee, S.; Abu, R. M.; Towfiqul, I.; Atiqur, R. M. Analyzing trend and forecasting of rainfall changes in India using non-parametrical and machine learning approaches. Scientifc Reports, v.10, p.1-21, 2020. https://doi.org/10.1038/s41598-020-67228-7
https://doi.org/10.1038/s41598-020-67228...
). This approach has been applied extensively to detect the discontinuity points in observed data series (Salehi et al., 2020Salehi, S.; Dehghani, M.; Mortazavi, S. M.; Singh, V. P. Trend analysis and change point detection of seasonal and annual precipitation in Iran. International Journal of Climatology , v.40, p.308-323, 2020. https://doi.org/10.1002/joc.6211
https://doi.org/10.1002/joc.6211...
; Lee & Dang, 2020Lee, S. K.; Dang, T. A. Extreme rainfall trends over the Mekong Delta under the impacts of climate change. International Journal of Climate Change Strategies and Management, v.12, p.639-652, 2020. https://doi.org/10.1108/IJCCSM-04-2020-0032
https://doi.org/10.1108/IJCCSM-04-2020-0...
). The SNHT is defined using Eq. 1:

Z i = Q i - Q ¯ σ Q (1)

where:

Qi - is the value of rainfall data series at time step i;

Q - is the mean value of rainfall data series; and,

σQ - is the standard deviation of rainfall data series at time step i.

Assuming the Zi is a normal distribution, N, the null hypothesis for all variants of the SNHT is Ho with Zi ? N(0, 1) i ? {1, …, n} i.e., the whole series is homogeneous.

When all values in the normalised series are normally distributed with a mean value equal to 0 and standard deviation equal to 1.

Pettitt’s test was developed to detect the interruption points in continuous data series (Lee & Dang, 2020Lee, S. K.; Dang, T. A. Extreme rainfall trends over the Mekong Delta under the impacts of climate change. International Journal of Climate Change Strategies and Management, v.12, p.639-652, 2020. https://doi.org/10.1108/IJCCSM-04-2020-0032
https://doi.org/10.1108/IJCCSM-04-2020-0...
; Dinh & Dang, 2022Dinh, T. K. H.; Dang, T. A. Potential risks of climate variability on rice cultivation regions in the Mekong Delta, Vietnam. Revista Brasileira de Engenharia Agrícola e Ambiental , v.26, p.348-355, 2022. https://doi.org/10.1590/1807-1929/agriambi.v26n5p348-355
https://doi.org/10.1590/1807-1929/agriam...
). Accordingly, the Pettitt test requires the data to be input independently, identically distributed random quantities, with the alternative being a stepwise shift in the mean, which is presented (Lee & Dang, 2019Lee, S. K.; Dang, T. A. Spatio-temporal variations in meteorology drought over the Mekong River Delta of Vietnam in the recent decades. Paddy and Water Environment, v.17, p.35-44, 2019. https://doi.org/10.1007/s10333-018-0681-8
https://doi.org/10.1007/s10333-018-0681-...
; Dinh & Dang, 2022Dinh, T. K. H.; Dang, T. A. Potential risks of climate variability on rice cultivation regions in the Mekong Delta, Vietnam. Revista Brasileira de Engenharia Agrícola e Ambiental , v.26, p.348-355, 2022. https://doi.org/10.1590/1807-1929/agriambi.v26n5p348-355
https://doi.org/10.1590/1807-1929/agriam...
). However, it cannot detect the change points in the distribution series if there is no change in the position of the applied points. Pettitt’s test is presented as follows:

By assuming the data series (X) are ranked from 1 to N, when the value of Vi,N is defined by Eq. 2:

V i = N + 1 - 2 R i   w i t h   I = 1 , 2 , 3 , . . . , N (2)

where Ri is the rank of Xi in the sample of N obtained data series.

And Ui is calculated by Eq. 3:

U i = U i - 1 + V i (3)

By assuming U1 = V1 and KN is calculated based on Eq. 4:

K N = m a x 1 i N U i (4)

Finally, the Pettitt test is established by Eq. 5:

P O A = 2 e - 6 K N 2 N 3 + N 2 (5)

If the value of POA is less than a critical value, the test is called the statistical significance when the H0 is rejected (Dang et al., 2021Dang, T. A.; Nguyen, V. H.; Mai, P. N. Utilizing rainfed supply and irrigation as a climate variability adaptation solution for coastal lowland areas in Vietnam. Agriculture and natural resources, v.55, p.485-495, 2021. https://doi.org/10.34044/j.anres.2021.55.3.19
https://doi.org/10.34044/j.anres.2021.55...
; Dinh & Dang, 2022Dinh, T. K. H.; Dang, T. A. Potential risks of climate variability on rice cultivation regions in the Mekong Delta, Vietnam. Revista Brasileira de Engenharia Agrícola e Ambiental , v.26, p.348-355, 2022. https://doi.org/10.1590/1807-1929/agriambi.v26n5p348-355
https://doi.org/10.1590/1807-1929/agriam...
).

To conduct this research, a non-parametric statistical method, namely the Mann-Kendall is applied to detect the CTRFs across the study area during the period 1984-2015. Accordingly, the Mann-Kendall test assumes the H0 of no trend is recorded while a monotonic trend (e.g., increasing or decreasing) is defined as the alternative hypothesis (H1) (Lee & Dang, 2019Lee, S. K.; Dang, T. A. Spatio-temporal variations in meteorology drought over the Mekong River Delta of Vietnam in the recent decades. Paddy and Water Environment, v.17, p.35-44, 2019. https://doi.org/10.1007/s10333-018-0681-8
https://doi.org/10.1007/s10333-018-0681-...
; Nguyen et al., 2022Nguyen, T. T. T.; Hoang, P. H. Y.; Dang, T. A. Climate variability induced drought across the coastal fringes of the Mekong Delta, Vietnam. MAUSAM, v.73, p.525-536, 2022. https://doi.org/10.54302/mausam.v73i3.5373
https://doi.org/10.54302/mausam.v73i3.53...
). The Mann-Kendall test is defined based on Eq. 6:

S = i = 1 n j = i + 1 n s g n X j - X i (6)

Where sgn(Xj - Xi) in Eq. 6 is defined by Eq. 7:

s g n X j - X i = + 1   i f   x j - x i > 0 0   i f   x j - x i = 0 - 1   i f   x j - x i < 0 (7)

In a time series, Xi, i = 1, 2, 3, ….., n, the value of S is assumed to be a normal distribution with a mean 0 while the discrepancy of statistics S has calculated applying Eq. 8:

v a r S = 1 18 n n - 1 2 n + 5 - j = 1 m t j t j - 1 2 t j + 5 (8)

The standard test (Zs) is used to define whether the time series information is a significant trend or not. When the Zs value is defined by applying Eq. 9:

Z s = S - 1 V a r S   i f   S > 0 0   i f   S = 0 S - 1 V a r S   i f   S < 0 (9)

Accordingly, the positive and negative values of Zs in a normalized test statistic reflect the increasing and decreasing trend, respectively, while the Zs having 0 values reflects a normally distributed data series (Banerjee et al., 2020Banerjee, A.; Chen, R. E.; Meadows, M.; Singh, R. B.; Mal, S.; Sengupta, D. An analysis of long-term rainfall trends and variability in the Uttarakhand Himalaya using google earth engine. Remote Sensing, v.12, p.1-24, 2020. https://doi.org/10.3390/rs12040709
https://doi.org/10.3390/rs12040709...
).

Sen’s slope estimator is also known as a non-parametric statistical test which is commonly used to define the true slope of an existing trend. In order hand, it was often used to define the magnitude of change trends in the observed data series (Dinh & Dang, 2022Dinh, T. K. H.; Dang, T. A. Potential risks of climate variability on rice cultivation regions in the Mekong Delta, Vietnam. Revista Brasileira de Engenharia Agrícola e Ambiental , v.26, p.348-355, 2022. https://doi.org/10.1590/1807-1929/agriambi.v26n5p348-355
https://doi.org/10.1590/1807-1929/agriam...
; Nguyen et al., 2022Nguyen, T. T. T.; Hoang, P. H. Y.; Dang, T. A. Climate variability induced drought across the coastal fringes of the Mekong Delta, Vietnam. MAUSAM, v.73, p.525-536, 2022. https://doi.org/10.54302/mausam.v73i3.5373
https://doi.org/10.54302/mausam.v73i3.53...
). The Sen’s slope is obtained using Eq. 10:

β = M e d i a n X i - X j i - j   w i t h   j < i (10)

where Xi and Xj are input data at time intervals ti and tj, respectively.

To constitute the spatial distribution map of rainfall, the Kriging method was adapted to interpolate the analyzed rainfall data points. This method is based on spatial autocorrelation data (Praveen et al., 2020Praveen, B.; Talukdar, S.; Shahfahad; Susanta, M.; Jayanta, M.; Pritee, S.; Abu, R. M.; Towfiqul, I.; Atiqur, R. M. Analyzing trend and forecasting of rainfall changes in India using non-parametrical and machine learning approaches. Scientifc Reports, v.10, p.1-21, 2020. https://doi.org/10.1038/s41598-020-67228-7
https://doi.org/10.1038/s41598-020-67228...
). It uses the spatial correlation between marked data points to interpolate the corresponding values in the spatial distribution field (Praveen et al., 2020Praveen, B.; Talukdar, S.; Shahfahad; Susanta, M.; Jayanta, M.; Pritee, S.; Abu, R. M.; Towfiqul, I.; Atiqur, R. M. Analyzing trend and forecasting of rainfall changes in India using non-parametrical and machine learning approaches. Scientifc Reports, v.10, p.1-21, 2020. https://doi.org/10.1038/s41598-020-67228-7
https://doi.org/10.1038/s41598-020-67228...
). In addition, the Kriging method allows for generating the estimation of the uncertainty surrounding each interpolated value point (Salehi et al., 2020Salehi, S.; Dehghani, M.; Mortazavi, S. M.; Singh, V. P. Trend analysis and change point detection of seasonal and annual precipitation in Iran. International Journal of Climatology , v.40, p.308-323, 2020. https://doi.org/10.1002/joc.6211
https://doi.org/10.1002/joc.6211...
). The Kriging method is commonly applied to conjecture the values Z*(x0) at the point x0 presented in Eq. 11.

Z * x 0 = i = 1 n λ i z x i (11)

Where λi interests to weights and the conjecture of error variance σk 2(x0) is defined as Eq. 12.

σ k 2 x 0 = μ + i = 1 n λ i γ x 0 - x i (12)

In Eq. 12, μ is the Lagrange constant, (x0-xi) presents the semi-variogram value equivalent to the distance between x0 and xi. While γ consists of the regionalized variable theory and intrinsic hypotheses and it is presented by Eq. 13:

γ h = 1 2 N h i = 1 N h Z x i - Z x i + h 2 (13)

In Eq. 13, γ(h) presents semi-variance, h presents lag distance, Z presents the rainfall-related factors.

Accordingly, the Kriging method was applied to draw the maps of rainfall spatial distribution across the Plain of Reeds based on ArcGIS (Version 10.8) software.

Results and Discussion

The daily rainfall data series from 16 observation stations in the Plain of Reeds during the period 1984-2015 was transformed into annual rainfall before conducting the breakpoint detection and applying the SNHT and Pettitt tests. The results indicated that the obtained critical values (p) of 16 rainfall observation stations are larger than the significance level (α = 0.05), which implies that the rainfall data series at all observation stations are of quality and meet the requirements for further studies (Table 1).

Table 1
Results of breakpoint detection for rainfall data series at all observation stations across the study area during the period 1984-2015

The basic statistical features of annual rainfall (AR) across the study area, including minimum, maximum, standard deviation (SD), and coefficient of variation (CV), are presented in Table 2. Statistical analysis indicates that the annual mean rainfall (AMR) across the study area is approximately 1,513.9 mm, while the maximum AMR is recorded in the northeastern part, Tan Thanh, Moc Hoa, Ben Luc and Thu Thua (1,658.7, 1,637.1, 1,620.1 and 1,567.8 mm, respectively). The minimum AMR appears in the western and southwestern parts, Tan Hung (1,231.6 mm), Hong Ngu (1,287.6 mm), Cay Lay (1,389.7 mm) and My Tho (1,396.6 mm) (Table 2).

Table 2
Features of annual rainfall at the observation stations in the period 1984-2015

The SD of annual rainfall varies from 217.4 to 330.6 mm with the highest SD detected in the northeastern part i.e., Moc Hoa (330.6 mm), followed by Tan Thanh (329.7 mm) and Thu Thua (308.6 mm). While the lowest SD was concentrated in the southern part, i.e., Chau Thanh (217.4 mm), followed by Cai Lay (241.5 mm) and Cai Be (234.8 mm) (Table 2). It means that the study area has a strong fluctuation in AR.

The spatial distribution maps of the CV of the AR and seasonal rainfall across the study area were established by applying the Kriging interpolation technique (Figure 3). The results pointed out that the southern part (Cai Be and Cay Lay stations) recorded the highest fluctuations of AR (e.g., 57.9 and 59.8%) while the southwest part (Cao Lanh station) recorded the lowest fluctuations of AR (51.9%) (Figure 3A). For RSR, the highest fluctuations were detected in the northeastern part, measuring 56.7% (Moc Hoa station), while the lowest fluctuations were recorded in the southern part (Cai Lay station), estimated at approximately 37.3% (Figure 3B). For DSR, the whole study area recorded a strong fluctuation in the CV, varying from 87.3 to 94.2% (Figure 3C). In general, a strong fluctuation in seasonal rainfall can lead to negative impacts on local agriculture activities, especially irrigation. For example, in rice cultivation paddies, rainfed agriculture plays an important role in irrigation in the dry season. The results of the Mann_Kendall test and Sen’s slope estimator for rainfall features across the study area during the period 1984-2015 are presented in Table 3. A value of α=0.05 was applied to the analysis for Zs-statistics of the Mann_Kendall test and β- statistics of Sen’s slope estimator.

Figure 3
The distribution of the coefficient of variation (CV) of annual rainfall-AR (A), rainy season rainfall-RSR (B), and dry season rainfall-DSR (C) across the Plain of Reeds

Table 3
Results of Mann_Kendall (Zs) test and Sen’ estimator slope (β) for rainfall features across the study area during the period 1984-2015

For AR, the Zs showed positive values varying from 0.37 to 1.83, which were detected in 14 observation stations (Figure 4A). It implies an uptrend of AR over most of the study area, while only 2 out of 16 observation stations (e.g., Hong Ngu and Moc Hoa) were detected with Zs values from -0.6 to -0.04 (Table 3). For the RSR, the Zs values varied from 0.21 to 1.37 at 6 out of the 16 observation stations. This indicates that a slight increasing trend in RSR was recorded in less than 40% of the total area of the study (Figure 4B). A decreasing trend in the Zs value, varying from -1.87 to -0.17 was recorded at 10 out of 16 observation stations for DSR (Figure 4C). A study on the spatiotemporal variations in drought events over the Mekong Delta of Vietnam by Lee & Dang (2019Lee, S. K.; Dang, T. A. Spatio-temporal variations in meteorology drought over the Mekong River Delta of Vietnam in the recent decades. Paddy and Water Environment, v.17, p.35-44, 2019. https://doi.org/10.1007/s10333-018-0681-8
https://doi.org/10.1007/s10333-018-0681-...
) stated that a decreasing trend in the RSR was detected in the eastern and northeastern parts of the Mekong Delta. A similar result of Nguyen et al. (2022Nguyen, T. T. T.; Hoang, P. H. Y.; Dang, T. A. Climate variability induced drought across the coastal fringes of the Mekong Delta, Vietnam. MAUSAM, v.73, p.525-536, 2022. https://doi.org/10.54302/mausam.v73i3.5373
https://doi.org/10.54302/mausam.v73i3.53...
) reported that a decreasing trend in the RSR is contributed to an increasing trend in the six-month drought timescales.

Figure 4
Spatial variation of rainfall features analyzed using the Mann_Kendall (Zs) test for annual rainfall - AR (A), rainy season rainfall - RSR (B), and dry season rainfall - DSR (C) across the study area during the period 1984-2015

In general, no observation stations for AR and RSR recorded significantly increasing or/decreasing trends. However, an increasing trend in DSR was detected at all observation stations, with a slight increase in trend (Zs value varying from 0.21 to 1.87) at 10 out of 16 observation stations, while a significant trend (Zs = 1.95÷2.97) was observed at 6 out of 16 observation stations (Table 3).

A study on the spatiotemporal variation trends in rainfall and drought over the Mekong Delta of Vietnam by Lee & Dang (2019Lee, S. K.; Dang, T. A. Spatio-temporal variations in meteorology drought over the Mekong River Delta of Vietnam in the recent decades. Paddy and Water Environment, v.17, p.35-44, 2019. https://doi.org/10.1007/s10333-018-0681-8
https://doi.org/10.1007/s10333-018-0681-...
) reported that an increasing trend in the DSR has occurred during the period 1984-2019. Accordingly, an increasing trend in DSR will bring a positive contribution to agricultural activities such as supplementing the scarce irrigation water from irrigation canals. A decreasing trend in RSR will contribute to reducing the damage caused by heavy rainfall events to sowing activities, as well as the crop harvest.

Conclusions

  1. No insignificant trends in increasing or decreasing annual and rainy seasonal rainfall were observed. However, an increasing trend was detected in dry season rainfall.

  2. In the eastern part of the Plain of Reeds (including Moc Hoa, Thu Thua, Tan Thanh, Cai Be, Cai Lay and Tan Phuoc stations), there apears to be a significant increasing trend in dry season rainfall.

  3. The findings indicated that the Plain of Reeds is experiencing a period of increased dry season rainfall, and it may contribute positively to agricultural production activities.

Literature Cited

  • Ahmadi, F.; Nazeri, T. M.; Mirabbasi, R.; Khalili, K.; Jhajharia, D. Spatiotemporal trend and abrupt change analysis of temperature in Iran. Meteorological Applications, v.25, p.314-321, 2018. https://doi.org/10.1002/met.1694
    » https://doi.org/10.1002/met.1694
  • Alhathloul, S. H.; Khan, A. A.; Mishra, A. K. Trend analysis and change point detection of annual and seasonal horizontal visibility trends in Saudi Arabia. Theoretical and Applied Climatology, v.144, p.127-146, 2021. https://doi.org/10.1007/s00704-021-03533-z
    » https://doi.org/10.1007/s00704-021-03533-z
  • AlSubih, M.; Kumar, M.; Mallick, J.; Ramakrishnan, R.; Islam, S.; Singh, C. K. Time series trend analysis of rainfall in last five decades and its quantification in Aseer Region of Saudi Arabia. Arabian Journal of Geosciences, v.14, p.1-15, 2021. https://doi.org/10.1007/s12517-021-06935-5
    » https://doi.org/10.1007/s12517-021-06935-5
  • Balcha, Y. A.; Malcherek, A.; Alamirew, T. Understanding future climate in the upper Awash Basin (UASB) with selected climate model outputs under CMIP. Climate, v.10, p.1-28, 2022. https://doi.org/10.3390/cli10120185
    » https://doi.org/10.3390/cli10120185
  • Banerjee, A.; Chen, R. E.; Meadows, M.; Singh, R. B.; Mal, S.; Sengupta, D. An analysis of long-term rainfall trends and variability in the Uttarakhand Himalaya using google earth engine. Remote Sensing, v.12, p.1-24, 2020. https://doi.org/10.3390/rs12040709
    » https://doi.org/10.3390/rs12040709
  • Bartels, R. J.; Black, A. W.; Keim, B. D. Trends in precipitation days in the United States. International Journal of Climatology, v.40, p.1038-1048, 2020. https://doi.org/10.1002/joc.6254
    » https://doi.org/10.1002/joc.6254
  • Dang, T. A. Grain yield optimization in the Plain of Reeds in the context of climate variability. Revista Brasileira de Engenharia Agrícola e Ambiental, v.25, p.591-596, 2021. https://doi.org/10.1590/1807-1929/agriambi.v25n9p591-596
    » https://doi.org/10.1590/1807-1929/agriambi.v25n9p591-596
  • Dang, T. A.; Nguyen, V. H.; Mai, P. N. Utilizing rainfed supply and irrigation as a climate variability adaptation solution for coastal lowland areas in Vietnam. Agriculture and natural resources, v.55, p.485-495, 2021. https://doi.org/10.34044/j.anres.2021.55.3.19
    » https://doi.org/10.34044/j.anres.2021.55.3.19
  • Dinh, T. K. H.; Dang, T. A. Potential risks of climate variability on rice cultivation regions in the Mekong Delta, Vietnam. Revista Brasileira de Engenharia Agrícola e Ambiental , v.26, p.348-355, 2022. https://doi.org/10.1590/1807-1929/agriambi.v26n5p348-355
    » https://doi.org/10.1590/1807-1929/agriambi.v26n5p348-355
  • Gershunov, A.; Shulgina, T.; Clemesha, R. E. S.; Guirguis, K.; Pierce, D. W.; Dettinger, M. D.; Lavers, D. A.; Cayan, D. R.; Polade, S. D.; Kalansky, J.; Ralph, F. M. Precipitation regime change in Western North America: the role of atmospheric rivers. Scientific Reports, v.9, p.1-11, 2019.
  • Gupta, N.; Banerjee, A.; Gupta, S. K. Spatio-temporal trend analysis of climatic variables over Jharkhand, India. Earth Systems and Environment, v.5, p.71-86, 2021. https://doi.org/10.1007/s41748-021-00204-x
    » https://doi.org/10.1007/s41748-021-00204-x
  • Iqbal, Z.; Shahid, S.; Ahmed, K.; Ismail, T.; Nawaz, N. Spatial distribution of the trends in precipitation and precipitation extremes in the sub-Himalayan region of Pakistan. Theoretical and Applied Climatology , v.137, p.2755-2769, 2019. https://doi.org/10.1007/s00704-019-02773-4
    » https://doi.org/10.1007/s00704-019-02773-4
  • Lee, S. K.; Dang, T. A. Extreme rainfall trends over the Mekong Delta under the impacts of climate change. International Journal of Climate Change Strategies and Management, v.12, p.639-652, 2020. https://doi.org/10.1108/IJCCSM-04-2020-0032
    » https://doi.org/10.1108/IJCCSM-04-2020-0032
  • Lee, S. K.; Dang, T. A. Spatio-temporal variations in meteorology drought over the Mekong River Delta of Vietnam in the recent decades. Paddy and Water Environment, v.17, p.35-44, 2019. https://doi.org/10.1007/s10333-018-0681-8
    » https://doi.org/10.1007/s10333-018-0681-8
  • Nguyen, T. T. T.; Hoang, P. H. Y.; Dang, T. A. Climate variability induced drought across the coastal fringes of the Mekong Delta, Vietnam. MAUSAM, v.73, p.525-536, 2022. https://doi.org/10.54302/mausam.v73i3.5373
    » https://doi.org/10.54302/mausam.v73i3.5373
  • Nikumbh, A. C.; Chakraborty, A.; Bhat, G. S. Recent spatial aggregation tendency of rainfall extremes over India. Scientifc Reports, v.9, p.1-7, 2019. https://doi.org/10.1038/s41598-019-46719-2
    » https://doi.org/10.1038/s41598-019-46719-2
  • Patakamuri, S. K.; Muthiah, K.; Sridhar, V. Long-term homogeneity, trend, and change-point analysis of rainfall in the arid district of Ananthapuramu, Andhra Pradesh State, India. Water, v.12, p.1-22, 2020. https://doi.org/10.3390/w12010211
    » https://doi.org/10.3390/w12010211
  • Praveen, B.; Talukdar, S.; Shahfahad; Susanta, M.; Jayanta, M.; Pritee, S.; Abu, R. M.; Towfiqul, I.; Atiqur, R. M. Analyzing trend and forecasting of rainfall changes in India using non-parametrical and machine learning approaches. Scientifc Reports, v.10, p.1-21, 2020. https://doi.org/10.1038/s41598-020-67228-7
    » https://doi.org/10.1038/s41598-020-67228-7
  • Risser, M. D.; Collins, W. D.; Wehner, M. F.; O’Brien, T. A.; Paciorek, C. J.; O’Brien, J. P.; Patricola, C. M.; Huang, H.; Ullrich, P. A.; Loring, B. A framework for detection and attribution of regional precipitation change: Application to the United States historical record. Climate Dynamics, v.60, p.705-741, 2023. https://doi.org/10.1007/s00382-022-06321-1
    » https://doi.org/10.1007/s00382-022-06321-1
  • Salehi, S.; Dehghani, M.; Mortazavi, S. M.; Singh, V. P. Trend analysis and change point detection of seasonal and annual precipitation in Iran. International Journal of Climatology , v.40, p.308-323, 2020. https://doi.org/10.1002/joc.6211
    » https://doi.org/10.1002/joc.6211
  • 1 Research developed at Faculty of Social Science Education, Vietnam

Edited by

Editors: Ítalo Herbet Lucena Cavalcante & Carlos Alberto Vieira de Azevedo

Publication Dates

  • Publication in this collection
    17 June 2024
  • Date of issue
    Aug 2024

History

  • Received
    02 Jan 2023
  • Accepted
    18 Mar 2024
  • Published
    29 Apr 2024
Unidade Acadêmica de Engenharia Agrícola Unidade Acadêmica de Engenharia Agrícola, UFCG, Av. Aprígio Veloso 882, Bodocongó, Bloco CM, 1º andar, CEP 58429-140, Campina Grande, PB, Brasil, Tel. +55 83 2101 1056 - Campina Grande - PB - Brazil
E-mail: revistagriambi@gmail.com