Next Article in Journal
Hydraulic Turbine Performance Assessment with Implementation of an Innovative Aeration System
Next Article in Special Issue
Using Adaptive Capacity to Shift Absorptive Capacity: A Framework of Water Reallocation in Highly Modified Rivers
Previous Article in Journal
Evaluating the Performance of Machine Learning Approaches to Predict the Microbial Quality of Surface Waters and to Optimize the Sampling Effort
Previous Article in Special Issue
Worldwide Research on Socio-Hydrology: A Bibliometric Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Socio-Hydrological Evolution Processes Based on a Modeling Approach in the Upper Reaches of the Han River in China

1
State Key Laboratory of Eco-Hydraulics in Northwest Arid Region, School of Water Resources and Hydropower, Xi’an University of Technology, Xi’an 710048, China
2
Hydrology and Water Resources Bureau of Henan, Yellow River Conservancy Commission, Zhengzhou 450000, China
3
School of Statistics and Mathematics, Central University of Finance and Economics, Beijing 100081, China
4
School of Environmental Studies, China University of Geosciences, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Water 2021, 13(18), 2458; https://doi.org/10.3390/w13182458
Submission received: 20 July 2021 / Revised: 5 September 2021 / Accepted: 6 September 2021 / Published: 7 September 2021

Abstract

:
The Han River is the water source of the South-to-North Water Diversion Project and the “Han River to Wei River Water Diversion Project” in China. In order to ensure that the water quality and quantity are sufficient for the water diversion project, the natural forest protection project, river chief system and other measures have been implemented in the Han River by the government. At the same time, several large reservoirs have been built in the Han River basin and perform the functions of water supply and hydropower generation, which is an important type of clean power. Under the influence of human activities, the coupling interaction between humans and water in the upper reach of the Han River drives the socio-hydrological evolution process. In this study, from the perspective of socio-hydrology, a model of socio-hydrological evolution (SHE) in the Han River in southern Shaanxi was built to simulate the potential evolution path of the socio-hydrological system and determine possible measures for the sustainable governance of the river basin. By adjusting the model parameters, four future scenarios were established: natural continuation, economic development, environmental protection and industrial adjustment scenarios. Taking 2018 as the base year, the evolution of socio-hydrology in the upper reaches of the Han River was predicted under the four scenarios from 2019 to 2045. The simulation results show that: (1) In the entire study area, except for domestic water, the water consumption of other departments show an upward or stable trend. There are differences in water consumption changes in the upper and lower sections, which are related to the different socio-economic conditions. (2) Comparing different development scenarios, the industrial adjustment and environmental protection scenarios are superior to the other scenarios. Natural continuation and economic development scenarios appear to be unfavorable for the sustainability of water resources and the economy. (3) In addition, based on the development scenarios, some policy suggestions are put forward, such as reducing the irrigation water quota, appropriately adjusting the industrial structure and promoting the growth of the urban population and the development of urbanization.

1. Introduction

The period of the industrial revolution is known as the Anthropocene, during which humans have changed the Earth on an unprecedented scale [1]. Human beings directly or indirectly change the land water cycle by intercepting rivers, transferring water across river basins and discharging wastewater into rivers [2,3]. In addition, changes in hydrological characteristics also cause changes in the ecosystem [4]. The connections between people and water systems are highly complex [5]. Several concepts, such as virtual water [6], the water footprint [7] and water value stream [8], attempt to combine humans and hydrology. Water is not only the core element of the ecological and environmental system but also a strategic resource for the development of human society and the economy. Fresh water supply and water quality are the key factors that affect and constrain human living, food production and economic growth. However, a growing number of signals, ranging from declining groundwater levels and deteriorating lake water quality to the disappearance of wetlands, indicate that current water system use is unsustainable [9,10], and the decline in water resources has been compounded by urbanization and population growth [11]. With the expanding human footprint in the hydrological cycle, it is clear that humans need to adapt their activities to water systems to achieve co-evolution [12,13]. The study of human–water co-evolution has long attracted the attention of researchers [14,15]. In 2012, after summarizing previous studies, Sivapalan first proposed the definition of socio-hydrology [16], which has since become an independent hydrologic discipline. At the same time, he encouraged a broader view of socio-hydrological systems in the context of socio-environmental systems.
Subsequently, the discussion of socio-hydrology has gradually intensified. Coupled systems are not limited to social and hydrological subsystems but often involve ecological subsystems as well. Firstly, the study of socio-hydrology focuses on building a theoretical framework. Elshiafei established a general conceptual framework for a socio-hydrological model applicable to agriculture and selected some general proxy indicators for different model modules, such as water quality change, vegetation change and bird species abundance, to represent ecological service functions [17]. In addition, human activities such as coal mining and reservoir operation have a significant impact on the hydrological process in basins [18,19]. Reservoir operation was modeled with the river ecohydrological process, and the demands of different types of hydropower stations on the river discharge were studied [20]. Subsequently, researchers attempted to establish a conceptual socio-hydrological model by describing regional socio-hydrological mechanisms using equations and by introducing socio-economic variables. Di Baldassarre explored the interaction between human settlements and flooding [21]. Srinivasan developed a human–hydrodynamic model in India [22]. Van Emmerik developed a conceptual model of the relationship between society and hydrology in river basins with extensive irrigation [23]. The co-evolution of humans and water in the Tarim River Basin of China was analyzed from a socio-hydrologic perspective using the Taiji-Tire Model, and Liu established a human–water co-evolutionary coupling model for the Tarim River Basin [24,25]. Lu used sentiment analysis to assess the cooperation needs of various countries and established a socio-hydrological model of trans-boundary river cooperation (TCSH model) to simulate the cooperation and conflict observed in the Lancang-Mekong River Basin [26]. Halder used social-hydrological approach to predict groundwater demand in different sectors of the Sundarbans [27]. Kumar analyzed the water scarcity on three very large river islands in Asia, and advocates that socio-hydrological research should focus on how to improve the adaptability of the vulnerable coastal island communities [28].
Studies on the socio-hydrological coupling system mainly adopt evolution models. The socio-hydrological models generally use dynamic equations to describe the process of evolution. Compared with the objective optimization model, which can only simulate the stable state of the system, the social hydrological model is based on process simulation. Thus, the socio-hydrological model can simulate the evolution process of the system under a changing environment and can be used to study the evolution of the socio-hydrological system, also known as the socio-hydrological evolution model. The models mainly include the system dynamics model (SDM), the agent-based model (ABM) and system of systems (SoS) [29]. The agent-based model (ABM) is mainly used to estimate the loss of human life and property in the event of environmental collapse [30,31]. Due to the need to predict the subject behavior of each person in the event of a disaster, the performance of the ABM is poor in the face of more complex socio-hydrological processes [32,33]. Some scholars established the SOS biofuel model to realize the interaction between humans and water in the context of the National Biofuel Policy [34,35]. The system model needs to consider not only the relationship between environment, policy and water quality but also the various subsystems involved in biofuel development (hydrology, economics, transport, refining and agriculture) [36]. Because the scale is large and the computational complexity is high, it is difficult to grasp and generalize. The SDM is based on the use of a feedback loop and time delay to simulate the co-evolutionary process of the system. The coupling simulation can be carried out using the MATLAB (MathWorks, Natick, MA, USA) and Vensim (Ventana Systems, Inc., Harvard, MA, USA) software packages as a visual platform, which makes the model easy to calculate and generalize [37].
The SDM is a qualitative and quantitative method for system integration simulation analysis, which can estimate changes in a complex system with time [38]. The system dynamics model, also known as the stylized model, can use a series of discrete or continuous equations to represent the change in state variables [39,40]. EI-Gafy established a model with a series of discrete equations to analyze the water footprint and virtual water balance of wheat production and consumption in Egypt [41]. Feng established several continuous forms of joint difference equations to analyze the co-evolution among water supply, power generation and the environment in the context of reservoir operation and interest preference. In this study, a water quantity equation was established based on the system dynamics method to simulate the evolution of the socio-hydrological system in the upper reach of the Han River in China [42].
During the last 10 years, there have been several studies on the socio-hydrological system [43]. By reviewing the research progress in social hydrology, several researchers have reported that the coupling scale, quantization method and system sensitivity are key factors that introduce challenges to the coupling of the social eco-hydrological system [29,44,45]. The spatial scale of ecological and hydrological systems is the river basin, while that of the social system is the administrative region; the time scale of socio-economic data is a year, while that of hydrological and meteorological data is usually a month or day. Therefore, conversion between different spatial and time scales is a problem that should be addressed. In this study, in order to overcome the problem of converting time and space scales, the watershed boundary was matched with the administrative unit scale when dividing the research area, which is more conducive to collecting accurate demographic and socio-economic data. Then, the collected social, hydrological and meteorological data were unified into the annual time scale. In addition, quantifying factors in the social system, such as socio-economic, population and environmental variables, is also a problem that needs to be addressed in the modeling process [46]. In the modeling of socio-hydrological models, the main limitation is that the models are often not accurate due to the uncertainty in input parameters and the process description [47]. On the other hand, the distribution of relevant data in socio-hydrology is very small, and it is difficult to completely reconstruct historical socio-hydrological processes from these data. The limitation of using data-driven methods is that causality cannot be easily determined through empirical data alone [48]. However, in many cases, despite the uncertainty in the data, hydrological and social data are generally more reliable and provide a better understanding of socio-hydrological processes than numerical models [49]. Therefore, considering the advantages and disadvantages of models and data, socio-hydrological modeling studies usually use statistical data to calibrate and validate models. However, in the face of complex system models, researchers need to strike the right balance between generality, accuracy and authenticity [50]. It is worth mentioning that Feng et al. established comprehensive equations and introduced many parameters [39], which are very useful for reshaping the socio-hydrological process in the past and for future development. The policy implications of the simulated results are discussed to provide suggestions for the future development of the study area. This study aims to establish a socio-hydrological evolution (SHE) model to simulate the potential evolution path of the socio-hydrological system in the upper reach of the Han River in China, demonstrate the capability of the system dynamics model in socio-hydrological modeling and determine possible measures for the sustainable governance of the river basin. The remainder of this paper is organized as follows. The study area, description of the input data and the model are introduced in Section 2. Section 3 and Section 4 present the construction and testing of the socio-hydrological model. Section 5 analyzes the future evolution of the system by setting different scenarios according to the current development situation. Section 6 presents the simulation results and some suggestions for the development of the study area. Finally, conclusions are drawn in Section 7.

2. Study Area and Data

2.1. Study Area

The Han River, the largest tributary of the Yangtze River, flows through Shaanxi and Hubei provinces. The main stream is 1577 km long, with a basin area of 159,000 km2 and average annual water resources of 56.6 billion m3. Over 2700 reservoirs have been built along the river, including Shiquan Reservoir, Ankang Reservoir and Danjiangkou Reservoir. The construction of seven cascade hydropower plants is planned in Shaanxi Province, and five of them have been built or are under construction.
The study area is located in the upper reach of the Han River in Shaanxi Province of China, mainly including Hanzhong City, Ankang City and some counties of Shangluo City. The study area is about 51,833 km2, accounting for 82.75% of the area of the Han River basin in Shaanxi Province. According to the distribution of the river network, the study area is divided into two sections, shown by red lines in Figure 1: the upper section is the upstream catchment of Shiquan hydrological station, and the lower section is the catchment between Shiquan hydrological station and Baihe hydrological station.
The social system is an important component of the SHE system. From 1990 to 2018, the total population of the study area increased from 6.28 million to 6.84 million, and the urbanization rate increased from 15.68% to 26.88%. With the rapid development of the social economy in the upper reach of the Han River, regional water consumption is constantly increasing. From 2001 to 2018, the water consumption of Ankang and Hanzhong increased from 2.083 billion m3 to 2.418 billion m3 (http://slt.shaanxi.gov.cn, accessed on 20 July 2021). Social systems affect the other two SHE systems. The development of the social economy is accompanied by pollution of the ecological environment, and the annual pollution emissions of the two cities increased from 78 million tons to 107 million tons (http://slt.shaanxi.gov.cn, accessed on 20 July 2021). The expansion of urbanization has led to a reduction in the area of cultivated land, with the study area decreasing from 520,000 hectares in 1990 to 370,000 hectares in 2018. Water pollution has led to the deterioration of river water quality, which, in turn, has triggered changes in the classification and structural function of the biome in the upper reach of the Han River [51,52]. In response to the deterioration of the ecological environment, the government implements macro-control policies. Moreover, rainfall and runoff in the upper reach of the Han River have been decreasing since 1980, while evaporation has been increasing [53]. Human activities and climate change are the main factors affecting runoff change in the upper reaches of the Han River [54,55]. The reduction in rainfall and runoff leads to less regional inflow and restricts social development. Therefore, the conflict between humans and water in the Han River basin is prominent, the interaction between humans and water is highly complex, and the social, ecological and hydrological systems form feedback loops. In the ever-changing environment, it is necessary to conduct coupling simulations among social, ecological and hydrological systems in the upper reaches of the Han River.

2.2. Data

The data include hydrological and meteorological data, land use data and population and socio-economic development data [53,56]. The potential evaporation and precipitation data (1980–2018) in the Upper section were obtained by using the weighted average of the four meteorological stations: Hanzhong, Liuba, Foping and Shiquan stations. The weighted average values of Shiquan, Ankang and Zhenan stations were used to obtain the precipitation and potential evaporation data (1980–2018) in the Lower section. Runoff data (1956–2018) were collected from Shiquan and Baihe hydrological stations, and the missing data from 2000 to 2005 at Baihe hydrological station were interpolated by the rainfall–runoff relation. Land use data (2001–2014) were extracted from MCD12Q1 remote sensing data from MODIS (https://ladsweb.modaps.eosdis.nasa.gov/search/, accessed on 20 July 2021). Demographic and socio-economic development data (1990–2018) were derived from the Shaanxi Statistical Yearbook, Shaanxi Regional Yearbook, Ankang Yearbook, Shangluo Yearbook, Hanzhong Yearbook, Xi’an Statistical Yearbook and Baoji Yearbook.

3. Socio-Hydrological Evolution Model Based on System Dynamics

We establish a socio-hydrological evolution model (SHE model) that can capture the evolutionary trajectory of coupled human-water systems. The SHE model based on system dynamics in the Han River describes the causal relationship among variables in the system, as shown in Figure 2, and then formulates the mathematical expression of system dynamics; that is, it constructs equations. The model framework consists of three main parts: the water consumption module, evaporation module and water balance module. There are three modules in the upper and lower sections, and each module contains the same variables, including state variables, rate variables, auxiliary variables and so on. There is no upstream inflow in the Upper section, and the outflow in the upper section is the inflow in the lower section. The outflow in the lower section is the runoff observed by the Baihe hydrological station. At the annual time scale, the socio-hydrological model constructed in this study uses the water balance equation to describe the key interactions in the complex system. Vensim PLE 8.2.1 software (https://vensim.com/, accessed on 20 July 2021) was used to build the model.

3.1. Governing Equations

The governing equations of the socio-hydrological system in the upper reach of the Han River are the water balance, population, socio-economic and water quota equations, which represent the hydrology, society, economy and policy subsystems, respectively. The subscript “U” represents the upper section, and the subscript “L” represents the lower section. See Appendix A Table A1 for the definitions and symbols of the variables in the model.

3.1.1. Water Balance Equation

With reference to Feng et al.’s research [39], we assumed that the study area is a closed system. The water storage is recharged by inflow, consumed by domestic use, agriculture, industry and so on, and finally released downstream. In addition, referring to Liu et al.’s method of calculating evaporation [25], evaporation was added to the water balance.
The water balance of each section is:
S U , t = S U , t 1 + I U O U E U W U S L , t = S L , t 1 + I L L O L E L W L
where S U , t oration in the Upper section, and the subscript t (year) represents time. The definitions of these variables in the Lower section are similar to those in the Upper section, so only the notation of the Upper I U :
I U = P U R A I N × A U I L = P L R A I N × A L + O U
where P U R A I N (mm) and P L R A I N (mm) are the rainfall in the Upper and Lower sections, and A U (km2) and A L (km2) are the areas of the Upper and Lower sections. The outflow ( O ) (m3) is related to the inflow ( I ), evaporation ( E ) and water storage ( S ). The formula of the inflow is shown in the latter section.
The actual evapotranspiration formulas are:
E U = E U G F + E U U + E U W E U G F = f A U G F , E U P   = A U G F × α U G F × E U P E U U = f A U U , E U P   = A U U × α U U × E U P E U W = f A U W , E U P   = A U W × α U W × E U P
where E U P (mm) is the annual potential evaporation in the Upper section. A U G F (km2), A U U (km2) and A U W (km2) are the grass and forest area, the urban area and the water area in the Upper section. E U G F (m3), E U U (m3) and E U W (m3) are the grass and forest evaporation, the urban evaporation and the water evaporation in the Upper section. α U G F , α U U and α U W are the empirical coefficients.
The water consumption formula is:
W U = W U A + W U D + W U I + W U E + W U T
where W U A (m3), W U D (m3), W U I (m3), W U E (m3) and W U T (m3) are the agricultural, domestic, industrial, environmental and tertiary industry water consumption in the Upper section. Water consumption is calculated by the quota method, and the equations are shown in Appendix A Table A2 and Table A3.

3.1.2. Social Equation

The Malthus population prediction model was used to predict the population in the upper reach of the Han River. This model regards the population growth rate as a constant. Based on the population in the base year, the annual population is predicted by calculating the statistical population growth trend:
P t = P 0 × 1 + r P n
where P t is total population, r P is the natural rate of population change, n is the year number, and P 0 is the base year population.

3.1.3. Economic Equation

Industrial added value, cultivated land area, livestock quantity, hospital beds and tertiary industry employees are used to represent the economic scale of the community. Cultivated land area, livestock quantity and hospital beds were predicted by the constant growth rate method or the constructed constitutive equation, which can be seen in Appendix A Table A2 and Table A3.
The logistic model was used to predict the industrial added value ( V A I , t ). The formula is:
V A I , t = V A I , m a x 1 + V A I , m a x V A I , 0 e r V A I × t t 0
where r V A I is the change rate of industrial added value, V A I , m a x is the maximum industrial added value that the environment can accommodate, and V A I , 0 is the initial value.

3.1.4. Water Quota Equation

In this study, the water quota method was used to calculate the water consumption of each industry. Since 2010, water consumption per 10,000 yuan (Chinese yuan) of GDP and water consumption per 10,000 yuan of industrial added value have been used as important indicators of China’s energy management. Water consumption per 10,000 yuan of industrial added value was selected to estimate industrial water consumption. The quotas of industrial added value water and irrigation water are from “Shaanxi Province Water Resources Bulletin” (http://slt.shaanxi.gov.cn/, accessed on 20 July 2021). The water quota per 10,000 yuan of industrial added value shows a decreasing trend. According to historical observation data, the DOSERESP model with the S-shaped curve can be appropriately adjusted to simulate the water quota for industrial added value ( q V I A ):
q V I A = A 1 + A 2 A 1 1 + 10 t 0 t p
where A 1 is the minimum value, A 2 is the maximum value, t 0 is the center point, and p is the growth rate.
The quota of irrigation water ( q W I ) has no obvious trend, and the measured values are used as input data. It is assumed that the quota of irrigation water from 1990 to 2005 is the same as that in 2006, and the irrigation water quota for future scenarios is the same as that in 2018. The other water quotas slightly change according to “Norm of water intake for industries of Shaanxi Province (DB61/T943–2020)”, as shown in Table 1.

3.2. Constitutive Relationship

The constitutive relationship, which is established through historical data, is an equation that quantifies the internal causality among variables and is also the decisive factor driving the co-evolution of the socio-hydrological system. A total of 42 constitutive equations were established for the model. At present, a method that can determine the functional relationships between variables of the system scientifically, universally and accurately has not been proposed. Sivapalan and Blöeschl point out that the functional relationships can be estimated using intuition, through data analysis (if the appropriate data exist), from related studies or based on consensus principles [40]. These methods were used in our study. According to Jia’s method [57], constitutive equations can be divided into two categories: in Class I, the simulation variables of the constitutive model are established based on the growth rate or time based on historical data (in Appendix A Table A2), and in Class II (in Appendix A Table A3), the constitutive relationship is estimated based on observation data or experts, combined with modeling purposes, model calibration and other relevant variables. In this study, the focus is on the estimation of the Class II constitutive relationship. For example, the outflow is controlled by the inflow, evaporation and water storage. When the inflow and water storage are large, the outflow released to the downstream area is large; when the evaporation is large, the outflow released to the downstream area is small. Domestic water consumption is determined by the population and its water quota. When one of these variables—the population or its water quota—increases while the other remains unchanged, the water consumption increases.

4. Calibration and Validation of Socio-Hydrological Evolution Model

Based on the current situation of the development and utilization of water resources in the upper reaches of the Han River, this study adopted the prototype framework proposed by Elshafei et al. [17] for socio-hydrological modeling and the general steps of socio-hydrological modeling proposed by Sivapalan and Blöschl [40]. The causality relationship among social, ecological and hydrological subsystems within the region is taken as the core to establish the equations, and the subsystems are coupled by the principle of system dynamics, which is consistent with some previous methods for establishing socio-hydrological models [42,57], to formulate a socio-hydrological model and simulate the evolutionary behavior of the system in the upper reaches of the Han River. The variables used as input data are ignored in the diagram in Figure 2. The time step of the model is 1 year. The model divides the study area into two parts (the upstream and downstream parts), adopts Shiquan hydrological station as the boundary and embeds the county boundary in the social system into the watershed scale for scale unification. The equations representing the relationships among the variables in the model quantify the interactions between humans and water in the socio-hydrological system [21,25,50,58].
The simulation period is 1990–2045, in which 1990–2018 is the historical data period, 2019–2045 is the model simulation period, and the simulation time step is 1 year. The historical period (1990–2018) is divided into the calibration period (1990–2009) and the validation period (2010–2018). The initial values of the main variables are shown in Table 2, which are estimated based on historical data.
According to the causality and equation of the model, the parameters that need to be calibrated in the model include the change rate of the population ( r P ), the change rate of the urbanization rate ( r R U ), the change rate of cultivated land area ( r C ), the change rate of cows ( r C O W ), the growth rate of industrial added value ( r V I A ) and the water quota change rate of industrial added value ( r q V I A ). The variables that are validated by historical data include total population ( P ), urban population ( P U ), cultivated land area ( A C ), the number of cows ( N C O W ), industrial added value ( V I A ) and industrial added value water quota ( q V I A ). In addition, the equations for the numbers of pigs and sheep ( N S H E E P , N P I G ) are fitted by historical data, and there are no parameters to calibrate. The potential evaporation ( E P ) and the annual precipitation are input data. The values of calibration parameters are shown in Table 3, and the evaluation results in the calibration and validation periods are shown in Table 4 and Table 5. Mean Absolute Percentage Error (MAPE), Nash–Sutcliffe efficiency coefficient (NSE) and the coefficient of determination (R2) are used to evaluate the variables simulated by the model. When NSE is 1, the simulated value equals the observed value.
The simulation results in the calibration period are good. Except for the number of cows, the simulation results are better in the calibration period. The simulation performance of livestock ( N C O W , N S H E E P , N P I G ) and cultivated land area ( A C ) in the validation period is poor, and the MAPE of industrial added value is more than 30%. This is because the rapid development of the economy and urbanization led to a sharp decline in the area of cultivated land around 2000; after that, under the influence of national policy regulation, the area of cultivated land was stable and showed a slight upward trend (Figure 3e,f). Cultivated land area was simulated in the validation period according to the change rate in the historical period, so the simulation results are poor. In Section 5, the growth rate of cultivated land area in the simulation period is adjusted according to the actual situation to better reflect the change law in the validation period. According to the irrigation water consumption data published in the Water Resources Bulletin (http://slt.shaanxi.gov.cn/zfxxgk/fdzdgknr/zdgz/, accessed date on 20 August 2021), the simulated irrigation water consumption is in line with the actual situation. The simulation results are reported in Section 6. N C O W , N S H E E P and N P I G are affected by the fluctuation of prices in the market. The number of cows dropped rapidly in 2007 before rising to average levels in 2010, and then it dropped sharply again in 2018 (Figure 3g,h). At the same time, the fluctuations of the number of pigs and sheep are larger (Figure 3i–l). These variables are greatly affected by external conditions (outside the catchment), so large errors in the validation period are acceptable. In addition, agricultural water consumption consists of livestock water consumption and irrigation water consumption; the proportion of livestock water consumption is small and is not the focus of the model simulation, so it has little impact on the results. The simulation performance of other variables is good. In general, the simulation results of the model are acceptable, and the system dynamics model can be used to analyze the possible trajectories of future co-evolution. The model is a process-based model, and the simulation is driven by external meteorological conditions, i.e., precipitation and potential evaporation. The system does not exhibit any tipping point behavior under average climate conditions.

5. Development Scenarios of the Socio-Hydrological System in the Han River

The study area is a water source protection area. Since 2014, the population in the study area has been declining. According to the Population Development Plan of Shaanxi Province (2016–2030), the urbanization rate of household registration in Shaanxi Province was 39% in 2015 (http://www.shaanxi.gov.cn, (accessed on 20 August 2021), while the urbanization rate in the study area was only 23.8% in 2014, which is low. In the region, agriculture consumes the most water, industrial development is rapid, and the proportion of the tertiary industry is low. This study mainly analyzes changes in regional water use with changes in the population, urbanization and industrial structure. According to the development characteristics and actual conditions of Hanzhong City and Ankang City, six sensitive parameters, namely, urbanization rate, cultivated land area change rate, irrigation water quota, change rate of forest and grassland area, growth rate of industrial added value and the proportion of tertiary industry employees, were selected in both the upper and lower sections in this study. Based on the current base year (2018), these six sensitive parameters were adjusted, and four development scenarios were established: natural continuation, environmental protection, industrial adjustment and economic development.
In the natural continuation scenario, all parameters are set to the values in the base year (2018). In the environment protection scenario, the main goal is to protect the ecology and environment, achieve high-quality and low-speed development of industry and increase the proportion of tertiary industry employees appropriately. In the industrial adjustment scenario, the structure of water use is changed, the urbanization level is maintained at a medium-high speed of development, the growth rate of industrial added value is reduced, and the proportion of the tertiary industry is increased. In the economic development scenario, the main goal is to increase the economic development, so the proportion of urbanization, the growth rate of industrial added value and the urban area increase. Since the hydro-meteorological conditions cannot be adjusted or analyzed through policies, scenario analysis was carried out under the assumption that the hydro-meteorological conditions remain stable. Taking 2018 as the base year, the evolution of the socio-hydrological system from 2019 to 2045 was simulated under the different development scenarios. The values of the sensitive parameters in the study area under the four scenarios are shown in Table 6 and Table 7.

6. Results and Discussion

6.1. Results

The simulation results are shown in Figure 4 and Figure 5. In both the upper and lower sections, the total water consumption shows an upward trend in all four scenarios. The greatest water consumption is by industrial and agricultural activities, accounting for more than 90%, followed by domestic water consumption, accounting for about 10%; the proportion of environmental water and tertiary industry water is the lowest. Because the parameters adopt the values in the scenarios after 2019, the simulated variables show jumps in Figure 4a,b,g,h and Figure 5a,b,g,h, respectively.
Figure 4 shows the future evolution path of the main variables in the upper section. As shown in Figure 4a, during the simulation period, the population declines significantly in the natural continuation scenario and slightly increases in the economic development scenario. The population is stable in the environmental protection and industrial adjustment scenarios. The urban population (Figure 4b) shows an increasing trend in the four scenarios. As shown in Figure 4e, the average total water consumption (2019–2045) in the region shows an upward trend. The order of the average total water consumption among scenarios is natural continuation (18.62 × 108 m3/year) > economic development (17.74 × 108 m3/year) > industrial adjustment (17.05 billion m3/year) > environmental protection (17.01 × 108 m3/year). Among all water uses, agricultural water consumption (Figure 4f) is the highest, followed by industrial and domestic water consumption, and the tertiary industry and environment water consumption are the lowest. Agricultural water consumption increases significantly under the natural continuity scenario, while it declines under the economic development scenario. The changing trend of agricultural water consumption under different scenarios is the same as that of cultivated land area (Figure 4c). Industrial water consumption increases rapidly, which is consistent with the increasing trend of industrial added value (Figure 4d). In all scenarios, the tertiary industry and environmental water consumption are below 0.3 × 108 m3.
In the upper section, comparing the values in 2018 and 2045 under different development scenarios shows that: the population drops by 14.17% in the natural continuation scenario (Figure 4a), domestic water consumption drops by 9.21%, cultivated land area increases by 9.83%, and industrial added value increases by about five times. Domestic and tertiary industry water consumption decrease. Agricultural and industrial water consumption increase the most, so the total water consumption is the highest, while the environmental water consumption is the lowest in all scenarios. Under the natural continuation scenario, the upper section will face several problems in the future, such as population loss, a low level of urbanization, reduced runoff and an imbalance between the supply and demand of water resources.
Under the economic development scenario, the proportion of industrial water rises from 4.3% in 2018 to 20.49% in 2045, and the cultivated land area decreases the most (by 15.34%). The water consumption increases by 24.05%, which is less than the total water use in the natural continuation scenario. Long-term development of this scenario may lead to a shortage of cultivated land resources, which will limit the development of agriculture.
In the environmental protection scenario, the water consumption in the upper section is the lowest among the four scenarios. In this scenario, the cultivated land shows a stable and slightly rising trend, and the level of economic development is low.
In the industrial adjustment scenario, the tertiary industry water consumption has a clear upward trend, while the proportion of agricultural water consumption drops by 10%; however, agriculture still consumes the most water. The urban population and industrial added value are higher than those in the environmental protection scenario, and the cultivated land area decreases by 8.44%. Water consumption in the industrial adjustment scenario is slightly higher than that in the environment protection scenario.
Figure 5 shows the future evolution path of the main variables in the lower section. In the four future scenarios, the total water consumption in the region shows an upward trend. The order of total water consumption among scenarios is economic development (11.94 × 108 m3/year) > environmental protection (10.57 × 108 m3/year) > industrial adjustment (10.51 × 108 m3/year) > natural continuation (9.71 × 108 m3/year). The ranking of scenarios according to total water consumption differs between the lower and upper sections. The current water consumption in the lower section is lower than that in the upper section. In the lower section, the changing trends of the population, urban population, industrial added value, tertiary industry personnel, domestic water consumption, tertiary industry water consumption and industrial water consumption are the same as those in the upper section. The cultivated land area is the lowest in the natural continuation scenario. Compared with the historical period, the cultivated land area in the simulation period is reduced by 27.5%, and the agricultural water consumption is the lowest. In all scenarios, the domestic water consumption in the upper section is higher than that in the lower section during the simulation period, because the population in the lower section is larger than that in the upper section. The rest of the water consumption in the upstream area is higher than that in the lower section.
In the lower section, comparing the values in 2018 and 2045 under different development scenarios indicates that: in the natural continuation scenario, the population and cultivated land area markedly decrease, and the urbanization and industrial levels are low, so the water consumption is the lowest among the four scenarios. In the economic development scenario, the most prominent features are the rapid increase in the population, industrial added value and urbanization rate. The proportion of industrial water consumption increases from 2.64% to 25.82%, and the proportion of agricultural water consumption drops from 82.36% to 58.85% in 2045. In the industrial adjustment and environmental protection scenarios, urbanization and industrialization development are moderate, and the cultivated land area increases by 1.20% only in the environmental protection scenario. The water consumption in the industrial adjustment scenario is lower than that in the economic development scenario and higher than that in the natural continuation scenario: this is because agriculture consumes the largest amount of water, and the cultivated land area drastically decreases in the natural continuation scenario, so the agricultural water consumption is much lower than in other scenarios. The tertiary industry and environmental water consumption proportions increase but remain at low levels.

6.2. Implications

Based on the analysis and comparison of the evolution paths of the socio-hydrological system in the four scenarios, and combined with research on regional policies and planning, the following suggestions are put forward to achieve sustainable development in the Han River basin in southern Shaanxi:
(1)
Improve agricultural irrigation technology and reduce the irrigation water quota. The study area is dominated by agricultural water consumption. In 2018, the proportion of agricultural water consumption in the upper section was 87.7%, and that in the lower section was 82.4%, so reducing agricultural water consumption can effectively reduce water consumption. The Chinese government and the Shaanxi provincial government have issued the “strictest farmland protection policy” and included cultivable land area protection as a target in the government assessment [59,60]. In order to control agricultural water consumption while protecting the cultivable land area, the irrigation water quota should be reduced. Under the natural continuation scenario, water consumption in the upper and lower sections in 2045 is predicted to be 21.61 × 108 m3 and 11.34 × 108 m3, respectively. As shown in Table 8, with other conditions remaining unchanged and the irrigation water quota reduced by 5%, 10%, 15% and 20%, the water consumption in the upper section will be reduced by 3.8%, 7.7%, 11.5% and 15.3%, respectively, and that in the lower section will be reduced by 3.4%, 6.9%, 10.3% and 13.8%, respectively, by 2045. The reduction in the irrigation water quota depends on the development of irrigation technology, which can also promote the further development of modern agriculture.
(2)
Accelerate the development of urbanization. According to the “13th Five-Year” Plan of Shaanxi Province, the target value of the urbanization rate in 2020 was 45%, and the actual value was 48.5% [61]. According to scenario analysis, the urbanization rate of the study area is the highest in the economic development scenario and would be about 30% in 2020. The urbanization development rate of southern Shaanxi in the Han River basin is relatively low. Therefore, the urbanization process should be accelerated. Furthermore, the transformation of the rural to urban population is conducive to the restoration of forest and grassland [62].
(3)
Adjust the industrial structure appropriately. In 2015, the proportion of GDP from the tertiary industry in the Han River basin in southern Shaanxi was 35.78%, and the simulation found that the tertiary industry water consumption in the study area accounted for 1.47% of the total water consumption. In the industrial adjustment scenario, if the industrial development speed is appropriately reduced and the proportion of tertiary industry personnel is increased, then the total water consumption will significantly decrease. Therefore, the adjustment of the industrial structure should be accelerated to encourage the development of the tertiary industry due to its low water consumption and high output value.

7. Conclusions

In this study, a socio-hydrological model for the upper reaches of the Han River in Shaanxi province of China was established by using the system dynamics method. The main variables representing the population, socio-economic development and hydrological changes were selected to establish the governing equations. The evolution of the water quota for industrial added value, which reflects technology development in industries and agriculture, was simulated using the model. The model was calibrated and validated by historical data of the upper reaches of the Han River. The simulated evolution processes reproduced the historical processes of the socio-hydrological system. After the calibration of the model, four scenarios, i.e., natural continuation, economic development, environment protection and industrial adjustment, were established, and the future evolution process of the socio-hydrological system in these four scenarios was analyzed.
According to the evolution analysis of the future scenarios, the population in the study area is declining, and the level of urbanization development is low. Water consumption rises in the four scenarios, and the proportion of agricultural and industrial water is more than 90%. The population loss is rapid in the natural continuation scenario. In the economic development scenario, industrial water consumption increases rapidly, but urban land occupies cultivable land, which is not conducive to agricultural development. If the river basin appropriately improves the level of urbanization and the speed of industrial development, then the industrial adjustment scenario is the best. If cultivated land resources are to be protected, then water consumption in the region should be strictly controlled, and the environmental protection scenario is the best. The simulation results of the four scenarios also provide a reference for water resource management in this region. Policymakers can formulate development policies to carry out macro-control and find a balance between economic development and environmental protection.
The socio-hydrological evolution model based on the system dynamics method was able to simulate the evolution process of the socio-hydrological system, and constitutive equations were established to couple factors representing social development and the hydrometeorology. The model provides a tool for understanding the potential development paths of socio-hydrology in the upper reaches of the Han River.
The socio-hydrological system is a complex system under the comprehensive influence of natural and social factors at different scales. In the present socio-hydrological model in the upper reaches of the Han River, the social and economic modules are only described by a few key variables, and dynamic descriptions of policy changes, environmental awareness and technological development level are lacking, so this is a preliminary simplified model. In future research, the model will be refined by applying a targeted approach, and the physical mechanism of the model will be improved.

Author Contributions

Conceptualization, methodology, and model developed by X.Z., D.L. and X.W.; validation and formal analysis conducted by X.Z., D.L. and X.W.; original draft prepared by X.Z.; review and editing carried out by D.L., L.M., M.L., and X.M.; and supervision supported by Q.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by the National Natural Science Foundation of China (NSFC) (Grant Nos. 51779203, 51609270, 51979252) and State Key Laboratory of Eco-hydraulics in Northwest Arid Region (Grant No. 2019KFKT-4). The authors thank the editor and anonymous reviewers for their constructive comments.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Land use data were extracted from MCD12Q1 of MODIS (https://ladsweb.modaps.eosdis.nasa.gov/search/, accessed on 20 July 2021). All of yearbooks are publicly published. The other data presented in this study are available on reasonable request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Variables and parameters used in the model.
Table A1. Variables and parameters used in the model.
DescriptionSymbol in Upper SectionSymbol in Lower SectionUnitType
Agricultural water W U A W L A m3Variable
Tertiary industry employees α U T P α L T P Variable
Cultivable land area A U C A L C haVariable
Domestic water consumption W U D W L D m3Variable
Environmental water consumption W U E W L E m3Variable
Forest and grass area A U G F A L G F km2Variable
Forest and grass evaporate E U G F E L G F m3Variable
Industrial added value V U I A V L I A 104/yuanVariable
Industrial water consumption W U I W L I m3Variable
Inflow I U I L m3Variable
Irrigation water consumption W U N W U N m3Variable
Livestock water consumption W U L S W L L S m3Variable
Number of cows N U C O W N L C O W perVariable
Number of sheep N U S H E E P N L S H E E P perVariable
Number of pigs N U P I G N L P I G perVariable
Outflow O U O L m3Variable
Population P U P L perVariable
Potential evaporation E U P E L P mmVariable
Precipitation P U R A I N P L R A I N mmVariable
Proportion of irrigated area α U I α L I Variable
Rural population P U R P L R perVariable
Surface water area A U W A L W km2Variable
Tertiary industry water W U T W L T m3Variable
Urban area A U U A L U km2Variable
Urban evaporation E U U E L U m3Variable
Water quota for industrial added value q U V I A q L V I A m3/104 yuanVariable
Catchment area A U A L km2Parameter
Conversion coefficient of forest and grassland evaporation α U G F α L G F Parameter
Conversion coefficient of urban evaporation α U U α L U Parameter
Conversion coefficient of water surface evaporation α U W α L W Parameter
Water quota for cows q U C O W q L C O W m3/perParameter
Water quota for hospital beds q U H B q L H B m3/bedParameter
Water quota for the tertiary industry employees q U T P q L T P m3/perParameter
Water quota for pigs q U P I G q L P I G m3/perParameter
Water quota for rural population q U R P q L R P m3/perParameter
Water quota for sheep q U S H E E P q L S H E E P m3/perParameter
Water quota for urban population q U U P q L U P m3/perParameter
Change rate of urbanization rate r U R U r L R U Parameter
Change rate of cultivated land r U C r L C Parameter
Change rate of industrial added value r U V I A r L V I A Parameter
Proportion of employees in the tertiary industry r U T P r L T P Parameter
Table A2. The constitutive equations of Class I.
Table A2. The constitutive equations of Class I.
VariableVariable NameEstimated Constitutive RelationCategory
N U C O W Number of cows in the upper section N U C O W , t = N U C O W , t 1 × 1 + r U C O W I
N L C O W Number of cows in the lower section N L C O W , t = N L C O W , t 1 × 1 + r L C O W I
N U P I G Number of pigs in the upper section N U P I G , t = 1252.82 t 2 + 5.11 × 10 6 t 5.21 × 10 9 I
N L P I G Number of pigs in the lower section N L P I G , t = 867.06 t 2 + 3.60 × 10 6 t 3.73 × 10 9 I
N U S H E E P Number of sheep in the upper section N U S H E E P , t = 160.46 t 2 + 6.52 × 10 5 t 6.61 × 10 8 I
N L S H E E P Number of sheep in the lower section N L S H E E P , t = 508.31 t 2 2.07 × 10 6 t + 2.11 × 10 9 I
A U C Cultivable land area in the upper section A U C , t = A U C , t 1 × 1 + r U C I
A L C Cultivable land area in the Lower section A L C , t = A L C , t 1 × 1 + r L C I
A U I Irrigation area in the Upper section A U I = A U C × α U I I
A L I Irrigation area in the Lower section A L I = A L C × α L I I
A U G F Grass and forest area in the Upper section A U G F , t = A U G F , t 1 × 1 + r U G F I
A L G F Grass and forest area in the Lower section A L G F , t = A L G F , t 1 × 1 + r L G F I
P U T Tertiary industry employees in the Upper section P U T = P U × α U T P I
P L T Tertiary industry employees in the Lower section P L T = P L × α L T P I
A U U Urban area in the Upper section A U U = 75 r U U P 0.25 r U P × 300 r U U P > 0.25 I
A L U Urban area in the Lower section A L U = 75 r L U P 0.25 r L P × 200 r L U P > 0.25 I
P U R Rural population in the Upper section P U R = P U P U U I
P L R Rural population in the Lower section P L R = P L P L U I
P U U Urban population in the Upper section P U U = P U × r U U P I
P L U Urban population in the Lower section P L U = P L × r L U P I
Table A3. The constitutive equations of Class II.
Table A3. The constitutive equations of Class II.
VariableVariable NameEstimated Constitutive RelationCategory
O U Outflow in the Upper section O U = f I U , E U , Δ S U = I U E U × 0.73 + Δ S U × 0.08 II
O L Outflow in the Lower section O L = f I L , E L , Δ S L = I L E L × 0.8 + Δ S U × 0.18 II
W U E Environmental water consumption in the Upper section W U E = f S U = 0.001 × 130 E X P ( α U E × S U II
W L E Environmental water consumption in the Lower section W L E = f S L = 0.001 × 120 E X P ( α L E × S L II
E U G F Grass and forest evaporation in the Upper section E U G F = f A U G F , E U P = A U G F × α U G F × E U P II
E L G F Grass and forest evaporation in the Lower section E L G F = f A L G F , E L P = A L G F × α L G F × E L P II
E U U Urban evaporation in the Upper section E U U = f A U U , E U P = A U U × α U U × E U P II
E L U Urban evaporation in the Lower section E L U = f A L U , E L P = A L U × α L U × E L P II
E U W Water surface evaporation in the Upper section E U W = f A U W , E U P = A U W × α U W × E U P II
E L W Water surface evaporation in the Lower section E L W = f A L W , E L P = A L W × α L W × E L P II
W U D Domestic water consumption in the Upper section W U D = f P U R , P U U = P U R × q U R + P U U × q U U II
W L D Domestic water consumption in the Lower section W L D = f P L R , P L U = P L R × q L R + P L U × q U U II
W U I Industrial water consumption in the Upper section W U I = f V U A I = V U A I × q U V A I II
W L I Industrial water consumption in the Lower section W L I = f V L A I = V L A I × q L V A I II
W U A Agricultural water consumption in the Upper section W U A = f W U N , W U L S = W U N + W U L S II
W L A Agricultural water consumption in the Lower section W L A = f W L N , W L L S = W L N + W L L S II
W U N Irrigation water consumption in the Upper section W U N = f A U I = A U I × q U I II
W L N Irrigation water consumption in the Lower section W L N = f A L I = A L I × q L I II
W U L S Livestock water consumption in the Upper section W U L S = f N U S H E E P , N U C O W , N U P I G = N U S H E E P × q U S H E E P + N U C O W × q U C O W + N U P I G × q U P I G II
W L L S Livestock water consumption in the Lower section W L L S = f N L S H E E P , N L C O W , N L P I G = N L S H E E P × q L S H E E P + N L C O W × q L C O W + N L P I G × q L P I G II
W U T Tertiary industry water consumption in the Upper section W U T = f P U T , N U H B = P U T × q U T P + N U H B × q U H B II
W L T Tertiary industry water consumption in the Lower section W L T = f P L T , N L H B = P L T × q L T P + N L H B × q L H B II

References

  1. Steffen, W.; Crutzen, P.J.; McNeill, J.R. The Anthropocene: Are Humans Now Overwhelming the Great Forces of Nature. AMBIO 2007, 36, 614–621. [Google Scholar] [CrossRef]
  2. Goudie, A.S.; Li, B.; Clark, W.C.; Kates, R.W.; Richards, J.F.; Mathews, J.T.; Meyer, W.B. The earth as transformed by human action; global and regional changes in the biosphere over the past 300 years. For. Ecol. Manag. 1990, 55, 341–342. [Google Scholar] [CrossRef]
  3. Luo, P.; Kang, S.; Zhou, M.M.; Nover, D. Water quality trend assessment in Jakarta: A rapidly growing Asian megacity. PLoS ONE 2019, 14, e0219009. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Zhu, Y.; Luo, P.; Zhang, S.; Sun, B. Spatiotemporal Analysis of Hydrological Variations and Their Impacts on Vegetation in Semiarid Areas from Multiple Satellite Data. Remote Sens. 2020, 12, 4177. [Google Scholar] [CrossRef]
  5. Steffen, W.; Grinevald, J.; Crutzen, P.; Mcneill, J. The Anthropocene: Conceptual and historical perspectives. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2011, 369, 842–867. [Google Scholar] [CrossRef]
  6. Emmett, C.F.; Rogers, P.; Lydon, P.; Wolf, A.T. Water in the Arab World: Perspectives and Prognoses. Geogr. Rev. 1996, 86, 152–154. [Google Scholar] [CrossRef]
  7. Hoekstra, A.Y. Virtual Water Trade: Proceedings of the International Expert Meeting on Virtual Water Trade; IHE Delft Institute for Water Education: Delft, The Netherlands, 2003; Volume 12, pp. 1–244. [Google Scholar]
  8. Jiang, W.L. Axiology of Water Resources; Science Press: Beijing, China, 1998. [Google Scholar]
  9. Viala, E. Water for Food, Water for Life: A Comprehensive Assessment of Water Management in Agriculture. Irrig. Drain. Syst. 2008, 22, 127–129. [Google Scholar] [CrossRef]
  10. Connor, R. Managing Water under Uncertainty and Risk: The United Nations World Water Development Report 4. United Nations World Water. 2012. Available online: https://www.issuelab.org/resource/managing-water-under-unce (accessed on 13 July 2021).
  11. Luo, P.; Sun, Y.T.; Wang, S.T.; Nover, D. Historical assessment and future sustainability challenges of Egyptian water resources management. J. Clean. Prod. 2020, 263, 121154. [Google Scholar] [CrossRef]
  12. Wagener, T.; Sivapalan, M.; Troch, P.A.; McGlynn, B.L.; Harman, C.J.; Gupta, H.V.; Kumar, P.; Rao, P.S.; Basu, N.B.; Wilson, J.S. The future of hydrology: An evolving science for a changing world. Water Resour. Res. 2010, 46, W05301. [Google Scholar] [CrossRef]
  13. Oreskes, N. How Earth Science Has Become a Social Science. Hist. Soc. Res. 2015, 2, 246–270. [Google Scholar]
  14. Pataki, D.E.; Carreiro, M.M.; Cherrier, J.; Grulke, N.E.; Jennings, V.; Pincetl, S.; Pouyat, R.V.; Whitlow, T.H.; Zipperer, W.C. Coupling biogeochemical cycles in urban environments: Ecosystem services, green solutions, and misconceptions. Front. Ecol. Environ. 2011, 9, 27–36. [Google Scholar] [CrossRef]
  15. Rasul, G.; Sharma, B. The nexus approach to water–energy–food security: An option for adaptation to climate change. Clim. Policy 2015, 16, 682–702. [Google Scholar] [CrossRef] [Green Version]
  16. Sivapalan, M.; Savenije, H.H.G.; Blöschl, G. Socio-hydrology: A new science of people and water. Hydrol. Process. 2012, 26, 1270–1276. [Google Scholar] [CrossRef]
  17. Elshafei, Y.; Sivapalan, M.; Tonts, M.; Hipsey, M.R. A prototype framework for models of socio-hydrology: Identification of key feedback loops with application to two Australian case-studies. Hydrol. Earth Syst. Sci. Discuss. 2014, 11, 629–689. [Google Scholar] [CrossRef] [Green Version]
  18. Luan, J.K.; Zhang, Y.Q.; Tian, J.; Meresa, H.; Liu, D.F. Coal mining impacts on catchment runoff. J. Hydrol. 2020, 589, 125101. [Google Scholar] [CrossRef]
  19. Zhang, L.P.; Huang, Q.; Liu, D.F.; Deng, M.J.; Zhang, H.X.; Pan, B.Z.; Zhang, H.G. Long-term and mid-term ecological operation of cascade hydropower plants considering ecological water demands in arid region. J. Clean. Prod. 2021, 279, 123599. [Google Scholar] [CrossRef]
  20. Chen, Q.; Chen, D.; Li, R.; Ma, J.; Blanckaert, K. Adapting the operation of two cascaded reservoirs for ecological flow requirement of a de-watered river channel due to diversion-type hydropower stations. Ecol. Modell. 2013, 252, 266–272. [Google Scholar] [CrossRef]
  21. Baldassarre, G.D.; Viglione, A.; Carr, G.; Kuil, L.; Salinas, J.L.; Blöschl, G. Socio-hydrology: Conceptualising human-flood interactions. Hydrol. Earth Syst. Sci. 2013, 17, 3295–3303. [Google Scholar] [CrossRef] [Green Version]
  22. Srinivasan, V. Reimagining the past—Use of counterfactual trajectories in socio-hydrological modelling: The case of Chennai, India. Hydrol. Earth Syst. Sci. 2015, 19, 785–801. [Google Scholar] [CrossRef] [Green Version]
  23. Van, E.; Li, Z.; Sivapalan, M.; Pande, S.; Kandasamy, J.; Savenije, H.; Chanan, A.; Vigneswaran, S. Socio-hydrologic modeling to understand and mediate the competition for water between agriculture development and environmental health: Murrumbidgee River basin, Australia. Hydrol. Earth Syst. Sci. 2014, 18, 4239–4259. [Google Scholar] [CrossRef] [Green Version]
  24. Liu, Y.; Tian, F.; Hu, H.; Sivapalan, M. Socio-hydrologic perspectives of the co-evolution of humans and water in the Tarim River Basin, Western China: The Taiji--Tire Model. Hydrol. Earth Syst. Sci. 2014, 18, 1289–1303. [Google Scholar] [CrossRef] [Green Version]
  25. Liu, D.F.; Tian, F.Q.; Lin, M.; Sivapalan, M. A conceptual socio-hydrological model of the co-evolution of humans and water: Case study of the Tarim River basin, western China. Hydrol. Earth Syst. Sci. 2015, 19, 1035–1054. [Google Scholar] [CrossRef] [Green Version]
  26. Lu, Y.; Tian, F.Q.; Guo, L.Y.; Borzi, I.; Patil, R.; Wei, J.; Liu, D.F.; Wei, Y.P.; Yu, D.J.; Sivapalan, M. Socio-Hydrologic Modeling of the Dynamics of Cooperation in the Transboundary Lancang-Mekong River. Hydrol. Earth Syst. Sci. 2020, 25, 1883–1903. [Google Scholar] [CrossRef]
  27. Halder, S.; Kumar, P.; Das, K.; Dasgupta, R.; Mukherjee, A. Socio-Hydrology Approach to Explore Groundwater–Human Wellbeing Nexus: Case Study from Sundarbans, India. Water 2021, 13, 1635. [Google Scholar] [CrossRef]
  28. Kumar, P.; Avtar, R.; Dasgupta, R.; Johnson, B.A.; Mukherjee, A.; Ahsan, M.N.; Nguyen, D.C.H.; Nguyen, H.Q.; Mishra, B.K. Socio-hydrology: A key approach for adaptation to water scarcity and achieving human well-being in large riverine islands. Prog. Disaster Sci. 2020, 8, 100134. [Google Scholar] [CrossRef]
  29. Liu, P.; Feng, M.Y.; Guo, S.L.; Li, Z.J. Methodologies and Challenges for Socio-Hydrology. J. Water Resour. Res. 2016, 5, 521–529. [Google Scholar] [CrossRef]
  30. Assaf, H.; Hartford, D. Physically-based modelling of life safety considerations in water resource decision-making. ASCE EWRI World Water and Environmental Resources Congress. Bridging the Gap: Meeting the World’s Water and Environmental Resources Challenges. In Proceedings of the World Water and Environmental Resources Congress, Salt Lake City, UT, USA, 27 June–1 July 2004. [Google Scholar] [CrossRef]
  31. Chen, X.; Zhang, F.B. Agent-based modelling and simulation of urban evacuation: Relative effectiveness of simultaneous and staged evacuation strategies. J. Oper. Res. Soc. 2008, 59, 25–33. [Google Scholar] [CrossRef]
  32. Dawson, R.J.; Peppe, R.; Miao, W. An agent-based model for risk-based flood incident management. Nat. Hazards 2011, 59, 167–189. [Google Scholar] [CrossRef]
  33. Evans, T.P.; Kelley, H. Multi-scale analysis of a household level agent-based model of landcover change. J. Environ. Manage. 2004, 72, 57–72. [Google Scholar] [CrossRef]
  34. Housh, M.; Cai, X.; Ng, T.L.; Mcisaac, G.F.; Ouyang, Y.; Khanna, M.; Sivapalan, M.; Jain, A.K.; Eckhoff, S.; Gasteyer, S. System of Systems Model for Analysis of Biofuel Development. J. Infrastruct. Syst. 2015, 21, 04014050. [Google Scholar] [CrossRef] [Green Version]
  35. Yaeger, M.A.; Sivapalan, M.; Mcisaac, G.F.; Cai, X. Comparative analysis of hydrologic signatures in two agricultural watersheds in east-central Illinois: Legacies of the past to inform the future. Hydrol. Earth Syst. Sci. 2013, 17, 4607–4623. [Google Scholar] [CrossRef] [Green Version]
  36. Yaeger, M.A.; Housh, M.; Cai, X.; Sivapalan, M. An integrated modeling framework for exploring flow regime and water quality changes with increasing biofuel crop production in the U.S. Corn Belt. Water Resour. Res. 2014, 50, 9385–9404. [Google Scholar] [CrossRef]
  37. Emmerik, T.V.; Elshafei, Y.; Mahendran, R.; Kandasamy, J.; Pande, S.; Sivapalan, M. Alternative socio-centric approach for model validation—A way forward for socio-hydrology. In Proceedings of the 19th EGU General Assembly Conference Abstracts, Vienna, Austria, 23–28 April 2017. [Google Scholar]
  38. Chen, Z.; Wei, S. Application of System Dynamics to Water Security Research. Water Resour. Manage 2014, 28, 287–300. [Google Scholar] [CrossRef]
  39. Feng, M.; Liu, P.; Li, Z.; Zhang, J.; Liu, D.; Xiong, L. Modeling the nexus across water supply, power generation and environment systems using the system dynamics approach: Hehuang Region, China. J. Hydro. 2016, 543, 344–359. [Google Scholar] [CrossRef]
  40. Sivapalan, M.; Blöeschl, G. Time scale interactions and the coevolution of humans and water. Water Resour. Res. 2015, 51, 6988–7022. [Google Scholar] [CrossRef] [Green Version]
  41. El-Gafy, I.K. System Dynamic Model for Crop Production, Water Footprint, and Virtual Water Nexus. Water Resour. Manag. 2014, 28, 4467–4490. [Google Scholar] [CrossRef]
  42. Feng, Z.K.; Niu, W.J.; Cheng, C.T. China’s large-scale hydropower system: Operation characteristics, modeling challenge and dimensionality reduction possibilities. Renew. Energ. 2019, 136, 805–818. [Google Scholar] [CrossRef]
  43. Herrera-Franco, G.; Montalván-Burbano, N.; Carrión-Mero, P.; Bravo-Montero, L. Worldwide Research on Socio-Hydrology: A Bibliometric Analysis. Water 2021, 13, 1283. [Google Scholar] [CrossRef]
  44. Tian, F.Q.; Cheng, T.; Lu, Y.; Xu, Z.X. A review on socio-hydrology and urban hydrology. Prog. Geo. 2018, 37, 46–56. [Google Scholar] [CrossRef] [Green Version]
  45. Ding, J.Y.; Zhao, W.W.; Fang, X.N. Advances in social hydrology. Chin. J. Appl. Ecol. 2015, 26, 1055–1063. [Google Scholar] [CrossRef]
  46. Xia, J.; Zhang, X.; Wei, F.L.; Wang, Q.; She, D.; Xu, J. Water system theory and its practices in China. South-to-North Water Transf. Water Sci. Technol. 2018, 16, 1–13. [Google Scholar] [CrossRef]
  47. Sawada, Y.; Hanazaki, R. Socio-hydrologic data assimilation: Analyzing human-flood interactions by model-data integration. Hydrol. Earth Syst. Sci. 2020, 1–54. [Google Scholar] [CrossRef] [Green Version]
  48. Dang, Q.; Konar, M. Trade Openness and Domestic Water Use. Water Resour. Res. 2018, 54, 4–18. [Google Scholar] [CrossRef] [Green Version]
  49. Mostert, E. An alternative approach for socio-hydrology: Case study research. Hydrol. Earth Syst. Sci. 2018, 22, 317–329. [Google Scholar] [CrossRef] [Green Version]
  50. Baldassarre, G.D.; Viglione, A.; Carr, G.; Kuil, L.; Yan, K.; Brandimarte, L.; Blschl, G. Debates-Perspectives on socio-hydrology: Capturing feedbacks between physical and social processes. Water Resour. Res. 2015, 51, 4770–4781. [Google Scholar] [CrossRef]
  51. Li, Z.; Wang, J.; Liu, Z.; Meng, X.; Heino, J.N.; Jiang, X.; Xiong, X.; Jiang, X.; Xie, Z. Different responses of taxonomic and functional structures of stream macroinvertebrate communities to local stressors and regional factors in a subtropical biodiversity hotspot. Sci. Total Environ. 2019, 655, 1288–1300. [Google Scholar] [CrossRef]
  52. Li, Z.; Liu, Z.; Heino, J.N.; Jiang, X.; Wang, J.; Tang, T.; Xie, Z. Discriminating the effects of local stressors from climatic factors and dispersal processes on multiple biodiversity dimensions of macroinvertebrate communities across subtropical drainage basins. Sci. Total Environ. 2020, 711, 134750. [Google Scholar] [CrossRef] [PubMed]
  53. Wei, X.; Liu, D.F.; Liu, H.; Li, M.; Lin, M.; Huang, Q. Analysis of socio-economic driving forces of cultivated land area and its forecast in the upstream of Han River of China. IOP Conf. Ser. Earth Environ. Sci. 2018, 189, 052037. [Google Scholar] [CrossRef]
  54. Xia, J.; Ma, X.Y.; Zou, L.; Wang, Y.L.; Jing, Z.X. Quantitative analysis of the effects of climate change and human activities on runoff in the Upper Hanjiang River basin. South-to-North Water Transf. Water Sci. Technol. 2017, 15, 1–6. (In Chinese) [Google Scholar] [CrossRef]
  55. Yang, Q.; Liu, D.F.; Meng, X.M.; Huang, Q.; Lin, M. Quantitative analysis of impacts of environment changes on runoff in upper Han River. J. Hydroelectr. Eng. 2019, 38, 73–84. (In Chinese) [Google Scholar] [CrossRef]
  56. Wei, X.; Liu, D.F.; Luan, J.K.; Li, M.; Lin, M. Analysis of the Changing Laws and Influencing Factors of Social and Economic Indicators in the Upper Reach of Han River Basin. Univers. J. Geosci. 2018, 6, 55–64. [Google Scholar] [CrossRef] [Green Version]
  57. Jia, B.; Zhou, J.; Zhang, Y.; Tian, M.; Ding, X. System Dynamics Model for the Coevolution of Coupled Water Supply-Power Generation-Environment Systems: Upper Yangtze River Basin, China. J. Hydrol. 2020, 593, 125892. [Google Scholar] [CrossRef]
  58. Li, B.; Sivapalan, M.; Xu, X. An Urban Sociohydrologic Model for Exploration of Beijing’s Water Sustainability Challenges and Solution Spaces. Water Resour. Res. 2019, 55, 5918–5940. [Google Scholar] [CrossRef]
  59. The Implementation Opinions of the General Office of Shaanxi Provincial People’s Government on Strictly Protecting Cultivated Land and Resolutely Stopping the “Non-Agricultural” Behavior of Cultivated Land. Available online: http://zrzyt.shaanxi.gov.cn/info/1150/51238.htm (accessed on 12 August 2021).
  60. Regulations for the Implementation of the Land Administration Law of the People’s Republic of China. Available online: http://zrzyt.shaanxi.gov.cn/info/1148/54120.htm (accessed on 12 August 2021).
  61. Notice of Shaanxi Provincial People’s Government on Printing and Distributing Population Development Plan (2016–2030). Available online: http://www.shaanxi.gov.cn/zfxxgk/zfgb/2018_3966/d5q_3971/201803/t20180320_1638241.html (accessed on 12 August 2021).
  62. Wei, J.; Zheng, K.; Zhang, F.; Fang, C.; Zhou, Y.; Li, F.; Jian, S. Migration of Rural Residents to Urban Areas Drives Grassland Vegetation Increase in China’s Loess Plateau. Sustainability 2019, 11, 6764. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Map of the upper reach of Han River Basin.
Figure 1. Map of the upper reach of Han River Basin.
Water 13 02458 g001
Figure 2. Causal loop diagram of socio-hydrological system: (a) the Upper section; (b) the Lower section.
Figure 2. Causal loop diagram of socio-hydrological system: (a) the Upper section; (b) the Lower section.
Water 13 02458 g002
Figure 3. Simulated and observed values of variables. (a) Population in the upper section, (b) population in the lower section, (c) urban population in the upper section, (d) urban population in the lower section, (e) cultivable land area in the upper section, (f) cultivable land area in the lower section, (g) number of cows in the upper section, (h) number of cows in the lower section, (i) number of sheep in the upper section, (j) number of sheep in the upper section, (k) number of pigs in the upper section, (l) number of pigs in the lower section, (m) industrial added value in the upper section, (n) industrial added value in the upper section, (o) Shiquan runoff in the upper section, (p) Baihe runoff in the lower section, (q) water quota for industrial added value in the upper section, (r) water quota for industrial added value in the lower section.
Figure 3. Simulated and observed values of variables. (a) Population in the upper section, (b) population in the lower section, (c) urban population in the upper section, (d) urban population in the lower section, (e) cultivable land area in the upper section, (f) cultivable land area in the lower section, (g) number of cows in the upper section, (h) number of cows in the lower section, (i) number of sheep in the upper section, (j) number of sheep in the upper section, (k) number of pigs in the upper section, (l) number of pigs in the lower section, (m) industrial added value in the upper section, (n) industrial added value in the upper section, (o) Shiquan runoff in the upper section, (p) Baihe runoff in the lower section, (q) water quota for industrial added value in the upper section, (r) water quota for industrial added value in the lower section.
Water 13 02458 g003aWater 13 02458 g003b
Figure 4. The future evolution of the main variables in the upper section. (a) Population in the upper section, (b) urban population in the upper section, (c) cultivable land area in the upper section, (d) industrial added value in the upper section, (e) water consumption in the upper section, (f) agricultural water consumption in the upper section, (g) tertiary industry water consumption in the upper section, (h) domestic water consumption in the upper section, (i) industrial water consumption in the upper section, (j) environmental water consumption in the upper section.
Figure 4. The future evolution of the main variables in the upper section. (a) Population in the upper section, (b) urban population in the upper section, (c) cultivable land area in the upper section, (d) industrial added value in the upper section, (e) water consumption in the upper section, (f) agricultural water consumption in the upper section, (g) tertiary industry water consumption in the upper section, (h) domestic water consumption in the upper section, (i) industrial water consumption in the upper section, (j) environmental water consumption in the upper section.
Water 13 02458 g004aWater 13 02458 g004b
Figure 5. The future evolution of the main variables in the lower section. (a) Population in the lower section, (b) urban population in the lower section, (c) cultivable land area in the lower section, (d) industrial added value in the lower section, (e) water consumption in the lower section, (f) agricultural water consumption in the lower section, (g) tertiary industry water consumption in the lower section, (h) domestic water consumption in the lower section, (i) industrial water consumption in the lower section, (j) environmental water consumption in the lower section.
Figure 5. The future evolution of the main variables in the lower section. (a) Population in the lower section, (b) urban population in the lower section, (c) cultivable land area in the lower section, (d) industrial added value in the lower section, (e) water consumption in the lower section, (f) agricultural water consumption in the lower section, (g) tertiary industry water consumption in the lower section, (h) domestic water consumption in the lower section, (i) industrial water consumption in the lower section, (j) environmental water consumption in the lower section.
Water 13 02458 g005aWater 13 02458 g005b
Table 1. The value of water quota.
Table 1. The value of water quota.
Water QuotaVariableUnitValue
Water quota for rural population q R P m3/per29.20
Water quota for urban population q U P m3/per47.45
Water quota for hospital beds q H B m3/bed60.23
Water quota for tertiary industry employees q T P m3/per36.50
Water quota for cows q C O W m3/per29.20
Water quota for sheep q S H E E P m3/per3.65
Water quota for pigs q P I G m3/per12.78
Table 2. The initial values of the main variables and calibrated parameters of the evolution model.
Table 2. The initial values of the main variables and calibrated parameters of the evolution model.
VariablesVariable DescriptionInitial Value in Upper SectionInitial Value in Lower
Section
Unit
P Population2,885,7473,211,930per
N P I G Number of pigs1,100,0001,520,090per
N S H E E P Number of sheep95,068345,078per
N C O W Number of cows277,192278,240per
V I A Industrial added value23,77465,753104 yuan
N H B Number of hospital beds11,0164500bed
A G F Forest and grass area20,00018,000km2
A W Water area16.2112.74km2
A U Urban area7525km2
α G F Empirical coefficient of forest and grass evaporation0.70.6
α E Empirical coefficient of environmental water consumption0.0080.008
α U Empirical coefficient of urban evaporation0.20.2
α W Empirical coefficient of water surface evaporation11
A C Cultivable land area219,871297,268km2
Δ S Water storage130120m3
E P Potential evaporation705.1813.3mm
α I Proportion of irrigated area0.600.70
V I A , m a x Maximum industrial added value1.83 × 1072.08 × 107104 yuan
r U P Urbanization rate0.15680.0904
Table 3. The parameter settings of the evolution model.
Table 3. The parameter settings of the evolution model.
VariablesVariable NameUpper SectionLower Section
r P Change rate of population0.00260.0036
r R U Change rate of urbanization rate0.00320.0041
r C Change rate of cultivated land−0.01−0.014
r C O W Change rate of cows−0.0046−0.003
r V I A Change rate of industrial added value0.170.22
r q V I A Change rate of the quota for industrial added value −0.128−0.149
Table 4. The evaluation results in the calibration period (1990–2009).
Table 4. The evaluation results in the calibration period (1990–2009).
Upper SectionLower Section
VariablesMAPE (%)NSER2MAPE (%)NSER2
P 0.400.920.970.630.880.95
P U 2.250.961.006.630.900.99
A C 5.700.640.945.760.960.93
N C O W 4.980.230.759.580.130.59
N S H E E P 27.410.770.8419.720.770.84
N P I G 12.100.550.6213.560.620.60
V I A 24.920.980.9917.991.000.99
q V I A 9.380.980.9914.900.980.99
R S h i q u a n 16.750.690.91
R B a i h e 6.820.910.96
Table 5. The evaluation results in the validation period (2010–2018).
Table 5. The evaluation results in the validation period (2010–2018).
Upper SectionLower Section
VariablesMAPE
(%)
NSER2MAPE
(%)
NSER2
P 0.890.870.800.430.950.11
P U 4.050.931.0011.640.860.97
A C 2.930.84−0.716.860.830.31
N C O W 3.680.290.586.990.320.18
N S H E E P 6.880.88−0.0413.570.75−0.21
N P I G 20.760.38−0.7931.640.67−0.19
V I A 32.490.540.9732.680.470.97
q V I A 7.911.000.9818.521.000.98
R S h i q u a n 25.870.400.93
R B a i h e 20.13−0.140.40
Table 6. Sensitive parameter values in the upper section.
Table 6. Sensitive parameter values in the upper section.
Scenario r U R U r U C r U P r U V I A α U T P
Natural continuation scenario0.0040.004−0.0030.110.160
Environmental protection scenario0.0060.0010.0030.060.162
Industrial adjustment scenario0.008−0.0030.0030.090.170
Economic development scenario0.010−0.0060.0060.120.165
Table 7. Sensitive parameter values in the lower section.
Table 7. Sensitive parameter values in the lower section.
Scenario r L R U r L C r L P r L V I A α L T P
Natural continuation scenario0.006−0.008−0.0060.110.160
Environmental protection scenario0.0080.0010.0030.080.162
Industrial adjustment scenario0.010−0.0030.0030.100.170
Economic development scenario0.012−0.0060.0060.160.165
Table 8. The water consumption in 2045 after changing water quota for irrigation.
Table 8. The water consumption in 2045 after changing water quota for irrigation.
ConditionsUpper Section (m3)Lower Section (m3)
0.95 × q U I 20.78 × 10810.94 × 108
0.90 × q U I 19.96 × 10810.55 × 108
0.85 × q U I 19.13 × 10810.16 × 108
0.80 × q U I 18.30 × 1089.77 × 108
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, X.; Liu, D.; Wei, X.; Ma, L.; Lin, M.; Meng, X.; Huang, Q. Analysis of Socio-Hydrological Evolution Processes Based on a Modeling Approach in the Upper Reaches of the Han River in China. Water 2021, 13, 2458. https://doi.org/10.3390/w13182458

AMA Style

Zhao X, Liu D, Wei X, Ma L, Lin M, Meng X, Huang Q. Analysis of Socio-Hydrological Evolution Processes Based on a Modeling Approach in the Upper Reaches of the Han River in China. Water. 2021; 13(18):2458. https://doi.org/10.3390/w13182458

Chicago/Turabian Style

Zhao, Xiaoyu, Dengfeng Liu, Xiu Wei, Lan Ma, Mu Lin, Xianmeng Meng, and Qiang Huang. 2021. "Analysis of Socio-Hydrological Evolution Processes Based on a Modeling Approach in the Upper Reaches of the Han River in China" Water 13, no. 18: 2458. https://doi.org/10.3390/w13182458

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