the Creative Commons Attribution 4.0 License.
Special issue: The changing water cycle of the IndoGangetic Plain
Research article 09 Aug 2018
Research article  09 Aug 2018
Modeling the changes in water balance components of the highly irrigated western part of Bangladesh
 ^{1}Hydrology Lab, Department of Earth and Environmental Sciences, Graduate School of Science and Technology, Kumamoto University, 2401 Kurokami, Kumamoto, Japan
 ^{2}Department of Statistics, University of Rajshahi, Rajshahi 6205, Bangladesh
 ^{3}Department of Geology and Mining, University of Rajshahi, Rajshahi 6205, Bangladesh
 ^{4}Institute of Bangladesh Studies, University of Rajshahi, Rajshahi 6205, Bangladesh
 ^{1}Hydrology Lab, Department of Earth and Environmental Sciences, Graduate School of Science and Technology, Kumamoto University, 2401 Kurokami, Kumamoto, Japan
 ^{2}Department of Statistics, University of Rajshahi, Rajshahi 6205, Bangladesh
 ^{3}Department of Geology and Mining, University of Rajshahi, Rajshahi 6205, Bangladesh
 ^{4}Institute of Bangladesh Studies, University of Rajshahi, Rajshahi 6205, Bangladesh
Correspondence: A. T. M. Sakiur Rahman (shakigeo@gmail.com)
Hide author detailsCorrespondence: A. T. M. Sakiur Rahman (shakigeo@gmail.com)
The objectives of the present study were to explore the changes in the water balance components (WBCs) by coutilizing the discrete wavelet transform (DWT) and different forms of the Mann–Kendall (MK) test and develop a wavelet denoise autoregressive integrated moving average (WDARIMA) model for forecasting the WBCs. The results revealed that most of the potential evapotranspiration (P_{ET}) trends (approximately 73 %) had a decreasing tendency from 1981–1982 to 2012–2013 in the western part of Bangladesh. However, most of the trends (approximately 82 %) were not statistically significant at a 5 % significance level. The actual evapotranspiration (A_{ET}), annual deficit, and annual surplus also exhibited a similar tendency. The rainfall and temperature exhibited increasing trends. However, the WBCs exhibited an inverse trend, which suggested that the P_{ET} changes associated with temperature changes could not explain the change in the WBCs. Moreover, the 8year (D3) and 16year (D4) periodic components were generally responsible for the trends found in the original WBC data for the study area. The actual data was affected by noise, which resulted in the ARIMA model exhibiting an unsatisfactory performance. Therefore, wavelet denoising of the WBC time series was conducted to improve the performance of the ARIMA model. The quality of the denoising time series data was ensured using relevant statistical analysis. The performance of the WDARIMA model was assessed using the Nash–Sutcliffe efficiency (NSE) coefficient and coefficient of determination (R^{2}). The WDARIMA model exhibited very good performance, which clearly demonstrated the advantages of denoising the time series data for forecasting the WBCs. The validation results of the model revealed that the forecasted values were very close to actual values, with an acceptable mean percentage error. The residuals also followed a normal distribution. The performance and validation results indicated that models can be used for the shortterm forecasting of WBCs. Further studies on different combinations of wavelet analysis are required to develop a superior model for the hydrological forecasting in the context of climate change. The findings of this study can be used to improve water resource management in the highly irrigated western part of Bangladesh.
 Article
(3402 KB) 
Supplement
(4673 KB)  BibTeX
 EndNote
The water balance model is considerably important for water resource management, irrigation scheduling, and crop pattern designing (Kang et al., 2003; Valipour, 2012). The model can also be used for the reconstruction of catchment hydrology, climate change impact assessment, and streamflow forecasting (Alley, 1985; Arnall, 1992; Xu and Halldin, 1996; Molden and Sakthivadivel, 1999; Boughton, 2004; Anderson et al., 2006; Healy et al., 2007; Moriarty et al., 2007; Karimi et al., 2013). Therefore, accurately forecasting the water balance components (WBCs) and detecting the changes in them is important for achieving sustainable water resource management. However, hydrometeorological time series are contaminated by noises from hydrophysical processes. This affects the accuracy of the analysis, simulation, and forecasting (Sang et al., 2013; Wang et al., 2014). Hence, denoising the time series is essential for improving the accuracy of the obtained results. In this study, the wavelet denoising technique was coupled with the ARIMA (autoregressive integrated moving average) model for forecasting the WBCs after detecting the changes in them by using different forms of the Mann–Kendall (MK) test. Moreover, the time period responsible for the trends in the WBC time series was identified using discrete wavelet transform (DWT) time series data.
Physicsbased numerical models are generally used for understanding a particular hydrological system and forecasting the water balance or water budget components (Fulton et al., 2015; Leta et al., 2016). To achieve reliable forecasting using numerical models, a large amount of hydrological data is required for assigning the physical properties of the grid and model parameters and calibrating the model simulation. However, numerical models have numerous limitations, such as the cost, time, and availability of the data (Yoon et al., 2011; Adamowski and Chan, 2011). Databased forecasting models and statistical models are suitable alternatives for overcoming these limitations. The most common statistical methods for hydrological forecasting are the ARIMA model and multiple linear regression (Young, 1999; Adamowski, 2007). Many studies have used the ARIMA model to predict water balance input parameters, such as rainfall (Rahman et al., 2016), temperature (Nury et al., 2016), and potential evapotranspiration (P_{ET}; Valipour, 2012). However, the ARIMA model cannot handle nonstationary hydrological data without preprocessing the input time series data (Tiwari and Chatterjee, 2010; Adamowski and Chan, 2011). Wavelet analysis, a new method in the area of hydrological research, can be used to effectively handle nonstationary data (Adamowski and Chan, 2011). Adamowski and Chan (2011) coupled wavelet analysis with artificial neural network (ANN) models for forecasting hydrological variables, such as the groundwater level, in Quebec, Canada. Kisi (2008), Partal (2010), and Santos and da Silva (2014) developed hybrid wavelet ANN models for monthly and daily streamflow forecasting. Rahman and Hasan (2014) found that the performance of waveletbased ARIMA models was superior to that of classical ARIMA models for forecasting the humidity of the Rajshahi meteorological station in Bangladesh. A comparative study of wavelet ARIMA models and wavelet ANN models was conducted by Nury et al. (2017). The study indicates that the wavelet ARIMA models are more effective than wavelet ANN models for temperature forecasting. Khalek and Ali (2016) developed the wavelet seasonal ARIMA (WSARIMA) and wavelet neural network autoregressive (WNNAR) models for forecasting the groundwater level. They observed that the WSARIMA model exhibited a superior performance to the WNNAR model. In all the aforementioned studies, the performance of the waveletaided model was better than that of the classical ARIMA and ANN models. Moreover, analyzing the periodicity using wavelettransformed details and using the approximation components of the hydrometeorological time series data can provide insight regarding the effects of the time period on the data trend (Nalley et al., 2013; Araghi et al., 2015; Pathak et al., 2016). As a result, detecting the periodicity through the wavelet transformation of hydrometeorological time series data has gained popularity in recent years (Partal and Küçük, 2006; Partal, 2009; Nalley et al., 2013; Araghi et al., 2015; Pathak et al., 2016). Studies have been conducted on the spatiotemporal characteristics of hydrometeorological variables, such as rainfall (Shahid and Khairulmaini, 2009; Ahasan et al., 2010; Kamruzzaman et al., 2016a; Rahman and Lateh, 2016; Rahman et al., 2016; Syed and Al Amin, 2016), temperature (Shahid, 2010; Nasher and Uddin, 2013; Rahman, 2016; Syed and Al Amin, 2016; Kamruzzaman et al., 2016a), and P_{ET} (Hasan et al., 2014; Acharjee, 2017), in Bangladesh. Karim et al. (2012) studied the WBCs, such as the P_{ET}, A_{ET}, deficit of water, and surplus of water, of 12 districts in Bangladesh. Kanoua and Merkel (2015) studied the water balance of Titas Upazila (subdistrict) in Bangladesh. Most of the studies conducted on hydrological variables in Bangladesh were limited to detecting trends and forecasting the rainfall and temperature. Therefore, this study was conducted to detect the trends and identify the periodicities in the WBCs, such as the potential evapotranspiration (P_{ET}), actual evapotranspiration (A_{ET}), and annual deficit and surplus of water, by coutilizing the DWT and different forms of the MK test in the western part of Bangladesh. Moreover, a wavelet denoise (WD)ARIMA model was developed for forecasting the WBCs. To date, no comprehensive study has coupled wavelet denoising methods with ARIMA models for forecasting the WBCs. Wavelet denoising methods are widely used in the engineering and scientific fields. However, these methods have been used to a limited extent in hydrology (Sang, 2013). The combination of wavelet denoising methods with ARIMA models is expected to provide insight regarding WBCs, which would ultimately help policymakers prepare sustainable water resource management plans.
2.1 Study area
The climate of Bangladesh is humid, warm, and tropical. The western part of Bangladesh covers approximately 41 % or 60 165 km^{2} of the country. The geographic coordinates of the study area extend between a latitude of 21^{∘}36^{′}–26^{∘}38^{′} N and longitude of 88^{∘}19^{′}–91^{∘}01^{′} E. The annual rainfall and average temperature in the study area vary from 1492 to 2766 mm, with an average of 1925 mm, and 24.18 to 26.17 ^{∘}C, with an average of 25.44 ^{∘}C, respectively (Kamruzzaman et al., 2016a). Bangladesh is the fourthlargest producer of rice in the world (Scott and Sharma, 2009), and the livelihood of a majority of the people (approximately 75 %) (Shahid and Behrawan, 2008; Kamruzzaman et al., 2016b) depends on agricultural practices. The crop calendar of Bangladesh is related to the climatic seasons. Rice is grown during three seasons (Aus, Aman, and Boro) in Bangladesh. Almost 73.94 % of the cultivable area in the country is used to cultivate Boro rice (Banglapedia, 2003). The Aus and Aman rice varieties are mainly rainfed crops. However, Boro rice is almost completely groundwaterfed (Ravenscroft et al., 2005) and requires approximately 1 m of water per square meter in Bangladesh (Harvey et al., 2006; Michael and Voss, 2009).
2.2 Data
The national climate database of Bangladesh prepared by the Bangladesh Agricultural Research Council (BARC) was used for this study. The database is available for research and can be obtained from the BARC website (http://climate.barcapps.gov.bd/, last access: 27 July 2018). The database has been prepared from the data recorded by the Bangladesh Meteorological Division and contains longterm monthly climate data, such as rainfall, minimum, maximum, and average temperatures, humidity, sunshine hours, wind speed, and cloud cover. The locations of the meteorological stations in the study area are displayed in Fig. 1. The data are rearranged according to the hydrological year for the period from 1981–1982 to 2012–2013. The hydrological year in Bangladesh begins in April and ends in March.
2.3 Methods
In this study, the WBCs were calculated and their trends were identified using the MK or Modified MK (MMK) test for evaluating the longterm water balance of the highly irrigated western part of Bangladesh. The DWT data of the WBC time series were analyzed for identifying the time period responsible for the trend in the data. The WBCs were forecasted using the ARIMA model, whose performance was statistically evaluated. If the performance of the model was unsatisfactory for forecasting the WBCs, denoising of the original time series was conducted using DWT techniques to improve the performance of the model. The descriptions of the methods are presented in the following sections.
2.3.1 Calculation of the potential evapotranspiration and water balance components
The potential evapotranspiration (P_{ET}) is a key parameter to estimate the WBCs. In this study, the potential evapotranspiration was calculated using the Penman–Monteith equation (Allen et al., 1998). The soil–water balance concept proposed by Thornthwaite and Mather (1955) is one of the most widely used methods for estimating the WBCs. This method is suitable for assessing the effectiveness of agricultural water resource management practices and regional water balance studies because it allows the actual evapotranspiration (A_{ET}), water deficit, and water surplus to be estimated (Chapman and Brown, 1966; Bakundukize et al., 2011; Karim et al., 2012; Viaroli et al., 2017). The actual evapotranspiration (A_{ET}) is the amount of water removed from the surface due to evaporation and transpiration. The amount by which the P_{ET} exceeds the A_{ET} is termed as the deficit. The surplus is the excess rainfall received after the soil has reached its waterholding capacity (de Jong and Bootsma, 1997). Calculating the field capacity of the soil is essential for estimating the WBCs. The field capacity of the soil in the study area was calculated using the soil texture map of Bangladesh prepared by the Soil Resource Development Institute, Bangladesh (SRDI, 1998), where the description of the soils was presented by Huq and Shoaib (2013). The values suggested by Thornthwaite and Mather (1957) for the waterholding capacity of the soil and rooting depth of the plants were used for estimating the WBCs in this study. The first step of the calculation involves subtracting 5 % rainfall from the monthly rainfall data because this amount of water is lost due to direct runoff (Wolock and McCabe, 1999; Karim et al., 2012; Kanoua and Merkel, 2015). The remaining rainfall amount is included in the calculation. The WBCs, such as the A_{ET}, surplus, and deficit, were estimated using the formulas presented in Table 1. The details of the WBC calculations are available in the Supplement.
2.3.2 Trend test
In this study, the trends in the WBCs were detected using the nonparametric MK test (Mann, 1945; Kendall, 1975) because it exhibits a better performance than the parametric test (Nalley et al., 2012) for identifying trends in hydrological variables, such as rainfall (Shahid, 2010), temperature (Kamruzzaman et al., 2016a), P_{ET} (Kumar et al., 2016), soil moisture (Tabari and Talaee, 2013), runoff (Pathak et al., 2016), groundwater level (Rahman et al., 2016), and water quality (Lutz et al., 2016). The MK test cannot be used to accurately calculate the test statistic (Z) if there exists a significant serial correlation at lag 1 in the time series data (Yue et al., 2002) because the variance is underestimated (Hamed and Rao, 1998). The autocorrelation at lag 1 was checked before analyzing the time series data. If there existed a significant lag1 autocorrelation at the 5 % level, the MMK test (Hamed and Rao, 1998) was applied instead of the MK test. The estimated Z statistic from the MK or MMK test was evaluated for the direction of the trend (a positive Z statistic indicated an increasing trend and vice versa). Moreover, the Z statistic indicated the level of significance of the obtained trend. If the calculated Z statistic is equal to or higher than the tabulated value of the Z statistic (+1.96), it indicates a significant positive trend at the 95 % confidence level. If the calculated Z statistic is equal to or less than −1.96, it indicates a significant decreasing trend. Moreover, the sequential values of the u(t) statistic derived from the sequential MK (SMK) test (Sneyers, 1990) are used for detecting the change point. The u(t) statistic is similar to the Z statistic (Partal and Küçük, 2006). The magnitude of the change was calculated using Sen's slope estimator (Sen, 1968). Numerous studies have already been conducted (notably Nalley et al., 2012) using the methods described in this section. Further details regarding these methods can be obtained from Mann (1945), Sen (1968), Kendall (1971), Hamed and Rao (1998), Sneyers (1990), and Yue et al. (2002).
2.3.3 Wavelet transform (WT) and periodicity
Wavelet analysis has been used in different parts of the world to identify the periodicity in hydroclimatic time series data (Smith et al., 1998; Azad et al., 2015; Nalley et al., 2012; Araghi et al., 2015; Pathak et al., 2016). WT, a multiresolution analytical approach, can be applied to analyze time series data because it offers flexible window functions that can be changed over time (Nievergelt, 2001; Percival and Walden, 2000). WT can be applied to detect the periodicity in hydroclimatic time series data (Smith et al., 1998; Pišoft et al., 2004; Sang, 2012; Torrence and Compo, 1998; Araghi et al., 2015; Pathak et al., 2016) and exhibits better a performance than traditional approaches (Sang, 2013). There exist two main types of WT, namely continuous WT (CWT) and DWT. Applying the CWT is complex because it produces numerous coefficients (Torrence and Compo, 1998; Araghi et al., 2015), whereas DWT is simple and useful for hydroclimatic analysis (Partal and Küçük, 2006; Nalley et al., 2012). The wavelet coefficients of the DWT with a dyadic format can be calculated as follows (Mallat, 1989):
where ψ is the mother wavelet, m is the wavelet dilation, and n is the wavelet translation. The specified fixed dilation step (s_{o}) is larger than 1, and τ_{o} is the location parameter. For practical application, the values of s_{o} and τ_{o} are considered as 2 and 1, respectively (Partal and Küçük, 2006; Pathak, 2016). After substituting these values in Eq. (1), the DWT for a time series x_{i} becomes the following:
where W indicates the wavelet coefficient at a scale s=2^{m} and location τ=2^{m}n
In the DWT, details (D) and approximations (A) of the time series can emerge from the original time series after passing through lowpass and highpass filters, respectively. When approximations are the highscale and lowfrequency components, details are the lowscale and highfrequency components. Successive iterations are performed to decompose the time series into its several lowresolution components (Mallat, 1989; Misiti et al., 1997). In this study, four levels (D1–D4) of decomposition were performed following the dyadic scales. The decompositions are referred to as D1, D2, D3, and D4, which correspond to a 2, 4, 8, and 16year periodicity, respectively. The Daubechies wavelet was used because of its superior performance in hydrometeorological studies (Nalley et al., 2012, 2013; Ramana et al., 2013; Araghi et al., 2015). To confirm the periodicity present in the time series, the correlation coefficient (Co) between u(t) of the original data, u(t) of the decomposition (D) time series data, and different models (D1 + A … D4 + D3 + A) of the time series data were calculated and the obtained results were compared (Partal and Küçük, 2006; Partal, 2009).
2.3.4 ARIMA models
ARIMA models (Box and Jenkins, 1976) are used in hydrological science to identify the complex patterns in data and project future scenarios (Adamowski and Chan, 2011; Valipour et al., 2013; Nury et al., 2017; Khalek and Ali, 2016). ARIMA models include (1) an autoregressive process (AR) represented by order p, (2) nonseasonal differences for nonstationary data termed as order d, and (3) a moving average (MA) process represented by order q. An ARIMA model of order $(p,d,q)$ can be written as follows:
where θ_{0} is the intercept with a mean of 0, U_{t} is the white process with constant variance, ∅_{p}(L) represents the AR term $(\mathrm{1}{\mathrm{\varnothing}}_{\mathrm{1}}L\phantom{\rule{0.25em}{0ex}}\mathrm{\dots}\phantom{\rule{0.25em}{0ex}}{\mathrm{\varnothing}}_{p}{L}^{p})$, and θ_{q}(L) represents the MA term $(\mathrm{1}{\mathit{\theta}}_{\mathrm{1}}L\phantom{\rule{0.25em}{0ex}}\mathrm{\dots}\phantom{\rule{0.25em}{0ex}}{\mathit{\theta}}_{p}{L}^{p})$.
2.3.5 Wavelet denoising
Wavelet denoising based on the thresholds introduced by Donoho et al. (1995) has been applied to hydrometeorological analysis (Wang et al., 2005, 2014; Chou, 2011). In this study, the following three analysis steps were performed for denoising the time series data.

Decomposing the time series data x(t) into M resolution levels for obtaining the detail coefficients (W_{j,k}) and approximation coefficients using the DWT.

The detail coefficients obtained from the DWT (1 to M levels) were treated using threshold (T_{j}) selection. A soft or hard threshold can be used to deal with detail coefficients and obtain the decomposed coefficient. In this study, a soft threshold was selected because it performed better than a hard threshold (Wang et al., 2014; Chou, 2011).
$$\begin{array}{ll}{\displaystyle}& {\displaystyle}\text{Softthresholdprocessing:}\\ {\displaystyle}& {\displaystyle}W{{}^{\prime}}_{j,k}=\left\{\begin{array}{ll}\mathrm{sgn}\left({W}_{j,k}\right)\left(\left{W}_{j,k}\right{T}_{j}\right)& \left{W}_{j,k}\right{T}_{j}\\ \mathrm{0}& \left{W}_{j,k}\right{T}_{j}\end{array}\right..\end{array}$$ 
Detail coefficients from levels 1 to M and approximate coefficients at level M were reconstructed to obtain denoising time series data.
Selecting the threshold value is essential for denoising the data. In this study, the universal threshold (UT) method (Donoho and Johnstone, 1994) was used for estimating the threshold value because it exhibited satisfactory performance in analyzing hydrometeorological data (Wang et al., 2005; Chou, 2011).
2.3.6 Assessment of model performance
There exist several indicators to assess the performance of the models. The Nash–Sutcliffe efficiency (NSE) (Nash and Sutcliffe, 1970) coefficient, a normalized goodnessoffit statistic, is the most powerful and popular method for measuring the performance of hydrological models (McCuen et al., 2006; Moussa, 2010; Ritter and MuñozCarpena, 2013). The NSE coefficient was used in this study to evaluate and compare the ARIMA and WDARIMA models. The NSE is calculated as follows (Nash and Sutcliffe, 1970):
where N, O_{i}, P_{i}, $\stackrel{\mathrm{\u203e}}{O}$, and SD are the sample size, number of observations, model estimates, mean, and standard deviation of the observed values, respectively. The performance of a model can be evaluated according to its NSE value as very good (NSE≥0.9), good (NSE = 0.8–0.9), acceptable (NSE≥0.65), and unsatisfactory (NSE < 0.65) (Ritter and MuñozCarpena, 2013). E_{RMS} is the rootmeansquare error and can be calculated as follows:
The coefficient of determination (R^{2}) is another goodnessoffit test to measure the performance of models. The perfect fit of the model draws a line between the actual values and fitted values, where R^{2} is 1. If y_{i} is the observation data, ${\widehat{y}}_{i}$ represents the modelforecasted values of y_{i} and N is the number of data points used. R^{2} is given as follows (Sreekanth et al., 2009):
Moreover, the mean percentage error (E_{MP}) and mean error (E_{M}) were also calculated to evaluate the validation of the model for forecasting. E_{MP} indicates the percentage of bias (large or small) between the forecasted and actual data (Khalek and Ali, 2016). E_{MP} and E_{M} can be calculated as follows:
3.1 Exploratory statistics of the water balance components
The mean annual P_{ET} in the study area between 1981–1982 and 2012–2013 varied from 1228 to 1460 mm (Fig. 2a), with an average of 1338 mm. High P_{ET} values were observed in the central part of the area, where the annual rainfall was low, but the temperature was high (Kamruzzaman et al., 2016a). The standard deviations of the P_{ET} varied from 205 (Jessore station) to 41 mm (Bhola station). The A_{ET} (Fig. 2b) (average = 925 mm) was almost 31 % less than the P_{ET} because during the dry months (December–May), the soil moisture condition reached a critical stage. The annual surplus of water varied from 515 to 1277 mm (Fig. 2c), with an average of 838 mm. According to Wolock and McCabe (1999), 50 % of the surplus water can be considered as runoff for the major parts of the world. A high amount of surplus water was found in the northern part of the study area and along the coastal area. The annual deficit of water, which mainly occurred during the dry season (December–May), varied from 329 to 556 mm, with an average of 416 mm (Fig. 2d). The highest annual deficit of water was observed in Rajshahi, which is located in the central–western part of the study area, where the depth of groundwater below the surface increases rapidly (Shamsudduha et al., 2009; Rahman et al., 2016).
3.2 Trend and periodicity of the water balance components
3.2.1 Potential evapotranspiration
The MK or MMK test based on lag1 autocorrelation was applied to detect the trend in the P_{ET}. Table 2 represents the Z statistic of the MK or MMK test for the original P_{ET} time series data and the Z statistic of the decomposition (D1–D4), approximation (A), and model (D1 + A … D3 + D4 + A) time series. The estimated Z statistic of the original data ranged from −2.07 (Satkhira station) to 2.37 (Bhola station). The Satkhira and Bhola stations exhibited significant P_{ET} trends. The plots of the sequential u(t) statistic obtained from the SMK test for these two stations are displayed in Fig. 3, where the dashed lines correspond to a 5 % significance level (±1.96). The decreasing P_{ET} trend for the Satkhira station began in 1985–1886, and a significant decreasing trend occurred in 1993–1994. The trend reversed after 2007–2008. However, the significant increasing P_{ET} trend of the Bhola station began very recently (2010–2011) after some fluctuation.
Most of the trends (73 %) observed in the P_{ET} time series data of the study area were negative and statistically insignificant at the 95 % confidence level or 5 % significance level. Moreover, the Z statistic of the approximation (A) time series obtained using the DWT indicated decreasing P_{ET} trends for all the stations. The calculated Z statistic of the approximation (A) time series was approximately −1.8 after rounding the figures for all the stations. The approximation time series data of all the stations exhibited a similar pattern (Fig. S1 of the Supplement) over time. The magnitude of P_{ET} changes ranged from −10.89 mm yr^{−1} for the Satkhira station to 1.67 mm yr^{−1} for the Bhola station (Fig. 4a). The MK or MMK test was also applied to the decomposition time series and model time series generated from the combination of the approximation and decomposition time series data. Table 2 represents the results for four stations arranged in alphabetical order, and the complete results can be found in Table S1 of the Supplement. To determine the dominant periodicity affecting the P_{ET} trends, a twostep analysis was performed. First, the value closest to the Z statistic of the original time series data was obtained from the Zstatistic values of different model and decomposition time series data. Second, the correlation coefficients (Co) of pairs of data (such as the Co between the u(t) statistics obtained from the SMK test for the original and decomposition time series data) were estimated, and the highest Co was determined from the estimated Co values for different pairs (Table 2). The Z statistic of the D4 time series data for the Barisal station was 0.76, which was the closest to the Z statistic (0.72) of the original time series data (Table 2). Moreover, the Z statistic of the model (D3 + D4 + A) time series data was 0.56, which is the secondnearest value to the Z statistic of the original time series and has the highest correlation coefficient (Co = 0.85). The D4 (16year) component was the dominant periodic component in the trend of the original data. However, D3 also affected the trend of the data. The Zstatistic value (2.47) of the original time series for the Bhola station was the closest to that (2.36) of the model (D2 + D4 + A) time series data. However, the Zstatistic values of the D2, D4, D2 + A, and D4 + A time series were 0.61, 1.2, 0.48, and 0.9, respectively. These values were not close to the Z statistic of the original time series data. Hence, in this case, the Z statistic was unable to determine which periodic component (D2∕D4) was the basic periodic component for the significant trend in the original data. To determine the dominant periodic component, the values of Co were analyzed. The correlation coefficient (Co) between the u(t) statistic of the SMK test for the original and D4 time series data was higher than the correlation coefficient between the u(t) statistic of the SMK test for the original and D2 time series data (Table 2). Moreover, the values of the Z statistic for time series with the D4 components, such as the D4 and D4 + A model time series, were higher than those for time series with the D2 component (D2 and D2 + A) (Table 2). Therefore, D4 was the main periodic component responsible for the P_{ET} trend of the Bhola station. However, the Zstatistic values of D4 and D4 + A were not close to the Z statistic of the original data (Table 2). Moreover, there existed a statistically significant positive trend in the original P_{ET} data of the Bhola station, whereas the trends of the D4 and D4 + A model time series data were not statistically significant. When the D2 time series was added to the D4 + A model time series data, the Z statistic of the resultant (D2 + D4 + A) model time series data was very close to that of the original time series data. The trend of the D2 + D4 + A model time series was statistically significant, similar to the trend in the original time series data (Table 2). Hence, D2 affected the trend of the original time series data. Stationwise analysis indicated that almost half of the stations exhibited harmony between the Zstatistic values of the D3 + D4 + A model and original time series data. Individual analysis of the D3 and D4 time series data indicated that a higher relationship existed between the D4 and original time series data. Three stations (Dinajpur, Ishurdi, and Jessore) exhibited similar Zstatistic values for the original and D1 + D4 + A model time series data, with higher Co values of the u(t) statistic for the SMK test on the D4 time series data than that for the SMK test on the original data (except for the Ishurdi station). Moreover, two stations (Bhola and Satkhira) exhibited significant trends in the original data. The closest Z statistic was found between the original and D2 + D4 + A time series data for both of the stations. D4 (16year periodicity) was the dominant periodic component according to the Co values for both these stations. Therefore, 16year periodicity was the main periodic component responsible for the trends in the P_{ET} data over the study area. Moreover, D3 (8year) periodicity also had an effect on the trends for some stations (Tables 2 and S1). D4 (16year) periodicity dominates the annual rainfall trend for the Marmara region in Turkey (Partal and Küçük, 2006). Araghi et al. (2015) determined that 8–16year (D3 to D4) periodicity is responsible for the trends in the annual temperature in Iran.
3.2.2 Actual evapotranspiration
All the stations except the Bogra station exhibited decreasing trends in the A_{ET}. The calculated Z statistic ranged from −2.90 for the Bogra station to 0.31 for the Ishurdi station. Similar to the P_{ET} trends, the A_{ET} trends were also insignificant at a 5 % significance level. However, the Ishurdi station exhibited a significant (at a 5 % significance level) decreasing trend. The magnitudes of the trends of the original A_{ET} data varied from −5 mm yr^{−1} for the Faridpur station to 0.75 mm yr^{−1} for the Bogra station. The distribution of the trend magnitude is displayed in Fig. 4b. The periodicity in the A_{ET} was marginally different from that in the P_{ET} (Table S2). For almost half of the stations (five), D2 (4year) was the main periodic component. D4 (16year) also affected the trend because the Z statistic of the D2 + D4 + A model time series was the nearest to that of the original series for the Khulna and Ishurdi stations. Moreover, D4 (16year) was the main periodic component for the Rangpur and Rajshahi stations. D1 (2year) was the dominant periodic component for the Barisal, Bhola, and Bogra stations. The A_{ET} value depends on climatic factors, such as the P_{ET}, rainfall, and soil moisture conditions. The variations in the periodicities of the A_{ET} and P_{ET} were mainly related to the soil moisture conditions of the area.
3.2.3 Surplus
Almost 82 % of the stations exhibited insignificant decreasing trends for the annual surplus of water. The magnitude of the trends of the original annual surplus data ranged from −11.63 to 6.71 mm yr^{−1} (Fig. 4c). The periodicity characteristics of the P_{ET} and surplus were similar (Table S3). D4 (16year) was the main periodic component present in seven stations. In most cases, D2 was also present (D2 + D4 + A), except in Rajshahi. D3 (8year) was mainly responsible for the surplus trend of three stations. Surplus mainly occurred during the rainy season (June–October) in the study area, when the soil pores were almost completely filled with water and the A_{ET} was equal to the P_{ET}. Surplus mainly depends on rainfall and hence provides insight regarding the periodicity in rainfall.
3.2.4 Deficit
Approximately 73 % of the stations exhibited increasing trends for the annual deficit of water. The increasing trends were significant for two stations at the 95 % confidence level (Table S4). However, the Satkhira station exhibited a significant decreasing trend (Z $=$2.08) in the annual deficit of water. The magnitude of the trends of the original annual deficit data ranged from −8.1 to 7.7 mm yr^{−1} (Fig. 4b). Periodicity analysis revealed that D4 was mainly responsible for the trends in the annual deficit of water. The Z statistic of the (D2 + D4 + A) model time series data was close to the Z statistic of the original time series data (Table S4). D3 (8year periodicity) was also responsible for the trends in the data of the two stations.
3.3 Model selection and forecasting ability
The ARIMA model was selected for forecasting the WBC time series. A fourstep analysis was performed during time series modeling. (1) First, the stationarity of the data was checked using the Augmented Dickey–Fuller (ADF) test. (2) Then, the autocorrelation function (ACF) was used for selecting the order of the MA process (Figs. S2–S5). (3) The partial autocorrelation function (PACF) was then used for selecting the order of the AR process (Figs. S2–S5). (4) Finally, the appropriate model was selected based on several trials and model selection criteria, such as Akaike information criterion (AIC) and Bayesian information criterion (BIC). In addition to the manual model selection based on the ACF, PACF, AIC, and BIC, the auto ARIMA function of the “forecast” package (Hyndman et al., 2017) of R (R 3.4.0 language developed by R Core Team, 2016) was used during the trails for model selection to obtain information regarding the nature of the data for modeling. The model with the lowest AIC and BIC values and highest R^{2} value was selected. The Q–Q plot was prepared to examine the normality of the residuals. The performance of the ARIMA model (parameters are given in Table S5) was evaluated using the NSE coefficient and R^{2} values (Table 3). The estimated NSE coefficient of the ARIMA model for the P_{ET} time series varied from −0.6 for the Bhola station to 0.81 for the Jessore station (Table 3). The ARIMA model exhibited an unsatisfactory performance for almost all the stations. The average NSE coefficient of the 11 stations was 0.38, and the R^{2} values ranged from 0.1 to 0.81, with an average of 0.38. Moreover, the NSE coefficient of the Bhola station indicated that the ARIMA model was unsuitable for forecasting the P_{ET}. The ARIMA model was also applied to the A_{ET}, surplus, and deficit time series data. There existed no significant spikes in the ACF and PACF of the A_{ET} (Fig. S3). Moreover, the results obtained from the auto ARIMA functions exhibited similar results. Therefore, the ARIMA model was unsatisfactory for forecasting the variability in the A_{ET}. For WBCs such as surplus and deficit, the performance of the ARIMA model was similar to that of the A_{ET}, except for a few cases. Because hydrometeorological data are affected by noises from different hydrophysical processes (Wang et al., 2014), the results obtained using the ARIMA models were unsatisfactory. To improve model performance, noise must be removed from the data. In this study, DWT denoising was applied to the WBC data and the quality of the denoising time series data was examined before further processing. When selecting a method for denoising the time series using WT, the mean of the original and denoising time series data should be close and the standard deviation of the denoising time series should be less than that of the original time series (Wang et al., 2014). Figure 5a displays the means of the actual and wavelet denoising P_{ET} time series. No visible difference was observed between the mean of the original and DWT wavelet denoising time series data. Moreover, the standard deviation of the P_{ET} for the wavelet denoising time series was lower than that for the original time series (Fig. 5b). The A_{ET}, surplus, and deficit time series also exhibited similar results. Furthermore, the lag1 autocorrelation of the wavelet denoising time series data must be higher than that of the original time series (Wang et al., 2014). Under this condition, the absolute lag1 value of autocorrelation for the wavelet denoising time series was higher than that for the original series (Figs. S2b, S3b, S4b, and S5b). The performance of the WDARIMA model is represented in Table 3. After denoising the data, the performance of the ARIMA model was satisfactory for all the WBC time series data (Table 3). The average NSE coefficient of the WDARIMA model for the P_{ET} time series of the 11 stations located in the western part of Bangladesh was 0.76, with an average R^{2} value of 0.67. The R^{2} and NSE coefficient values indicated that the performance of the WDARIMA model was better than that of the classical ARIMA model for the modeling of P_{ET} (Table 3). Moreover, the average NSE value of the WDARIMA model for the A_{ET} time series of the 11 stations was 0.92, which indicated that the performance of the model was very good. The average R^{2} value was 0.89, which indicated that the model could explain almost 89 % of the variance in the data (Table 3). The WDARIMA model also exhibited a very good forecasting performance for the annual surplus and deficit (Table 3). The average NSE coefficient of the WDARIMA model for the annual surplus of the 11 stations was approximately 0.92, and the average R^{2} value was 0.9. The WDARIMA model exhibited a good performance in forecasting the annual deficit (average NSE = 0.88). The performance of the WDARIMA model was good or very good for forecasting the A_{ET}, annual surplus, and annual deficit. However, the performance was acceptable for forecasting the P_{ET}. This deviation may have arisen because the variability of the P_{ET} was higher than that of the other WBCs or the deviation may be related to the variability of climatic variables.
The WDARIMA models were validated to explore their forecasting ability. The mean percentage error (E_{MP}) of the forecasted values for the 4year period from 2008–2009 to 2012–2013 was calculated to determine the percentage bias of the forecasted data (Table 4). The average E_{MP} of the WDARIMA model for the P_{ET} values of the 11 stations was −0.6 (ranging from 0.75 to −3.34), which indicated that the forecasted values were marginally lower than the actual values. The typical plots of the actual time series data versus the fitted model data, normal Q–Q plots of the residuals of the models, and actual and observed values for the WBCs (plots for all the stations are displayed in Figs. S6–S9) are illustrated in Fig. 6. The plot of the actual values versus the forecasted values (Fig. 6) indicates that the actual and forecasted values were very close for the hydrologic years 2009–2010 and 2010–2011. The normal Q–Q plots revealed that the residuals of the models were near normal. However, the differences in the values increased after these two hydrologic years for all the WBCs (Figs. S6–S9). The E_{MP} values of WDARIMA models for the A_{ET} ranged from −0.7 to 0.2, with an average of −0.09, which indicated that the forecasted A_{ET} values were marginally lower than the actual A_{ET} values. The E_{MP} values for the annual surplus (average $=$0.75) and annual deficit (average $=$0.12) were similar to that for the A_{ET} and P_{ET}. The average E_{MP} values for all the WBCs were negative, which indicated that the forecasted values for the WBCs were marginally lower than the actual values for most of the stations.
3.4 Discussion
This study indicated that a decreasing P_{ET} trend dominated the study area. However, positive trends in the rainfall and temperature dominated the western part of Bangladesh (Shahid and Khairulmaini, 2009; Kamruzzaman et al., 2016a). Moreover, a recent study found a negative trend in the evapotranspiration for four stations located in northwest Bangladesh (Acharjee et al., 2017). Although the annual rainfall and temperature of the Satkhira station exhibited positive trends (Kamruzzaman et al., 2016a), its P_{ET} exhibited a significant decreasing trend. Increasing temperature and decreasing P_{ET} trends were observed in the Yunnan Province of South China (Fan and Thomas, 2012). McVicar et al. (2012) also found decreasing P_{ET} trends in different parts of the world. Therefore, although the temperature is the primary factor driving changes in the P_{ET} (IPCC, 2007), temperaturebased models cannot suitably explain the causes of P_{ET} changes. To obtain a detailed insight regarding the mechanisms underlying the P_{ET} changes, a detailed analysis must be conducted of all climatic variables, such as rainfall, temperature, sunshine hours, wind speed, and humidity, and climatecontrolling phenomena, such as El Niño–Southern Oscillation.
The WDARIMA model was used in this study for forecasting the WBCs. The performance of the model indicated the benefit of denoising hydrological time series data, such as the P_{ET}, A_{ET}, surplus, and deficit. However, the NSE coefficient indicated that the performance of the model was acceptable for P_{ET} forecasting (NSE ≥ 0.65). The deviation between the forecasted values and actual values increased with increasing time steps. Therefore, the WDARIMA model was unsuitable for longterm forecasting. The WDARIMA model was developed by coupling the discrete wavelet denoising time series data and ARIMA model. The soft threshold method was selected for denoising the time series data, and the UT method was used for determining the threshold value. However, there exist other approaches, such as SURE (Stein, 1981) and MINMAX (Donoho and Johnstone, 1998), for determining the threshold value. Moreover, Wang et al. (2014) developed a hybrid method called the adaptive wavelet denoising approach using sample entropy (AWDASE) for denoising hydrometeorological time series data, such as rainfall and streamflow data. The study (Wang et al., 2014) indicated that the performance of the developed denoising method was better than that of conventional methods for denoising rainfall and streamflow data. The aforementioned approaches may be used to increase the performance of the ARIMA model for forecasting hydrological variables, such as the P_{ET}. Moreover, there exist several mother wavelet families, such as Daubechies, Haar, Coiflets, Morlet, and Mexican hat (Sang, 2013). In this study, only Daubechies 6 from the Daubechies wavelet family was applied as the mother wavelet for the DWT. The WDARIMA model exhibited very good performance for forecasting the A_{ET}, surplus, and deficit, whereas the classical ARIMA model exhibited poor performance or was unable to forecast the WBCs. Moreover, studies (Chou, 2011; Kisi, 2008; Partal, 2009; Santos and da Silva, 2014; Rahman and Hasan, 2014; Nury et al., 2016; Adamowski and Chan, 2011; Khalek and Ali, 2016) have indicated that the performance of waveletaided models is better than that of the classical ARIMA and ANN models for forecasting nonstationary hydrometeorological variables. Because traditional methods such as Wiener filtering, Kalman filtering, and Fourier transform are unsuitable for nonstationary hydrological time series data (Adamowski and Chan, 2011; Sang, 2013), wavelet denoising can be used to improve the performance of the classical ARIMA model for forecasting hydrological variables.
In this study, the changes in the WBCs were explored using various forms of the waveletaided MK test. Moreover, a waveletaided ARIMA model was used for forecasting the WBCs. The results obtained from trend analysis indicated that decreasing trends were dominant in all the WBCs in the western part of Bangladesh during the period from 1982–1983 to 2012–2013. However, most of the trends were insignificant at the 95 % confidence level. One significant positive and one significant negative P_{ET} trend was found for the Satkhira and Bhola stations, respectively. Different combinations of the D and A (i.e., D + A and D + A + A) components of the DWT were analyzed using the Co value of the u(t) statistic from the SMK test, which provides detailed information regarding the dominant periodicity and time period affecting the trend of the original data (see the Trend and periodicity section or the example of the Bhola station). The findings of this study revealed that to obtain details regarding the time period responsible for the trends in the data, different combinations of components (D + A and D + A + A) must be analyzed rather than only the details (D) or approximation (A) components of the WT data. Moreover, this study indicated that the changes in temperature and rainfall were not only associated with the changes in the P_{ET}. To determine the attributes of P_{ET} changes, a detailed analysis must be conducted of all the relevant climatic variables. In the western part of Bangladesh, the D3 (8year) and D4 (16year) components had a dominant effect on the trends in the original WBC time series data. D2 (4year) periodicity was also present in some cases, especially for the A_{ET}. Because surplus occurs during the monsoon season and most of the rainfall occurs during this season, the rainfall pattern may have a similar periodicity (D3 to D4).
Modeling of the study revealed that the WBC time series data was affected by noises from different hydrophysical interactions. As a result, the classic ARIMA model exhibited unsatisfactory performance in most of the cases (e.g., P_{ET}) or was unable to model the variability and changes in the A_{ET}, surplus, and deficit. This study indicated that the ARIMA model can be used to model the time series data of WBCs after denoising the data using DWT with a UT. The quality of the wavelet denoising time series data was evaluated, and satisfactory results were obtained for WBC data denoising. The performance of the fitted WDARIMA model was evaluated using the NSE and R^{2} values. The average NSE and R^{2} values of the 11 stations located in the western part of Bangladesh were 0.76 and 0.67, respectively, for the P_{ET}; 0.92 and 0.89, respectively, for the A_{ET}; 0.92 and 0.9, respectively, for the annual surplus; and 0.88 each for the annual deficit. The validation of the WDARIMA model for the period of 2009–2010 to 2012–2013 provided an acceptable E_{MP} value. Thus, the WDARIMA model had an acceptable to very good performance for the shortterm forecasting of WBCs. However, the gap between the actual and forecasted data increased with increasing time. The obtained results encourage further studies to determine a realistic model for realworld application under changing climate. The results of this study can be incorporated into water resource management plans for the highly irrigated western part of Bangladesh, where the groundwater resource is at a critical stage. Further studies regarding the denoising of hydrological time series data using different mother wavelets, such as Haar and Coiflet, and the determination of thresholds by using the MINMAX, SURE, or entropybased adaptive denoising approaches would enable the development of superior models for forecasting hydroclimatic time series in the context of climate change and be beneficial for sustainably managing water resources.
The national meteorological database of Bangladesh prepared by the Bangladesh Agricultural Research Council (BARC) was used to accomplish this study. Data are available for research and can be obtained from the BARC website (http://climate.barcapps.gov.bd/).
The supplement related to this article is available online at: https://doi.org/10.5194/hess2242132018supplement.
ATMSR designed and wrote the manuscript with input from all coauthors. MSA, MK, MAK, and ATMSR prepared the R code and ATMSR, MK, MAK, and MSA performed the statistical analysis. HMA and ATMSR performed the water balance analysis. QHM and CSJ supervised the whole work.
The authors declare that they have no conflict of interest.
This article is part of the special issue “The changing water cycle of the IndoGangetic Plain”. It is not associated with a conference.
We thank the two anonymous reviewers for their constructive comments that
greatly improved the manuscript. We would like to thank editor Ana Mijic of
the special issue for her comments and support for publication in Hydrology and Earth System Sciences.
Edited by: Ana Mijic
Reviewed by: two anonymous referees
Acharjee, T. K., Halsema, G., Ludwig, F., and Hellegers, P.: Declining trends of water requirements of dry season Boro rice in the northwest Bangladesh, Agr. Water Manage., 180, 148–159, 2017.
Adamowski, J.: Development of a shortterm river flood forecasting method based on wavelet analysis, Warsaw Polish Academy of Sciences, Monograph, 172, 2007.
Adamowski, J. and Chan, H. F.: A wavelet neural network conjunction model for groundwater level forecasting, J. Hydrol., 407, 28–40, 2011.
Ahasan, M. N., Chowdhary, M. A. M., and Quadir, D. A.: Variability and trends of summer monsoon rainfall over Bangladesh, J. Hydrometeorol., 7, 1–17, 2010.
Allen, R. G., Pereira, L. S., Raes, D., and Smith, M.: Crop evapotranspiration: guidelines for computing crop water requirements, FAO Irrigation and Drainage Paper, No. 56, Rome, Italy, 328 pp., 1998.
Alley, W. M.: Water balance models in onemonthahead streamflow forecasting, Water Resour. Res., 21, 597–606, 1985.
Anderson, R., Hansen, J., Kukuk, K., and Powell, B. Development of watershedbased water balance tool for water supply alternative evaluations, Proceedings of the Water Environment Federation, (WEF`06), Water Environment Federation, 2817–2830, 2006.
Araghi, A., Baygi, M. M., Adamowski, J., Malard, J., Nalley, D., and Hasheminia, S. M. Using wavelet transforms to estimate surface temperature trends and dominant periodicities in Iran based on gridded reanalysis data, Atmos. Res., 155, 52–72, https://doi.org/10.1016/j.atmosres.2014.11.016, 2015.
Arnall, N. W.: Factors controlling the effects of climate change on river flow regimes in a humid temperate environment, Journal Hydrogeology, 132, 321–342, 1992.
Azad, S., Debnath, S., and Rajeevan, M. Analysing predictability in Indian monsoon rainfall: a data analytic approach, Environ. Process., 2, 717–727, 2015.
Bakundukize, C., Camp, M. V., and Walraevens, K.: Estimation of groundwater recharge in Bugesera region (Burundi) using soil moisture budget approach, Geol. Belg., 14, 85–102, 2011.
Banglapedia: National Encyclopedia of Bangladesh, Asiatic Society of Bangladesh, Dhaka, 2003.
Boughton, W.: Catchment water balance modelling in Australia 1960–2004, Agr. Water Manage., 71, 91–116, 2004.
Box, G. E. P. and Jenkins, G. M.: Time Series Analysis: Forecasting and Control (Revised edition), San Francisco, Holden Day, 1976.
Chapman, L. C. and Brown, D. M.: The climates of Canada for agriculture. Canada Land Inventory Report No. 3, Environment Canada, Lands Directorate, 24, 1966.
Chou, C.: A threshold based wavelet denoising method for hydrological data modelling, Water Resour. Manag., 25, 1809–1830, https://doi.org/10.1007/s1126901197763, 2011.
de Jong, R. and Bootsma, A.: Estimates of water deficits and surpluses during the growing season in Ontario using the SWATRE model, Can. J. Soil Sci., 77, 285–294, 1997.
Donoho, D. L.: Denoising by softthresholding, IEEE T. Inform. Theory, 41, 613–627, 1995.
Donoho, D. L. and Johnstone I. M.: Ideal Denoising in an Orthonormal Basis Chosen from a Library of Bases, CR Acad. Sci. IMath., 319, 1317–1322, 1994.
Donoho, D. L. and Johnstone, I. M.: Minimax estimation via wavelet shrinkage, Ann. Stat., 26, 879–921, 1998.
Fan, Z. and Thomas, A.: Spatiotemporal variability of reference evapotranspiration and its contributing climatic factors in Yunnan Province, SW China, 1961–2004, Climatic Change, 116, 309–325, https://doi.org/10.1007/s1058401204794, 2012.
Fulton, J. W., Risser, D. W., Regan, R. S., Walker, J. F., Hunt, R. J., Niswonger, R. G., Hoffman, S. A., and Markstrom, S. L.: Waterbudgets and rechargearea simulations for the Spring Creek and Nittany Creek Basins and parts of the Spruce Creek Basin, Centre and Huntingdon Counties, Pennsylvania, Water Years 2000–06: U.S. Geological Scientific Investigations Report 2015–5073, 86 pp., https://doi.org/10.3133/sir20155073, 2015.
Hamed, K. H. and Rao, A. R.: A modified Mann–Kendall trend test for autocorrelated data, J. Hydrol., 204, 182–196, 1998.
Harvey, C. F., Ashfaque, K. N., Yu, W., Badruzzaman, A. B. M., Ali, M. A., Oates, P. M., Michael, H. A., Neumann, R. B., Beckie, R., Islam, S., and Ahmed, M. F.: Groundwater dynamics and arsenic contamination in Bangladesh, Chem. Geol., 228, 112–136, 2006.
Hasan, M. A., Islam, A. K. M. S., and Bokhtiar, S. M.: Changes of reference evapotranspiration ETo in recent decades over Bangladesh, 2nd International Conference on Advances in Civil Engineering, 26–28 December 2014 CUET, Chittagong, Bangladesh, 2014.
Healy, R. W., Winter, T. C., LaBaugh, J. W., and Franke, O. L.: Water Budgets: Foundations for Effective Water Resources and Environmental Management, United States Geological Survey, Reston, Virginia, 2007.
Huq, S. M. I. and Shoaib, J. U.: The Soils of Bangladesh, 1–172, https://doi.org/10.1007/9789400711280, 2013.
Hyndman, R., Mitchell O'Hara, Wild, M., Bergmeir, C., Razbash, S., and Wang, E. R.: Language Forecast Package Development team, 131 pp., 2017.
IPCC (Intergovernmental Panel on Climate Change): Technical summary of climate change 2007: the physical science basis, Contribution of working group I to the fourth assessment report of the intergovernmental panel on climate change, edited by: Solomon, S., Qin, D., Manning, M., Marquis, M., Averyt, K., Tignor, M. M. B., Miller Jr., H. L., and Chen, Z., Cambridge University Press, Cambridge, 2007.
Kamruzzaman, M., Rahman, A. T. M. S., Kabir, M. E., Jahan, C. S., Mazumder, Q. H., and Rahman, M. S.: Spatiotemporal Analysis of Climatic Variables in the Western Part of Bangladesh, Environ. Dev. Sustain., 18, 89–108, https://doi.org/10.1007/s106680169872x, 2016a.
Kamruzzaman, M., Kabir, M. E., Rahman, A. T. M. S., Mazumder, Q. H., Rahman, M. S., and Jahan, C. S.: Modeling of Agricultural Drought Risk Pattern using Markov Chain and GIS in the Western Part of Bangladesh, Environ. Dev. Sustain., 18, 569–588, doi 10.1007/s1066801698980, 2016b.
Kang, S., Gu, B., Du, T., and Zhang, J.: Crop coefficient and ratio of transpiration to evapotranspiration of winter wheat and maize in a semihumid region, Agr. Water Manage., 59, 239–254, 2003.
Kanoua, W. and Merkel, B. J.: Groundwater recharge in Titas Upazila in Bangladesh, Arab. J. Geosci., 8, 1361–1371, https://doi.org/10.1007/s1251701413052, 2015.
Karim, M. R., Ishikawa, M., and Ikeda, M.: Modeling of seasonal water balance for crop production in Bangladesh with implications for future projection, Ital. J. Agron., 7, e21, https://doi.org/10.4081/ija.2012.e21 2012.
Karimi, P., Bastiaanssen, W. G. M., and Molden, D.: Water Accounting Plus (WA+) – a water accounting procedure for complex river basins based on satellite measurements, Hydrol. Earth Syst. Sci., 17, 2459–2472, https://doi.org/10.5194/hess1724592013, 2013.
Kendall, M. G.: Rank Correlation Methods, Griffin, London, 1975.
Khalek, M. A. and Ali, M. A.: Comparative Study of WaveletSARIMA and WaveletNNAR Models for Groundwater Level in Rajshahi District, IOSR Journal of Environmental Science, Toxicology and Food Technology (IOSRJESTFT), 10, 1–15, 2016.
Kisi, O.: Stream flow forecasting using neurowavelet technique, Hydrol. Process., 22, 4142–4152, 2008.
Kumar, M., Denis, D. M., and Suryavanshi, S.: Longterm climatic trend analysis of Giridih district, Jharkhand (India) using statistical approach, Modeling Earth Systems and Environment, 2, 116, https://doi.org/10.1007/s4080801601622, 2016.
Leta, O. T., ElKadi, A. I., Dulai, H., and Ghazal, K. A.: Assessment of climate change impacts on water balance components of Heeia watershed in Hawaii, J. Hydrol.Regional Studies 8, 182–197. https://doi.org/10.1016/j.ejrh.2016.09.006, 2016.
Lutz, S. R., Mallucci, S., Diamantini, E., Majone, B., Bellin, A., and Merz, R.: Hydroclimatic and water quality trends across three Mediterranean river basins, Sci. Total Environ., 571, 1392–1406, 2016.
Mallat, S. G.: A theory for multiresolution signal decomposition: the wavelet representation, IEEE T. Pattern Anal., 11, 674–693, 1989.
Mann, H. B.: Nonparametric tests against trend, Econometrica, 13, 245–259, 1945.
McCuen, R. H., Knight, Z., and Cutter, A. G.: Evaluation of the Nash–Sutcliffe Efficiency Index, J. Hydrol. Eng., 11, 597–602, 2006.
McVicar, T. R., Roderick, M. L., Donohue, R. J., Li, L. T., Van Niel, T. G., Thomas, A., Grieser, J., Jhajharia, D., Himri, Y., Mahowald, N. M., Mescherskaya, A. V., Kruger, A. C., Rehman, S., and Dinpashoh, Y.: Global review and synthesis of trends in observed terrestrial nearsurface wind speeds: Implications for evaporation, J. Hydrol., 416–417, 182–205, https://doi.org/10.1016/j.jhydrol.2011.10.024, 2012.
Michael, H. A. and Voss, C. I.: Controls on groundwater flow in the Bengal Basin of India and Bangladesh: regional modeling analysis, Hydrogeol. J., 17, 1561–577, 2009.
Misiti, M., Misiti, Y., Oppenheim, G., and Poggi, J.: Wavelet Toolbox User's Guide, The Math Works, Inc., 1997.
Molden, D. and Sakthivadivel, R.: Water accounting to assess use and productivity of water, Int. J. Water Resour. D., 15, 55–71, https://doi.org/10.1080/07900629948934, 1999.
Moriarty, P., Batchelor, C., AbdAlhadi, F., Laban, P., and Fahmy, H.: The Empowers Approach to Water Governance Guidelines, Methods and Tools. Jordan: InterIslamic Network on Water Resources Development and Management (INWRDAM), 2007.
Moussa, R.: When monstrosity can be beautiful while normality can be ugly: assessing the performance of event–based flood models, J. Hydrolog. Sci., 55, 1074–1084, 2010.
Nalley, D., Adamowski, J., and Khalil, B.: Using discrete wavelet transforms to analyze trends in streamflow and precipitation in Quebec and Ontario (1954–2008), J. Hydrol., 475, 204–228, 2012.
Nalley, D., Adamowski, J., Khalil, B., and OzgaZielinski, B.: Trend detection in surface air temperature in Ontario and Quebec, Canada during 1967–2006 using the discrete wavelet transform, Atmos. Res., 132–133, 375–398, 2013.
Nash, J. E. and Sutcliffe, J. V.: River flow forecasting through conceptual models, part I: a discussion of principles, J. Hydrol., 10, 282–290. https://doi.org/10.1016/00221694(70)902556, 1970.
Nasher, N. M. R. and Uddin, M. N.: Maximum and Minimum Temperature Trends Variation over Northern and Southern Part of Bangladesh, Journal of Environmental Science and Natural Resources, 6, 83–88, 2013.
Nievergelt, Y.: Wavelets Made Easy, Birkhäuser, 297 pp., 2001.
Nury, A. H., Hasan, K., and Alam, J. B.: Comparative study of waveletARIMA and waveletANN models for temperature time series data in northeastern Bangladesh, Journal of King Saud University – Science, 29, 47–61, 2017.
Nury, A. H., Hasan, K., Erfan, K. M., and Dey, D. C.: Analysis of Spatially and Temporally Varying Precipitation in Bangladesh, Asian Journal of Water, Environment and Pollution, 13, 15–27, https://doi.org/10.3233/AJW160023, 2016.
Partal, T.: Modelling evapotranspiration using discrete wavelet transform and neural networks, Hydrol. Process., 23, 3545–3555, 2009.
Partal, T.: Wavelet transformbased analysis of periodicities and trends of Sakarya basin (Turkey) streamflow data, River Res. Appl., 26, 695–711, 2010.
Partal, T. and Küçük, M.: Longterm trend analysis using discrete wavelet components of annual precipitations measurements in Marmara region (Turkey), Phys. Chem. Earth, 31, 1189–1200, 2006.
Pathak, P., Kalra, A., and Ahmed, S.: WaveletAided Analysis to Estimate Seasonal Variability and Dominant Periodicities in Temperature, Precipitation, and Streamflow in the Midwestern United States, Water Resour. Manag., 30, 4649–4665, https://doi.org/10.1007/s1126901614450, 2016.
Percival, D. B. and Walden, A. T.: Wavelet Methods for Time Series Analysis, 1106 Cambridge University Press, New York, 594 pp., 2000.
Pišoft, P., Kalvová, J., and Brázdil, R.: Cycles and trends in the Czech temperature series using wavelet transforms, Int. J. Climatol., 24, 1661–1670, 2004.
Rahman, A. T. M. S.: Sustainable Groundwater Management in the Context of Climate Change in Drought Prone Barind Area, NW Bangladesh, Unpublished Master of Philosophy Thesis, Institute of Environmental Science, University of Rajshahi, Bangladesh, 2016.
Rahman, A. T. M. S., Jahan, C. S., Mazumder, Q. H., Kamruzzaman, M., and Hossain, A. Evaluation of spatiotemporal dynamics of water table in NW Bangladesh: An integrated approach of GIS and Statistics, Sustainable Water Resource Management, 2, 297–312, https://doi.org/10.1007/s4089901600574, 2016.
Rahman M. A., Yunsheng, L., and Sultana, N.: Analysis and prediction of rainfall trends over Bangladesh using Mann–Kendall, Spearman's rho tests and ARIMA model, Meteorol. Atmos. Phys., 129, 409–424, https://doi.org/10.1007/s0070301604794, 2016.
Rahman, M. J. and Hasan, M. A. M.: Performance of Wavelet Transform on Models in Forecasting Climatic Variables, in: Computational Intelligence Techniques in Earth and Environmental Sciences, edited by: Islam, T., Srivastava, P., Gupta, M., Zhu, X., and Mukherjee, S., Springer, Dordrecht, 2014.
Rahman, M. R. and Lateh, H.: Spatiotemporal analysis of warming in Bangladesh using recent observed temperature data and GIS, Clim. Dynam., 46, 2943–2960, 2016.
Ramana, R. V., Krishna, B., Kumar, S. R., and Pandey, N. G.: Monthly rainfall prediction using wavelet neural network analysis, Water Resour. Manag., 27, 3697–3711, 2013.
Ravenscroft, P., Burgess, W. G., Ahmed, K. M., Burren, M., and Perrin, J.: Arsenic in groundwater of the Bengal Basin, Bangladesh: Distribution, field relations, and hydrogeological setting, Hydrogeol. J., 13, 727–751, 2005.
R Core Team: A language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, Austria, 2016.
Ritter, A. and MuñozCarpena, R.: Performance evaluation of hydrological models: Statistical significance for reducing subjectivity in goodnessoffit assessments, J. Hydrol., 480, 33–45, https://doi.org/10.1016/j.jhydrol.2012.12.004, 2013.
Sang, Y.F.: A practical guide to discrete wavelet decomposition of hydrologic time series, Water Resour. Manag., 26, 3345–3365, 2012.
Sang, Y. F.: A review on the applications of wavelet transform in hydrology time series analysis, Atmos. Res., 122, 8–15, https://doi.org/10.1016/j.atmosres.2012.11.003, 2013.
Sang, Y.F., Wang, Z., and Liu, C.: Discrete waveletbased trend identification in hydrologic time series, Hydrol. Process., 27, 2021–2031, 2013.
Santos, C. A. G and da Silva, G. B. L.: Daily streamflow forecasting using a wavelet transform and artificial neural network hybrid models, Hydrolog. Sci. J., 59, 312–324, 2014.
Scott, C. A. and Sharma, B.: Energy supply and the expansion of groundwater irrigation in the IndusGanges Basin, International Journal of River Basin Management, 7, 1–6, 2009.
Sen, P. K.: Estimates of the regression coefficient based on Kendall's tau, J. Am. Stat. Assoc., 63, 1379–1389, 1968.
Shahid, S. and Behrawan, H.: Drought risk assessment in the western part of Bangladesh, Nat. Hazards, 46, 391–413, 2008.
Shahid, S. and Khairulmaini, O. S.: Spatial and temporal variability of rainfall in Bangladesh, AsiaPac. J. Atmos. Sci., 45, 375–389, 2009.
Shahid, S.: Recent trends in the climate of Bangladesh, Clim. Res., 42, 185–193, 2010.
Shamsudduha, M., Chandler, R. E., Taylor, R. G., and Ahmed, K. M.: Recent trends in groundwater levels in a highly seasonal hydrological system: the GangesBrahmaputraMeghna Delta, Hydrol. Earth Syst. Sci., 13, 2373–2385, https://doi.org/10.5194/hess1323732009, 2009.
Smith, L. C., Turcotte, D. L., and Isacks, B. L.: Stream flow characterization and feature detection using a discrete wavelet transform, Hydrol. Process., 12, 233–249, 1998.
Sneyers, R.: On the Statistical Analysis of Series of Observations, Secretariat of the World Meteorological Organization, 192 pp., 1990.
SRDI (Soil Resources Development Institute): Soil map of Bangladesh, Soil Resources Development Institute, 1998.
Sreekanth, P., Geethanjali, D. N., Sreedevi, P. D., Ahmed, S., Kumar, N. R., and Jayanthi, P. D. K.: Forecasting groundwater level using artificial neural networks, Curr. Sci., 96, 933–939, 2009.
Stein, C. M.: Estimation of the Mean of a Multivariate NormalDistribution, Ann. Stat., 9, 1317–1322, 1981.
Syed, M. A. and Al Amin, M.: Geospatial Modeling for Investigating Spatial Pattern and Change Trend of Temperature and Rainfall, Climate, 4, 21, https://doi.org/10.3390/cli4020021, 2016.
Tabari, H. and Talaee, P. H.: Moisture index for Iran: spatial and temporal analyses, Global Planet. Change, 100, 11–19, 2013.
Thornthwaite, C. W. and Mather, J. R.: The Water Balance, Publications in Climatology, VIII(1), 1–104, Drexel Institute of Climatology, Centerton, New Jersey, 1955.
Thornthwaite, C. W. and Mather, J. R.: Instructions and tables for computing potential evapotranspiration and the water balance, Publications in Climatology, Laboratory of Climatology, Drexel Institute of Technology, Centerton, New Jersey, USA, 10, 183–311, 1957.
Tiwari, M. K. and Chatterjee, C.: Development of an accurate and reliable hourly flood forecasting model using wavelet–bootstrap–ANN (WBANN) hybrid approach, J. Hydrol., 1, 458–470, 2010.
Torrence, C. and Compo, G. P.: A practical guide to wavelet analysis, B. Am. Meteorol. Soc., 79, 61–78, 1998.
Valipour, M.: Ability of BoxJenkins Models to Estimate of Reference Potential Evapotranspiration (A Case Study: Mehrabad Synoptic Station, Tehran, Iran), IOSR Journal of Agriculture and Veterinary Science, 1, 1–11, 2012.
Valipour, M., Banihabib, M. B., and Behbahani, M. R.: Comparison of the ARMA, ARIMA, and the autoregressive artificial neural network models in forecasting the monthly inflow of Dez dam reservoir, J. Hydrol., 476, 433–441, 2013.
Viaroli, S., Mastrorillo, L., Lotti, F., Paolucci, V., and Mazza, R.: The groundwater budget: a tool for preliminary estimation of the hydraulic connection between neighboring aquifers, J. Hydrol., 556, 72–86, doi.org/10.1016/j.jhydrol.2017.10.066, 2017.
Wang, W., Ding, J., and Li, Y.: Hydrologic Wavelet Analysis, Chemical Industry Press, Beijing, China, 2005 (in Chinese).
Wang, D., Singh, V. P., Shang, X., Ding, H., Wu, J., Wang, L., Zou, X. Chen, Y., Chen, X., Wang, S., and Wang, Z.: Sample entropy based adaptive wavelet denoising approach for meteorological and hydrologic time series, J. Geophys. Res.Atmos., 119, 8726–8740, https://doi.org/10.1002/2014JD021869, 2014.
Wolock, D. M. and McCabe, G. J.: Effects of potential climatic change on annual runoff in the conterminous United States, J. Am. Water Resour. As., 35, 1341–1350, 1999.
Xu, C.Y. and Halldin, S.: The effect of climate change in river flow and snow cover in the NOPEX area simulated by a simple water balance model, Proc. of Nordic Hydrological Conference, Alkureyri, Iceland, 1, 436–445, 1996.
Yoon, H., Jun, S. C., Hyun, Y., Bae, G., and Lee, K. K.: A comparative study of artificial neural networks and support vector machines for predicting groundwater levels in a coastal aquifer, J. Hydrol., 396, 128–138, 2011.
Young, P. C.: Nonstationary time series analysis and forecasting, Progress in Environmental Science, 1, 3–48, 1999.
Yue, S., Pilon, P., Phinney, B., and Cavadias, G.: The influence of autocorrelation on the ability to detect trend in hydrological series, Hydrol. Process., 16, 1807–1829, 2002.