1. Introduction
Excluding glacial water, groundwater is the largest fresh water resource in the world. Due to characteristics such as well water quality, strong regulation ability and steady storage, groundwater plays a vital role in agriculture, industry, public supply and ecosystems in many countries [
1]. The demand for groundwater has been intensely increased with economic development and population growth. Overexploitation of groundwater resources will distort the natural recharge–discharge equilibrium of the groundwater system and thereby result in the decline of groundwater levels, leading to freshwater scarcity, saltwater intrusion and land subsidence in different parts of the world [
2,
3]. For example, there are many areas throughout the world where aquifers face depletion, such as California’s Central Valley in the United States [
4], the North China Plains Aquifers [
5], and the Neogene and Dammam Aquifers in eastern Saudi Arabia [
6]. Many studies have been undertaken, from both technical and managerial perspectives, to address groundwater-related issues [
7,
8,
9,
10]. The need for improved groundwater management is becoming increasingly obvious, as it is becoming clear that technical changes alone cannot solve the problems [
11]. Effective groundwater resource management is a favorable guarantee for the rational allocation of groundwater resources and the coordinated development of resources, the environment, the economy and society [
3,
12].
Effective groundwater management needs the support of national laws [
13]. Many countries, such as USA [
12], Australia [
14] and Europe [
15,
16], have developed a number of laws and regulations for the management of groundwater resources [
17]. The healthy development of groundwater resources has benefited from strong legal support and good administrative decisions about groundwater management in these countries. In China, with more and more attention being paid to the protection of groundwater resources, the Chinese government has issued a series of water resource management systems and policies, wherein a strict water resource management system is a new policy provided by the Ministry of Water Resources of China in 2012, and has become a fundamental measure to solve the problem of water shortage in China. This policy clearly notes that the management and the dynamic monitoring of groundwater must be strengthened, both the total groundwater quantity (GWQ) and groundwater level (GWL) must be controlled, and the overexploitation of groundwater must be cut down, especially in the overexploitation region, to realize the balance of groundwater recharge and recovery. The water-related government emphasizes the establishment of water resource management responsibilities and investigative systems, which can offer an important basis for comprehensive assessment by the local government managers. Previously in China, groundwater resource management was mainly based on groundwater use. However, the quantity of groundwater use is difficult to supervise and control under actual operation. GWL is the comprehensive reflection of the hydrogeological condition; not only can it reflect the variation law of groundwater recharge, runoff and discharge but it is also an important indicator of whether there exists a potential environmental problem, as well as its severity.
From the perspective of sustainable utilization of groundwater resources, the determination of a reasonable GWL is a fundamental basis for the evaluation of groundwater resources and control of the total amount of groundwater intake. The observation of the GWL is an important and direct means of understanding the status of groundwater. To relieve the stress of groundwater overdraft, the Chinese government has invested over one billion Yuan in recent years to implement the national groundwater-monitoring project, which is jointly carried out by the Ministry of Land and Resources and the Ministry of Water Resources, China. At the local level in China, groundwater management legislation has been issued. Some provinces have carried out practical studies on the management of the groundwater level. In Shanxi Province, the GWL monitoring network has been developed, and the fluctuation of the GWL has been incorporated into the appraisal index system for the government of local districts, including 11 cities and 119 counties since 2008. A percentage grading system is used to assess the effectiveness of groundwater management, which classifies grades into good, qualified and unqualified, according to the degree of controlled GWL. The government whose assessment is unqualified must make an improvement action, and the relevant leaders must shoulder the responsibility. Eight second-level indicators are used to evaluate deep groundwater management efficiency from three aspects, which are law enforcement capability, management ability and management effectiveness in Tianjin City [
18]. According to the targets of groundwater management, the variation index criteria of the GWL are put forward as the basis for assessment in each district and county. The city of Shanghai, Hebei and Shandong Provinces have also taken active measures on the prohibited value of the total GWQ and completed an assessment according to the change in groundwater in each district. Overall, groundwater management legal support is still relatively weak in China, so it is urgent to introduce relevant management regulations, especially on groundwater “double control” (GWQ and GWL), as soon as possible to strengthen the management of GWL [
11].
Numerical simulation is an important and common method to analyze groundwater dynamics [
19,
20,
21]. An accurate numerical model depends on ample data related to the groundwater system. However, it is difficult to obtain these important and sufficient data due to the limitation of human and material resources, as well as acquisition methods and techniques in a short time period. Moreover, there is a strong requirement for high-performance computers that have high division accuracy in space and time. Therefore, it is necessary to find a surrogate model that can make an accurate prediction, while mainly depending on a limited amount of data. This model not only can reduce the cost of modeling but also may be more suitable and available for managers to act in groundwater management due to its speed, convenience and friendlier operability. Therefore, data-driven models are a good choice when data are limited and the groundwater system is very complex [
22]. Data-driven models such as multiple linear regression (MLR), Back-Propagation Artificial Neural Network (BP-ANN), Support Vector Machines (SVM) and Power Function Models (PFM), are gradually used for groundwater level prediction, and thus, comparative studies on the above prediction techniques have been widely used and applied by some researchers [
23,
24,
25].
Beijing, the political, economic and cultural center of China, should be guaranteed the safety of its water supply. However, Beijing is one of the most water-scarce cities, with a population of more than 20 million. The annual available water is less than 100 m
3 per person per year, which is one-twentieth of the average in China [
26]. Groundwater plays an importance role in urban water supply, which accounts for approximately 2/3 of the total water supply. Due to a lack of water resources, groundwater in the Beijing Plain area has experienced excessive exploitation that has continued for over 30 years, leading to problems such as land subsidence and deterioration of water quality [
27]. Considering the above problems, a series of laws and regulations have been implemented in recent years, especially the implementation of the assessment system for groundwater resources in all districts and counties. These policies mainly strengthen the management and protection of groundwater resources from the two aspects: the total amount of groundwater and the water level.
The objectives of this paper aim at developing surrogate models, which can dually control the GWL and GWQ to evaluate groundwater dynamics quickly and accurately in each district of the Beijing Plain. These models use the data of groundwater use, precipitation and a series of depths to groundwater (GWDs) by multiple linear regression (MLR) and the back propagation artificial neural network (BP-ANN) method. The surrogate model uses modern science and technology to offer a more convenient way for managers to quickly estimate the amplitude of fluctuation of GWD in each administrative region, and can provide effective, timely, and scientifically based decision support for the management of water resources in the Beijing Plain.
3. Background of the Study Area
The study area is the Beijing Plain, located in the northwest region of the North China Plain (NCP) at longitude 115°25′~117°30′ E and latitude 39°28′~41°05′ N. The study area covers an area of approximately 6400 km
2. The Beijing Plain consists of 15 districts that are shown in
Figure 4. Beijing belongs to typical warm, semi-humid continental monsoonal climate zone with distinguishing seasons that have dry, cold winters and rainy, high temperature summers [
27]. According to meteorological data, districts 1, 2, 5 and 9 have the largest amount of precipitation, where annual precipitation is generally greater than 600 mm. In general, the precipitation reduces gradually from east to west and from northwest to southeast in Beijing. The annual potential evaporation is approximately 1800 mm. Observation wells are generally uniformly distributed in this area.
3.1. Hydrogeology Conditions
In the study area, the aquifer systems can be divided into porous aquifers, which are the main sources for groundwater use and the focus of this paper; karst aquifers; and fissure aquifers. Porous aquifers are mainly distributed in the sedimentary strata formed by fluvial processes. The sedimentary characteristics of alluvial and flood fans are basically uniformed in the whole plain. The transmissivity is relatively higher at the top of the alluvial fan, and the transmissivity can reach more than 10,000 m2/day in districts 1, 2, 5, 12 and 13. The transmissivity is lower at the center of the plain, which is less than 500 m2/day southeast of districts 4 and 8.
Precipitation is the main source of groundwater recharge and accounts for approximately 50% of the total recharge. The average annual precipitation was approximately 585 mm in the Beijing Plain from 1956 to 2000 [
26]. However, Beijing suffered the longest drought in its history during 1999–2011 [
39].
Figure 5 shows the average annual precipitation and groundwater use in the Beijing Plain during 2000–2014 according to the observations offered by the central hydrological station of Beijing. As the map shows, the average annual rainfall was greater than 585 mm only in 2008, 2011 and 2012. Persistent drought and increasing area of hardening road have caused and intensified the decrease in precipitation infiltration. The discharge consists of natural discharge and artificial extraction, which is the main discharge outflow. To meet the water supply demand in Beijing, several emergency water sources were developed in the plain [
40]. During the continuous drought period since 1999, the demand for groundwater in the plain area has been intensely increased compared to before. Because of the continuous excessive groundwater use, the GWD is increasing, which makes the phreatic water evaporation and groundwater overflow decrease gradually.
The variations in the GWL and flow field in the Beijing Plain have a close relation with the changes in groundwater recharge and water use. Due to the reduced groundwater recharge caused by continuous drought and an increasing water supply demand for sharply expanded population [
26], the groundwater use has been increased in each district in last ten years. The use of emergency water sources has been routinized in districts 2, 5 and 6. Thus, the GWL has declined significantly in most areas in recent years. The rate of fall is largest at the top of the alluvial fan and the GWL has been rapidly declining around the emergency water sources.
Figure 6 describes the dynamic changes in the GWL at location “a” (see
Figure 4) in district 3 during 2000–2014.
3.2. Cones of Depression
The overall trend of groundwater flow is from the northwestern alluvial fan to the top southeastern front into the plain area in the Beijing Plain. According to the official
Water Resources Bulletin of Beijing (
http://zfxxgk.beijing.gov.cn/110030/), the cone of depression was mainly located at districts 3 and 7, and the total area of all cones in the whole Beijing Plain was approximately 796 km
2, and the area, where the average GWD was usually greater than 10 m, was approximately 2292.5 km
2 in 2000. In 2014, the two cones of depression in districts 3 and 7 were connected, and at the same time, a large cone of depression was produced in district 4. The total area of all cones was increased to 1058 km
2, during the period from 2000 to 2014. And the area, where the GWD was greater than 10 m, was increased to 5470.0 km
2 at the end of 2014.
3.3. Data Sources Preparation
The available data of observation wells is a period of 15 years, from 2000 to 2014, in 13 districts of the Beijing Plain. For each district, the data from the first fourteen years (2000–2013) were used for model calibration (or training), and the data from 2014 were used for model testing.
Table 1 shows the data sources of precipitation, groundwater use and depth to groundwater. The monthly records of rainfall amount for several locations within the study area for the 15-year period (2000–2014) were obtained from 14 rainfall stations in each district of the study area.
Figure 4 presents the locations of the rainfall stations and observation wells in the Beijing Plain. Groundwater use data were obtained from the hydrologic station and office for water conservation in each district of the Beijing Plain. To calculate the average depth to groundwater in each district, monthly observed depth to groundwater data from 2000 to 2014 were obtained from the central hydrological station of Beijing.
4. Results and Analysis
In this study, both MLR and BP-ANN are employed to develop monthly dynamic models for quantifying the relation of depth to groundwater (GWD) and groundwater quantity (GWQ) based on GWD, groundwater use and precipitation of each district, respectively.
4.1. MLR Model
Taking the average GWD of the administrative area as the output variable, the groundwater use and rainfall of each month, as well as the GWD, groundwater use and rainfall of the previous month are used as input variables. Taking one month as a time unit of all these variables, the correlation among GWD, rainfall and groundwater use during the years from 2000 to 2013 is analyzed and the regression coefficients between the variables are determined. Then, the statistical regression equation is established. Thus, the comparison of predicted and observed GWD in each district can be obtained. The developed models can then be used to predict the regional average GWD of 2014 according to the estimation of precipitation and the plan of groundwater use. Finally, by comparing the predicted GWD and observed GWD, the reliability of the models can be tested and verified.
4.2. BP-ANN Model
Taking the average GWD of each administrative region as the output variable, the precipitation, and groundwater use in each district and its adjacent districts are regarded as the input variable. Each variable is based on the monthly time series from 2000 to 2014. Each district has 168 records of training data, and the number of input variables is determined by the number of neighboring districts. In this study, the BP-ANN model is constructed based on the C++ program. Based on many experiments and experiences, the number of hidden layers is set as 5. After training 50,000,000 times, the comparison of GWD between calculated and observed data in each district can be obtained. Then, using the data of 2014, we can predict the average GWD of each district in 2014 by the developed model. By comparing the predicted GWD with the observed data, the accuracy of the models can be tested.
4.3. Comparison of Results from Two Models
To demonstrate the applicability and effectiveness of the developed surrogate models, the comparisons of predicted and observed GWDs from developed models in the typical districts are shown in
Figure 7. Comparing the predicted GWD by MLR models with observed GWD, it is evident that the predicted results effectively captured the trends in the inter-monthly variability of observed GWD, and both the curves dynamically increased from Jan 2000 to Dec 2013. By and large, the predicted GWD has a perfect match with the observed data in most of the districts. Likewise, comparing the predicted GWD to the ANN models with the observed GWD, the predicted results captured the trends in the inter-monthly variability of the observed GWD quite effectively, and both curves dynamically increased from January 2000 to December 2013. Thus, the predicted GWD is a perfect match with the observed data in nearly all districts. Overall, both of the models have nearly the same well behavior for GWD prediction in each district of the Beijing Plain.
4.4. Evaluation of Model Efficiency
To provide information about the overall match between the observed and simulated depths to groundwater in the typical districts, the Nash–Sutcliffe efficiency (NSE) coefficient, which is calculated using Equation (3) and which is used to evaluate the model efficiency, is used.
where
E is the NSE coefficient,
n is the total number of data,
i is the serial number,
O indicates observed data,
is the average of all observed data, and
P indicates simulated values. The closer the
E value is to the maximum value 1, the higher the simulation accuracy.
Table 2 summarizes the results of the model efficiency evaluation. The NSE coefficients for both the MLR and ANN models are almost over 0.96. In addition, the correlation coefficient, usually denoted as
R, provides a measure of how well-observed outcomes are replicated by the model. The performance of most MLR and ANN models yielded an R value of over 0.97. These coefficients demonstrate that the simulations of monthly GWD using two methods match the curve of the observed data very well. Although the model efficiency is not very satisfactory in district 4, the model efficiency of these two models on the monthly GWD meets the requirement for groundwater management in most districts.
To further validate the efficiency of the results predicted by the developed models, the comparison between the observed and predicted GWDs in each district of the Beijing Plain is listed in
Table 3. The absolute errors between the predicted and observed annual average GWD are generally less than 0.23 m, and the relative errors are approximately less than 0.23% for MLR models. The absolute errors are basically less than 1.38 m, and the relative errors are approximately less than 0.05% for ANN models. The simulated GWD and the measured GWD are lower than 4.65% in all districts. The better fitting and validation indicate that both the developed models are feasible and reliable.
Comparing BP-ANN models with MLR models, the results are consistent with the observed results. The fitting from the BP-ANN model is better than that from the MLR, while the predicted GWD of 2014 is worse, especially in districts 2 and 3. Generally, the better the fitting is, the more accurate the predicted GWD. From
Table 3, the more coefficients there are, the less the relative deviations predicted by the MLR are in most districts, except districts 4 and 7. However, the prediction of the BP-ANN does not observe this. For example, the observed GWD and predicted GWD do not match very well in districts 3, 10 and 11, where NSEs are 0.96, 0.99, 0.99, respectively, and R values are 0.98, 0.99, 0.99, respectively, while its relative deviations of prediction are nearly 0.90 m, 1.12 m and 0.96 m, respectively. Therefore, the MLR models are more appropriate and available as surrogate models for predictions of the Beijing Plain. Due to the good fitting of the prediction using the MLR models with the observed GWD in the test phase and the low deviation between the predicted and measured GWDs, the developed models can be used as surrogate models of regular and complex numerical simulation models for GWD prediction of the Beijing Plain. Using the weights of the developed models after training, the average GWD of each district can be predicted by precipitation, groundwater use and the historical GWD of each district, which is more convenient for groundwater management. Moreover, we can compare the predicted annual average GWD by the reliable surrogate model, namely, the MLR model, with the measured annual average GWD of 2014 in each district of Beijing by the rainfall, groundwater use and GWD before 2014. Thus, the grade of groundwater of each district can be evaluated, and the corresponding management strategies can be made.
5. Model Application and Discussion
To test the applicability of the developed surrogate models in groundwater management, different scenarios of groundwater use are set to discuss groundwater management. Five groundwater use scenarios of groundwater consist of 80%, 90%, 100%, 110%, and 120% compared to that of 2014. Three scenarios of precipitation are set by the dry, normal and wet seasons. Precipitation of the three scenarios and the groundwater use in 2014 are shown in
Figure 8, and based on this figure, we set fifteen scenarios of groundwater management for five kinds of groundwater use and three precipitation conditions in three typical districts. Taking 2014 as the evaluation year, we use the developed MLR models to predict the yearly average GWD of the next year to analyze the changes in the groundwater level under different scenarios. To achieve the groundwater assessment for management, combining the above method as shown in
Figure 3, we also show the monthly precipitation of wet, normal and dry seasons, respectively in
Figure 8. Thus, the monthly regional GWD can be predicted one by one under different water periods. The predicted results in three typical districts are shown in
Figure 9. From
Figure 9, we can see that the predicted monthly GWDs in three typical districts have the same patterns. Usually, the low peak of GWD will emerge in June and December, while there is a high peak of GWD in September and October. Variations in groundwater level changes are less than approximately 0.7 m, 1.0 m and 2.0 m for three precipitation scenarios, respectively, which is conformed with the observed data in this year. Therefore, these predictions can offer a reference for groundwater management.
Taking district 2 as an example, the amplitude of GWD ranges from 0.19~0 m, −0.13~−0.5 m and −0.02~−1.0 m, under wet, normal and dry seasons, respectively. Then, at the end of 2015, the measured amplitude in district 2 can be obtained by subtracting the GWD at the end of 2014 from the GWD at the end of 2015. If the measured amplitude is equal to 0.11 m and the actual precipitation is a wet season in 2015 by comparing the predicted amplitude with the measured amplitude, the grading level of groundwater management in 2015 is good in district 2. However, if the measured amplitude is equal to −0.1 m, the grading level of groundwater management is bad. Based on the management grade, corresponding management strategies would be determined, and the reasonable water level and mining schemes can be strengthened in the next year. It is worthwhile to note that the developed surrogate models need to be calibrated anew using actual related data at the end of each year for better prediction and management.
6. Conclusions
China faces a considerable challenge in groundwater management due to the shortage of groundwater resources, especially in North China. Efficient management of groundwater resources is conducive to solving these groundwater problems. With the deterioration of the groundwater situation, the Chinese government has strengthened groundwater regulations. Groundwater management legislation has been issued, which puts GWLs into the assessment system for related managers. The GWL is an important indicator to judge whether a potential environmental problem exists. Compared with groundwater use, GWL is more appropriate as an indicator for groundwater resource management. Therefore, it is significant and necessary to predict the GWL accurately. Compared with numerical models, surrogate models are much easier to operate and more efficient for quick groundwater resource management. The multiple linear regression and back propagation artificial neural network models have been used as the most common approaches for GWL prediction.
For the convenience of water administrative management, the Beijing Plain can be divided into 13 different sub-regions by administrative boundaries. Using two frequently used methods, MLR and BP-ANN, we developed surrogate models in each district. Groundwater use data for 180 months, including precipitation and average depth to groundwater from 2000 to 2014 of the Beijing Plain was used for training (2000–2013) and testing (2014) the models. The results of the developed models showed that the performance of the MLR models is better than that of the BP-ANN models and that the MLR models are a better choice for the surrogate models than the BP-ANN models. Fifteen scenarios under different groundwater use and precipitation conditions are set to show the applicability of the developed model in groundwater management. The monthly predictions show the changes of groundwater level under different scenarios match well with the observed results. Finally, the established surrogate models are used for groundwater management in the sub-regions, and then, the grade of groundwater management can be made.
The surrogate models can provide effective, timely, and scientifically based decision support that can be used by decision-makers for the management of water resources in Beijing, as well as offer some reference and demonstration for other Chinese cities and regions. This presented study mainly talks about how the surrogate can be used for groundwater management, but the relationship between the total groundwater quantity and groundwater level is important and ongoing work and not discussed in detail. The model can be used for prediction and management of the groundwater level in any region to be interested. For the purpose of actual requirements of local governments in China, the study only discusses the change of groundwater level in different districts. However, since this model relies heavily on the availability of related data, especially the groundwater observations, the prediction depends on the accuracy and locations of input variables. Uncertainty of the model prediction is caused by the quality of observed GWL. Effective improvement on the groundwater monitoring network used for groundwater management of each region will need to be enhanced in the future.