The analysis of groundwater table variations in Sylhet region, Bangladesh
Article information
Abstract
The trend analysis of the study was acquired by selecting multiyear monthly groundwater table data and monitors the wells in each sub-district under the study area. The intention of this research was to analyze the outcome of the non-parametric Mann-Kendall test at greater than the significance level which is 95% of groundwater level in Sylhet. The aptitude is effective at two conjunctures where the confidence bounds are 95% and it meets the estimate line of Sen’s. To calculate and assess the spatial differences in the inanition of groundwater table, geostatistical methods was applied based on data from 27 groundwater wells during the period from January 1975 to December 2011 which were obtained from a secondary source, Bangladesh Water Development Board. The geographic information system was used to assess the spatial change in order to find the level of groundwater. Cross-validation errors were found within an advisable level in estimating the groundwater depth with different interpolation models of ordinary kriging methods. Finally, surface maps were generated with the best-fitted model. The southeast region was found highly vulnerable from groundwater level point of view. Northern region was detected highest hazard prone area for diverge groundwater using kriging method.
1. Introduction
Groundwater is considered as the most significant source of freshwater consumption in the world. Around 33% of world’s population relies on groundwater for the daily, industrial and agricultural purpose [1]. As world’s population is increasing day by day, the extraction of groundwater is also rising significantly. Groundwater extraction is the primary cause of increasing groundwater table in many areas [2–5]. Nowadays attenuation of groundwater level is ascertained as a global problem as it causes a couple of adverse effects including water pollution, saltwater intrusion, groundwater contamination, etc. [6–7]. Groundwater depletion has been estimated about 4,500 sq. km in total for the past century where maximum rates of depletion occurred in the period of 2000–2008 [8]. Groundwater consumption for irrigation purpose has risen significantly throughout the world, especially in South Asia, irrigation related to the producing rice directly rely on groundwater, and, Bangladesh is the world’s fourth-biggest rice producing country [9–10]. Recent studies showed that increasing trend in water table depth or declining trend in groundwater level which implies an unsuitable condition of groundwater in Bangladesh. In Dhaka city, dynamics of water table has been studied by Sarkar and Ali [11] and showed a steeply increasing trend in the water table. The trend analysis was carried out by nonparametric Mann-Kendall (MK) Test and Sen’s slope estimator. The similar procedure was also used in order to study the sustainability of groundwater resources by Ali et al. [12]. The study showed a similar output; rising trend in water table depth or declining trend in groundwater level in the North-Eastern region of Bangladesh. Also, the study revealed that if the similar trend continues, water table depth will increase significantly and will be double in most cases by 2060. Parametric regression approach has also been done in North-Western Bangladesh which displayed a drop in groundwater level in Barind area [13].
There are several statistical methods for trend analysis which vary from simple linear regression to more advanced parametric and non-parametric methods [14–15]. The most popular non-parametric method for analyzing the trend in the time series is the MK test [16–17]. MK test is a popular technique for analyzing the trends in climatic parameters, evaporation, reference evapotranspiration, stream flows and groundwater fluctuations [18–21]. The result was based on the MK trends analysis using the Sen’s template application [22]. In order to determine the magnitude of change per unit time of the trends detected, the Sen’s estimator was adopted [23]. This non-parametric aligned rank test was considered adequate because no underlying frequency distribution of data could be obtained. Salmi et al. [22] ruminate, in the cases where important significant trends are observed, statistical computation gives an eminent level of sharpness with narrow angles between the confidence lines. In this paper, the aptitude was substantial at two conjunctures where the confidence bounds were 95% and it met the Sen’s line of estimation.
Moreover, Sylhet district is the capital of Sylhet division and one of the most developed regions of Bangladesh. As most of the people primarily depend on groundwater, it is under tremendous pressure. Although recent studies have proved groundwater level is rapidly depleting over the last decade, no study was conducted for this area addressing trend and variation of groundwater level fluctuation. Therefore, it has been aimed in the study to display trend and spatial analysis of groundwater level fluctuation.
2. Materials and Methods
2.1. Study Area
Sylhet district is the center of the Sylhet division and consists of twelve sub-districts that lie on the banks of Surma river in north-east of Bangladesh with an area of 3,452.07 km2. The average elevation of Sylhet district is 35 m. The climate of Sylhet is torrid monsoon with predominantly thermal and moist summer and a relatively cold winter. The annual average highest temperatures are 23°C (Aug-Oct) and the average lowest temperature is 7°C (Jan) because the settlement is in the monsoon climate zone. In between May and September, approximately 80% of the annual average precipitation (3,334 mm) occurs.
In the study area, there is limited information that is available or accessible in the aquifer systems. The main aquifer in the north-eastern region varies from semi-confined to confined types. A prospective aquifer in the north-eastern hills is the highly weathered alluvial sands of the Dupi Tila formation [24]. These sands are fine to medium grained and crop out in small hillocks in Sylhet and Moulvibazar districts and in some parts of Habiganj district. However, the permeability of these sands is lower than that of the alluvial deposits. The young gravelly sands also form a potential aquifer, although they are poorly sorted and contain large amounts of gravel and pebbles, making it difficult to use low-cost drilling techniques.
At present Groundwater is considered as the most important source of water supply in Bangladesh [25]. As a part of regular monitoring work by the Bangladesh Water Development Board (BWDB), the depths to groundwater are measured in piezometric observation wells situated in different parts of the study area. There are 27 observation wells in the study area which are considered in this study in order to determine the trend of groundwater depth. The observed variations are due to regional groundwater flow, head difference between the aquifers, and confining clay. Also, the population density of Sylhet district is 990/ km2 which indicates a tremendous pressure as it is the major freshwater source.
2.2. Data Collection
There are 27 piezometric observation wells in Sylhet district. Data were collected from secondary sources. Data has been collected from BWDB. Monthly data for a period of 1975–2011 were obtained in the study for trend analysis. It was observed that groundwater table depth is increasing over the period. In 1975, the groundwater depth was found from 0.35–5.66 m where 22% increase in water table depth was observed by the end of 1985. Finally, in 2011, the depth to groundwater varied from 15–1.52 m which clearly indicates the groundwater level is declining over the period. In this study, trend analysis was conducted by collecting secondary data and using MK test and Sen’s estimator of the slope.
2.3. Methods of Analysis
The MK [16, 26] and Sen’s slope estimator test were employed for trend analysis and the slope of the trend line [23].
2.3.1. Mann-Kendall (MK) test
This method is significantly tested by applying time series data if there is a trend exists. The preliminary value of the MK statistic, S, is assumed to be 0 (for examples, no trend). S is incremented by 1 if the later time period of a date value is higher than another previous time period of the data value. Contrarily, S is incremented by 1, if the data value from an additional time period is lower than the previous time period. The procedure to compute this probability is stated in [27–28].
Where,
Where, S is eminent and equitable, it signifies that the trend is rising and a very small value suggests that the aptitude is declining. Nevertheless, it is necessary to compute this probability to statistically quantify the significance of the trend [27–28].
2.3.2. Sen’s estimator of slope
The method of calculating the Sen’s slope estimator requires a time series of equally spaced data. Sen’s method proceeds by calculating the slope as a change in measurement per change in time, as shown here in equation [23]:
Where, Q is the slope between data points xj and xk, xj is the data measurement at time j, xk is the data measurement at time k and j is the time after time k.
2.4. Geostatistical Method
Geostatistics is a section of statistics which concentrate on spatial or spatiotemporal datasets. It is basically manifested to forecast probability distributions of ore grade mining operations. Nowadays it is applied in several systems including petroleum geology, hydrogeology, hydrology, meteorology, oceanography, geochemistry, geo-metallurgy, geography, forestry, environmental control, landscape ecology, soil science, and agriculture [29–30].
2.5. Spatial Prediction Method
2.5.1. Kriging
Kriging is a branch of geostatistical methods to interpolate the value of a random field (e.g., the elevation, z, of the landscape as a function of the geographic location) at unsighted whereabouts from the regard of its value at nearby locations [31]. This method is best suitable for normally distributed data. The data need to be resolved into normally distributed data using the transformation methods; if they are not normally distributed. The most usual transformation type is a logarithmic method because of its simplicity. The log transformation is as follows:
For Z(s) > 0 Where Z(s) is observed data, Y(s) is transformed normal data and in is the natural logarithm. Detailed discussions of kriging methods and their descriptions can be found in [30]. In kriging method, the most commonly used variogram models are spherical, exponential and Gaussian.
Where, γ(h) is semi variance, h is lag, a is the range, Co is Nugget variance, c0 + c1 = sill.
2.6. Generation of Best Fitted Models
The best-fitted models are generated by comparing the Mean Error (ME), Root Mean Square Error (RMSE), Average Standard Error (ASE) and Root Mean Square Standardized Effect (RMSSE). For the best prediction, the models must satisfy the following criteria [32].
The ME should be zero.
RMSE and ASE should be small as possible.
RMSSE should be close to 1.
For comparison of these models and data transformation, the following formula is used.
2.6.1. Root Mean Square Error (RMSE)
It is probably the most easily interpreted statistic since it has the same units as the parameter estimated. The RMSE is thus the difference, on average, of an observed data and the estimated data.
2.6.2. Mean Square Error (MSE)
For every data point, take the difference of the observed to the corresponding estimated values, and square the values. Then add up all those values for all data points, and divide by the number of points. The squaring is done so negative values do not cancel positive values. Smaller MSE indicates a better prediction of the data. The MSE has the units squared of the parameter estimated [33].
2.6.3. Root Mean Square Standardized Effect (RMSSE)
This standardized measure of effect size is used in the analysis of Variance to characterize the overall level of population effects. It is the square root of the sum of squared standardized effects divided by the number of degrees of freedom for the effect [33].
Where, σ̃2 (Xi) is the Kriging variance for location Xi. So, using ordinary kriging method each groundwater quality parameters can be generated.
3. Results and Discussion
Below Table 1 and Table 2 are provided in order to display cross-validation results and trend analysis of groundwater level in Sylhet district.
3.1. Result of Mann-Kendall Test and Sen’s Estimator of Slope
Variation in Groundwater Level Trends in some selected stations in Sylhet district. The results for groundwater level trend analysis are presented in Table 2 of the twenty-seven stations in the Sylhet District, thirteen stations had a significant trend at greater than 95% with test statistic Z ranging from −0.50 to 6.85. The trend is said to be decreasing if Z is negative and the computed probability is greater than the level of significance and increasing if Z is positive and the computed probability is greater than the level of significance. If the computed probability is less than the level of significance, there is no trend. From Table 2, it is evident that some of the stations in the Sylhet districts show no significant trend. However, Balaganj (GT9108002) show a very pronounced positive trend (increasing groundwater table depth) at Z equals to 6.85. Gowainghat and Jaintapur show relatively weaker negative trends at −3.53 and −1.65, respectively. Although there are many stations in Table 2 that showed relatively large upwards trends (increasing water table depth), only four are statistically significant. These Balaganj (GT9108001) (0.062 m/y), Balaganj (GT9108002) (0.063 m/y), Bishwhanath (GT9120006) (0.041 m/y), Golabganj (GT9138010) (0.088 m/y), Gowainghat (GT9141015) (0.0250 m/y).
Kanaighat (GT9159018) (0.023 m/y), Sylhet Sadar (GT9162021) (0.050 m/y), Sylhet Sadar (GT9162022) (0.068 m/y), (0.041 m/y), and Sylhet Sadar (GT9162026) (0.359 m/y). The magnitude of Sen’s slope estimate is shown in the bracket for each station. There are a few stations which showed a statistically significant downward trend. In Fig. 4 twenty-seven stations with statistical significant trend are presented.
Below Figs. 5–7 display the rainfall trend (MK-Z test) for three periods of time, post-monsoon, rainy and winter, respectively.
Seasonal and annual trends for the period 2001–2012 is provided below in Fig. 8 that shows Sen’s slopes estimated in seasonal and annual time scales. The median of slopes in winter is lower compared with the other seasons. As expected, rainfall trends show large variability in magnitude and direction of the trend from one station to another. All the stations have a positive value of the Sen’s estimator for the post monsoon and summer seasons. As expected, rainfall trends show large variability in magnitude and direction of the trend from one station to another. All the stations have a positive value of the Sen’s estimator for the summer season.
Figs. 9–11 explain the Sen’s slope variations in rainfall trends in pre-monsoon, rainy and winter.
4. Conclusions
By applying trend detection on homogenous zones groundwater level data (below PWD datum) in Sylhet district in Bangladesh, has resulted in identifying some significant trends. The study also represents the spatiotemporal variation of groundwater level in the study area. The direction of groundwater level trend was, in general, downward (i.e. groundwater table was upward) and statically significant across the study area. More than half of the stations that have shown highly significant trends, Sen.’s slope estimates varied between 0.359 and 0.023 mm/y, indicating groundwater level is lowering day by day. There were a few stations which reported a negative Sen’s estimator value in the whole of stations. Generally, a significant positive trend was observed to start in the mid-1975s in almost all the zones. These results demonstrate the impacts of anthropogenic developments in the study.
This research has clearly demonstrated that climate studies focusing on localized areas are possible and important since local-scale climate has direct relevance to large populations dependent on subsistence agriculture. A continuous decrease in groundwater level can be explained by the lack of strong local influences on rainfall formation like vegetation on mountains, also, increase in population that use more land for agriculture and settlement (urban sprawling). Moreover, Climate change effects at global to local scale might also be playing a role. Another concluding remark is that; continuous significant decline (at the > 95% level) was observed in the groundwater level. The scarcity of groundwater can be explained by the over-exploitation of groundwater resources for irrigation and domestic purpose, also, lacks of surface water is responsible for the growing demand of groundwater. The study suggests proper management so as to increase groundwater recharge and to find other alternative water resources in order to fulfill the water demand. However, if over-exploitation of forests will be avoided; for example, through assistance to meet communities’ basic needs within come generation opportunities, the vegetative cover can be restored. Poverty mitigation measures can include planting trees (for their products and services) in major a forestation schemes, woodlots, linear plantations, windbreaks, and agroforestry.
Acknowledgments
I wish expressing my profound thanks and sincerest appreciation to Bangladesh Water Development Board (BWDB) for providing and allowing the necessary data related to the research work.