Next Article in Journal
Multi-Residue Screening of Pesticides in Aquaculture Waters through Ultra-High-Performance Liquid Chromatography-Q/Orbitrap Mass Spectrometry
Next Article in Special Issue
Flood Risk Evaluation in Ungauged Coastal Areas: The Case Study of Ippocampo (Southern Italy)
Previous Article in Journal
Using Steel Slag for Dissolved Phosphorus Removal: Insights from a Designed Flow-Through Laboratory Experimental Structure
Previous Article in Special Issue
New Empirical Model Using Landscape Hydric Potential Method to Estimate Median Peak Discharges in Mountain Ungauged Catchments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Response of LUCC on Runoff Generation Process in Middle Yellow River Basin: The Gushanchuan Basin

School of Water Conservancy Science and Engineering, Zhengzhou University, Zhengzhou 450001, China
*
Author to whom correspondence should be addressed.
Water 2020, 12(5), 1237; https://doi.org/10.3390/w12051237
Submission received: 20 March 2020 / Revised: 10 April 2020 / Accepted: 21 April 2020 / Published: 26 April 2020

Abstract

:
Runoff reduction in most river basins in China has become a hotpot in recent years. The Gushanchuan river, a primary tributary of the middle Yellow river, Northern China, showed a significant downward trend in the last century. Little is known regarding the relative contributions of changing environment to the observed hydrological trends and response on the runoff generation process in its watershed. On the basis of observed hydrological and meteorological data from 1965–2010, the Mann-Kendall trend test and climate elasticity method were used to distinguish the effects of climate change and human activities on runoff in the Gushanchuan basin. The results indicate that the runoff in the Gushanchuan Basin has experienced significant declines as large as 77% from 1965 to 2010, and a mutation point occurred around 1997; the contribution rate of climate change to runoff change is 12.9–15.1%, and the contribution rate of human activities to runoff change is 84.9–87.1%. Then we divided long-term data sequence into two stages around the mutation point, and analyzed runoff generation mechanisms based on land use and cover changes (LUCC). We found that the floods in the Gushanchuan Basin were still dominated by Excess-infiltration runoff, but the proportion in 1965–1997 and 1998–2010 decreased gradually (68.46% and 45.83% in turn). The proportion of Excess-storage runoff and Mixed runoff has increased, which means that the runoff is made up of more runoff components. The variation law of the LUCC indicates that the forest area increased by 49.61%, the confluence time increased by 50.42%, and the water storage capacity of the watershed increased by 30.35%.

1. Introduction

The world’s population is increasing rapidly so that the scale of the economy is expanding, and humans are over-interfering with the natural ecosystems. It has resulted in a series of ecological and environmental problems in the Loess Plateau located in the arid and semi-arid regions of the north in China. That is the main problem of soil erosion [1,2]. For a long time, water shortage and soil erosion have become the main sources that restrict the economic development in the Loess Plateau. It is known as “Frequent drought and hardship rank first among the world”, and 1/5 of the poverty counties in China are distributed in the Loess Plateau [3]. Vegetation is an important part of soil and water conservation, and it is the core of the ecological environment construction of the Loess Plateau. It is also one of the main ways to restore the ecological decline of the Loess Plateau [4,5]. Over the past 50 years, the country has successively implemented a series of ecological construction projects such as key construction of soil and water conservation, construction of the Three-North Shelterbelt Program (TNSP), protection of natural forest resources, returning farmland to forests and grasslands, and construction of silt dams on the Loess Plateau. Most of these major initiatives are based on vegetation restoration and reconstruction. All these initiatives are in order to enhance the benign development of regional water cycle and promote the balance of the socio-economic-ecological system. The melioration of soil erosion in the Loess Plateau through vegetation restoration can be divided into two aspects: (1) changing the water circulation path to prevent soil erosion, perennial vegetation can reduce surface runoff by increasing effective vegetation coverage; (2) increasing soil water content and increasing vegetation productivity.
The Loess Plateau in China is the main runoff area of the Yellow River Basin. After the 1980s, the Yellow River water volume showed a notable downward trend [6,7] and the cause of the sharp decline of runoff in the Yellow River has always been the focus of academic debate. Feng [8] pointed out that due to land use and cover changes (LUCC) on the Loess Plateau, its water yield may decline. Besides, with the large-scale vegetation restoration on the Loess Plateau, the flow process is significantly affected, and its impact depends on the vegetation type [9,10]. However, Shi [11] found that the contribution rate of climate change to the reduction of runoff in the middle reaches of the Yellow River accounted for more than 40%. The impact of changing environment on runoff is still controversial. Currently, most of the research focuses on quantitative analysis of land use change on the basin’s hydrological cycle. This study intends to explore changes in runoff process in typical basins in the middle reaches of the Yellow River from the perspective of runoff generation patterns. Runoff generation, referring to the process of generating various runoff components in the basin. It is essentially the movement and development process of water under the combined effects of various factors on the underlying surface, it is also the redistribution process of rainfall on the underlying surface (ground and aeration zone). Different underlying surface conditions have different runoff generation mechanisms, and different runoff generation mechanisms affect the development of the entire runoff process. It presenting different runoff characteristics. At present, the natural runoff generation mechanism can be divided into four types: Surface runoff (Rs), Saturated ground runoff (Rsat), Interflow (Rint) and Underground runoff (Rg) [12,13,14]. These four classes of runoff mechanisms comprehensively show three phenomena: Excess-infiltration runoff, Excess-storage runoff and Mixed runoff. The Loess Plateau has sparse rainfall, high rainfall intensity, low vegetation coverage, and the thickness of the aerated zone up to tens of meters or even hundreds of meters, which is a typical Excess infiltration runoff area [15].
The Gushanchuan Basin is located in the first sub region of Loess Hilly and gully area, with extremely serious soil erosion and fragile ecological environment. However, the hydrological cycle in this special area still have not been investigated thoroughly. This paper aims to identify the runoff process, and gains insight into the runoff generation mechanism variation to the response of LUCC, based on the quantitative analysis of the impacts of climate change and human activities on runoff in the Gushanchuan Basin. It also tries to seek the cause of the decrease in the amount of runoff in the Yellow River from the perspective of the runoff process, and provide the necessary scientific basis for the effectiveness of ecological construction of the Loess Plateau and the formulation of regional sustainable development countermeasures.

2. Materials and Methods

2.1. Study Area

Gushanchuan Basin is a primary tributary on the right bank of the middle reaches of the Yellow River in Shaanxi and Inner Mongolia province, China (110°32′24″–111°05′24″E, 39°00′00″–39°27′36″N). Originating from Urigol Township, Inner Mongolia, through the Junggar Banner and Fugu County of Shaanxi Province, and finally flows into the Yellow River near Fugu County, with a drainage area of 1272 km2 (Figure 1). Mean annual precipitation of this basin is about 420 mm, and the interannual variation of rainfall is unevenly distributed. The flood season (June-September) accounts for 80% of total precipitation. Affected by the semi-arid continental climate in the mid-temperate zone, the basin has long sunshine duration and large evapotranspiration. Short-term heavy rainfall is the main cause of runoff generation in the basin, often causing mountain torrents and severe soil erosion. The geomorphic types in the basin are mainly loess hills and gullies, and there are a few loess sand covered areas in the upper reaches. The main soil type is loessial soil, followed by chestnut soil. The soil structure is loose, with great permeability, and it is easy to form the surface runoff. Since the 1970s, a large number of soil and water conservation measures were implemented in the basin, and the land use changed significantly which has affected the characteristics of regional runoff. Therefore, detecting and analyzing the impact of changing environment of Gushanchuan basin can be important to manage water resources in this watershed.

2.2. Data Collection

For the purposed research, we collected annual and hourly data series between 1965 and 2010 of runoff, from Gaoshiya hydrological station in Gushanchuan Basin; hourly rainfall data was obtained from four gauges (Xinmiao, Xinmin, Gushan and Gaoshiya) with a time series of 46 years, from 1965 to 2010. Besides, there is no national level meteorological station in the basin, the evaporation data of Fugu meteorological station in 1965–2010 near the basin is selected, and the evapotranspiration of the basin was calculated by inverse distance weighted interpolation (IDW). All the above data are obtained from the Yellow River Hydrological Yearbook [16]. The spatial resolution of the basin’s DEM (Digital Elevation Model) data is 90 m. The land use data comes from Landsat series of remote sensing images, which are acquired by human-computer interactive interpretation based on ENVI (The Environment for Visualizing Images) and ArcGIS, the resolution is 1 km, provided by Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Sciences (http://www.gscloud.cn).

2.3. Methodology

2.3.1. Analysis of Meteo-Hydrological Changes

In order to identify meteo-hydrological trends, we first apply the widely used Mann-Kendall trend test [17] to the time series of precipitation, runoff and potential evapotranspiration and point out significant changes in Gushanchuan basin. For a time series X (X = x1, x2,…, xn) which has n samples, the Mann- Kendall statistics S and UF were given as:
  S = n 1 i = 1 n j = i + 1   s g n ( x i x j )
U F = S k E [ S k ] V a r [ S k ] ( k = 1 , 2 , ,   n )
where xi and xj represent the values in years i and j, i > j, n is the record length of the series, and sgn(x − xj) is a characterization of the function. E(SK) and Var(SK) represent the mean and variance of SK, respectively.
Then the time series order is reversed (i.e., xn, xn-1,…, x1) and the above process is repeated to compute statistical variables UB. Finally, we draw curves of UF and UB, and if the two curves have an intersection point in the confidence level, then it can be considered as a mutation point. Researchers believe that after the hydrological sequence is significantly affected by climate change or human activities, the stability of its distribution sequence will be disturbed or destroyed, showing a certain degree of phase or trend. Therefore, according to the detected runoff mutation point, the research period can be divided into base period and mutation period. The base period watershed is in a natural state, means that the runoff is only affected by climate change. However, in the mutation period, human activities may have a significant impact on the runoff change.

2.3.2. Estimating the Impact of Climate Variability and Human Activities on Runoff

For a certain watershed, the total changes (ΔR, mm) of runoff is a synthetic production of climate change (ΔRC, mm) and human activities (ΔRH, mm), it can be calculated from the difference between the measured runoff before and after the mutation point. Therefore, to quantitatively distinguish the impact of climate change and human activities on runoff, only one factor should be quantified, and the effect of the other factor can be calculated by formula (3). Based on long-term water balance assumptions and the Budyko hypothesis [18], the evapotranspiration (Ea) is a function of the dryness index (ϕ = E0/P). Then the contribution of climate change to runoff can be derived by calculating the elastic coefficient of runoff on precipitation and potential evaporation, which is expressed as:
Δ R = Δ R C + Δ R H
Δ R C = Δ R P + Δ R E = ε P R ¯ P ¯ Δ P + ε E 0 R ¯ E 0 ¯ Δ E 0
ε P + ε E 0 = 1
ε P = 1 + ϕ F ( ϕ ) 1 F ( ϕ )
where ΔRP, ΔRE are contributions of rainfall and potential evapotranspiration to runoff, respectively. ε P ,   ε E 0 are dimensionless parameters which represent elastic coefficients of runoff on precipitation and potential evapotranspiration, respectively; P ¯ , R ¯ and E 0 ¯ are multi-year average of precipitation (mm), runoff depth (mm) and potential evapotranspiration (mm) of the study basin; F(ϕ) and F′(ϕ) can be calculated by the following Budyko forms [18,19,20,21] (Table 1):

2.3.3. Runoff Generation Pattern

Due to the uneven spatial distribution of rainfall and underlying surface conditions, the runoff mechanism is changing during the rainfall process, and nine types of different runoff mechanisms are combined under specific aerated zone structures and rainfall characteristics [22] (Table 2). However, according to the impact factors, it can be divided into two kind of patterns, namely the Excess saturation runoff and the Excess infiltration runoff patterns. Generally, the impact factors of Excess saturation runoff are P, E, W0, and the Excess infiltration runoff patterns are P, E, W0, i [23], which is expressed as:
R = f s t o r a g e ( P ,   E ,   W 0 )
R = f e x c e s s ( P ,   E ,   W 0 ,   i )
where R is runoff depth generated by rainfall (mm); P is the total amount of precipitation (mm); E is watershed evapotranspiration (mm); W0 is initial water storage capacity (mm); i is rainfall intensity (mm/h).
In addition, owing to the complexity of the spatial and temporal distribution of rainfall, the above two runoff patterns are intertwined which leads to the so-called Mixed runoff. This usually occurs in semi-arid and semi-humid areas. Therefore, the runoff pattern is the embodiment of different runoff mechanism combinations of each unit in the watershed. However, for a specific watershed, its runoff process is not static and difficult to describe quantitatively with simple indicators. Table 3 shows the integrated analysis of this study along with the runoff generation pattern to deduce the runoff process qualitatively, quantitatively, and comprehensively [24], on the basis of flood generation mechanism in the Gushanchuan Basin.

2.3.4. Antecedent Precipitation Index

Antecedent precipitation index (Pa) reveal the amount of rainfall that was retained in the soil during the previous rainfall [25]. In this study, the recursive formula was used to estimate the amount of Pa in the early stage of the flood, computed as:
P a , t = k P t 1 + k 2 P t 2 + + k n ( P a , t n + P t n )
And from this, the formula for calculating the amount of rainfall in the early stage is derived:
P a , t + 1 = k ( P a , t + P t )
where Pa, t is the antecedent precipitation index (mm) on day t; n is the number of previous rain days affecting the runoff, usually takes about 15 days, Pt is the t-day rainfall (mm); Pt − 1, Pt − 2… is the rainfall (mm) on the 1 day before t day, 2 days before. The daily decrease coefficient of soil water content k is calculated by the following formula:
k = 1 E m W m
where: Em is the watershed evapotranspiration capacity (mm); Wm is the water storage capacity (mm) of the basin

2.3.5. Watershed Water Storage Capacity

The watershed water storage capacity (Wm) is derived from multiple flood events [26], which selects the rainfall floods that can make the water content of the basin saturated and generate a large flood process. Then let Wm, 0 = Im, where the loss Im = PR. The calculation process of Wm can be estimated as:
K j = 1 E m W m , j
W 0 , j = n 0 P i × k j i
W m , j = W m , j 1 + W 0 , j 1
where Kj is daily decrease coefficient of soil water content under different iterations; Pi is the rainfall for the i day before the flood (mm); i is the power; j is the number of iterations; Wm, 0 is the assumed value for the first iteration (mm). After several iterations, Wm tends to a stable value, which is the water storage capacity of the basin.

3. Results and Discussion

3.1. Analysis of Meteo-Hydrologic Changes

Over the past 50 years, the annual runoff in the Gushanchuan Basin has shown a downward trend. From 1965 to 1980, the annual runoff fluctuated greatly, and after 1980, the runoff decreased sharply, the annual average runoff was 9.9 million m3 (1965–1980), it plummeted to 4.0 million m3 (1981–2010), the latter being merely 40.3% of the former (Figure 2A). The Mann-Kendall trend test was used to analyze the precipitation, runoff, and potential evapotranspiration in the Gushanchuan Basin. The results show that the statistical ZC values of the above meteo-hydrological factors are 0.08, −0.29, and −4.53, respectively. However, only the absolute value of runoff statistical value ZC > 1.96, indicates that the runoff change in the study area has reached a significant level. According to the Figure 2B, two curves intersected in 1997, and the intersection point was laid in two 95% confidence lines, showing that the annual runoff was abruptly changed around 1997, considering the vegetation restoration on the Loess Plateau since 1997, the vegetation adaptation to flooding events in terms of physiological, morphological, and biomechanical properties, may lead to the decrease of runoff [27]. Compared with the pre-mutation period, the annual runoff and precipitation in the catchment during the post-period decreased by 2.5% and 77%, but only the annual runoff series shows a significant downward trend within the confidence level of 0.5 (Table 4). The annual potential evapotranspiration changes are too small to attract attention. Hence, significant changes in runoff are not just caused by climate change, and human activities may have played an important role during this period.
According to the principle of water balance, the change of annual runoff for many years is the product of the combined effect of precipitation and potential evapotranspiration. In order to deeply understand the reasons for the decrease of runoff in Gushanchuan basin, the spatial distribution of changes in the average annual precipitation and potential evapotranspiration are shown in Figure 3. From the base period to the mutation period, the multi-year average precipitation of each rainfall station has decreased to varying degrees, especially in the northwest and southwest regions of the basin. The largest decrease in rainfall occurs at Xinmiao Station, which is −48–−44 mm (Figure 3A). The decrease in precipitation at each rainfall station is the direct cause of the runoff reduction in Gushanchuan basin. The spatial patterns of changes in potential evapotranspiration were generally consistent, which characterized by the more decrease in the upstream region than that in the downstream region (Figure 3B).
The impact of climate change and human activities on runoff change based on Budyko hypothesis is represented in Table 5. The results manifest that the elasticity coefficients of runoff to precipitation estimated by using four functions based on the Budyko hypothesis are very close, with an average value of 2.5. The contribution rate of climate change to runoff change is between 12.9% and 15.1%, and the contribution rate of human activities to runoff change is between 84.9% and 87.1%, and the average values are 13.9% and 86.1%, respectively. This result is similar to that of Zhan [28], which investigated the relative contributions of human activity and climate change to the hydrological response of the Weihe River Basin, the largest tributary of the Yellow River. However, all the methods adopted by them belong to the category of statistical methods and lack of physical mechanism. Wang [29] used data from 1958–2008 to explore the factors influencing the reduction of runoff in Huangfuchuan Basin, based on the methods of accumulative anomaly analysis. The results show that the contribution rates of human activities to the runoff change in Huangfuchuan Basin are 63.57% and 83.19%, during the period of 1980–1997 and 1998–2008, respectively. Overall, the analysis seems to suggest that the human activities play a critical role in the reduction of runoff in the Yellow River basin.

3.2. Characteristics of LUCC (Land Use and Cover Changes)

In the past 50 years, many soil and water conservation measures such as afforestation, check dams and terracing have been implemented. As of 2010, a total of 49 small and medium check dams were built, with a dam area of 8.9 km2, accounting for 0.70%. The scale of terracing has gradually increased, with a terracing ratio of 2.97%, and the number of reservoirs has also increased [30]. All these human activities greatly altered the land use and land cover, and further affected the land surface hydrological processes [31,32]. Then we analyzed the relationship of land use conversion based on three periods of remote sensing data from 1976, 1996, and 2006. Showing that the proportion of grassland area in the Gushanchuan Basin is the largest (60.89–62.92%), followed by farmland (30.59–32.52%), and the urban construction land area is the smallest (0.35–0.55%). Between 1976 and 2006, the area of farmland decreased year by year, and the area of forest, grassland and urban construction land increased in different extents, and the water area remained basically stable. In addition, farmland has the greatest change in area among all land use types. The total cultivated area of 24.52 km2 has been converted to other land use types. The main conversion form is conversion of farmland into forest land, and grassland into forest land. The forest area is increased by 18.22 km2, the grass area increased by 5.39 km2, demonstrating that the effect of the policy “reverting farmland to forests, returning farmland to grassland” was significant (Figure 4, Table 6).

3.3. Correlation of Rainfall-Runoff

Based on the mutation point of the Gushanchuan Basin, the time sequence data was divided into two stages for the research, and in accordance with the principle that the flood peak flow is greater than 100 m3/s, the first stage (1965–1997) collected 111 flood events, and the second stage (1998–2010) collected 24 times, with a total of 135 rainfall and flood data. Then we analyzed the rainfall-runoff characteristics of each flood event. Due to the influence of climate and land use changes, the rainfall-runoff relationship is different at each stage. The rainfall-runoff correlation coefficients in two stages are 0.48 and 0.65, respectively. In comparison, the rainfall of flood events in has decreased significantly associated with the second stage (Figure 5A). The correlation between rainfall intensity and runoff shows a slight downward trend, namely 0.56 and 0.51. For the flood events of the same magnitude of rainfall intensity in the two stages, the runoff is greater in the second stage (Figure 5B). Different from the humid area, the correlation coefficient of rainfall and runoff in this basin is relatively low, which can be analyzed from two aspects: (1) The effect of land cover changes. The grassland in Gushanchuan watershed can alleviate the impact of raindrops on the surface soil, increase the roughness of the surface soil, and reduce the impact of rainfall on runoff generation. With the increase of rainfall intensity, the rainfall that can penetrate the grass layer increases, namely, the infiltration capacity of the surface soil increases, so the runoff yield increases apparently. (2) The influence of antecedent precipitation index (Pa). The amount of Pa in the soil not only affects how much water the surface soil can store at the beginning of rainfall, but also affects the infiltration rate of soil water. Therefore, the Pa is also a very important factor affecting runoff. The lower the soil water content, the later the runoff generation begins [33].
The rainfall duration of the Gushanchuan Basin has increased, and most of the rainfall lasted between 20–40 h. The mean rainfall at each stage were 47.9 and 38.0 mm. Based on the criteria of meteorological classification, the flood events are divided into four levels: light rain (P < 10 mm), moderate rain (10–25 mm), heavy rain (25–50 mm) and rainstorm (P > 50 mm). Among the whole flood events, there were 21 light rains, 45 moderate rains, 39 heavy rains and 30 rainstorms. Therefore, the Gushanchuan Basin mainly suffered from moderate rains and heavy rains. The flood duration showed a slightly increase trend, as well as the peak flow showed a significant decrease over the years (Table 7). Under the influence of land use and cover changes, a similar phenomenon also appeared in the Chabagou basin, a tributary of the Yellow River [34].

3.4. Diagnostic Analysis of the Runoff Generation Process

We collected the characteristics of three typical floods in the Gushanchuan Basin. In the 19680821 flood (Figure 6A), the total rainfall was 31.36 mm, the rainfall lasted for 12 h, the average rainfall intensity was 2.6 mm/h, and the maximum rainfall intensity was 8.7 mm/h. The total rainfall during this rainfall process was small and short in duration, but the short-term rainfall intensity was heavy. According to the law of runoff generation, the main components of the runoff of this flood were surface runoff and interflow. The antecedent precipitation index (Pa) was only 2.3 mm, with dry and compact soil, combined with the Table 3, indicated that the flood is mainly dominated by Excess infiltration runoff. During the 15 days before the 19950728 flood (Figure 6B), light rain continued to occur and the surface soil layer was relatively humid, the Pa was calculated to be 23.2 mm in early stage. The rainfall intensity of this flood was not great, but it was more evenly distributed in space, and the runoff coefficient was 0.42. It can be considered that the surface layer of the basin is full, and the type of runoff generation mechanism was Rs + Rint + Rg, namely, Excess saturation runoff pattern; The 19940802 flood (Figure 6C) was a multi-peak flood. There was almost no rain in the first half month of the basin, which confirms that the runoff generation pattern in the early phase was mainly Excess infiltration runoff. With continuous rainfall leading to an increase in soil water content, the runoff generation type has changed from a single surface runoff mechanism (Rs) to Rs + Rint + Rg, and it is comprehensively judged that this flood is a Mixed runoff pattern. We can see that the rainfall and Pa in the 19680821 flood was small, while the rainfall and Pa in the 19950728 flood were both large. Although the rainfall in the Mixed runoff flood was larger, reaching 72.53 mm, the Pa in the early stage was extremely small, which leads to the Excess infiltration runoff pattern in the early stage of the basin. Hence, the antecedent precipitation index in the early stage is also a very important factor for basin runoff generation process.
From 1965 to 2010, the Gushanchuan Basin was dominated by Excess infiltration runoff, accounting for 64.44% of the flood events, and the Mixed runoff accounted for 32.60%, only 4 Excess saturation runoffs happened during this period. In these two stages, the floods with Excess infiltration runoff accounted for 68.46% and 45.83% of the total floods, respectively, as well as the proportion of Excess infiltration runoff decreased significantly. The proportion of Mixed runoff increased apparently by the stage, namely 30.63% and 41.67%. Nevertheless, the various stages are still dominated by Excess infiltration runoff, but the proportion of its occurrence is decreasing, while the proportion of another two kinds of runoff generation patterns is increasing (Table 8).
The area of forest in the Gushanchuan Basin increases rapidly, and vegetation can increase rainfall infiltration by improving soil organic matter and porosity. Therefore, land use cover change is an important reason for the change of rainfall runoff [35], which affects the process of runoff generation by affecting the outflow capacity of groundwater, infiltration capacity and water storage capacity [36]. The annual minimum daily flow and the minimum daily flow in flood season can reflect the outflow capacity of groundwater, which has increased gradually in the Gushanchuan Basin, especially in 1998–2010, the groundwater outflow capacity increased significantly, resulting in increased confluence time and reduced peak flow. The area of forest increased by 49.61%, the discharge capacity of groundwater increased by 44.04%, the average precipitation before generate flow increased by 90.72%, and the confluence time increased by 50.42% (Table 9).
Generally, with the increase of rainfall amount and intensity, runoff coefficient will gradually increase. The process of infiltration is also an important reason for the change of the runoff generation mechanism [37]. The stable infiltration rate of forest and grassland is greater than that of farmland, and the vegetation affects soil structure, chemical composition and bulk density. Hence, changes in vegetation play an important role in the redistribution of rainfall. Normally, deep root vegetation can alter soil properties to accommodate its water requirements, which may change the soil channels, increase soil porosity and permeability [38,39]. Studies have shown that the water storage capacity of forest land is higher than that of other land-use types, and the increase of water storage capacity will lead to an increase in underground runoff, runoff coefficient and confluence time, a decrease in peak flow [40]. The forest area in Gushanchuan Basin increased by 49.61%, the catchment water storage capacity increased by 30.35% (Table 9), and other studies have reached similar conclusions [41].

4. Conclusions

In this research, we selected Gushanchuan basin as the research object, to detect annual runoff trend and analyze the changes of the runoff generation process based on the flood events data from 1965 to 2010. Showing that the Gushanchuan Basin has experienced runoff decline as large as 77% in the last century, which is mainly caused by human activities. The runoff generation pattern of the basin is dominated by Excess infiltration runoff, but it has also changed to some extent, which is manifested in the increase of runoff components, the proportion of Excess infiltration runoff decreased, and the proportion of Excess saturation runoff and Mixed runoff increased stage by stage. The variety of land use conditions by human activities is the main factor affecting the change of the runoff generation process of the basin. Land use and cover change, as important factors that cause runoff change, especially in the flood season, cannot be ignored. In the present paper, we only take the Gushanchuan Basin as an example. How much of the land use type of the Loess Plateau has changed? Whether its runoff generation process has co-varied? It can also be expanded in the later period in order to further understand the law of the runoff generation of the Yellow River Basin. In addition, this study considers the effects of various human activities as a whole, but its comprehensive effects are extremely sophisticated. We also expect to further refine the research on the response of runoff changes to different human activities in future research.

Author Contributions

All authors have read and agree to the published version of the manuscript. Conceptualization, C.H. and S.J.; methodology, L.Z.; software, Q.W.; validation, C.H., Q.W. and L.Z.; formal analysis, L.Z.; investigation, L.Z.; resources, C.H.; data curation, L.Z.; writing—original draft preparation, L.Z.; writing—review and editing, C.H.; supervision, S.-e-h.S.; project administration, C.H.; funding acquisition, C.H.

Funding

This research was funded by [National Key Research Priorities Program of China], grant number [2016YFC04024023], [National Natural Science Foundation of China], grant number [51979250], [Key projects of National Natural Science Foundation of China], grant number [51739009], [National Natural Science Foundation of China] grant number [31700370].

Acknowledgments

We would like to thank the potential reviewer very much for their valuable comments and suggestions. We also thank my other colleagues’ valuable comments and suggestions that have helped improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhao, G.; Mu, X.; Wen, Z.; Wang, F.; Gao, P. Soil erosion, conservation, and eco-environment changes in the loess plateau of china. J. Land Degrad. Dev. 2013, 24, 499–510. [Google Scholar] [CrossRef]
  2. Zhang, B.; He, C.; Burnham, M.; Zhang, L. Evaluating the coupling effects of climate aridity and vegetation restoration on soil erosion over the loess plateau in china. J. Sci. Total Environ. 2016, 539, 436–449. [Google Scholar] [CrossRef] [PubMed]
  3. Wei, J.; Zhou, J.; Tian, J.; He, X.; Tang, K. Decoupling soil erosion and human activities on the Chinese loess plateau in the 20th century. J. Catena 2006, 68, 10–15. [Google Scholar] [CrossRef]
  4. Feng, X.; Fu, B.; Piao, S.; Wang, S.; Ciais, P.; Zeng, Z. Revegetation in china’s loess plateau is approaching sustainable water resource limits. J. Nat. Clim. Chang. 2016, 6, 1019. [Google Scholar] [CrossRef]
  5. Zhou, J.; Fu, B.; Gao, G.; Lü, Y.; Liu, Y.; Lü, N.; Wang, S. Effects of precipitation and restoration vegetation on soil erosion in a semi-arid environment in the loess plateau, china. J. Catena 2016, 137, 1–11. [Google Scholar] [CrossRef]
  6. Mu, X.; Zhang, X.; Shao, H.; Gao, P.; Wang, F.; Jiao, J.; Zhu, J. Dynamic changes of sediment discharge and the influencing factors in the yellow river, china, for the recent 90 years. J. Clean-Soil Air Water 2012, 40, 303–309. [Google Scholar] [CrossRef]
  7. Wang, S.; Fu, B.J.; Piao, S.; Lü, Y.; Ciais, P.; Feng, X.; Wang, Y. Reduced sediment transport in the yellow river due to anthropogenic changes. J. Nat. Geosci. 2015, 9, 38. [Google Scholar] [CrossRef]
  8. Feng, X.; Sun, G.; Fu, B.; Su, C.; Liu, Y.; Lamparski, H. Regional effects of vegetation restoration on water yield across the loess plateau, china. J. Hydrol. Earth Syst. Sci. Discuss. 2012, 16, 2617–2628. [Google Scholar] [CrossRef] [Green Version]
  9. Pasquino, V.; Gualtieri, P.; Doria, G. On evaluating flow resistance of rigid vegetation using classic hydraulic roughness at high submergence levels: An experimental work. In Hydrodynamic and Mass Transport at Freshwater Aquatic Interfaces; GeoPlanet: Earth and Planetary Sciences; Rowiński, P., Marion, A., Eds.; Springer: Cham, Switzerland, 2016; pp. 269–277. [Google Scholar]
  10. Aberle, J.; Järvelä, J. Flow resistance of emergent rigid and flexible floodplain vegetation. J. Hydraul. Res. 2013, 51, 33–45. [Google Scholar] [CrossRef]
  11. Shi, C.; Zhou, Y.; Fan, X.; Shao, W. A study on the annual runoff change and its relationship with water and soil conservation practices and climate change in the middle yellow river basin. J. Catena 2013, 100, 31–41. [Google Scholar] [CrossRef]
  12. Horton, R. The role of infiltration in the hydrologic cycle. J. Eos Trans. Am. Geophys. Union 1933, 14, 446. [Google Scholar] [CrossRef]
  13. Hewlett, J.; Hibbert, A. Moisture and energy conditions within a sloping soil mass during drainage. J. Geophys. Res. 1963, 68, 1081–1087. [Google Scholar] [CrossRef] [Green Version]
  14. Dunne, T.; Black, R. Partial area contributions to storm runoff in a small new England watershed. J. Water Resour. Res. 1970, 6, 1296–1311. [Google Scholar] [CrossRef] [Green Version]
  15. Liu, D.; Tian, F.; Hu, H.; Hu, H. The role of run-on for overland flow and the characteristics of runoff generation in the loess plateau, china. Hydrol. Sci. J. 2012, 57, 1107–1117. [Google Scholar] [CrossRef]
  16. Hydrological Yearbook of People’s Republic of China-Hydrological Data of Yellow River Basin; Hydrological Bureau of the Ministry of Water Resources of the People’s Republic of China: Zhengzhou, China, 2013; pp. 1965–2010.
  17. Mann, H.B. Nonparametric tests against trend. J. Econom. J. Econom. Soc. 1945, 13, 245–259. [Google Scholar] [CrossRef]
  18. Budyko, M.I.; David, H. Climate and life; Academic Press: New York, NY, USA, 1974; Volume 508, pp. 23–24. [Google Scholar]
  19. Pike, J. The estimation of annual run-off from meteorological data in a tropical climate. J. Hydrol. 1964, 2, 116–123. [Google Scholar] [CrossRef]
  20. Fu, B. On the calculation of the evaporation from land surface. Chin. J. Atmosph. Sci. 1981, 5, 23–31. [Google Scholar]
  21. Zhang, L.; Dawes, W.; Walker, G. Response of mean annual evapotranspiration to vegetation changes at catchment scale. Water Resour. Res. 2001, 37, 701–708. [Google Scholar] [CrossRef]
  22. Dunne, T.; Richard, D. An experimental investigation of runoff production in permeable soils. J. Water Resour. Res. 1970, 6, 478–490. [Google Scholar] [CrossRef]
  23. Rui, X.; Gong, X.; Zhang, C. Formation and calculation of watershed runoff yield. J. Hydroelectr. Eng. 2009, 28, 146–156. [Google Scholar]
  24. Barthold, F.; Woods, R. Stormflow generation: A meta-analysis of field evidence from small, forested catchments. J. Water Resour. Res. 2015, 51, 3730–3753. [Google Scholar] [CrossRef]
  25. Kohler, M.A.; Linsley, R.K. Predicting the Runoff from Storm Rainfall; US Department of Commerce, Weather Bureau: Washington, DC, USA, 1951; p. 34.
  26. Hu, C.; Li, S.; Zhang, W.; Chen, X.; Wang, J. Establishing and Testing Runoff Model Based on Runoff Coefficient. J. China Hydrol. 2016, 36, 8–13. [Google Scholar]
  27. Pasquino, V.; Saulino, L.; Pelosi, A.; Allevato, E.; Rita, A.; Todaro, L.; Saracino, A.; Chirico, G. Hydrodynamic behaviour of European black poplar (Populus nigra L.) under coppice management along Mediterranean river ecosystems. J. River Res. Appl. 2018, 34, 586–594. [Google Scholar] [CrossRef]
  28. Zhan, C.; Jiang, S.; Sun, F. Quantitative contribution of climate change and human activities to runoff changes in the Wei River basin, China. Hydrol. Earth Syst. Sci. 2014, 18, 8. [Google Scholar] [CrossRef] [Green Version]
  29. Wang, S.; Yan, Y.; Yan, M.; Zhao, X. Quantitative estimation of the impact of precipitation and human activities on runoff change of the Huangfuchuan River Basin. J. Geogr. Sci. 2012, 22, 906–918. [Google Scholar] [CrossRef]
  30. Liu, X. The Cause of the Sharp Decline of Water and Sediment in the Yellow River in Recent Years; Science Press of China: Beijing, China, 2016; pp. 45–69. [Google Scholar]
  31. Meerkerk, A.; Wesemael, B.; Bellin, N. Application of connectivity theory to model the impact of terrace failure on runoff in semi-arid catchments. J. Hydrol. Process. 2009, 23, 2792–2803. [Google Scholar] [CrossRef]
  32. Li, E.; Mu, X.; Zhao, G.; Gao, P.; Sun, W. Effects of check dams on runoff and sediment load in a semi-arid river basin of the Yellow River. Stoch. Environ. Res. Risk Assess. 2016, 31, 1791–1803. [Google Scholar] [CrossRef]
  33. Lee, K.; Huang, J. Runoff simulation considering time-varying partial contributing area based on current precipitation index. J. Hydrol. 2013, 486, 443–454. [Google Scholar] [CrossRef]
  34. Mo, S.; Gong, Y.; Li, J. Impact of human activities on flood process in Chabagou basin of the Loess Plateau, China. J. Basic Sci. Eng. 2019, 27, 492–508. [Google Scholar]
  35. Zuo, D.; Xu, Z.; Yao, W.; Jin, S.; Xiao, P.; Ran, D. Assessing the effects of changes in land use and climate on runoff and sediment yields from a watershed in the loess plateau of china. J. Sci. Total Environ. 2016, 544, 238–250. [Google Scholar] [CrossRef]
  36. Bergkamp, G. A hierarchical view of the interactions of runoff and infiltration with vegetation and microtopography in semiarid shrublands. J. Catena 1998, 33, 201–220. [Google Scholar] [CrossRef]
  37. Guo, Y.; Liu, S.; Baetz, B. Probabilistic rainfall-runoff transformation considering both infiltration and saturation excess runoff generation process. J. Water Resour. Res. 2012, 48, 6513. [Google Scholar]
  38. Kou, M.; Garcia, P.; Hu, S.; Jiao, J. The effect of Robinia pseudoacacia afforestation on soil and vegetation properties in the loess plateau (China): A chronosequence approach. J. For. Ecol. Manag. 2016, 375, 146–158. [Google Scholar] [CrossRef]
  39. Jian, S.; Zhao, C.; Fang, S.; Yu, K.; Peng, S. Characteristics of rainfall interception by Caragana korshinskii and Hippophae rhamnoides in Loess Plateau of Northwest China. J. Appl. Ecol. 2012, 23, 2383–2389. [Google Scholar]
  40. Li, N. Study on the Mechanism of Runoff Production and Confluence in the Loess Plateau under the Change of Underlying Surface. Ph.D. Thesis, Zhengzhou University, Zhengzhou, China, 2018. [Google Scholar]
  41. Mascha, J.; Nadine, W.; Christian, P.; Matthias, S.; Christoph, L.; Frank, M.; Thomas, F. Nutrient release from decomposing leaf litter of temperate deciduous forest trees along a gradient of increasing tree species diversity. J. Soil Biol. Biochem. 2009, 41, 2122–2130. [Google Scholar]
Figure 1. Map of the Gushanchuan basin. (A) Inset map showing the location of the Yellow River Basin in China; (B) Location of the Gushanchuan watershed along the Yellow River basin; (C) The control catchment of Gushanchuan basin and distribution of rainfall-runoff gauge stations.
Figure 1. Map of the Gushanchuan basin. (A) Inset map showing the location of the Yellow River Basin in China; (B) Location of the Gushanchuan watershed along the Yellow River basin; (C) The control catchment of Gushanchuan basin and distribution of rainfall-runoff gauge stations.
Water 12 01237 g001
Figure 2. Trend analysis (A) and mutation point of annual runoff (B) in Gushanchuan basin.
Figure 2. Trend analysis (A) and mutation point of annual runoff (B) in Gushanchuan basin.
Water 12 01237 g002
Figure 3. Changes of precipitation (A) and potential evapotranspiration (B) in Gushanchuan basin.
Figure 3. Changes of precipitation (A) and potential evapotranspiration (B) in Gushanchuan basin.
Water 12 01237 g003
Figure 4. Land use changes of Gushanchuan basin.
Figure 4. Land use changes of Gushanchuan basin.
Water 12 01237 g004
Figure 5. Correlation of runoff depth against precipitation (A) and rainfall intensity (B) in Gushanchuan Basin.
Figure 5. Correlation of runoff depth against precipitation (A) and rainfall intensity (B) in Gushanchuan Basin.
Water 12 01237 g005
Figure 6. 19680821 (A), 19950728 (B) and 19940802 (C) flood process line of Gushanchuan Basin.
Figure 6. 19680821 (A), 19950728 (B) and 19940802 (C) flood process line of Gushanchuan Basin.
Water 12 01237 g006
Table 1. Four commonly used forms of F(ϕ) and F′(ϕ) based on the Budyko hypothesis.
Table 1. Four commonly used forms of F(ϕ) and F′(ϕ) based on the Budyko hypothesis.
FormsF(ϕ)F′(ϕ)
Budyko[ϕtanh(1/ϕ)(1 − e−ϕ)]0.50.5[ϕtanh(1/ϕ)(1 − e)]0.5 × [(tanh(1/ϕ) − sech2(1/ϕ)/ϕ)(1 − e) + ϕtanh(1/ϕ) e−ϕ]
Turc-Pike(1 + ϕ−2)−0.51/ϕ3 [(1 + (1/ϕ)2)1.5]
Fubaopu1 + ϕ − (1 + ϕα)1/α1 − (1 + ϕα)1/α − 1ϕα − 1
Zhang(1 + ωϕ)/(1 + ωϕ + 1/ϕ)(ω + 2ω/ϕ−1 + 1/ ϕ2)/(1 + ωϕ + 1/ϕ)2
Note: α and ω are parameters of each formula. According to optimization method, α = 2.5 and ω = 0.9.
Table 2. Combination of basin runoff generation mechanism.
Table 2. Combination of basin runoff generation mechanism.
TypeComposition
of Runoff
Impact FactorsTypeComposition
of Runoff
Impact Factors
1R = RsP, E, W0, i6R = Rsat + RintP, E, W0, i
2R = Rs + RgP, E, W07R = Rsat + Rint + RgP, E, W0
3R = RgP, E, W08R = RintP, E, W0, i
4R = Rs + RintP, E, W0, i9R = Rint + RgP, E, W0
5R = Rs + Rint + RgP, E, W0
Table 3. Comprehensive analysis table of runoff generation patterns.
Table 3. Comprehensive analysis table of runoff generation patterns.
ObjectsExcess Saturation RunoffExcess Infiltration Runoff
Average annual rainfallP > 1000 mmP < 400 mm
Runoff coefficient>0.4<0.2
Symmetry of Flow Process LineHighLow
Impacts of rainfall intensityLowHigh
Influencing factors of flowPa, PPa, i
Surface soil structureLooseCompact
Water shortage capacitySmallBig
Proportion of RgHighLow
RelativityRelated to PRelated to i
Table 4. Variation of meteo-hydrological factors around the mutation point of Gushanchuan basin.
Table 4. Variation of meteo-hydrological factors around the mutation point of Gushanchuan basin.
VariablesZCSignificantMean Value (mm)Change (%)
1965–19971998–2010
P−0.08ns427.1402.4−2.5
E0−0.29ns2426.32410.7−0.0
R−4.53**60.513.9−77.0
Note: ns indicates significance level exceeds 0.05; ** indicates significance levels of 0.0.5.
Table 5. Separating the effects of climate change and human activities by climate elasticity method.
Table 5. Separating the effects of climate change and human activities by climate elasticity method.
VariablesBudykoTurc-PikeFubaopuZhangAverage
εP2.62.72.42.32.5
ΔRC (mm)−6.8−7.0−6.3−6.0−6.5
ΔRH (mm)−39.8−39.6−40.3−40.6−40.1
ΔRC (%)14.615.113.512.913.9
ΔRH (%)85.484.986.587.186.1
Table 6. Land use classification statistics of Gushanchuan basin.
Table 6. Land use classification statistics of Gushanchuan basin.
TypeFarmland
(km2)
Forest
(km2)
Grass
(km2)
Water
(km2)
Urban Construction
(km2)
Sum
(km2)
Year
1976413.6865.47774.5512.286.021272
1996405.9448.41800.3212.824.511272
2006389.1683.69779.9412.167.051272
Average402.9365.85784.9412.425.861272
Proportion (%)31.685.1861.710.970.46100
Table 7. Characteristics of flood events in different stages of Gushanchuan basin.
Table 7. Characteristics of flood events in different stages of Gushanchuan basin.
Characteristics1965–19971998–2010Change (%)
Mean rainfall duration (h)12.317.340.7
Mean rainfall (mm)47.938.0−20.7
Mean flood duration (h)47.653.813.0
Mean peak flow (m3/s)1480.6314.7−78.7
Table 8. Statistics of runoff generation patterns in different stages.
Table 8. Statistics of runoff generation patterns in different stages.
PatternExcess Infiltration RunoffExcess Saturation RunoffMixed RunoffSum
Year
First stage
1965–1997
76 (68.46%)1 (0.10%)34 (30.63%)111 (100%)
Second stage
1998–2010
11 (45.83%)3 (12.50%)10 (41.67%)24 (100%)
Sum87 (64.44%)4 (2.96%)44 (32.60%)135 (100%)
Table 9. Statistics on the convergence characteristics of Gushanchuan Basin.
Table 9. Statistics on the convergence characteristics of Gushanchuan Basin.
Items1965–19971998–2010Change
(%)
Forest areas (km2)55.9483.6949.61
Annual minimum flow (m3/s)1.011.4032.67
Minimum daily flow in flood season (m3/s)1.932.7844.04
Mean precipitation before generate flow (mm)9.718.590.72
Duration of flood recession (h)14.3221.5450.42
Catchment water storage capacity (mm)73.1795.3830.35

Share and Cite

MDPI and ACS Style

Hu, C.; Zhang, L.; Wu, Q.; Soomro, S.-e.-h.; Jian, S. Response of LUCC on Runoff Generation Process in Middle Yellow River Basin: The Gushanchuan Basin. Water 2020, 12, 1237. https://doi.org/10.3390/w12051237

AMA Style

Hu C, Zhang L, Wu Q, Soomro S-e-h, Jian S. Response of LUCC on Runoff Generation Process in Middle Yellow River Basin: The Gushanchuan Basin. Water. 2020; 12(5):1237. https://doi.org/10.3390/w12051237

Chicago/Turabian Style

Hu, Caihong, Li Zhang, Qiang Wu, Shan-e-hyder Soomro, and Shengqi Jian. 2020. "Response of LUCC on Runoff Generation Process in Middle Yellow River Basin: The Gushanchuan Basin" Water 12, no. 5: 1237. https://doi.org/10.3390/w12051237

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