Next Article in Journal
A Review of Non-Contact Water Level Measurement Based on Computer Vision and Radar Technology
Next Article in Special Issue
River Ecological Restoration and Groundwater Artificial Recharge II
Previous Article in Journal
Life Cycle, PESTLE, and Multi-Criteria Decision Analyses of Novel Process for Nitrogen Recovery from Reject Water: Combining Electroconcentration and Stripping Methods
Previous Article in Special Issue
Hydro-Geochemical Characteristics of the Shallow Alluvial Aquifer and Its Potential Artificial Recharge to Sustain the Low Flow of the Garonne River
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Coupling Simulation and Prediction of Sustainable Utilization of Water Resources in an Arid Inland River Basin under Climate Change

1
China Institute of Geo-Environment Monitoring, Beijing 100081, China
2
Center for Hydrogeology and Environmental Geology, CGS, Baoding 071051, China
3
Department of Geology, Northwest University, Xi’an 710069, China
4
School of Water Resources and Environment, Hebei GEO University, Shijiazhuang 050031, China
5
Zhejiang Monitoring Center of Geology Environment, Hangzhou 310007, China
*
Author to whom correspondence should be addressed.
Water 2023, 15(18), 3232; https://doi.org/10.3390/w15183232
Submission received: 24 May 2023 / Revised: 15 August 2023 / Accepted: 15 August 2023 / Published: 11 September 2023
(This article belongs to the Special Issue River Ecological Restoration and Groundwater Artificial Recharge II)

Abstract

:
The arid endorheic basin of northwest China is characterized by rich land resources, water shortage, and a fragile ecological environment. The establishment of a credible coupling model of groundwater and surface water based on multi-source observation data is an effective means to study the change in basin water cycles and the sustainable utilization of water resources in the past and future. Based on the latest understanding of hydrogeological conditions, hydrology and water resource utilization data in the middle reaches of the Heihe River Basin (HRB), this paper constructs an up-to-date coupling model of surface water and groundwater to study the water balance change of the basin. The water resources data series under historical replay and CMIP5 climate model prediction are constructed to predict future changes in water resources. The study shows that, under the joint influence of natural conditions and human activities, the average annual recharge of groundwater in the study area from 1990 to 2020 is 17.98 × 108 m3/a, the average annual discharge is 18.62 × 108 m3/a, and the difference between recharge and discharge is −0.64 × 108 m3/a. The total groundwater storage is −19.99 × 108 m3, of which the groundwater storage from 1990 to 2001 was −17.52 × 108 m3 and from 2002 to 2020 was −2.47 × 108 m3. Abundant water from 2002 to 2020 in the basin significantly improved the loss of groundwater storage. Under the prediction of historical reappearance and the CMIP5 CNRM-CM5 model RCP4.5 and RCP8.5 pathways, the groundwater level of the Heihe River–Liyuanhe River inclined plain falls first because the HRB has just experienced a wet season and then rises according to future climate change. The groundwater level of the inclined plain east of the Heihe River and Yanchi basin decreases continuously because of the change in water cycle caused by human activities. The erosion accumulation plain is located in the groundwater discharge zone, and the water level is basically stable. Under the conditions of water resource development and utilization, the runoff of Zhengyixia hydrological station cannot meet the requirements of the “97 Water Dividing Plan” of the State Council in most years in the future, and the ecological and production water in the lower reaches of HRB cannot be effectively guaranteed. With the implementation of water-saving irrigation under the RCP4.5 and RCP8.5 scenarios, the runoff of Zhengyixia can meet the “97 Water Diversion Plan”. It is suggested to further improve the level of agricultural water savings in the middle reaches of the HRB and control the reasonable scale of cultivated land in order to reduce water consumption in the middle reaches of the HRB and implement sustainable utilization of water resources in the HRB.

1. Introduction

The arid region of northwest China is short on water resources and has a fragile ecological environment [1]. As an active natural resource, water resources play an important role in supporting the economic, social, and ecological civilization construction of arid areas. Water resources in arid areas have their own unique formation and evolution modes, mainly forming in mountainous areas, dissipating in plains, and existing in various forms such as precipitation, glacial meltwater, surface water, groundwater, etc. The groundwater and surface water in the basin show a huge amount of mutual transformation [2,3,4]. It is of great significance to scientifically understand the water cycle mechanism of the basin, make full use of the function of large-thickness aquifer regulation and storage, and realize the joint storage and utilization of groundwater and surface water for the sustainable utilization of water resources in arid areas [1,2,3,4,5].
The Heihe River Basin (HRB) is the second-largest arid inland river basin in China, and its water cycle and water resource characteristics have always been the research hotspots of domestic and foreign scholars [6,7]. Numerical simulation is an effective means to carry out relevant research, and groundwater models, including groundwater models, soil water and groundwater coupling models, and surface water and groundwater coupling models in different regions of the HRB, have been established. Based on the groundwater model, the relationship between base flow and recharge in mountainous areas in the upper reaches of the Heihe River [8], groundwater equilibrium and water cycle evolution under water resource management scenarios such as the history and future of agricultural irrigation and water conservation in the middle reaches of Zhangye Basin [9,10], and the eco-hydrological effect of Ejina Banner under environmental flow management of the Heihe River downstream were studied [11]. Based on the SWAT model and intelligent optimization, the optimal allocation of agricultural water and soil resources in the middle reaches of the HRB and Zhangye areas was studied [12]. The coupling model of soil water and groundwater was used to study the water cycle process in Zhangye area under historical and future water-saving scenarios [13], and the soil water balance model and numerical model of farmland gas wrapping zone based on empirical assumptions were used to study the soil water balance in Zhangye agricultural irrigation area [14]. The eco-hydrological model is used to evaluate the water cycle flux in the middle and lower basins of the HRB [15], and the role of groundwater in the eco-hydrological system of the middle and lower basins [16]. Based on the coupling model of surface water and groundwater, the water use contradiction and water cycle flux of agriculture and ecosystems in the middle reaches of the Heihe River Basin and Yanchi Basin were studied [17], and the relationship between surface water and groundwater transformation was studied [18]. The coupling model of eco-hydrology and socio-economy at the basin scale was used to study multiple eco-hydrological processes in the basin and deepen the understanding of the coupling relationship between water, soil, atmosphere, vegetation, and humans in the arid basin [19].
The accuracy of water recycling resources is key to the effectiveness of numerical models. In the middle reaches of the HRB, the water cycle mainly includes river runoff, water diversion, and extraction, among which river runoff can be used to calculate river seepage and groundwater discharge [3,4]. The amount of river leakage depends on the amount of water coming from the river, the length of the riverbed, the leakage capacity, etc. Through the preliminary geological survey work, empirical functions such as the power function regression equation of the leakage rate per unit length of the Heihe River channel and the power function of the incoming water have been established, and the curves of the cumulative discharge of groundwater and the discharge per unit length of the midstream basin have also been obtained [3,4,20,21,22]. Water diversion and extraction quantities can be used to calculate irrigation channel leakage and field infiltration quantities. The preliminary work summarized the relationship between irrigation channel leakage and the lithology of the unsaturated zone, canal lining, groundwater depth, and the relationship between field infiltration and stratum lithology, irrigation water quantity, and burial depth of water level.
Climate change is a global issue of universal concern in today’s society. In the past century, major changes have taken place in the global climate and environment, whose main feature is global warming [23], which has had an important impact on the global water cycle [24]. Due to the fragile ecological environment in northwest China, climate change prediction is extremely important for maintaining water security, ecological security, and food security in the region. Global climate models (GCMs) are currently an effective tool for predicting climate change [25]. GCMs can better simulate the characteristics of large-scale meteorological elements, but they are difficult to directly apply to meteorological and water resource prediction at the basin scale because of their low spatial resolution [26]. The downscaling method can transform large-scale, low-resolution global climate models into small-scale, high-resolution regional climate information [25], which has the advantages of a simple method, a small calculation amount, and no need to consider boundary conditions. It is widely used in China and abroad [27].
In the past 33 years, the HRB has experienced the most fluctuating hydrological change process in the dry period and the most abundant period since observation data was obtained. The cultivated land area has continued to expand, the water-saving level of agricultural irrigation has increased, the water consumption has changed significantly, and the water cycle has new characteristics. In order to deepen the understanding of the mechanism of the long hydrological cycle and future changes in water resources in the basin of the middle reaches of the HRB under the conditions of climate change and human activities, based on the latest 1:50,000 hydrogeological survey, National Natural Science Foundation projects, groundwater monitoring, and other multi-source data, this paper further analyzes the characteristics of the river–aquifer system in the middle reaches of the HRB. Based on the new observation data at the turning point of river leakage and groundwater discharge to the river, the law of river leakage in the middle reaches of the Heihe River is recognized, and a coupled numerical model of surface water and groundwater on a long time scale is constructed based on multi-source data. The process and trend of the water cycle in the study area are studied, and the reasons are analyzed. The statistical downscaling model is used to predict the future precipitation in both the upper and middle reaches of the HRB, and then the mountain runoff is estimated. Three future climate change scenarios are constructed, and the future changes in water resources in the study area are predicted by numerical models. The dynamic changes in groundwater levels, groundwater storage, and Zhengyixia station runoff are analyzed. The rationality of the development and utilization of water resources in the basin is judged, and the problems that should be paid attention to in the sustainable utilization of water resources are put forward so as to provide reference for the sustainable utilization of water resources in the arid inland river basins.

2. Study Area

Located in the middle section of the Qilian Mountains and Hexi Corridor, The Heihe River Basin is the second-largest basin of the most arid inland river basin in China, bounded by the Dahuangshan Mountain in Shandan County of Zhangye City in the east and the Black Mountain in Jiayuguan City in the west, bordering the Shiyang River Basin and the Shule River Basin. It covers the Qilian Mountain Watershed in the south to the terminal Juyanhai Lake in the north, with a basin area of 22.25 × 104 km2 (Figure 1).
From 1944 to 2022, the total runoff of the Heihe River Basin was about 37.87 × 108 m3/a, of which the main stream of the Heihe River was 16.86 × 108 m3/a and the Liyuanhe River was 2.48 × 108 m3/a. The water in the basin is mainly used for agricultural irrigation, production, life, and natural oases. Water resources are mainly formed in the Qilian Mountains in the Southern part of the basin and dissipated in the Southern and Northern basins. Water resources generally flow from south to north, eventually converging in the East Yanhai Sea and the West Yanhai Sea in the north. The main ecological and geological environmental problems in the Heihe River Basin are ecological degradation and land desertification [2].
The study area is located in the middle reaches of the Heihe River Basin, including Yanchi Basin and Zhangye Basin. Zhangye Basin consists of the Heihe River–Liyuanhe River inclined plain, the inclined plain east of the Heihe River, the erosion accumulation plain in the lower middle reaches of the HRB, and Yanchi Basin, with a study area of 0.93 × 104 km2 (Figure 1).

3. Materials and Methods

3.1. Data Source

Over the years, the natural resources department, the water resources department, the meteorological department, and the Chinese Academy of Sciences have carried out a lot of work in the Heihe River Basin. Among them, the representative research projects in recent years include the 2014–2018 “Hydrogeological Survey of Key Areas in the Heihe River Basin” and “1:50,000 Hydrogeological Survey of the Heihe River Basin in the Hexi Corridor” by the China Geological Survey, the 2010–2017 Major Research Plan “Integrated Study of Ecological-Hydrological Processes in the Heihe River Basin” funded by the National Natural Science Foundation and the National Groundwater Monitoring Project launched in 2015. Based on the above-mentioned work, a groundwater, water conservancy, and remote sensing monitoring network for the Heihe River Basin has been established, constituting rich groundwater, surface water, water resource development and utilization, land use, and evapotranspiration data. On the basis of various types of hydrogeological drilling and geophysical prospecting results, the relationship between the mountains and the basin, as well as the hydrogeological conditions of the basin, are deeply investigated. A multi-layer, layered pumping system a layered monitoring scientific test base, and a test site for hydrological parameters such as infiltration coefficient are built to update the key hydrological parameters of the basin [28,29]. The above resources provide abundant basic data for this study.
The monthly hydrological data from 1990 to 2022 on the runoff of the Yingluoxia hydrological station, the Zhengyixia station of the Heihe River, and the Yinggezui station of the Liyuanhe River comes from Zhangye Hydrological Station (Zhangye, China). The monthly data from 1990 to 2020 on the development and utilization of water resources, such as surface water diversion and groundwater extraction quantity in each irrigation area, are from the Zhangye Water Bureau (Zhangye, China). Groundwater level dynamics and groundwater flow fields are from the China Institute of Geo-Environment Monitoring (Beijing, China).
The water resources forecast period is 2023–2055, including three future climate change scenarios. The first one is a historical reappearance scenario, that is, natural conditions such as mountain runoff and rainfall use data from 1990 to 2022 to reproduce the dry-to-flood season process in the HRB. The other two are climate prediction scenarios, namely the CNRM-CM5 model RCP4.5 and RCP8.5 pathways scenarios of CMIP5. The development and utilization of water resources, such as water diversion and groundwater exploitation in irrigation areas, are all reappearances of current data in 2020.
The climate prediction data includes three kinds of data: meteorological station measurement data (China Meteorological Data Network, Beijing, China), global atmospheric NCEP reanalysis data (NOAA, Washington, DC, USA) and CMIP5 climate model prediction data (Lawrence Livermore National Laboratory, Livermore, CA, USA).

3.1.1. Heihe Mountain Runoff

The observation section of the runoff from the Heihe River is the Yingluoxia hydrological station (Figure 2), and the multi-year runoff dynamics of this station show the characteristics of alternating dry-flat-dry-flat-dry-abundant (Figure 2). The period from 1990 to 2022 is a dry period up to 33 years, of which 1990–2001 is a dry period and 2002–2022 is a flood period, which is very typical. The multi-year dynamic changes in the difference between the runoff of Yingluoxia and the runoff of Zhengyixia in the Heihe River roughly reflect the change in water consumption scale in the middle basin, dominated by agricultural irrigation.

3.1.2. Heihe River Leakage

In the middle basin of the Heihe River Basin, there is a huge amount of groundwater and surface water conversion, of which the leakage of the Heihe River and the discharge of groundwater through the river channel are the main manifestations. The infiltration section and discharge section of the main stream of the Heihe River in the Zhangye Basin are roughly bounded by the G312 bridge, also known as the turning point of groundwater and surface water conversion, which is located about 1 km upstream of the G312 bridge. The Yingluoxia–G312 bridge section is the seepage section of the hanging river, including the Yingluoxia–Caotanzhuang section and the Caotanzhuang–G312 bridge section, and the river infiltration and replenishment are mainly controlled by the actual water flow into the river. The G312 bridge–Zhengyixia section is a groundwater discharge section. Based on such observational data as hydrological station runoff, irrigation channel diversion, G312 section runoff, and evaporation in the Heihe River–Zhangye Basin, the water leakage law of the Heihe River is analyzed by the river water balance method in order to obtain the river water leakage amount over the years [30].
In 1990–2020, the runoff from Yingluoxia was 17.96 × 108 m3/a. Under the current conditions of water resource development and utilization, the annual average river water leakage of the Yingluoxia–G312 Heihe bridge section was 7.11 × 108 m3/a, accounting for 39.59% of the total runoff of Yingluoxia. The annual average infiltration of Yingluoxia–Caotanzhuang was 2.48 × 108 m3/a, and the annual average infiltration of Caotanzhuang–G312 bridge was 4.63 × 108 m3/a. The annual average water intake was 6.68 × 108 m3/a, accounting for 37.19% of the runoff from Yingluoxia. The discharge under the transformation section was 4.02 × 108 m3/a, accounting for 22.38% of the runoff from Yingluoxia. Compared with the natural conditions, the amount of water flowing through the river is reduced under human activity conditions, and the leakage and recharge rate of the river channel increase.

3.1.3. Groundwater Lateral Runoff Discharge to the Plain of the Southern Qilian Mountain

The amount of lateral runoff recharge to plain areas by mountain groundwater in the arid inland basin of northwest China has long been controversial, and this amount has been estimated to be small or large. Combined with the latest geo-hydrogeological survey, geophysical prospecting, drilling, pumping test, groundwater dynamic observation, water chemistry test, basin groundwater level measurement, and comprehensive research, the multi-year average lateral runoff recharge of groundwater in the plain area of Qilian Mountain is estimated to be 0.77 × 108 m3/a, the recharge of runoff in the shallow mountain zone is 0.40 × 108 m3/a, the lateral recharge of bedrock in the valley section is 0.07 × 108 m3/a, and the recharge of the fourth series in the valley section is 0.30 × 108 m3/a. The total recharge accounts for about 3.0% of the basin’s groundwater resources [31].

3.1.4. Infiltration Replenishment Coefficient

Irrigation water and rainfall infiltration recharge are important sources of groundwater recharge in arid endorheic basins. For a long time, infiltration recharge was calculated using empirical parameter methods, lacking mature monitoring methods and empirical data. Based on the artificial bromine tracing method, the infiltration and recharge coefficients of three irrigation conditions and natural precipitation conditions in the middle reaches of the Heihe River were determined. The infiltration recharge coefficients under large-flood irrigation, small-flood irrigation, and drip irrigation were 0.16, 0.07, and 0.11, respectively, and that is 0.1 under natural precipitation conditions [32].

3.1.5. Hydrogeological Structure and Parameters

Among the 314 boreholes in the middle reaches of the HRB, 98 core boreholes were selected, and the total drilling advance was 17,632.27 m, of which the deepest Quaternary borehole was 570.34 m. The hydrogeological structure of the basin was identified based on boreholes and geophysical data [28,33]. The hydrogeological parameters of the basin were updated based on the pumping test data [23,24]. At the top of the alluvial fan of the Qilian Mountain and the Northern Longshou Mountain, etc., the aquifer is a single phreatic aquifer, and the fine soil plain area is a multi-layer aquifer. The main lithology of the basin is pebbles and gravel, with a little sand, sub-sandy, and sub-clay. The groundwater storage space in the midstream basin is huge.

3.1.6. Distribution of Cultivated Land Area

There are a total of 35 irrigated areas in the study area. In the inclined plain north of the Qilian Mountain, the groundwater is buried deeply, and irrigation mainly relies on surface water; the buried depth of groundwater decreases along the north direction, and the amount of groundwater extraction increases, which constitutes the agricultural irrigation water source together with surface water. At the spring discharge area, the proportion of spring irrigation increased, and the Wujiang irrigation area was totally a spring irrigation area.

3.1.7. Climate Prediction Data

The meteorological data comes from the China Meteorological Data Network, and the precipitation data from Qilian, Yeniugou, Shandan, Zhangye, and Gaotai stations is selected. The global atmospheric NCEP reanalysis data are from the global atmospheric NCEP reanalysis jointly released by the American Center for Environmental Forecast (NCEP) and the National Center for Atmospheric Research (NCAR). The alternative predictors of NCEP reanalysis data are land surface temperature (air), surface relative humidity (rhum), surface pressure (pres), sea level pressure (slp), surface wind speed (wspd), surface zonal wind speed (uwnd), surface meridional wind speed (vwnd), as well as temperature, relative humidity, wind speed, zonal wind speed, meridional wind speed, vertical wind speed (omega), and geopotential height (hgt) at 500/850 hPa. On the one hand, these factors affect the longitude and latitude of climate zonality by affecting the meridional and zonal distribution of atmospheric circulation. On the other hand, meteorological circulation factors such as wind speed, sea level pressure, convective temperature field, potential height field and surface temperature affect the diffusion and transport processes of water and heat and then affect the change in precipitation. The CNRM-CM5 model, which has less absolute error in precipitation simulation in the arid region of northwest China, was chosen [34].

3.2. Hydrogeological Conceptual Model

The middle basin of the HRB is a fault basin, and the north and south sides of the model are bounded by the piedmont reverse fault; the Eastern boundary of the model is bounded by faults; and the Western boundary of the model is the Western boundary of the Yanchi Basin. Vertically, the top boundary is the groundwater surface, and the bottom boundary is the Quartenary basal boundary. The model boundaries are all flow boundaries. Horizontally, the model area is split using a 500 × 500 m finite difference grid. Vertically, the model is divided into 20 layers (Figure 3). The hydrogeological parameters, such as hydraulic conductivity, are given initial values first and then adjusted during model identification and validation.
The model recharge items included 8 items, namely rainfall infiltration recharge quantity (RCH), river seepage recharge quantity (RIV), irrigation canal seepage recharge quantity (CHA), canal irrigation (field) infiltration recharge quantity (CIR), well irrigation (field) recharge quantity (WIR), undercurrent lateral recharge quantity of the piedmont river valley (LRP), lateral recharge quantity of the piedmont non-river valley (LR), and spring water irrigation recharge quantity (SIR). The discharge items include three items, namely diving evaporation quantity (EVT), groundwater extraction quantity (WEL), and groundwater discharge quantity (GD).
The initial flow field of the numerical model is the flow field in January 1990. The simulation period is from 1990 to 2020, and the stress period is 1 month. The sources and sinks, water resources, and other elements of the water cycle are constant during each stress period.

3.3. Mathematical Model

The numerical model of groundwater flow in this study is a heterogeneous anisotropic transient three-dimensional flow model that is simulated using the MODFLOW-2005 simulation program issued by the United States Geological Survey [35]. The mathematical model of groundwater flow is:
x K h H x + y K h H y + z K z H z p = S s H t , x , y , z Ω , t > 0 H x , y , z , t t = 0 = H 0 x , y , z , x , y , z Ω , t = 0 K n H n x , y , z Γ 2 = q ( x , y , z ) , x , y , z Γ 2 , t > 0 K h H x 2 + K h H y 2 K z H z + w = μ H t , x , y , z S 1 K z H z = 0 , x , y , z S 2
In the equation: Ω —seepage area: H —aquifer head (m); S 1 —diving surface; S 2 —bulkhead; μ —yield of water (−); S s —water storage rate (1/m); K h —horizontal permeability coefficient; K z —vertical permeability coefficient (m/d); w —diving surface source sink term (m/d); p —groundwater pumping source sink term (1/d); H 0 —initial water level of the aquifer (m); Γ 2 —second boundary; n —the outer normal direction of the second boundary; and q —specific discharge of the second boundary (m/d).
The mathematical model of river flow in the study area is:
Q R L , t L = E ( L , Q R ) + j Q C j t δ L L j L > 0 , t > 0 Q R L , t L = 0 = Q R 0 t t > 0 Q R L , t 0 L > 0 , t > 0
In the equation: Q R ( L , t ) —the flow of the Heihe River, which varies with time and distance; Q R 0 ( t ) —inbound flow of Heihe River (Yingluoxia); E ( L , Q R ) —the exchange intensity between Heihe River and aquifer (gain from groundwater is positive, loss to groundwater is negative); L —distance along the Heihe River; Q C j —inflows of tributaries along the Heihe River (inflow is positive, outflow is negative); L j —the distance of the tributaries along the Heihe River.

3.4. Statistical Downscaling Model

Firstly, the statistical downscaling model is used to predict the future precipitation of the basin, including the upstream precipitation that determines the mountain runoff, as well as the precipitation in the study area. Secondly, mountain runoff is predicted. The essence of the statistical downscaling model is to establish a significant relationship between climate forecasting factors and regional meteorological elements. The process of establishing the future climate change prediction of the CMIP5 model through the statistical downscaling model is shown in Figure 4. The data period used to establish the model is from 1960 to 2014, including 1960–2005 for model calibration and 2006–2014 for model verification.

3.4.1. Prediction Factor Selection

The steps of selecting forecasting factors are as follows: first, all candidate factors are screened by the stepwise regression method, and the collinearity of the included factors in the regression equation is checked. Secondly, for similar factors with significant collinearity, according to their correlation and partial correlation with the predicted quantity, the factors with greater correlation are selected. Thirdly, after collinear exclusion, if there are still too many prediction factors in the regression equation, the standardized coefficient in the equation is small, i.e., the factors that have little influence on the prediction.

3.4.2. Multiple Regression Downscaling Model

The multiple regression downscale model of forecast factor and forecast quantity is established, and the sub-rate is predicted periodically and in the verification period, in which the prediction threshold of precipitation is set at 0.2 mm, i.e., no precipitation occurs in that month when the monthly precipitation is less than 0.2 mm. The degree of fitting of the monthly value and the relative error of the annual value are selected as indicators to evaluate the simulation effect.

3.4.3. Supplementary Multiple Regression Downscaling Model

The residuals of simulated and measured precipitation at periodic stations are calculated, and it is found that each residual can be well fitted by linear regression. Therefore, the deviation correction is carried out, and a supplementary multiple regression model based on residual analysis is established to improve the simulation effect of the downscale model of the HRB [34].

3.4.4. Evaluation of Regional Applicability of CMIP5 Model

The simulation process of the Global Climate Model (GCM) is uncertain, so it is necessary to evaluate its regional applicability before using the GCM to predict the future climate. The starting time of CMIP5′s CNRM-CM5 model prediction is January 2006, so the prediction factors of the period from January 2006 to December 2014 are used to simulate the precipitation and temperature of each station according to the two models established above. Before simulation, the spatial resolution of the CNRM-CM5 model is adjusted to 2.5 × 2.5 by using the distance reciprocal weight interpolation method, which is consistent with the NCEP reanalysis data. Regional applicability evaluation is divided into two paths: RCP4.5 and RCP8.5.

3.4.5. Prediction of Future Precipitation

According to the downscaling model determined by regional applicability evaluation, the future precipitation of the basin is predicted under the two paths of RCP4.5 and RCP8.5.

3.4.6. Prediction of Future Runoff

According to the regression model of mountain runoff and precipitation in the upper reaches of the basin, the mountain runoff in the study area is estimated.

4. Results and Discussion

4.1. Calibration and Verification of Numerical Model

Using a variety of data, such as groundwater level dynamics, groundwater flow field, groundwater discharge along the river, annual groundwater discharge to the river, go runoff, etc., the model is calibrated and verified by manual and automatic parameter adjustments [36]. The model calibration period is from 1990 to 2015, and the model verification period is from 2016 to 2020. Using Nash efficiency coefficients [4], Taylor diagrams [37], Box plots [38], Radar charts [39], and other statistics, the numerical model is considered effective in simulation [30].

4.2. Groundwater Budget in the Midstream Basin

From 1990 to 2020, the average annual groundwater recharge in the study area was 17.98 × 108 m3/a, the average annual discharge was 18.62 × 108 m3/a, and the storage change was −0.58 × 108 m3/a. From 1990 to 2001, the average annual groundwater recharge in the basin was 16.39 × 108 m3/a, the average annual discharge was 17.85 × 108 m3/a, and the storage change was −1.46 × 108 m3/a. From 2002 to 2020, the average annual groundwater recharge in the basin was 18.98 × 108 m3/a, the average annual discharge was 19.11 × 108 m3/a, and the storage change was −0.13 × 108 m3/a (Table 1). The total groundwater storage change during the simulation period was −19.99 × 108 m3, of which the groundwater storage change was −17.52 × 108 m3 from 1990 to 2001 and −2.47 × 108 m3 from 2002 to 2020.
The interannual fluctuation of groundwater recharge in the middle basin of the Heihe River Basin is large, while the interannual fluctuation of discharge is small (Figure 5). This reflects the regulation and storage effect of the large-thickness aquifer system in the basin, and the interannual fluctuation of discharge is more stable than that of recharge after the adjustment of the aquifer. At the same time, groundwater discharge has a self-adaptation feature. For example, when the amount of groundwater extraction in the shallow buried area on both sides of the Heihe River increases, the groundwater discharge decreases accordingly, and the sum of the two reflects the macroscopic change in groundwater discharge in the midstream basin.

4.3. Influence Area of the Heihe River

Based on the groundwater flow field in typical years, the influence area of the Heihe River was identified (Figure 6). From 1990 to 2001, due to the continuous drought in the HRB and the increase in groundwater extraction, the area affected by the recharge of the Heihe River decreased. From 2001 to 2014, the river basin turned from dry to wet, and the influence area increased. From 2014 to 2020, the HRB continued to be abundant in precipitation, and under the combined effects of mountain runoff, irrigation water diversion, and groundwater extraction, the area affected by the leakage of the Heihe River decreased slightly.

4.4. Proportions of Groundwater Discharge to the Heihe River

Four typical years were selected, and the multi-year variation in the proportions of discharge of the Heihe River from recharge of the Heihe River was analyzed by streamlining tracking (Figure 7, Table 2). Compared with the median water year in 1990, the proportion of leakage recharge of the Heihe River to the total recharge in the dry year of 2001 decreased from 36.92% to 34.22%, which was mainly affected by the decrease in mountain runoff. The ratio of discharge of the Heihe River to recharge of the Heihe River increased from 63.37% to 67.80%. The reason may be that, with the significant increase in groundwater exploitation in the irrigation area, more groundwater comes from storage nearby since the recharge area of the Heihe River is relatively far away. Compared with the dry year in 2001, the proportion of leakage recharge of the Heihe River in the 2014 wet year increased significantly, the discharge of the Heihe River increased, and the ratio of discharge in the Heihe River from recharge of the Heihe River also increased significantly. Compared with 2014, the mountain runoff in 2020 decreased, and the proportion of leakage recharge of the Heihe River to the total recharge also fell. In general, the larger the mountain runoff, the greater the proportion of leakage recharge of the Heihe River to the total recharge, and the greater the ratio of discharge of the Heihe River from recharge of the Heihe River.

4.5. Prediction of Mountain Runoff in the Watershed

According to the principle of prediction factor selection, the precipitation prediction factors of the Qilian station are air850, rhum500, and omega850, and the prediction factors of the Yeniugou station are air, rhum500, and omega850.
Based on the multiple regression downscaling model, the fitting degree of the monthly value of the precipitation rate is 0.536 to 0.837, and the verification period is 0.454 to 0.837. The relative error of the annual value of the precipitation rate is −22.17 percent, and the verification period is −43.47 percent. The overall fitting rate of the verification period is poor. For different stations, with the decrease in average annual precipitation, the fitting degree of the monthly value of the rate periodicity and the verification period becomes smaller, and the relative error of the annual value becomes larger. Generally speaking, there is a big error in the simulated precipitation values of each meteorological station during the verification period, and the simulated values are all smaller than the measured values, which cannot meet the prediction requirements.
The simulation effect of the precipitation supplementary regression model is obviously better than that of the multiple regression model, and the annual mean relative error of the simulation is obviously reduced. The annual mean relative error of the supplementary regression model rate of precipitation is −1.63~0.01%, and the verification period is −18.88~2.30%, which are significantly smaller than those of the multiple regression model.
The evaluation of regional applicability shows that under the two pathways, the simulation annual mean relative error of the supplementary regression models of precipitation is generally smaller than the multiple regression models. Under the two pathways, the simulated precipitation is higher than the measured value from January to April and from October to December and lower than the measured value from May to September. The simulation effect of the model for the annual distribution of precipitation is slightly poor, but the simulation effect for annual precipitation is better.
The mountain runoff of the Heihe River and other rivers is predicted according to the prediction of precipitation in the upper reaches of the HRB. Under the RCP4.5 scenario, the predicted annual runoff of Yingluoxia from 2023 to 2055 is 21.17 × 108 m3/a, and that of Yingluoxia from 2023 to 2055 is 21.46 × 108 m3/a under the RCP8.5 scenario (Figure 8).

4.6. Groundwater Budget from 2023 to 2055

Under the predicted historical reappearance scenario, the average annual recharge of groundwater in the study area from 2023 to 2055 is 18.93 × 108 m3/a, the average annual discharge is 19.59 × 108 m3/a, and the difference between recharge and discharge is −0.66 × 108 m3/a. Under the RCP4.5 scenario, the average annual recharge of groundwater in the study area from 2023 to 2055 is 19.67 × 108 m3/a, the average annual discharge is 20.23 × 108 m3/a, and the difference between recharge and discharge is −0.56 × 108 m3/a. Under the RCP8.5 scenario, the average annual recharge of groundwater in the study area from 2023 to 2055 is 19.77 × 108 m3/a, the average annual discharge is 20.30 × 108 m3/a, and the difference between recharge and discharge is −0.53 × 108 m3/a (Figure 9).

4.7. Prediction of Groundwater Levels in the Middle Reaches of the HRB

Under the three scenarios, the groundwater level in the inclined plain east of the Heihe River is in a state of continuous decline from 2023 to 2055, such as in the H65 observation hole, which is due to excessive surface water diversion and reduced groundwater recharge. Under the RCP4.5 and RCP8.5 scenarios, the water resources in the future basin increased, and the decline in groundwater level was significantly lower than that of the historical reappearance scenario. The groundwater level dynamics in the Heihe River–Liyuanhe River inclined plain are mainly affected by mountain runoff, and there is a certain trend of falling first and then rising in the future (H13-0). The falling process reflects the rise of the groundwater level during the strongest flood period from 2002 to 2022, and the rising process reflects the future increase in water resources in the basin. The erosion accumulation plain is located in the groundwater discharge zone of the midstream basin, and the groundwater table is basically stable, such as in the H77-1 observation hole. The groundwater level in the Yanchi Basin continued to decline, most typically in the Southern mountainous area, mainly affected by overexploitation in the irrigation area. For example, the continuous decline of the groundwater level in the H11-0 observation hole is due to the continuous exploitation of the Luotuocheng County irrigation area. Over time, the downward trend of H11-0 water level slowed down, which was due to the increase in lateral runoff recharge of groundwater caused by the decline of the groundwater table and the increase in the hydraulic gradient around the irrigation area. The downward trends of the groundwater levels in the Yanchi basin under the RCP4.5 and RCP8.5 scenarios are slower than those in the historical reappearance scenario (Figure 10 and Figure 11).

4.8. Groundwater Budget in the Midstream Basin

In the future, the strong groundwater runoff zone in the middle reaches of the HRB will change with the amount of water resources and human activity in the basin. In the flood year, the groundwater runoff is large, and the recharge characteristics of Heihe River leakage to groundwater are significant; in the dry year, the groundwater runoff is small, and the groundwater recharge caused by Heihe River leakage is relatively weak (Figure 12).

4.9. Future Runoff of Zhengyixia

The “97 Water Diversion Plan” of the State Council stipulates the functional relationship between the runoff of Zhengyixia and Yingluoxia, so as to ensure a certain amount of runoff from Zhengyixia to the lower reaches of the HRB and safeguard the downstream ecology. Under the historical reappearance scenario, the annual average difference between the runoff of Zhengyixia and the dividing runoff amount in the future is predicted to be −1.35 × 108 m3/a. Under the RCP4.5 and RCP8.5 scenarios, the differences between the two are −3.74 × 108 m3/a and −3.48 × 108 m3/a, respectively (Figure 13). In the future, the basin will be in a wet period and will not be able to meet the “97 Water Diversion Plan” for many more years. In the previous study, with the implementation of the water-saving irrigation under the RCP4.5 and RCP8.5 scenarios, the runoff of Zhengyixia can meet the “97 Water Diversion Plan” from 2023 to 2033 continuously [4]. It is suggested that the middle reaches of the HRB further improve the level of agricultural water saving and control the reasonable scale of cultivated land to reduce the water consumption of agricultural irrigation and ensure ecological water use in the lower reaches, aiming to realize the sustainable utilization of water resources in the basin.

5. Conclusions

Based on the latest hydrogeological survey and other multi-source data, this paper constructs a coupled numerical model of surface water and groundwater in the middle reaches of the HRB and analyzes the changing process of the water cycle. The future precipitation in the upper and middle reaches of the HRB is predicted by the statistical downscaling model, and the mountain runoff is predicted. Three future climate change scenarios are constructed, and the numerical model is used to predict the future change in water resources and judge the rationality of the development and utilization of water resources in the basin.
Under the combined influence of natural conditions and human activities, the average annual recharge of groundwater in the study area from 1990 to 2020 was 17.98 × 108 m3/a, the average annual discharge was 18.62 × 108 m3/a, and the difference between recharge and discharge was −0.64 × 108 m3/a. The total groundwater storage change was −19.99 × 108 m3, of which the groundwater storage change from 1990 to 2001 was −17.52 × 108 m3 and from 2002 to 2020 was −2.47 × 108 m3.
Under the scenarios of historical reappearance and RCP4.5 and RCP8.5 pathways of the CMIP5 CNRM-CM5 model, the groundwater level in the Heihe River–Liyuanhe River inclined plain falls first and then rises; the groundwater levels in the inclined plain east of the Heihe River and the Yanchi basin continue to decline; and the water level in the erosion accumulation plain is basically stable. In the future, the wet period will significantly improve the groundwater deficit in the basin. The future runoff of Zhengyixia cannot meet the requirements of the State Council’s “97 Water Diversion Plan” in most years, and the ecological and production water usage in the lower reaches of the HRB cannot be effectively guaranteed. With the implementation of the water-saving irrigation under the RCP4.5 and RCP8.5 scenarios, the runoff of Zhengyixia can meet the “97 Water Diversion Plan”.
The water resources in the inland river basin have multiple functions, such as resources, environment, ecology, etc. The middle reaches of the HRB should maintain a reasonable scale of land development, develop water-saving agriculture, appropriately return farmland to forest and grassland, protect wetlands and natural vegetation, reduce water consumption in the basin, and achieve water diversion and exploitation in an orderly manner to promote the coordinated and healthy development of hydrology, ecology, economy, and society in the inland river basins.

Author Contributions

Conceptualization, W.L.; Investigation, H.C. and X.Q.; Writing—original draft, X.Q.; Writing—review & editing, X.Q. and W.L.; Supervision, W.L., Y.Z., W.K., Z.L. and X.S.; Project administration, W.L. and Y.Z.; Funding acquisition, W.L. and Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Projects of the National Natural Science Foundation of China (Grant Nos. 91025019; 91325107); the China Geological Survey Project (DD20160292; DD20190322; DD20230075); the China Geological Survey Administration Project (102202220180000009060); and the Xinjiang Uygur Autonomous Region Water Conservancy Department Financial Special Project (2022.B-004).

Data Availability Statement

The data presented in this study are available within the article.

Acknowledgments

The authors are very grateful to all anonymous reviewers, chief editors and associate editors for their valuable comments and suggestions, which improved the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, W.; Zhou, H.; Zhou, Y.; Jiao, P.; Hao, A. Groundwater Flow Systems of Typical Arid Areas in Northwest China; Earthquake Press: Beijing, China, 1995. (In Chinese) [Google Scholar]
  2. Li, W.; Kang, W.; Liu, Z.; Hao, A.; Li, Z.; Li, S.; Zhao, C.; Zheng, Y.; Yang, X.; Li, J.; et al. Water Resources Regulation and Optimal Utilization Models of Typical Inland Basins in Northwest China—A Case Study of Heihe River Basin; Geology Press: Beijing, China, 2010. (In Chinese) [Google Scholar]
  3. Cui, H. Water Cycle Mechanism and Co-Evolution of Water Cycle and Oasis in Heihe River Basin. Ph.D. Thesis, Northwest University, Xian, China, 2016; p. 143. (In Chinese). [Google Scholar]
  4. Qi, X. Study on Sustainable Utilization of Water Resources in an Arid Inland Basin under Conditions of Climate Change; China University of Geosciences: Beijing, China, 2018. (In Chinese) [Google Scholar]
  5. Zhou, Y.; Nonner, J.C.; Meekes, S.; Griffieon, J.; Su, Z.; Li, W.; Hao, A.; Zhou, Y.; Li, X.; Shao, X.; et al. Strategies and Techniques for Groundwater Resources Development in Northwest China; China Land Press: Beijing, China, 2007. [Google Scholar]
  6. Zhang, M.; Wang, S.; Fu, B.; Gao, G.; Shen, Q. Ecological effects and potential risks of the water diversion project in the Heihe River Basin. Sci. Total Environ. 2017, 619–620, 794–803. [Google Scholar] [CrossRef]
  7. Cai, W.; Jiang, X.; Sun, H.; Lei, Y.; Nie, T.; Li, L. Spatial scale effect of irrigation efficiency paradox based on water accounting framework in Heihe River Basin, Northwest China. Agric. Water Manag. 2023, 277, 108118. [Google Scholar] [CrossRef]
  8. Yao, Y.; Zheng, C.; Andrews, C.; Zheng, Y.; Zhang, A.; Liu, J. What controls the partitioning between baseflow and mountain block recharge in the Qinghai-Tibet Plateau? Geophys. Res. Lett. 2017, 44, 8352–8358. [Google Scholar] [CrossRef]
  9. Wen, X.; Wu, Y.; Lee, L.; Su, J.; Wu, J. Groundwater flow modeling in the Zhangye Basin, Northwestern China. Environ. Geol. 2007, 53, 77–84. [Google Scholar] [CrossRef]
  10. Chen, S.; Yang, W.; Huo, Z.; Huang, G. Groundwater simulation for efficient water resources management in Zhangye Oasis, Northwest China. Environ. Earth Sci. 2016, 75, 647. [Google Scholar] [CrossRef]
  11. Yao, Y.; Zheng, C.; Tian, Y.; Li, X.; Liu, J. Eco-hydrological effects associated with environmental flow management: A case study from the arid desert region of China. Ecohydrology 2017, 11, e1914. [Google Scholar] [CrossRef]
  12. Zhang, Z.; Wang, Q.; Guan, Q.; Xiao, X.; Mi, J.; Lv, S. Research on the optimal allocation of agricultural water and soil resourcesin the Heihe River Basin based on SWAT and intelligent optimization. Agric. Water Manag. 2023, 279, 108177. [Google Scholar] [CrossRef]
  13. Li, J.; Mao, X.; Li, M. Modeling hydrological processes in oasis of Heihe River Basin by landscape unit-based conceptual models integrated with FEFLOW and GIS. Agric. Water Manag. 2017, 179, 338–351. [Google Scholar] [CrossRef]
  14. Li, J.; Mao, X.; Shang, S.; Steenhuis, T.S. Modeling Regional Soil Water Balance in Farmland of the Middle Reaches of Heihe River Basin. Water 2017, 9, 847. [Google Scholar] [CrossRef]
  15. Yao, Y.; Zheng, C.; Liu, J.; Cao, G.; Xiao, H.; Li, H.; Li, W. Conceptual and numerical models for groundwater flow in an arid inland river basin. Hydrol. Process. 2015, 29, 1480–1492. [Google Scholar] [CrossRef]
  16. Yao, Y.; Tian, Y.; Andrews, C.; Li, X.; Zheng, Y.; Zheng, C. Role of Groundwater in the Dryland Ecohydrological System: A Case Study of the Heihe River Basin. J. Geophys. Res. Atmos. 2018, 123, 6760–6776. [Google Scholar] [CrossRef]
  17. Wu, X.; Zheng, Y.; Wu, B.; Tian, Y.; Han, F.; Zheng, C. Optimizing conjunctive use of surface water and groundwater for irrigation to address human-nature water conflicts: A surrogate modeling approach. Agric. Water Manag. 2015, 163, 380–392. [Google Scholar] [CrossRef]
  18. Hu, L.; Xu, Z.; Huang, W. Development of a river-groundwater interaction model and its application to a catchment in Northwestern China. J. Hydrol. 2016, 543, 483–500. [Google Scholar] [CrossRef]
  19. Li, X.; Zhang, L.; Zheng, Y.; Yang, D.; Wu, F.; Tian, Y.; Han, F.; Gao, B.; Li, H.; Zhang, Y.; et al. Novel hybrid coupling of ecohydrology and socioeconomy at river basin scale: A watershed system model for the Heihe River basin. Environ. Model. Softw. 2021, 141, 105058. [Google Scholar] [CrossRef]
  20. Fan, X.; Ao, S.; Jin, L.; Jiang, Z.; Hou, X.; Wan, J.; Chen, S.; Zuo, M.; Xu, Z.; Dang, Y.; et al. Study on Groundwater Distribution and Rational Development and Utilization in Hexi Corridor of Gansu Province; Geological Science Institute of Gansu Geological and Mineral Bureau: Lanzhou, China, 1984. (In Chinese) [Google Scholar]
  21. Zhou, X.; Zhao, J.; Zhang, A.; Zhao, Y.; Chen, X. Two-Water Conversion and Comprehensive Development and Utilization of Water Resources in the Middle and Lower Reaches of Heihe River in Western Gansu and Inner Mongolia; Second Hydrogeological Engineering Geological Team of Gansu Geological and Mineral Bureau: Zhangye, China, 1990. (In Chinese) [Google Scholar]
  22. Zhou, X.; Zhao, J.; Wang, Z.; Zhao, Y.; Zhang, Y.; Hua, J.; Zhang, A.; Wang, F.; Ding, H. Investigation on Groundwater Resources and Their Rational Development and Utilization in the Middle Reaches of Heihe River in Gansu Province; Second Hydrogeological Engineering Geological Team of Gansu Geological and Mineral Bureau: Zhangye, China, 1990. (In Chinese) [Google Scholar]
  23. Mora, C.; Frazier, A.G.; Longman, R.J.; Dacks, R.S.; Walton, M.M.; Tong, E.J.; Sanchez, J.J.; Kaiser, L.R.; Stender, Y.O.; Anderson, J.M.; et al. The projected timing of climate departure from recent variability. Nature 2013, 502, 183–187. [Google Scholar] [CrossRef] [PubMed]
  24. Core Writing Team; Pachauri, R.K.; Meyer, L.A. Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; IPCC: Geneva, Switzerland, 2014. [Google Scholar]
  25. Chang-Ming, L.; Wen-Bin, L.; Guo-Bin, F.; Ru-Lin, O. A discussion of some aspects of statistical downscaling climate impacts assessment. Adv. Water Sci. 2012, 23, 427–437. [Google Scholar]
  26. Wilby, R.L.; Harris, I. A framework for assessing uncertainties in climate change impacts: Low-flow scenarios for the River Thames, UK. Water Resour. Res. 2006, 42, W02419. [Google Scholar] [CrossRef]
  27. Acharya, N.; Kar, S.C.; Mohanty, U.C.; Kulkarni, M.A.; Dash, S.K. Performance of GCMs for seasonal prediction over India—A case study for 2009 monsoon. Theor. Appl. Clim. 2011, 105, 505–520. [Google Scholar] [CrossRef]
  28. An, Y.; Wang, W.; Wu, X.; Xie, W.; Gong, L.; Yin, D.; Sun, Y.; Wang, X.; Zhang, J.; Yang, J.; et al. Report on Hydrogeological Survey Results in Key Areas of Heihe River Basin; Center for Hydrogeology and Environmental Geology, CGS: Baoding, China, 2016. (In Chinese) [Google Scholar]
  29. Wu, X.; An, Y.; Yin, D.; Cui, H.; Wang, X.; Wu, B.; Wang, X.; Wang, W.; Xie, W.; Sun, Y.; et al. Report of 1:50,000 Hydrogeological Survey Secondary Project in Heihe River Basin, Hexi Corridor; Center for Hydrogeology and Environmental Geology, CGS: Baoding, China, 2020. (In Chinese) [Google Scholar]
  30. Qi, X.; Li, W.; Cui, H.; Kang, W.; Liu, Z.; Shao, X. Study on the conversion mechanism of surface water and groundwater in the middle reaches of the Heihe River Basin. Hydrogeol. Eng. Geol. 2022, 49, 29–43. (In Chinese) [Google Scholar]
  31. Li, W.; Shao, X.; Qi, X. Estimation of groundwater lateral flow in the southern mountainous area of the middle Heihe River Basin. Hydrogeol. Eng. Geol. 2022, 49, 1–10. (In Chinese) [Google Scholar]
  32. Cui, H.; Li, W.; Kang, W. A study of groundwater recharge under different irrigation conditions in the middle reaches of the Heihe River. Hydrogeol. Eng. Geol. 2022, 49, 22–28. (In Chinese) [Google Scholar]
  33. Ma, N. Construction of 3D Hydrogeological Structure Model for the Middle Reaches of Heihe River Basin; China University of Geosciences: Beijing, China, 2019. (In Chinese) [Google Scholar]
  34. Qi, X.; Li, W.; Li, H.; Liu, H. Future climate change prediction of arid inland riverbasin based on CMIP5 model. Arid. Land Geol. 2017, 40, 987–996. (In Chinese) [Google Scholar]
  35. Harbaugh, A.W. MODFLOW-2005, the U.S. Geological Survey modular ground-water model—The Ground-Water Flow Process: U.S. Geological Survey Techniques and Methods 6-A16. Available online: https://pubs.usgs.gov/publication/tm6A16 (accessed on 6 March 2015).
  36. Doherty, J. PEST: Model-Independent Parameter Estimation, User Manual, 5th ed.; Watermark Computing: Brisbane, QL, Australia, 2005. [Google Scholar]
  37. Elbeltagi, A.; Di Nunno, F.; Kushwaha, N.L.; de Marinis, G.; Granata, F. River flow rate prediction in the Des Moines watershed (Iowa, USA): A machine learning approach. Stoch. Environ. Res. Risk Assess. 2022, 36, 3835–3855. [Google Scholar] [CrossRef]
  38. Di Nunno, F.; Granata, F. Future trends of reference evapotranspiration in Sicily based on CORDEX data and Machine Learning algorithms. Agric. Water Manag. 2023, 280, 108232. [Google Scholar] [CrossRef]
  39. Jiang, L.; Wu, Y.; Wang, J.; Wang, Z.; Jiang, Q. Scenario simulation and comprehensive assessment of water footprint sustainability system in Heilongjiang Province, China: A model combining system dynamics with improved radar chart approach. Ecol. Indic. 2023, 154, 110527. [Google Scholar] [CrossRef]
Figure 1. Sketch map of the study area.
Figure 1. Sketch map of the study area.
Water 15 03232 g001aWater 15 03232 g001b
Figure 2. Annual runoff of Yingluoxia.
Figure 2. Annual runoff of Yingluoxia.
Water 15 03232 g002
Figure 3. Numerical model space discretization.
Figure 3. Numerical model space discretization.
Water 15 03232 g003
Figure 4. Prediction process of future climate change based on statistical downscaling and the CMIP5 model.
Figure 4. Prediction process of future climate change based on statistical downscaling and the CMIP5 model.
Water 15 03232 g004
Figure 5. Difference between groundwater recharge and recharge discharge in the midstream basin over the years.
Figure 5. Difference between groundwater recharge and recharge discharge in the midstream basin over the years.
Water 15 03232 g005
Figure 6. Influence area of the Heihe River in typical years.
Figure 6. Influence area of the Heihe River in typical years.
Water 15 03232 g006
Figure 7. Streamline diagram of Heihe River overflow in a typical year.
Figure 7. Streamline diagram of Heihe River overflow in a typical year.
Water 15 03232 g007
Figure 8. Predicted Yingluoxia runoff from 2023 to 2055.
Figure 8. Predicted Yingluoxia runoff from 2023 to 2055.
Water 15 03232 g008
Figure 9. The groundwater recharge, discharge, and storage changes in the study area.
Figure 9. The groundwater recharge, discharge, and storage changes in the study area.
Water 15 03232 g009aWater 15 03232 g009b
Figure 10. Future groundwater regime variations.
Figure 10. Future groundwater regime variations.
Water 15 03232 g010
Figure 11. Groundwater level variations.
Figure 11. Groundwater level variations.
Water 15 03232 g011
Figure 12. Strong groundwater runoff zones of typical years in the middle reaches of HRB.
Figure 12. Strong groundwater runoff zones of typical years in the middle reaches of HRB.
Water 15 03232 g012aWater 15 03232 g012bWater 15 03232 g012cWater 15 03232 g012d
Figure 13. Annual differences between the amount of the predicted Zhengyixia runoff and the water dividing amount in the future.
Figure 13. Annual differences between the amount of the predicted Zhengyixia runoff and the water dividing amount in the future.
Water 15 03232 g013
Table 1. Groundwater budget in the middle reaches of the HRB.
Table 1. Groundwater budget in the middle reaches of the HRB.
Recharge ItemRecharge/(108 m3·a−1) Discharge ItemDischarge/(108 m3·a−1)
199020011990–2001201420202002–2020199020011990–2001201420202002–2020
RCH1.20.981.091.021.581.21EVT2.211.922.172.072.092.02
RIV8.296.367.6611.3010.4510.08WEL0.913.571.954.365.314.45
CHA6.285.035.865.915.855.64GD15.1112.2713.7313.2312.9612.64
CIR1.191.011.041.101.091.00
WIR0.100.400.210.450.550.47
LRP0.210.180.200.260.240.23
LR0.260.220.260.320.300.29
SIR0.080.070.070.060.060.06
Total in.17.6114.2516.3920.4220.1218.98Total out.18.2317.7617.8519.6620.3619.11
Storage change −0.62−3.51−1.460.76−0.24−0.13
Table 2. Statistical values of leakage and discharge characteristics of the Heihe River.
Table 2. Statistical values of leakage and discharge characteristics of the Heihe River.
YearHeihe River Leakage Recharge /(108 m3·a−1)The Proportion of Leakage Recharge of the Heihe River to the Total Recharge (%)Discharge of the Heihe River /(108 m3·a−1)Discharge of the Heihe River from Recharge of the Heihe River/(108 m3·a−1)Ratio of Discharge of the Heihe River from Recharge of the Heihe River(%)
19906.5036.9211.157.0763.37
20014.8834.229.786.6367.80
20148.8043.1010.427.6373.24
20208.1540.5110.267.2070.13
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Qi, X.; Li, W.; Zheng, Y.; Cui, H.; Kang, W.; Liu, Z.; Shao, X. Coupling Simulation and Prediction of Sustainable Utilization of Water Resources in an Arid Inland River Basin under Climate Change. Water 2023, 15, 3232. https://doi.org/10.3390/w15183232

AMA Style

Qi X, Li W, Zheng Y, Cui H, Kang W, Liu Z, Shao X. Coupling Simulation and Prediction of Sustainable Utilization of Water Resources in an Arid Inland River Basin under Climate Change. Water. 2023; 15(18):3232. https://doi.org/10.3390/w15183232

Chicago/Turabian Style

Qi, Xiaofan, Wenpeng Li, Yuejun Zheng, Huqun Cui, Weidong Kang, Zhenying Liu, and Xinmin Shao. 2023. "Coupling Simulation and Prediction of Sustainable Utilization of Water Resources in an Arid Inland River Basin under Climate Change" Water 15, no. 18: 3232. https://doi.org/10.3390/w15183232

APA Style

Qi, X., Li, W., Zheng, Y., Cui, H., Kang, W., Liu, Z., & Shao, X. (2023). Coupling Simulation and Prediction of Sustainable Utilization of Water Resources in an Arid Inland River Basin under Climate Change. Water, 15(18), 3232. https://doi.org/10.3390/w15183232

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop