Climate change and Land Use/Cover Change (LUCC) have been identified as two primary factors affecting watershed hydrological regime. This study analyzed the trends of streamflow, precipitation, air temperature and potential evapotranspiration (PET) from 1962 to 2008 in the Jihe watershed in northwestern Loess Plateau of China using the Mann-Kendall test. The streamflow responses to climate change and LUCC were quantified independently by the elasticity method. The results show that the streamflow presented a dramatic decline with a turning point occurred in 1971, while the precipitation and PET did not change significantly. The results also show that the temperature rose markedly especially since 1990s with an approximate increase of 1.74°C over the entire research period (1962-2008). Using land use transition matrix, we found that slope cropland was significantly converted to terrace between 1970s and 1990s and that forest cover increased relatively significantly because of the Grain for Green Project after 2000. The streamflow reduction was predominantly caused by LUCC and its contribution reached up to 90.2%, while the contribution of climate change to streamflow decline was only 9.8%. Although the analytical results between the elasticity method and linear regression model were not satisfactorily consistent, they both indicated that LUCC (human activity) was the major factor causing streamflow decline in the Jihe watershed from 1962 to 2008.
Changes in streamflow reflect the combined effects of climate, vegetation and soil (Rodriguez-Iturbe et al., 2001). Climate and hydrology are intricately linked and the response of hydrological cycle to climate change has been widely investigated by researchers, e.g., Xu et al. (2013), Perazzoli et al. (2013) and Ashraf Vaghefi et al. (2014). Researches consistently unearthed basin-wide rising of temperature, spatial variations of precipitation, and basin-wide decrease in water resources in the Yellow River Basin (Xu et al., 2007; Hu et al., 2011). Li et al. (2011) used four climate models with the prescribed three emission scenarios to assess the impact of future climate changes on hydrology in a tableland region of Loess Plateau, and the modeling predicted -2.6%-17.4% changes for precipitation and 10%-130% changes for streamflow under the three emission scenarios. In brief, examination of precipitation-runoff relationships under various climatic change scenarios can provide the needed scientific references for sustainable water resources management.
Apart from the climate variability, Land Use/Cover Change (LUCC) induced by human activities (e.g., agriculture, industrialization, urbanization, ecosystem conservation and reservoir building) also influence regional hydrology regime (Cuo et al., 2013) and could profoundly alter the hydrological characteristics. For example, planting trees and grasses can alter the water balance of a catchment through increasing precipitation interception and enhancing evapotranspiration (Brown et al., 2005). Nunes et al. (2011) studied hydro-geomorphic responses between different land use and land cover types and found that vegetation dynamic was the key factor in modulating the hydrological responses to LUCC in a marginal region of Portugal. In addition, land reclamation and management changes, especially extensive channel engineering, were found to be more important than climatic variation in causing the observed increasing river runoff in some regions (Schilling et al., 2010). On the Chinese Loess Plateau, the comprehensive practices have been carried out to control severe soil erosion in the past 30 years, including afforestation, building terraces and check dams. Reservoirs were built for both irrigation and drinking water supply with an additional aim to mitigate soil and water losses through reducing surface runoff.
Interactions between LUCC and regional climate variations are important in controlling watershed hydrological processes. The impacts of climate change and LUCC on hydrological regime are actually theoretically interlinked and therefore cannot be completely separated. Liu et al. (2003) concluded that climate change was the major reason for the total runoff decline while the land cover change was a minor reason in the Yellow River source region in 1980s and 1990s. The Loess Plateau, located in northwest China, has attracted much attention because of highly erodible loessial soils, steeply-sloped landscape, frequent rainfall storms, and low vegetation cover (Li et al., 2010; Fan et al., 2015). A study in Wudinghe watershed found that the contributions of climate change and human activities to streamflow decline were 35% and 65%, respectively. Channel (e.g., check-dam) and slope measures (e.g., terrace) and the Grain for Green Project were identified as the dominant factors responsible for the streamflow decline (Liang et al., 2015). Another study on the Loess Plateau in Shanxi Province attributed the streamflow reduction to increased forest cover and concluded that the effects of LUCC and climate change on streamflow reduction were 53.2% and 46.8%, respectively (Tang et al., 2010). These conclusions were drawn from the different studies that covered different periods using different methods. Further researches are thus necessary for improving the understanding of combined effects of climate change and LUCC on streamflows.
Numerous methods have been developed for estimating the impacts of climatic change and LUCC on streamflow and these methods include climate models coupled with spatially-distributed hydrological models (Chen et al., 2005; Li et al., 2013; Al-Mukhtar et al., 2014; Fiseha et al., 2014). However, climate models are always a result of downscaling regional climate models or globe climate model and the coarse resolution cannot capture geographical details precisely. Therefore, they may not provide full representation of variables to the scale required by hydrological models. Consequently, coupling climate models with hydrological models bring great uncertainties. Alternatively, methods that use the observed meteorological and hydrological data and that use the interpolator of an appropriate form based on simple theories (e.g., water balance and Budyko hypothesis) can be more realistic for estimating the impacts of climatic change and LUCC on streamflow. Zhang et al. (2001) developed a conceptual model based on catchment water balance, which can assess the long-term effects of vegetation cover changes on evapotranspiration and other hydrological factors. Xu et al. (2014) and Yang et al. (2014), both based on the Budyko framework, estimated the streamflow responses to climate change using an elasticity method. Studies repeatedly showed that this simple model performed well and was potentially useful for investigating the watershed hydrological responses to environmental changes.
This study investigated the trends of climate variables and streamflow in the Jihe watershed for the period of 1962-2008 and the objectives were (1) to detect the statistical variability trends in meteorological and hydrological regimes; (2) to identify the abrupt turning points; and (3) to search for the dominant factors (climate change or human activity). The research results would be of significance for water resources management.
The Jihe watershed is located in Tianshui, a city of Gansu Province in China (34° 20′ 19″-34° 38′ 59″N, 105° 07′ 50″-106° 00′ 45″E; 1158 to 2699 m a.s.l.; Fig. 1). The catchment area is 1019 km2 and the branch stream is 61 km long. The mean annual precipitation is about 556 mm with the most falling in summer. The air temperature has large diurnal but small seasonal variations with an annual mean of 10.8° C. The land uses are predominantly croplands and grasslands. The watershed is mainly covered by cinnamon soil, loessal soil and alluvial soil, and terracing is the main soil conservation measure. Revegetation and pasture establishment have been implemented in the Jihe watershed to improve the eco-environment.
The daily streamflow data of Tianshui II hydrological station and the daily precipitation observation data of the four raingauge stations (Huangjizhai, Guanzi, Xujiadian and Tianshui; Fig. 1) for the study period were obtained from Chinese Hydrologic Data Yearbook. The air temperature, sunshine hours, relative humidity and wind speed were downloaded from the National Meteorological Information Centre of China Meteorological Administration (http://data.cma.cn/site/index.html). The watershed mean annual precipitation was computed by Thiessen polygon method.
Four satellite images, Landsat MSS (1975), Landsat TM (1990, 2000), and Landsat ETM+ (2008), were acquired from the Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences. We processed the images based on geospatial analysis platform in ArcGIS 9.2 to identify the land cover categories. Land use maps of the four periods (1975, 1990, 2000, 2008) were used to assess the impact of LUCC on hydrologic regime. According to the field investigation on geomorphologic features from April to September in 2009, combining with the measured data on soil and vegetation cover types in the hilly region of Loess Plateau, seven land use types were identified as forest, grassland, slope cropland, terrace, settlement, unused land, and water body.
3.2.1 Trend test of hydroclimatic variables
Based on the major (or abrupt) change point of the streamflow, we divided the meteorological and hydrological data series into the baseline period (before the change point) and the post-change period (after the change point). We assumed that there was no human interference in the baseline period. The potential evapotranspiration (PET) was estimated by Hargreaves equation that used the easily obtainable data including daily average, maximum and minimum air temperatures, and extraterrestrial radiation (Berti et al., 2014). The non-parametric Mann-Kendall (M-K) test and the moving T-test were combined to detect trends in hydroclimatic time series. The test statistic sequence UF is estimated by the sequential data series while the backward sequence UB is estimated by the inversed data series. The intersection point of UF and UB located between the confidence interval is the time when an abrupt change or a change point occurred in the M-K test (Zhang et al., 2013). The details about the two tests can be seen in Chen et al. (2010), Myronidis et al. (2012), and Yang et al. (2014).
3.2.2 Climate elasticity method
The elasticity coefficient, which was firstly proposed by Schaake (1990) and subsequently improved by others, is defined as the ratio of the streamflow variation rate to the variation rate of a certain climatic factor (precipitation or PET; Eq. 1).
where \(\varphi _X\) is the streamflow elasticity coefficient to climatic factors, Qis the annual average streamflow (mm), and X is a certain climatic factor. An analytical method, deriving climate elasticity of streamflow based on Budyko hypothesis, is well defined in theory and does not depend on a large amount of historical data, and consequently it has been applied in many studies (Yang and Yang, 2011). Equation 2 formulated by Zhang et al. (2001) has been used most.
where E, P and E0 are the mean annual evapotranspiration, precipitation, and PET, respectively. The parameter w describes the watershed characteristics of the combined underlying surface such as land cover, soil type and vegetation. Given P, Q and E0, the parameter w is estimated by Equation 2 through minimizing the difference between modelled E and Ecalculated from watershed water balance. Assuming the response of the streamflow to the climatic factors is principally caused by precipitation and PET, the changes in mean annual streamflow due to climate variability (∆ Qclim) can be calculated by Equation 3 through the complete differential theory (Milly and Dunne, 2002; Zhan et al., 2014).
where ∆ Qclim is the changes in annual streamflow due to climate variability, ∆ P and ∆ E0 are the changes in annual precipitation and PET, respectively; φ p is the precipitation elasticity of streamflow (Eq. 4) and \(\varphi _{E_0}\)is the PET elasticity of streamflow (Eq. 5).
where E0/P is the aridity index.
3.2.3 Separating the effect of climate change from the effect of LUCC on streamflow
For a given catchment, changes in streamflow can be attributed by the coupling effects of all human activities and climate parameters (precipitation and potential evapotranspiration induced). We assume that the total streamflow change (∆ Qtot) could be partitioned into these two factors (human and climate) and can be calculated by Equation 6.
$$\Delta Q^{tot}=\Delta Q^{clim}+\Delta Q^{LUCC}, \ \ (6)$$
where ∆ QLUCC is the changes in annual streamflow due to LUCC, ∆ Qtot can be calculated by the following equation:
$$\Delta Q^{tot}= \overline{Q}_2 -\overline{Q}_1 \ \ (7)$$
where \(\overline{Q}_1\)and \(\overline{Q}_2\)are the observed mean annual streamflow before and after the change point respectively, which can be determined and identified by the M-K test after the change point. Then the ∆ Qtot subtracts the result of Equation 3 and the response of streamflow to LUCC (∆ QLUCC) can be assessed through Equation 6.
Figure 2 shows the hydrograph for the period 1962-2008. It exhibited a significant and monotonous decreasing trend at a slope of -2.48 mm/a. The coefficient of variation reached about 75%, presenting rather large fluctuations over the 47-year study period. The average annual streamflow depth was 76 mm, and the maximum and minimum levels were 247.3 mm (1967) and 1.7 mm (1996), respectively. The average annual streamflow of the 1990s declined dramatically, only accounted for 43.1% of the mean annual of the study period. This was caused by the reduction in precipitation and the impact of LUCC. As seen in Figure 3a, annual cumulative anomaly streamflow fluctuated frequently since the middle 1970s and there were two major (or abrupt) change points occurred: one in 1974 and another in 1982. According to the M-K test (Fig. 3b), the UF statistic showed a downward trend since 1968. The mean UF series was -2.48 after 1990, apparently lower than the 95% confidence level of -1.96.
There were three intersections appeared in 1971, 1975 and 1976 for UF and UB (Fig. 3b) with the intersection in 1975 dropped out of the lower limit confidence level. Based on the principle that the intersection must be within the confidence lines, it indicated the beginning of an abrupt change. Therefore, the change points took place in 1971 and 1976, signifying that the streamflow exhibited two abrupt changes during the research period. The UF values in the four years between 1972 and 1975 (-2.34, -2.47, -2.62 and -2.14, respectively) had all surpassed the 95% confidence level, and so it is proper to define 1971 as the change point after which the steamflow significantly declined. Furthermore the statistic value had completely surpassed the 95% significance level since 1991, demonstrating that the streamflow presented a markedly sharp decline since then.
Figure 4 shows the variations of annual precipitation and PET in the Jihe watershed during 1962-2008. The decadal mean precipitation was lower in the 1990s than in the other periods. There was no significant fluctuation in precipitation with a linear slope at -3.0 mm/decade and PET exhibited a slightly rising linear trend at 11.7 mm/decade. Besides, as shown in Figure 5, the M-K test and the moving T-test demonstrated that the trends in annual precipitation and PET were not at 95% significant confidence level. The annual PET presented an increasing trend from the middle 1980s and the trend was also collaborated by the UF statistical value (Fig. 5a) and the moving T-test (Fig. 5b). However, as for annual precipitation, there were differences about the two tests (Figs. 5c and d). The two change points occurred in 1993 and 1994, meaning that the monotonic trend in annual precipitation was statistically significant after 1994.
Annual average temperature rose evidently since 1962 with a linear warming trend at 0.37° C/decade (Fig. 6a). Overall, the annual average temperature was increased by 1.74° C during the entire research period (1962-2008). A change point occurred in 1993 according to the M-K test (Fig. 6b), indicating the long-term monotonic trend in temperature is statistically significant. Furthermore, the M-K statistic values exceeded the confidence level of 95% since 1994, meaning that the annual mean temperature in the Jihe watershed presented a distinctive increasing trend during the 15-year period from 1994 to 2008. As for seasonal mean temperature (Fig. 6c), it rose more in winter and spring (the linear increasing trend at 0.49° C/decade and 0.37° C/decade, respectively) than in summer and autumn (the linear increasing trend at 0.27° C/decade and 0.34° C/decade, respectively).
Land use types of the Jihe watershed for 1975, 1990, 2000 and 2008 are presented in Table 1. In 1975, the slope cropland and grassland occupied nearly 80% of the total watershed, whereas the terrace only accounted for 6.5% (Fig. 7). From 1975 to 1990, the terrace was rapidly increased by 397.1% from 66.3 to 329.6 km2 and the slope cropland was reduced dramatically by 65.7%. The forest land, water body and settlement all declined slightly. With the implement of the Grain for Green Project in 1998, the forest cover increased steadily. The terrace, forest and grassland accounted for 87.8% of the watershed, and the three types (terrace, forest and grassland) were the predominant land use types in the Jihe watershed. It should be noted that the area of slope cropland was reduced by 77.2% over the research period and the area of terrace cropland was increased by more than four times with the largest changes occurred in 1975-1990. The areas of the other five land use types had relatively small variations (Fig. 7).
The land use transition matrix showed the inter-conversions among the land use types and reflected the dynamic processes of the LUCC. The four-period land use maps were fused using the dissolve tool of ArcMap to merge into one land use record to enhance the computation efficiency. The unchanged area and converted area for a land use type were calculated using the superposition function of the ArcGIS spatial analysis, and the land use transition matrixes in 1975-1990, 1990-2000, 2000-2008 and 1975-2008 were then obtained. The results were listed in Tables 2 and 3 and also shown in Figures 8 and 9. Table 2 shows that the main transition was the conversion of slope cropland to terrace cropland and grassland, the flow-out probability of slope cropland reached 78.8%, with 71.6% to terrace, 6.2% to grassland, and 1.0% to forest and settlement. Both forest and grassland showed very small conversions during 1975-1990, with the flow-out probabilities of 6.0%. The inter-conversions of land use types have been relatively small since 1990s (Table 3). Along with the Grain for Green Project, some terraces and slope croplands were transferred to forest and grassland. The period of 2000-2008 was the time period when the land use types had the smallest changes (Fig. 8).
Figure 9 further demonstrated the inter-conversions of land use types over the entire study period 1975-2008. The conversions mainly occurred in the middle and lower reaches of the Jihe watershed. The slope cropland declined continuously with 350 km2 being converted to the other land use types. The flow-out probability reached 85.7% with 74.5% to terrace and 8.0% to grassland. In short, slope cropland, terrace and forest were the land use types that had large portions of inter-conversions, and grassland remained relatively stable.
The mean, standard deviation (SD), coefficient of variation (CV) and extreme ratio (ER) of climatic factors and streamflow were summarized for the baseline period and also for the post-change period (Table 4). Average annual stramflow was reduced by 52.2% in the post-change period than that in the baseline period. The CV and ER in the post-change period were also much higher than those in the baseline period. ER in the post-change period (1972-2008) achieved 103.3, indicating that the streamflow during this period was more discrete and underwent more dramatic changes. Annual precipitation did not present evident changes since the above four statistical parameters were similar in the two periods. Table 5 shows the contributions of climatic change and LUCC to streamflow variations estimated using the elasticity method. The climatic change and LUCC accounted for 9.8% and 90.2%, respectively, demonstrating that the effect of human activities on hydrological regime was far more important than the climate change in the watershed. Wang et al. (2006) concluded that 35.9% of total annual streamflow reduction occurred in 1970-1999 was contributed by the climate change and 64.1% by human activity in the middle reaches of Yellow River, being consistent with the conclusions of this study.
Linear regression model was adopted in the study to compare with the elasticity method. The period division of hydrological series was consistent with the M-K test: the baseline period of 1962-1971 and post-change period of 1972-2008. The relationship between annual precipitation and annual streamflow in the baseline period was built to predict the natural streamflow in the post-change period (Eq. 8).
Y=0.5063X-139.66 (r=0.8055), (8)
where Y is the annual streamflow (mm), X is the annual precipitation (mm), and 0.5063 and -39.66 are regression slope and intercept, respectively (Li et al., 2011).
The simulated value in the post-change period minus the observed value in baseline period was the change in streamflow due to climatic change, and the observed streamflow in the post-change period minus the simulated value in the post-change period was the change in streamflow due to LUCC. Table 6 shows the responses of streamflow to climate change and LUCC. The contribution of LUCC to streamflow reduction reached 119.2 % in the 1970s while climate change imposed a counteracting impact on the streamflow decline (i.e., -19.2%). However, climatic factors contributed 33.1 % to the observed streamflow decline in the 1990s and LUCC contributed 66.9%. Overall, during the period from 1972 to 2008, LUCC contributed 104.8% to the streamline decline and climate change imposed a counteracting impact on the streamflow decline (i.e., -4.8%). The relationships among the hydrological elements (precipitation, streamflow, evapotranspiration, etc.) are complex and inconstant, it is thus difficult to adopt a precisely-defined equation to describe the relationships. Consequently the error unavoidably exists when we use regression model to analyze the linearity of precipitation-streamflow relationships.
Precipitation and PET are supposed to be the major meteorological factors that affect streamflow in river basins (Zhan et al., 2014). However, in the Jihe watershed, they changed relatively insignificantly compared to the streamflow, the latter (i.e., streamflow) displaying a remarkable decreasing tendency in 1962-2008. Quantitative assessment proved that LUCC was the major anthropogenic driving force for the streamflow change. Uncertainties of this study were mainly resulted from two aspects. Firstly, change points were identified by the M-K test to have occurred in 1971 and 1976. However, it should be admitted that the period divisions according to the different change points might lead to a slight difference in the assessment results of the contributions of climate change and LUCC to the streamflow. Secondly, although LUCC is the most relevant human activities, the widely-practiced comprehensive soil erosion control measures are also relevant human activities. However, due to the difficulties in data collection, the soil erosion control measures were neglected in this study. It should be emphasized that the widely-practiced comprehensive soil erosion control measures such as check dam should be investigated to reduce the uncertainties of the estimation in the future.
This research identified the changes of hydro-meteorological factors in 1962-2008 and quantified the effects of climate change and LUCC on the streamflow in the Jihe watershed. We concluded that the streamflow presented a significant decreasing trend during 1962-2008 with a change point occurred in 1971. The streamflow dropped sharply in the 1990s and the M-K statistic value exceeded the 95% significance level since 1991. The annual streamflow for the period (1990s) was about 43% of the mean annual streamflow of the entire study period (1962-2008). Corresponding to the streamflow decline in 1990s was a marked rise in temperature and neither decline in precipitation nor decline in PET was observed for the same time. The 1990s was an important period of climatic change when the temperature showed a trend of a continuous increase with a change point occurred in 1994 while the precipitation decreased. Based on the land use type analysis, the watershed was dominated by slope cropland and grassland in the1970s while the area of terrace was very small. Terrace and grassland became the main land use types and forest cover to some extent increased because of the Grain for Green Project in the 1990s. The land use transition matrix shows that the greatest inter-conversion was the flow-out of the slope cropland to terrace. LUCC was the main reason for the rapid streamflow reduction with the contribution reaching 90.2% while the contribution of climate change to streamflow decline only accounted for 9.8%. Although the analytical results between the elasticity method and linear regression model were not satisfactorily consistent, they both indicated that LUCC (human activity) was the major factor causing streamflow decline in the Jihe watershed from 1962 to 2008.
Acknowledgements
This study was funded by the National Natural Science Foundation of China (41501025, 51609083, 41401038, 51509089) and the 2016 Key Scientific Research Projects for Universities of Henan Province (16A170014). We sincerely thank the anonymous reviewers and editors for their valuable comments.
The authors have declared that no competing interests exist.
1 |
|
2 |
|
3 |
|
4 |
|
5 |
|
6 |
|
7 |
|
8 |
|
9 |
|
10 |
|
11 |
|
12 |
|
13 |
|
14 |
|
15 |
|
16 |
|
17 |
|
18 |
|
19 |
|
20 |
|
21 |
|
22 |
|
23 |
|
24 |
|
25 |
|
26 |
|
27 |
|
28 |
|
29 |
|
30 |
|
31 |
|
32 |
|
33 |
|
34 |
|