The Impact of the Changes in Climate, Land Use and Direct Human Activity on the Discharge in Qingshui River Basin, China

Climate, land use and human activity have an impact on the Qingshui River in Chongli County. The Soil and Water Assessment Tool (SWAT) was used to separately analyze the contributions of climate, land use and direct human activity on the discharge variations. The results indicated that human activity had been the dominant factor for the discharge decrease, while climate and land use change had a positive inﬂuence on the discharge increase. The contributions of these three factors were − 56.24%, 38.59% and 5.17%, respectively. Moreover, on the seasonal scale, the impact of those factors was consistent with their impact on the annual scale. Human activity was the main factor for discharge decrease in the summer, the contribution accounting for − 77.13%. Due to the over-extraction of groundwater for irrigation and use in the mining industry, the discharge showed a decreasing tendency, which has the potential to place stress on sustainable water use in the future. The result of the study may contribute to the optimization of water resource allocation and management.


Introduction
Understanding the watershed hydrological process and its main driving factor is essential for water management. Due to climate change and human activity, watershed runoff has changed significantly in recent years [1][2][3][4][5][6][7]. Through the redistribution of precipitation and changes in temperature, climate variation may have a direct influence on runoff [8][9][10]. The human activities factor concerns not only land cover but also direct human activity [11]. As regards the impact of human activities on runoff, many studies simply focus on land cover and land change [4,12]. Generally, the decreasing use of land for cultivation and the expansion of urbanization are the main reasons for land use change. Land use cover and change impact the characteristics of the underlying surface in basins, which can affect the water cycle [13]. With the development of urbanization and population growth, increasing numbers of hydropower construction and agriculture irrigation projects have been built in some areas. These kinds of water projects have a direct impact on the amount of runoff. Therefore, it is essential to consider both land use change and direct human activity in the human activities factor [11,14].
According to a water resource vulnerability assessment conducted for the city of Zhangjiakou [15], water resources in Chongli County were vulnerable, especially for the hydrological and eco-environmental system. Groundwater in the study area was extracted for irrigation annually, taking up 54% of the whole water resource. Apart from that, groundwater is also extracted for industrial purposes. In 2022, skiing games will be Water 2021, 13, 3147 2 of 13 carried out in the study area. An event of this type threatens to seriously exacerbate water shortages [16,17]. Considering the importance of water resources within the study area, it is crucial to implement water management. The analysis of discharge change indicates valuable evidence-based information of watershed response to past changes in climate, land use change and human activity separately, which is important for water resources management [18].
The approaches to be applied in this field can be classified into three parts: a statistical analysis, elasticity-based methods and hydrological simulation. Through analyzing the linear relationship between meteorological data and runoff [4,19], the statistical method does not consider the physical processes of the hydrological cycle [20,21]. The Budyko framework is widely used as an elasticity-based method for considering the changes in precipitation (P), potential evapotranspiration (E 0 ) and the n (a parameter presenting catchment characteristics) on runoff (Q) [22,23]. However, it neglects the difference on the spatial scale. Watershed models are powerful tools for simulating the effect of watershed processes on soil and water resources [13]. Moreover, hydrological models can simulate water cycles under different scenarios, such as climate change and human activity, on both time and spatial scales [24]. Among watershed models, the Soil and Water Assessment Tool (SWAT) can be used to simulate the impacts of alternative data items, such as changes in land use, land management practice and climate [25][26][27]. The model is highly rated for its computational efficiency and long term continuous simulations [28].
The objectives of this study are as follows: (1) to analyze the precipitation trend and then to compare the relationship between precipitation and discharge change at the Qingshui station (the whole basin outlet); (2) to separate the contributions of climate change, land use change and direct human activity on the discharge variation on annual and seasonal scales; (3) to discuss how these three factors bring a change in the runoff.

Study Area
Chongli County, with a total area of 2334 km 2 and a population of 0.126 million, is located in the northwest of Hebei Province as a transition zone between the Inner Mongolia plateau and the North China Plain. Generally, the precipitation annual average is 466.05 mm and it varies seasonally, specifically in the flood season, which is concentrated from June to September. The rainfall in this period amounts to approximately 91%. The upper reaches of Qingshui hydrological station, covering 2100 km 2 , were chosen to perform the study. The climate station (Chongli station) is located in the eastern sub-basin (Donggou) ( Figure 1). The main land uses were cultivated land, forest and grassland.
Water 2021, 13, x FOR PEER REVIEW 3 of 12 (HWSD). The data regarding water supply and demand were available from the water resources bureau in Chongli County.

The Mann-Kendall Test Method
This analysis involved the application of the nonparametric Mann-Kendall test, which detects trends in temporal sequences. The test is widely used to determine tem-

Materials
Climate data at Chongli meteorological stations from 1980 to 2016 were provided by the National Climate Center of China. These include daily precipitation, temperature, solar radiation and relative humidity data. Monthly flow data from the Qingshui station (1980-1989, 2001, 2003, 2007-2016) were gathered from the Hydrological Yearbook of China, but data were missing in some years (1990-2000, 2002, 2004-2006). The digital elevation model (DEM) was provided by the Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Sciences (http://www.gscloud.cn (accessed on 5 November 2021)), with a resolution ratio of 30 m. The land use data with 30 m resolution from 1990 to 2015 were provided by Data Center for Resources and Environmental Sciences, Chinese Academy of Science (RESDC). The soil data (1:1,000,000) was provided by National Tibetan Plateau Data Center (http://data.tpdc.ac.cn (accessed on 6 November 2021)), obtained from the Harmonized World Soil Database version 1.1 (HWSD). The data regarding water supply and demand were available from the water resources bureau in Chongli County.

The Mann-Kendall Test Method
This analysis involved the application of the nonparametric Mann-Kendall test, which detects trends in temporal sequences. The test is widely used to determine temporal trends of hydrological and meteorological data [29,30]. Statistic Z of the Mann-Kendall rank trend test is based on Equations (1) and (2).
If the sample size n is large enough (n > 10), the standard normal variable Z can be used to evaluate the trend of the time series: where x j , x i are the sequential data values, n is the length of the data set and Σ denotes the summation over all ties. Var (S) is the expected value and variance of S. This time series has increasing trends if the value of Z is greater than zero; otherwise, the time series has decreasing trends. If |Z| ≥ Z_(1 − (α/2)), the null hypothesis of no trend is rejected at the significance level α. The null hypothesis can be tested at the 5% significance level. The value of Z_(1 − (α/2)) at the 5% significance level is 1.96 [31].

The Double Mass Curve Method
The double mass curve method was used for investigating the homogeneity of the data. For this purpose, cumulative values of annual precipitation were plotted against cumulative values of annual discharge. The double mass curve is a simple, visual and practical method, widely applied to study the consistency of hydro-climatic data [32]. This method is based on the direct relationship between precipitation and discharge, in which a straight line is plotted on the cumulative values of these two parameters.
The lack of deviation from the line shows that human activities have not made noticeable contributions to the changes of stream flow and it is only influenced by precip-itation [33]. A significant deviation from the straight line plotted on cumulative values of precipitation and discharge indicates that human activities, such as land use change, irrigation or the mining industry, could change the runoff.

Soil and Water Assessment Tool (SWAT) Model
SWAT is a semi-distributed river basin scale model developed by the United States Development of Agriculture (USDA) to assess the impact of management on water supplies and nonpoint source pollution [34]. The hydrologic cycle of the SWAT model is based on the water balance equation related to land and water routing sections of the hydrological cycle [28,35]: where SW t is the final soil water content (mm), SW 0 is the initial soil water content on day i (mm), t is the time (days), R day is the amount of precipitation on day i (mm), Q surf is the amount of surface runoff on day i (mm), E a is the amount of evapotranspiration on day i (mm), W seep is the amount of water entering the vadose zone from the soil profile on day i (mm), and Q gw is the amount of return flow on day i (mm). The processes simulated by the SWAT include precipitation, runoff, infiltration, evaporation and percolation to shallow and deep aquifers. The input data include spatial data and attribute data, comprised of the Digital Elevation Model, land use, soil parameters for hydrological characteristics, daily climate attributes and monthly runoff. Based on the topography, the whole basin was divided into different sub-basins and then into different hydrologic response units (HRU) with distinct combinations of land use, soil and slope qualities. The runoff was predicted separately for each HRU and routed to obtain the total runoff based on the watershed.

Model Calibration and Validation
Parameter sensitivity analysis determines which inputs contribute the most to output variations, which parameters are most highly correlated with the output [36], and also helps simulation data to best match observation data. In addition, sensitivity analysis, model calibration and validation were performed using the algorithm for Sequential Uncertainty Fitting (SUFI-2) [37,38] and it was implemented in the SWAT-CUP [38]. In addition to automated calibration, SWAT-CUP also ranks the parameters based on their sensitivity to the flow regime [27]. The suitable ranges for Nash-Sutcliffe efficiency coefficient (NSE) (Nash and Sutcliffe, 1970) and coefficient of determination (R 2 ) (Equations (6) and (7)) were chosen to measure the matching between the simulated runoff and the observation. NSE changes from −∞ to 0 and it is acceptable when the value varies from 0 to 1, especially when it is showing a satisfactory rating up to 0.5 [39]. R 2 is a complementary statistical criterion for efficiency statistics [37]. R 2 means the correlation between the two variables. The performance of the SWAT model on the monthly scale can be considered satisfactory if NSE > 0.5, R 2 > 0.6 [14].
where Q obs i is the observed value, Q sim i is the simulated value, Q mean is the mean of the observed data, Q smean is the mean of the simulated data, n is the total number of observations. To make a comparison, there were two periods used for the whole time of the study: a benchmark and a study period. The model was simulated and calibrated by the benchmark period. Through keeping one factor stable and changing another, the simulated runoff in the study period can be used to compare with the runoff in the benchmark period [40]. The change scenarios using the SWAT simulation are as follows: S 1 : climate data in the study period and land use in the baseline; S 2 : climate data in the baseline and land use in the study period; S 3 : both climate data and land use in the study period.
The simulated runoff after applying the SWAT model was affected by the climate and land use. The formulas of the impact of the climate change and land use are as follows: where ∆Q C , ∆Q l represent the relative discharge change due to the climate and land use change under scenario 1 and scenario 2 respectively, ∆Q h means the relative discharge change due to the direct human activity, Q S1 , Q S2 and Q s3 mean the runoff from scenarios 1, 2 and 3, respectively, Q b , Q o mean the runoff in the baseline and the study area, respectively. The relative contributions of climate change, land use and the direct human activity, respectively, can be calculated as follows: where η c , η l and η h are the impact of climate change, land use change and direct human activity, respectively.

Trend Analysis of Discharge
The M-K test Z-index is −2.01 for annual runoff, which indicates that the annual runoff showed a decreasing trend, occurring from 1980 to 2016. The M-K test Z-index results for annual precipitation and temperature are 0.03 and 2.67 respectively, indicating an increasing trend. The value of cumulative precipitation and discharge during the period of 1980-2016 is shown in Figure 2. This deviation strongly suggests that, in addition to precipitation, human activity has also contributed to runoff changes since 2002. The results show that the data were not homogeneous. It is evident that runoff deviation in the years

Model Calibration and Validation
Having considered the results of statistical analysis and the limitation data, the baseline was selected from 1980 to 1989 and the study period was 2016. Furthermore, the baseline was divided into three parts: first, a one (1980) was selected to warm up the model, next, calibration was performed period of 1981 to 1985 and finally, the period of validation was from 1986 data from 1982 were removed because of a monitoring problem. The tim monthly. Figure 3 shows a close match between the simulated and observed di during the period, which indicates the validity of the model. What is more, the rainfall is consistent with the runoff in the simulation. The NSE and R calibration period were 0.65 and 0.67, respectively. Moreover, the two ind separately 0.59 and 0.72 during the validation period. In general, the monthl for the Qingshui station works well for the two periods, which means tha model can be used in this study area.
In the sensitivity analysis, eight parameters were found to be the mos the stream runoff changes (Table 1). GWQMN is the most sensitive param means that the groundwater flow to the reach is allowed only if the depth of shallow aquifer is equal to or greater than GWQMN, pointing to the impo relationship between the shallow aquifer and the main channel. What is mor and third parameters are SOIL_BD and SOIL_AWC, which are related to th water stored in the soil layer, affecting the process of water transpiration and amount of the precipitation flowing to the runoff.

Model Calibration and Validation
Having considered the results of statistical analysis and the limitation of available data, the baseline was selected from 1980 to 1989 and the study period was from 2008 to 2016. Furthermore, the baseline was divided into three parts: first, a one-year period (1980) was selected to warm up the model, next, calibration was performed during the period of 1981 to 1985 and finally, the period of validation was from 1986 to 1989. The data from 1982 were removed because of a monitoring problem. The time step was monthly. Figure 3 shows a close match between the simulated and observed discharge data during the period, which indicates the validity of the model. What is more, the trend of the rainfall is consistent with the runoff in the simulation. The NSE and R 2 during the calibration period were 0.65 and 0.67, respectively. Moreover, the two indicators were separately 0.59 and 0.72 during the validation period. In general, the monthly simulation for the Qingshui station works well for the two periods, which means that the SWAT model can be used in this study area.
In the sensitivity analysis, eight parameters were found to be the most sensitive to the stream runoff changes (Table 1). GWQMN is the most sensitive parameter, and it means that the groundwater flow to the reach is allowed only if the depth of water in the shallow aquifer is equal to or greater than GWQMN, pointing to the importance of the relationship between the shallow aquifer and the main channel. What is more, the second and third parameters are SOIL_BD and SOIL_AWC, which are related to the amount of water stored in the soil layer, affecting the process of water transpiration and altering the amount of the precipitation flowing to the runoff.   Notes: "v_": parameter should be replaced by a given value; "r_": parameter should be multiplied by (1+ given value).

Contribution of Different Factors to Runoff Change
Through model simulation, discharges in different scenarios at the annual and seasonal scale were noted and are shown in Table 2, and the contributions of three factors during the study period are shown in Figure 4. Compared with the runoff in the base period, it may be noted that climate and land use changes increase the discharge, while human activity has the opposite impact at the annual scale. Human activity is the main factor decreasing the discharge, accounting for 56.24%. Climate contributes to discharge increasing, accounting for 38.59%, and land use makes a slight impact, accounting for just 5.17%. Calculating the contribution at the seasonal scale shows that climate makes the main positive contribution to the increasing runoff in spring, autumn and winter, at 52.26%, 61.02% and 55.92%, respectively. Human activity is the dominant factor for the decreasing runoff in summer with a contribution of 77.13%. The contributions of land use are −7.12%, 11.97%, 12.44% and −16.05%, respectively, in different seasons.  Notes: "v_": parameter should be replaced by a given value; "r_": parameter should be multiplied by (1+ given value).

Contribution of Different Factors to Runoff Change
Through model simulation, discharges in different scenarios at the annual and seasonal scale were noted and are shown in Table 2, and the contributions of three factors during the study period are shown in Figure 4. Compared with the runoff in the base period, it may be noted that climate and land use changes increase the discharge, while human activity has the opposite impact at the annual scale. Human activity is the main factor decreasing the discharge, accounting for 56.24%. Climate contributes to discharge increasing, accounting for 38.59%, and land use makes a slight impact, accounting for just 5.17%. Calculating the contribution at the seasonal scale shows that climate makes the main positive contribution to the increasing runoff in spring, autumn and winter, at 52.26%, 61.02% and 55.92%, respectively. Human activity is the dominant factor for the decreasing runoff in summer with a contribution of 77.13%. The contributions of land use are −7.12%, 11.97%, 12.44% and −16.05%, respectively, in different seasons.

Limitations and Uncertainties
The results have a low representativeness as they are based on from only one measuring station in such a large catchment area a catchment. However, this catchment was chosen to assess the hydr the water balance in the location of the Olympic Games arena, hen recognize the size of the runoff and the impact of human activity o Olympic Games will be held in the east part of the basin, which m from the whole Qingshui river in the Chongli District. Besides this will improve its social utilities and economic development by mak Olympic Games in the future.
The analyzed catchment area, despite its considerable size, is sented by meteorological stations. However, this is not an isolate

Limitations and Uncertainties
The results have a low representativeness as they are based on rainfall and flow data from only one measuring station in such a large catchment area as the Qingshui River catchment. However, this catchment was chosen to assess the hydrological conditions of the water balance in the location of the Olympic Games arena, hence, it is necessary to recognize the size of the runoff and the impact of human activity on its size. The Winter Olympic Games will be held in the east part of the basin, which may make use of water from the whole Qingshui river in the Chongli District. Besides this, the Chongli District will improve its social utilities and economic development by making use of the Winter Olympic Games in the future.
The analyzed catchment area, despite its considerable size, is not properly represented by meteorological stations. However, this is not an isolated situation with the problem of insufficient spatial coverage of the hydro meteorological monitoring network. This is a typical situation in many areas (especially mountain areas). Based on the observation sequence from the one meteorological station present in this area, this article tries to ensure the proper quality of the presented data, using for this purpose an M-K analysis to examine the homogeneity of the observations. It can also be assumed that the performed calibration and verification of the results of hydrological modeling allow, to some extent, the influence of the uncertainty of the input data on the results to be limited.

Assessment of the Impacts of the Three Factors
The hydrological components under different scenarios were simulated by the SWAT model (Table 3). Under the climate scenario and combined with the M-K analysis, the precipitation and temperature both show an increasing trend ( Figure 5). What is more, according to Table 3, the precipitation and evaporation increased by 2.87% and 3.32%, respectively, while the ration about E/P decreased from 87% to 83%. The surface and lateral runoffs both show an increasing trend of 2.51 mm and 2.12 mm, respectively.

Assessment of the Impacts of the Three Factors
The hydrological components under different scenarios were simulated by th SWAT model (Table 3). Under the climate scenario and combined with the M-K analysi the precipitation and temperature both show an increasing trend ( Figure 5). What more, according to Table 3, the precipitation and evaporation increased by 2.87% an 3.32%, respectively, while the ration about E/P decreased from 87% to 83%. The surfac and lateral runoffs both show an increasing trend of 2.51 mm and 2.12 mm, respectively  The land use from 1990 to 2015 changed slightly (Table 4). Under the land us change scenario, precipitation and ET are nearly the same as the data in the baseline pe riod. In addition, the hydrological components were nearly the same as those during th base period. Human activity was the main factor for the discharge change. Owing to intensiv The land use from 1990 to 2015 changed slightly (Table 4). Under the land use change scenario, precipitation and ET are nearly the same as the data in the baseline period. In addition, the hydrological components were nearly the same as those during the base period.
Human activity was the main factor for the discharge change. Owing to intensive human activity, natural single water circulation was disrupted and transformed into a dual combination of natural and social water circulation, which clearly increased water consumption. As shown in Figure 6, from 2005 to 2016, the main use of water was for the purpose of irrigation. The structure of water supply is single, with nearly 80% of the water supply coming from ground water. It is mostly used for agricultural irrigation and in the mining industry, which are the main vehicles for economic development. purpose of irrigation. The structure of water supply is single, with nearly 80% of the water supply coming from ground water. It is mostly used for agricultural irrigation and in the mining industry, which are the main vehicles for economic development.
Because of the increasing water demand and the low-efficiency irrigation method (furrow irrigation), there is an overexploitation of the groundwater, which decreased the ground water level in agricultural and mining areas. As shown in Figure 7, according to the available monitoring groundwater data from March 2015 to June 2018, the groundwater level in the vegetable planting area showed a tendency to decline. Thus, further research regarding water resource allocation should be carried out to increase water use efficiency and decrease the strain on water resources.

Conclusions
The purpose of this paper is to quantify the contribution of climate, land use and direct human activities on discharge variations on annual and seasonal scales, as simulated by the SWAT model. It may be helpful in identifying the driving factor for runoff Because of the increasing water demand and the low-efficiency irrigation method (furrow irrigation), there is an overexploitation of the groundwater, which decreased the ground water level in agricultural and mining areas. As shown in Figure 7, according to the available monitoring groundwater data from March 2015 to June 2018, the groundwater level in the vegetable planting area showed a tendency to decline. Thus, further research regarding water resource allocation should be carried out to increase water use efficiency and decrease the strain on water resources.
Water 2021, 13, x FOR PEER REVIEW 10 of 12 purpose of irrigation. The structure of water supply is single, with nearly 80% of the water supply coming from ground water. It is mostly used for agricultural irrigation and in the mining industry, which are the main vehicles for economic development.
Because of the increasing water demand and the low-efficiency irrigation method (furrow irrigation), there is an overexploitation of the groundwater, which decreased the ground water level in agricultural and mining areas. As shown in Figure 7, according to the available monitoring groundwater data from March 2015 to June 2018, the groundwater level in the vegetable planting area showed a tendency to decline. Thus, further research regarding water resource allocation should be carried out to increase water use efficiency and decrease the strain on water resources.

Conclusions
The purpose of this paper is to quantify the contribution of climate, land use and direct human activities on discharge variations on annual and seasonal scales, as simu-

Conclusions
The purpose of this paper is to quantify the contribution of climate, land use and direct human activities on discharge variations on annual and seasonal scales, as simulated by the SWAT model. It may be helpful in identifying the driving factor for runoff change and optimizing water allocation and management throughout different seasons. The main conclusions are as follows: In the study period, the surface runoff showed a decreasing tendency. Direct human activity was the dominant factor for the discharge decrease, while climate and land use change made a positive contribution to the discharge increase. On the annual scale, the contributions of climate, land use and human activity were 38.59%, 5.17% and −56.24%, respectively. On the seasonal scale, human activity was the main factor for the reduction in the discharge in summer, accounting for −77.13%. The climate had an important influence on the discharge increase, except for during the summer. The contribution of land use was slight.
Human activity, such as the construction of irrigation, impacted the discharge decrease. The climate showed a positive contribution to the discharge variation due to increasing precipitation and temperature during the study period. The reason that land use had a slight impact on the discharge increase was due to a minor land cover change.