Separation of the Climatic and Land Cover Impacts on the Flow Regime Changes in Two Watersheds of Northeastern Tibetan Plateau

1Key Laboratory of Ecohydrology of Inland River Basin, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou, Gansu 730000, China 2University of Chinese Academy of Sciences, Beijing 100049, China 3School of Agricultural, Computational and Environmental Sciences, Institute of Agriculture and Environment, University of Southern Queensland, Springfield, QLD 4300, Australia 4College of Earth Environmental Sciences, Lanzhou University, Lanzhou, Gansu 730000, China


Introduction
Climate change and land use/cover change (LUCC) are overlapped factors and function as two critical drivers that directly influence hydrological processes in various watersheds [1,2].This is especially true in arid and semiarid areas, where the ecosystems are fragile and more sensitive to climate change due to the restricted availability of water supply.Researchers have thus dedicated significant efforts towards broadening the understanding of the influence of these two variables on hydrological processes, as well as on the flow regimes [3].As it is well-known, the increase in the concentration of the greenhouse gases emission is causing global climatic warming [4].Consequently, the temporal spatial distribution and the magnitude of precipitation, frequency, and intensity of floods and recharge and discharge of streamflow are likely to be altered accordingly [5].LUCC can alter the associations of streamflow generation, transpiration, interception, and ecosystem conservation [6,7].Deforestation, for instance, can produce higher streamflow generation and reduce lag time between precipitation and runoff leading to increase in peak flow, flood damage risks, and severe soil erosion [8,9].Conversely, reforestation and intensive agriculture rehabilitation can dramatically decrease the generation of streamflow and trigger an increase in the flow concentration time.Although several sources in literature have reported the interactions of forest change, climate, and hydrological activities [10], how these factors are incorporated into hydrological processes as well as the different flow regimes still remain unclear.It is thus important to gather further insights into the impacts of the climate and LUCC factors on the flow regimes in order to accurately quantify their respective roles in different watershed regions.
Available hydrological impact assessment methods commonly used in literature include the catchment paired experiments, statistics analysis, and measurements with hydrological models [11], but seldom have these been involved in streamflow regime studies [12,13].Moreover, the experimental and statistical methods treat the study basin as a black box and rarely examined the complexities of precipitation changes, underlying surface conditions, and the interactions between climate change and the respective hydrological processes [2].Hydrological models, however, offer a viable framework for conceptualizing and investigating the relationships between climate, underlying surface, and hydrological processes in various categories of time and space [14,15] and the approach of one-factor-at-a-time (OFAT) has been widely applied [16,17].For instance, Karlsson et al. [18] modeled the combined effects of land use and climate changes on the hydrology for a catchment located in Denmark.Zhao et al. [19] evaluated the climate variability and land use influences on green and blue water resource in Weihe River basin, NW China.These methods considered the hydrological processes and their interactions with the environment.It should be noted that hydrological models are more effective as they relate the model parameters directly to the physically observable land surface characteristics [20].However, the assumption within OFAT method is that in the course of evaluating the influence of a given factor on hydrological processes, the other factors' effects are not considered.In fact, the other factors not considered are changing concurrently over the entire period of observation which can contribute as an apparent bias in the separation results.Separation of the individual influence of the two factors on hydrological processes thus warrants the need for inputting the applicable status of two factors from the baseline period to the entire period.Yang et al. [21] recommended that a combined statistical and modeling method could be used to resolve this bias and further identified the climate and land cover change impacts on the hydrological processes in the Heihe River.In that study, the results were more reasonable and offered greater accuracy than the traditional OFAT and other statistical methods.
The Yingluoxia catchment (YLC) and Minxian catchment (MXC), utilized in this study, are situated in the northeastern Tibetan Plateau, which are primarily influenced by the northern hemispheric middle latitude Westerlies and the Asia summer monsoon, respectively [22,23].The two watersheds are nearly located at the same longitude, but quietly different latitude.It means the further north the watersheds, the less the water vapor content received.Meanwhile, these regions are among the most sensitive zones in terms of climate change [24].More than 1.3 million people and 266,000 ha of irrigated agricultural land rely on the streamflow from YLC region [25].This flow is crucial for maintaining the oasis and agricultural ecosystems [26,27].The MXC is considered as an important water resource hub for the Yellow River, but it is also a source region addressing serious water shortage area in the Gansu Loess Plateau [28].Therefore, the quantification of the impacts of climate change and LUCC on flow regimes in these regions is crucial in order to clarify the hydrological response to different climate types and understand the flow regime behaviors with different climate and land cover combinations.The objectives of this study are, thus, (1) to define a conceptual framework and to propose a modified OFAT method for separating the individual contribution of climate and LUCC to flow regimes, (2) to improve the separation accuracy and to compare the separation results with climatic conditions, and (3) to compare and gain further insights and understand the effects of the two major factors on streamflow.The results can provide a beneficial reference for an assessment and proper management of the water resources for the associated policymakers and stakeholders within the YLC and MXC regions.

The Yingluoxia Catchment (YLC)
. YLC, the headwater region of the Heihe River basin, is located in the middle of the Hexi Corridor, China, with a total drainage area of 10,009 km 2 lying between 99 ∘ and 101 ∘ E and 38 ∘ and 39 ∘ N. The elevation varies from 5058 m to 1668 m (Figure 1).YLC is the main region for the generation of streamflow in the entire basin.Approximately 90% of the water resources in the middle and lower reaches are therefore recharged by the streamflow from the YLC region.The climate, mainly controlled by the northern hemispheric middle latitude Westerlies, is characterized as cold and moist with large spatial and temporal heterogeneity [29].The average annual precipitation is more than 400 mm and increases by 15.5-16.4mm for every 100 m increase in the elevation.Nearly 70% of annual precipitation occurs from June to September and less than 10% from November to March.The annual average air temperature varies from −5 ∘ C to 4 ∘ C, and decreases 0.5 ∘ C along with the elevation increase for every 100 m [30].The total annual runoff is 16.05 × 10 8 m 3 with significant interannual variability [31].Conditions within the Heihe River have attracted great attention in China due to increasing pressures in the nation, particularly in water resources and the degradation of natural environmental sectors [32,33] 2.1.2.The Minxian Catchment (MXC).MXC, the headwater of the Taohe River, is covered by a total drainage area of 14,912 km 2 with an annual streamflow of 26.8 × 10 8 m 3 , lying between 101 ∘ and 105 ∘ E and 34 ∘ and 35 ∘ N, dominated by the Asia summer monsoon.The land cover and topography change from forests and mountains to grasslands and open valleys [28].The annual average air temperature increases from 1 to 6 ∘ C at elevations ranging from 4562 to 2,411 m.The mean annual precipitation decreases from 640 to 560 mm from west to east of MXC region.The dominant climate
It is noteworthy that the base flow index (BFI) is the fraction of the base flow over the total streamflow [16], which is contributed by the groundwater.We have used BFI to reflect the low-flow characteristic of the streamflow component.Due to the restrictions by the simulation performance of the long term daily flow regimes in our study area, it was impossible to describe the fine-detail characteristics of daily flow regimes.Monthly flow metrics would thus be a valid alternative, and this has been widely used in the present water management strategies.However, the monthly flow regime could still capture several characteristics of daily flow regimes [36].

SWAT Model.
SWAT is a physically and mathematically based hydrological model developed by the research centers at US Department of Agriculture [37].The model requires significant quantities of specialized information, including weather, soil properties, topography, vegetation, and land management practices.The resulting data is used to determine parameters controlling hydrological processes in a given watershed.A detailed description of the model and its input/output variables can be found in the model's public domain website at http://www.brc.tamus.edu/swat.SWAT processes hydrological simulation by applying the principle of water balance in each of the hydrologic response units (HRUs) [38], which define areas with identical combinations of surface slope, land use, and soil type [18].The water balance equation is given as where SW  is the final soil moisture (mm), SW 0 is the initial soil moisture at the time of  (mm),  day is the precipitation at the time of  (mm),  surf is the streamflow at the time of  (mm),   is evapotranspiration at the time of  (mm),  seep is the water flow to the unsaturated zone from the soil profile at the time of  (mm), and  gw is the water flow from the watershed from underground at the time of  (mm).In this study the SCS (Soil Conservation Service) curve number method is applied for computing streamflow volume.The Muskingum and Penman-Monteith methods are used for flow routing and estimating the potential evapotranspiration.

Framework for Separating the Effects of Climate Change and LUCC.
In a given river basin, the changes in hydrological processes are caused by climate change and LUCC, both of which are assumed to be independent factors [39].In this study, the hydrological regimes are treated as a function of these two factors.To separate the discrete effects of climate change on the flow regime, we have created a schematic in order to illustrate the variation in hydrological processes along with climate change under differing LUCC scenarios (Figure 2).Due to the different response relationships between different land cover and hydrological processes, we show the schematic of flow regimes response to land covers  1 and  2 (Figure 2).Note that Δ 1 and Δ 2 are the responses of flow regimes to the climate change (Δ) under land cover conditions  1 and  2 .Normally, Δ 1 and Δ 2 are different, when the land cover change (Δ) is smaller, Δ 1 is closer to Δ 2 .Thus, we use the average of Δ High: 4562 to denote separate impacts of climate change on hydrological elements (Δ  ): where  and  2 under the climate condition of  1 and  2 , denoted as Δ 1 and Δ 2 , respectively.The smaller the change in climate change (Δ), the closer Δ 1 is to Δ 2 .Thus, we accept the arithmetic mean of Δ 1 and Δ 2 as the separate impact of LUCC to hydrological elements (Δ  ) as The total change in hydrological element Δ equals the sum of the two factors contributions.The total change also can be estimated by the differences of the hydrological elements in the two periods, as where Δ represents the total change of a given hydrological element; hydrological elements can serve as statistic mean values over annual, seasonal, or monthly time scales [21].
In this study, we select the period 1982-1995 as baseline situation and 1996-2010 as affection period, along with land use for each period reflected on a separate map.The land use map for the 1980s and 2000s is used to represent the land use patterns for the two periods.A calibrated SWAT model is then applied to each of the four permutations derived from the two time periods and the two land use maps; hereinafter these are referred to as four scenarios.The influences of the LUCC and climate change are quantified by comparing the SWAT outputs of the four scenarios by using the separation equations ((2) and ( 3

Model Calibration and Validation.
The SWAT model is calibrated using historical data from January 1982 to December 1995.The validation period is based on the observed runoff from January 1996 to December 2010.And the warming up period is January 1980 to December 1981.Correlation coefficient (), Nash-Sutcliffe coefficient (NSE), and Relative Bias (BIAS) are used to evaluate the performance of SWAT during calibration and validation [40,41].The equations are defined as where  is the number of data points,  obs is the observed runoff,  sim is the simulated runoff, and  obs and  sim are the mean values of the observed and simulated runoff, respectively.Model performance is considered satisfactory if NSE is greater than 0.5 and BIAS is less than 25% [42].

Calibration and Validation of SWAT Model.
The simulated monthly streamflow has good correspondence to the observed series during calibration and validation (Figure 3).In accordance with Table 2 and Figure 4, there exist significant correlations between the observed and simulated streamflow at monthly and annual scales.That is, the BIAS is found to be approximately 1.14% and 6.80% for the calibration and −0.15% and 2.7% for the validation period in the YLC and MXC regions, respectively.The correlation coefficient in the YLC region is 0.95 and 0.96 for the monthly streamflow and approximately 0.94 and 0.96 for the yearly streamflow, in the calibration and the validation period, respectively.In the MXC region, the correlation coefficient is 0.90 and 0.94 for the monthly streamflow and approximately 0.91 and 0.96 for the yearly streamflow, in the calibration and the validation periods, respectively.The NSE in the YLC region is 0.91 and 0.90 for the monthly streamflow and approximately 0.88 and 0.76 for the yearly streamflow, in the calibration and the validation periods, respectively.In the MXC region, the NSE are 0.85 and 0.74 for the monthly streamflow and approximately 0.80 and 0.78 for the yearly streamflow during the calibration and the validation periods, respectively.Although for the both watersheds the performance of the SWAT model during the validation period is found to be less accurate compared to the calibration period and the monthly performance generally better than the yearly results, the results of the performance are still considered   satisfactory, which indicates that the fundamental rainfallrunoff relationship and water balances including the intraannual distribution are well captured.We thus aver that these results are within the range of "the very good performance" benchmarks [42].

Characteristics of Land Use/Cover Changes.
In order to investigate the change in land use from 1980 to 2010, we integrated the different sources of land use conducts in 1980s and 2000s and compared the fraction of main land use in two time slices from 1980s to 2000s.Figure 5 shows the fraction of land use in YLC and MXC regions.The land use structure shows slight changes from 1980s to 2000s.Pasture is the dominate land use pattern, occupying about 52% and 61% of total area of YLC and MXC, respectively, followed by barren land and forest land, occupying 23% and 22% of total YLC, and forest land and farm land, occupying 27% and 7% of total MXC, respectively.The other land use types occupy less than 2.5%, which have limited effects on the runoff regimes.
Figure 6 shows the spatial distribution of land use conversion in the YLC and MXC regions.There are approximately 331 km 2 in YLC and 176 km 2 in MXC region that have changed their land use types occupying about 3.3% and 1.2% of their entire catchment that underwent type conversions, respectively.We mainly focus on the land use conversions that occupied large areas and may lead to significant change in runoff regimes.In the YLC region, the forest increased mainly due to the conversion from the pasture (63.65 km 2 ), which is primarily attributed to the afforestation programs and ecological environment restoration programs in Qilian mountain launched in early 2000s [43].Pasture area increased results from the conversion from barren land (226.03km 2 ), which is due to the global warming that led to pasture extending to higher altitude area [33].The other two conversion types in YLC are pasture to farm land (13.3 km 2 ) and glacier to barren land (28.37 km 2 ), which occupies about 12% of total land use conversion area.The conversions of land use indicate the increase of vegetation dynamics in the YLC region, which may lead to increase in canopy interception and transpiration through vegetation, meanwhile decreasing streamflow.In the MXC region, the three main conversion types are forest to pasture (34.44 km 2 ), pasture to farmland (92.37 km 2 ), and farmland to pasture (37.23 km 2 ) (Figure 6(b)).The forest conversion to pasture mainly occurred due to deforestation at three forest farms.The conversion of pasture to farmland is the largest land use conversion type in the MXC region, which, due to herdsmen, has settled down and grasslands are reclaimed [44].The conversion of farmland to pasture mainly results from grain for green program implemented in early 2000s in loess area of low stream [34].The other three conversion types occupy small area: forest to farmland (3.26 km 2 ), pasture to residential (3.59 km 2 ), and farmland to residential (5.89 km 2 ).These conversions are relevant to urbanization and reclamation.From the analysis above, we can conclude that the climate change influences the land use conversion in the YLC region, and human activities influence the land use conversion in the MXC region over the 30 years.

Observed Changes in Precipitation, Air Temperature, and
Streamflow. Figure 7 shows the yearly variation of precipitation, streamflow, and air temperature in the YLC region and MXC for the entire modeling period as well as a linear fitted to the observed values.Similar changes of these elements occurred in the two watersheds.The trends of precipitation and runoff decrease from 1980 to 1995 and increase from 1996 to 2010 in both watersheds.However these decrease or increase trends in YLC are slighter than that in MXC.This indicates that variation of water vapor sources due to climate change is more obvious in MXC than in YLC.This is mainly because the Asia summer monsoon tends to affect the MXC region more significantly than the YLC region resulting in more water vapor received in MXC than YLC region [22].However, the results of precipitation and streamflow trend magnitudes did not pass the significance test of 0.05, except for the streamflow in the MXC region during the period 1980-1995.The increasing trend of air temperature in the YLC region is more obvious than in the MXC region with rate of 0.0381 ∘ C/a ( < 0.05) in the YLC compared to 0.0265 ∘ C/a ( < 0.05) in the MXC region from 1980 to 1995.It is important to note that the significant increasing of air temperature magnitude means significant effects on the watershed evapotranspiration.From the difference of the change rates between precipitation and runoff, we can find that evapotranspiration in YLC is not intensive compared to that in MXC.Based on the above analysis, it is evident that the precipitation and streamflow have changed to a notable extent around 1995 and 1996.In accordance with this, we divided the period 1980-1995 as baseline phase and the period 1996-2010 as the affection phase in our study.

Effects of Climate Change and LUCC on Streamflow
Regimes.Table 3 shows the results of streamflow simulated by the SWAT model under the four hypothetical scenarios as presented in Section 3.3.Compared to (1), the simulated average annual runoff in 4 appeared to increase by approximately 5.6 mm in the YLC region and appeared to decrease by approximately 8.0 mm in the MXC region.These changes consequently result in, respectively, 3.1% and 4.3% of their total streamflow changes.It represents the combination of climate and land use change effects during the period 1980 to 2000 in the two watersheds.The difference between (1) and (2) and between (3) and (4) shows the influence of land use change during the two subject periods.The separation equation ( 3) is used to assess the impact of land use change.The results also suggest that the land use change leads to a decrease in average streamflow by approximately 0.4 mm and 0.5 mm, which then accounts for approximately 7.30% and 6.25% of the annual average streamflow changes change shows opposite effects on streamflow and the impact of climate exceeds by far that of LUCC in two studied watersheds.

Effects of LUCC on Flow
Regimes.In this paper, the simulated results instead of the measured data are used to compare hydrological model outcomes.The flow duration curve (FDC) is applied to the monthly stream flow to further quantify the temporal variation of the hydrological regimes.Flow duration curves are considered as a simple yet comprehensive approach to stream flow analysis that gives a graphical view of the overall variability associated with stream flow over a period of time [16].Flow duration curves for the YLC and the MXC region are shown in Figure 8  This means the LUCC impacts on the stream flow in MXC are almost negligible.Based on the above analysis, it is evident that the LUCC has limited impacts on flow regimes in both watersheds, primarily attributable to the area of land use changes that are small compared to entire watershed area.4, we can see that climate change leads to all streamflow metrics increasing in the YLC region, except for the decreasing of CvMF.This implies that climate change increases both the high-flow and low-flow, while the variation of streamflow is decreasing.In contrast, climate factors have an opposite influence on streamflow metrics in the MXC region.That is, climate has led to a decrease in all streamflow metrics, except the BFI, indicating that the climate change was able to induce a decrease in the highflow, low-flow, and variation of streamflow.LUCC cause the MMF, Lw75, Hg25, and CvMF decrease and DMF and BFI increase in both two watersheds.These changes imply that LUCC can cause a decrease in the magnitudes and variations of streamflow while increasing the base flow.Meanwhile, climate change factor has larger contributions to all streamflow metrics compared to LUCC in both watersheds.In accordance with this, we can conclude that all streamflow metrics have a larger contribution from climate change than LUCC, indicating that climate change is more relevant to a shift in the flow regimes.It can thus be averred that the contribution of climate change acts to dominate the change in flow regimes and the contribution of land use deduction is relatively small.This study attempted to disentangle the effects of climate change and LUCC on streamflow in two selected watersheds and to compare their different effects towards the dominant climate types and land cover conditions.In published literature, there have been some previous research works that have focused on the investigation of climate and land cover changes impacts [21,26,[45][46][47][48].The results of these studies agree with our study with the primary finding that climate change has influenced the streamflow regime more significantly than land use/cover change in the study regions [21,48,49].It is worthwhile to note that our study also emphasize the impacts of climate and land use/cover change on flow regimes and the latitudinal effect of climate on streamflow regimes.The results of this study further confirmed that the climate change offers a dominative effect on the magnitude, variability, duration, and component of streamflow regimes.

Conclusions
Deficit of water resources in China has raised the importance of investigating the causality of hydrological processes changes.It is widely accepted that these processes are controlled mainly by climate conditions.Notwithstanding this, it is also believed that they are strongly influenced by LUCC.The latter effect is produced by anthropogenic activities, which, in part, are the causes of climate change.In the present study, the SWAT method is successfully applied to separate the climatic and land cover impacts on the flow regime changes in two watersheds of northeastern Tibetan Plateau.This study has confirmed the usefulness of SWAT model for separating the impacts of climate change compared to those from LUCC on flow regimes.
The integration of the effects of LUCC and climate change was seen to increase annual streamflow by 5.6 mm in the YLC region and a subsequent decrease of annual streamflow by 8.0 mm in the MXC region, of approximately 3.2% and 3.4% of the total streamflow, respectively.It was evident that climate change shows an opposite effect in two watersheds.The results showed that climate change is confirmed to outweigh the impacts of LUCC on streamflow, leading to an increase in streamflow in the YLC and a decline in the MXC region, measured as approximately 107.3% and 93.75% of the total streamflow changes, respectively.The LUCC influence showed changes of approximately −7.3% and 6.25% of total runoff.These results outcomes were based on analysis of data collected over the last 30 years, showing that the climate change has significant impact compared to LUCC on flow regimes in the study region.Climate change is thus more relevant to shifting the flow regimes and can be seen to dominate the hydrological processes changes in northeast Tibetan Plateau.
Evident from study, we conclude that it is critical to apply an integrated approach, which combines the effects of climate and LUCC, for an accurate assessment of flow regimes.Future investigations should thus address the projection of different environmental scenarios for the study area and also take into account the uncertainty inherent in predicting both climate change and LUCC.

Figure 1 :
Figure 1: Location of Yingluoxia catchment (YLC) and Minxian catchment (MXC) together with distribution of digital elevation model and hydrometeorology stations.
and  2 2 donate the flow regimes value at points A, B, C, and D, respectively.Similarly, the responses of flow regime to LUCC are calculated by applying difference of hydrological component values of land use conditions  1

Figure 2 :
Figure 2: Schematic for separating the contribution of climate change and LUCC to flow regimes.

Figure 3 :
Figure 3: Simulated and observed monthly runoff variation in YLC during (a) calibration and (b) validation and in MXC during (c) calibration and (d) validation.

Figure 4 :
Figure 4: Scatter plots of observed and simulated runoff in YLC at (a) monthly and (b) annual scales and in MXC at (c) monthly and (d) annual scales.

Figure 5 :
Figure 5: Fraction of land uses in YLC (a) and MXC (b) from 1980s to 2000s.

Figure 6 :Figure 7 :
Figure 6: Spatial distribution of dominated land use conversion in YLC (a) and MXC (b).
with different time slices.The contrast between (1) and (2) and between (3) and (4) shows the influence of land use change in the two study periods.From Figures8(a) and 8(b), we can see the differences between the two land use situations in the YLC region.The curves for the land use situation within 2000s are almost under the below of the curves for the land use situation of 1980s.This indicates that the LUCC leads to a decrease in monthly stream flow in the YLC region although the decreasing effects are very small.Figures 8(c) and 8(d) present the differences of LUCC impacts on the flow duration curve in the MXC region.Evidently, the two curves for the land use situation of 1980s and 2000s are almost overlapped.
Figure 9 presents the climate change impacts on stream flow under different land use situation in two catchments.The curves for two slices are intertwined in each land use situations and respective catchments.This reflects the uncertainty and complexity of the climate variation within these regions.Compared to the impacts of LUCC, climate change impacts on the stream flow are more clear, indicating that the influences of climate factor on streamflow are predominant in these two study catchments.However, the impacts of climate change and LUCC on stream flow are various in different phases of flow duration curve.Thus, we selected six representative streamflow regime metrics to identify the separate impacts of climate change and LUCC on stream flow in two watersheds.Due to the combined climate and land use change influences, the six flow regimes appear to have changed accordingly.From Table

Figure 9 :
Figure 9: Impacts of climate change on flow duration curve in YLC for (a) LU1980s and (b) LU2000s and in MXC for (c) LU1980s and (d) LU2000s.

Table 1 :
The flow regime metrics selected to describe streamflow paradigm.

Table 2 :
Performance assessments of the SWAT model in the two watersheds.

Table 4 :
Separate contribution of climate change and LUCC on runoff regimes in YLC and MXC.