Hydrologic Responses to Land Use Change in the Loess Plateau: Case Study in the Upper Fenhe River Watershed

We applied an integrated approach to investigate the impacts of land use and land cover (LULC) changes on hydrology at different scales in the Loess Plateau of China. Hydrological modeling was conducted for the LULC maps from remote sensing images at two times in the Upper Fenhe River watershed using the SWAT model. The main LULC changes in this watershed from 1995 to 2010 were the transformation of farmland into forests, grassland, and built-up land. The simulation results showed that forested land contributed more than any other LULC class to water yield, but built-up land had most impact due to small initial loss and infiltration. At basin scale, a comparison of the simulated hydrological components of two LULC maps showed that there were slight increases in average annual potential evapotranspiration, actual evapotranspiration, and water yield, but soil water decreased, between the two intervals. In subbasins, obvious LULC changes did not have clear impacts on hydrology, and the impacts may be affected by precipitation conditions. By linking a hydrological model to remote sensing image analysis, our approach of quantifying the impacts of LULC changes on hydrology at different scales provide quantitative information for stakeholders in making decisions for land and water resource management.


Introduction
Changing land use represents the greatest human impact on the surface of the earth, and between one-third and onehalf has been altered over the last 300 years [1][2][3].Land use directly impacts the hydrologic components of a watershed, such as evapotranspiration, surface runoff, groundwater recharge, and streamflow [4][5][6][7].Many previous studies from around the world have demonstrated that land use significantly impacts hydrological processes [8][9][10][11][12][13][14][15].For example, Ghaffari et al. (2010) found that a 34.5% decrease of grassland with increases of shrubland (13.9%) and rain-fed agriculture (12.1%) led to a 33% increase in the amount of surface runoff and a 22% decrease in the groundwater recharge in northwest Iran [8].Nie et al. (2011) revealed the effects of urbanization and agricultural development on the runoff in North America [9], and Wang et al. (2014), Li et al. (2015), and Yao et al. (2015) revealed the effects of urbanization, land use patterns, and changes on the runoff, temperature, evapotranspiration, and surface climate in Chinese landscapes [10][11][12].López-Vicente et al. (2013) and Ohana-Levi et al. (2015) modeled the effects of land use change on the runoff and sediment yield in Europe's Pyrenees and rainfall-runoff relationships in an Eastern Mediterranean watershed, respectively [13,14].
2 Advances in Meteorology Wang et al. (2014) assessed the impacts of land use changes on evapotranspiration, quickflow, infiltration, and baseflow in the upstream regions of the Three Gorges reservoir in China [15].Assessing impacts of land use and land cover (LULC) changes on hydrology is essential for watershed management and ecological restoration [7,15,16].
In general, regional impacts of land use change on hydrology vary from place to place and need to be considered for specific circumstances and environments [15,16].The Loess Plateau of western China is infamous for serious soil erosion and is the main source of sediment in the middle reaches of the Yellow River and has become the focus of world attention.As early as the 1950s, China started soil and water conservation works on the Loess Plateau in order to improve the local environment and mitigate sediment accessions to the middle reaches of the Yellow River, and works were expanded after the 1970s [17,18].As a consequence, land use in the Loess Plateau has undergone great changes.Likewise, in the past century, a range of management measures were carried out in the Fenhe River Basin on the eastern margin of the Loess Plateau.Some research has been done on the impacts of land use change on basin hydrology [19][20][21][22], but research on the contribution of individual LULC to the total runoff and the impacts of LULC changes on watershed hydrology is lacking.In particular, there is a lack of information for evaluating the benefits of soil and water conservation in the Fenhe River basin and the whole of the Loess Plateau, where it is difficult to distinguish the impacts of LULC changes on hydrology.Greater understanding of the contribution of individual LULC change to runoff and the impacts of LULC changes on the hydrology at different scales is needed to guide comprehensive natural resources management in this region.
Methods for assessing hydrological impacts of land use changes in watersheds include comparisons of paired catchments, statistical analysis, and hydrological modeling [5,23,24].Among these approaches, hydrological modeling is the most suitable one for use in scenario studies at different scales.Widely used hydrological models in studies on the impacts on watershed hydrology include the Hydrologic Simulation Program, FORTRAN, the Soil and Water Assessment Tool (SWAT), and WaTEM/SEDEM [4,5,9,15].The SWAT model is widely used to assess hydrology and water quality in agricultural catchments around the world (https://www.card.iastate.edu/swatarticles/).It is readily available and user-friendly for data input [25].
The overall objectives of this study are to investigate the contributions of individual LULC change to runoff and to determine the impacts of LULC changes on the hydrology of the Upper Fenhe River watershed by an integrated approach that combines hydrological modeling and remotely sensed digital maps of LULC changes.To avoid the impacts of largescale water resource development projects on hydrological processes, particularly the water diversion project from the Yellow River to Fenhe River at Wanjiazhai in the study area after 2002, the period from 1990 to 2000 was selected as the study period, during which sufficient data are available.The specific objectives are (1)

Study Area
The Fenhe River is located in the eastern part of the Loess Plateau, North China, and runs north to south through almost all of Shanxi Province [26].The study area is the Upper Fenhe River watershed above the Hecha hydrological control station, which measures the amount of water and sediment entering Fenhe reservoir (Figure 1).The Upper Fenhe River watershed has an area of about 3242 km 2 and lies between latitudes 38 ∘ 00  and 39 ∘ 00  and longitudes 111 ∘ 40  and 112 ∘ 40  .Elevations in the basin range from 1,124 to 2,762 m and it has an arid to semiarid climate with a mean annual precipitation of 430 mm.The precipitation is concentrated between the months of June and September and accounts for approximately 70% of total annual precipitation, coming in the form of high-intensity storms that often cause extreme soil erosion.The average annual temperature is 9.0 ∘ C, with a minimum average temperature of −5 ∘ C in winter and a maximum average temperature of 20 ∘ C in summer.The land use classes of the watershed include farmland, forest, grassland, water or wetland, built-up land, and unused land; their percentages in 2010 were 30%, 28.3%, 39.57%, 1.39%, 0.73%, and 0.01%, respectively.The region is covered with highly erodible loess-derived soils.Major soil series include loess soil, cinnamon soil, brown soil, and calcium skeletal soil.

Methodology
A conceptual framework consisting of three stages was developed (Figure 2): (1) transition matrix analysis to describe the transformations of land use from 1995 to 2010; (2) hydrological modeling to simulate hydrological components for two LULC maps under the same climatic conditions; and (3) comparison of the contribution of each land use class to the watershed's water yield and their hydrological effects on subbasins.

Transition Matrix. The transition matrix of land use change was calculated as
where  is the area of a certain kind of land use,  is the number of land use types in the study area,  and  (,  = 1, 2, . . ., ) are the types of the land use before and after  transition, respectively, and   is area of the land use  turned into the land use .The study period was defined as 1995 to 2010 because soil and water conservation practices were first implemented from 1988 to 1997, and a second tenyear conservation project was carried out from 1998 to 2007 [19,20,27].The LULC maps in 1995 and 2010 represent the surface cover conditions in the latter part of the first conservation project and after the projects, respectively.The LULC map for 1995 was used as input data to calibrate the model.

SWAT Model Description.
The SWAT 2005 model [28] was applied to data of the Upper Fenhe River watershed to assess the impacts of LULC changes on hydrological components.The SWAT model is a continuous, long-term, physically based distributed model developed to assess impacts of climate and land management on hydrological components, sediment loading, and pollution transport in watersheds [25].In the SWAT model, a watershed is divided into subbasins.Subbasins are further divided into a series of uniform hydrological response units (HRUs) based on slope, soil, and LULC changes.Hydrological components, sediment yield, and nutrient cycles are simulated for each HRU and then aggregated for the subbasins.Hydrological components simulated by the SWAT model include evapotranspiration (ET), surface runoff, percolation, lateral flow, groundwater flow (return flow), and transmission losses [25].Evaporation and transpiration are simulated using exponential functions of soil depth and water content and a linear function of potential evapotranspiration (PET) and  1).
The outputs of the SWAT model used in this study were hydrological components, including precipitation (P), PET, actual evapotranspiration (AET), water yield (WY), and soil water (SW).

SWAT Model Calibration and Validation.
Simulations using the 1995 LULC map were used to calibrate monthly streamflows from January 1992 to December 1995 at two hydrological stations (Jingle and Hecha, Figure 1).After model calibration, simulations with the same model parameters and LULC map data were used to validate monthly streamflow from January 1996 to December 2000 (Table 2).

Performance Evaluation Criteria.
The model was calibrated manually and three criteria were used to evaluate ) describes the degree of colinearity between simulated and measured data and the proportion of the variance in the measured data explained by the model.Percent bias (PBIAS) measures the average tendency of the simulated data to be larger or smaller than their observed counterparts [29].The calibration and validation performance for the SWAT model were considered acceptable when  2 and NSE were greater than 0.5.The SWAT model performance is satisfactory when NSE is larger than 0.5, adequate when NSE ranges from 0.54 to 0.65, and very good when NSE is larger than 0.65.The SWAT model is rated as satisfactory, when the absolute value of PBIAS ranges from 15 to 25, rated good when from 10 to 15, and very good when smaller than 10 [30].It was transformed into grassland in the northern part of the watershed due to the implementation of the "grain for green" government conservation program that compensated farmers for converting arable land to grassland.Meanwhile, due to the social development, the area of built-up land increased markedly at the expense of farmland.Finally, a tiny proportion of it was changed into water or wetland and unused land.In the southeast, a part of the watershed grassland was turned into farmland.Overall, the vegetation cover increased from 1995 to 2010, during periods of increased soil and water conservation measures.

Model Calibration and Validation.
Parameter sensitivity analysis was carried out using the model parameters identified from the instructions for the calibration of the SWAT model, as given in the user's manual [29].The model was calibrated with a comprehensive model calibration method [30].The results of sensitivity analysis for the Upper Fenhe River watershed and the value of the parameters are shown in Table 4.
Comparisons between simulated and observed monthly streamflow values in the periods of calibration and validation are shown in Figure 4. Good matches were shown between simulated and observed values, except in 1997 and 2000.In these two years, the simulated runoff was overestimated for both Jingle and Hecha.For 1997, this may be because the water stored in the watershed in two consecutive wet years (1995 and 1996) was released in the subsequent year.For 2000, the overestimation may have happened because 1999 was an extremely dry year, and water bodies such as small swamps and lakes were dry, and the soil was also very dry.This was followed by heavy rain in 2000 that firstly wetted the soil, filled the dry water bodies, and did not form surface  runoff.This study did not take into account impounding of water in small water bodies which could have resulted in the overestimation of the actual runoff in 2000.
All NSE and  2 values for the monthly calibration and validation were above 0.80, and PBIAS values were in the range of −15% to 0 (most PBIAS values were in the range of −10% to 0), suggesting very good model performance [30].The overall performance of the model was very good as shown in Figure 4.

Contribution of Individual LULC Changes to Watershed
Water Yield.The proportional contributions of five LULC changes (i.e., farmland, forest, grassland, water or wetland, and built-up land) and simulated water yield (WY) are shown in Figure 5. Grassland, farmland, and forest were the main LULC changes in the Upper Fenhe River watershed, and the sum of their areas accounted for 98% of the total area.Their contributions to water yield were up to about 98% but there was no obvious positive correlation between the area of individual LULC changes and their proportional contributions to watershed WY.The grassland class includes mainly moderately dense and sparse grassland and very little dense grassland, so its water-yielding capacity was large.Although the area of forested land was much less than that of the grassland, its contribution to watershed WY was slightly more than the grassland's, because forests in this watershed are mainly shrubs and sparse trees.In addition, the rain mostly came in the form of high-intensity storms, which could reduce the amount of canopy interception compared with light rain.This is similar to the results obtained by Hao et al. (2004) in the upper Luohe River watershed and Song and Ma (2008) in the Nanhe River basin in the Longxi Loess Plateau [6,7].The area of farmland was between that for forests and grasslands, but its water-yielding capacity was much less.The farmland in the Upper Fenhe River watershed is terraced and used to grow close and straight row-seeded wheat and maize which greatly reduce the water yield.
By comparing the percentages of the area of each land use type and their contribution to water yield, it was found that an increase of built-up land was the main cause of increases  of runoff, although the area of built-up land was very small.This is because its water-yielding capacity was the highest due to little initial loss and infiltration and rapid rainfall-runoff processes [9,31].
As LULCs changed from 1995 to 2010, the contribution of individual LULC change to watershed WY changed too, but the change process was complex, because the WY is an integrated result of LULCs, soil, topography, and climate.  5. Compared to the LULC change baseline year of 1995, the average annual WY over the watershed was 0.13 mm higher in 2010.Similar to WY, average annual PET with LULC in 1995 was 1026 mm; it increased to 1028 mm by 2010.The AET increased from 391 mm for LULC in 1995 to 393 mm for LULC in 2010.On the contrary, the average annual SW for LULC in 2010 was 0.71 mm lower than that in 1995.Furthermore, compared to the LULC baseline year 1995, the annual PET, AET, and SW for each year from 1992 to 2000 for LULC change in 2010 changed in a similar manner as the average annual PET, AET, and SW.But in contrast to the size of the average annual WY for LULC in 1995 and 2010, the annual WY for each year in this period did not increase consistently.For example, in years of abundant precipitation, such as 1995 and 1996, the annual WY for the LULC change in 1995 was higher than that in 2010.This indicated that the precipitation could affect the impacts of the LULC changes on the hydrology in this region, which was also found by Hao et al. (2004) in the upper Luohe River watershed [6].
The comparison of variations of PET and AET and changes in LULCs suggested that the increase of annual PET could be mainly attributed to returning cropland to forest and grassland and to urban expansion from 1995 to 2010.Further comparison between changes in WY and changes of LULCs indicated that an increase of WY was mainly due to increases of urbanization, which can increase the area of impervious surfaces, increase runoff, and decrease  infiltration [9,31].An association between the decreases of SW and forest, grassland, and built-up land expansion from 1995 to 2010 could be indicated from the comparison between variations of average annual SW and changes in LULC from 1995 to 2010.Expansion of forest, grassland, and built-up land by replacing them with farmland could promote water infiltration and drainage because of well-developed root systems and prevention of infiltration due to increases in the areas of impervious surfaces.

Impacts of LULC Changes on Hydrology at the Subbasin
Scale.Three subbasins (numbers 7, 17, and 35 in Figure 1) with obvious LULC changes were selected for analysis of the impacts of LULC changes on hydrology at the subbasin scale, and the results are shown in Table 6.In subbasin number 7, the area of farmland decreased sharply, and the areas of forest, grassland, water, and built-up land increased from 1995 to 2010.The PET, AET, and WY increased, and, in contrast, the SW decreased.In subbasin number 17 only the forested area decreased and grassland area increased a little, so the impacts on hydrology were that AET increased clearly and WY and SW decreased.In subbasin number 35, areas of farmland, forest, and built-up land increased and grassland decreased, and WY increased and SW decreased.By comparing the LULC changes and the hydrological responses to them in these three sun-basins, obvious LULC changes may not show obvious impacts on hydrology, which means different combinations of LULCs may produce similar hydrological effects.Furthermore, the impacts of the same LULC changes on hydrology may be different under the conditions of different precipitation intensity and distribution.

Conclusions
In this study we aimed to investigate the impacts of LULC changes on hydrology in China's Loess Plateau.Major research findings and their implication on practices and future research are as follows.
From 1995 to 2010, the main land use changes were the transformations of farmland to forests and from farmland to grassland, and the area of built-up land increased at the expense of farmland.These changes were due to the implementation of watershed management measures and social and economic development.For the contribution of each LULC to the total WY of the watershed, the forest, grassland, and farmland areas were the largest contributors, with up to about 39%, 38%, and 21%, respectively.The land use that produced the greatest water yield was shown to be built-up land, which was higher than any other land use, followed by forest, grassland, and farmland.
Accompanying the LULC changes in the Upper Fenhe River watershed, increases in PET, AET, and WY indicated soil and water conservation practices increased runoff, while the expansion of the forest and grassland increased the PET and AET.In addition, the consumption of SW was increased by vegetation combined with decreased infiltration because increases of the built-up land resulted in decreases in SW.Furthermore, the precipitation could affect the impacts of LULC changes on the hydrologic components in this arid and semiarid area.
The approach used in this study simply determined contributions of individual LULC classes to the total discharge, providing quantitative information for decision-makers to make better choices for land and water resource planning and management.This approach also provides a solid example of the potential of hydrologic modeling using remotely sensed digital LULCs in understanding the impacts of landscape change on water provisioning, a vital ecosystem service in the Loess Plateau of China.It can be widely applied to a variety of watersheds, where time-sequenced digital land cover data are available, and to predict hydrological consequences to LULC changes.

Figure 1 :
Figure 1: Locations of study area, hydrologic stations, and meteorological stations.

Figure 2 :
Figure 2: The conceptual framework for analyzing the impacts of LULC changes on the hydrology.
Land Use Classes.The calibrated and validated hydrological model was run for each of the LULC maps (1995 and 2010) with constant DEM and soil data, from January 1990 to December 2000 (11 years), to quantify contributions of individual LULC classes to the total discharge and to evaluate impacts of LULC changes on hydrology at the basin and subbasin scales.

4. 1 .FarmlandFigure 3 :
Figure 3: Maps for land uses and land cover in 1995 and 2010 in the Upper Fenhe River watershed.

Figure 5 :
Figure 5: Percentage of each LULC change and its contribution to the total average annual water yields (WY) for the Upper Fenhe River watershed from 1992 to 2000.

4. 4 .
Impacts of LULC Changes on Hydrology at the Basin Scale.The average basin values of P, PET, AET, WY, and SW simulated from each LULC map in 1995 and 2010 are shown in Table to investigate the LULC changes of the watershed from 1995 to 2010 with LULC maps determined by remote sensing images in 1995 and 2010; (2) to calibrate and validate the SWAT model in terms of monthly streamflow based on the DEM, soil map, land use in 1995, and climate data from 1990 to 2000; (3) to assess the hydrological effects of individual land uses, and (4) to simulate responses of hydrologic components to land use changes at basin and subbasin scales under the same climatic conditions.The results should assist decision-makers in target water resources planning and vegetation restoration on the Loess Plateau.

Table 1 :
Hydrometeorological stations in the study area.
dominant land use type in this watershed.The surface runoff is estimated using a modification of the SCS curve number method with daily rainfall amounts.Percolation is calculated using the combination of a storage routing technique and a crack flow model.The lateral flow is estimated simultaneously with percolation using a kinematic storage model.The groundwater flow (baseflow) into the channel is calculated from hydraulic conductivity of the shallow aquifer, distance from subbasin to main channel, and water table height.3.2.2.SWAT Model Inputs and Outputs.The input data used in the SWAT model included a digital elevation model (DEM), soil data, digital LULC maps, and climate data.The DEM was derived from the USGS National Elevation Dataset with a resolution of about 90 m, and the soil data were obtained from the Resources and Environment Science Data Center, Chinese Academy of Sciences, with the scale of 1 : 1,000,000.The LULC data for two times (1995 and 2010) were used to assess the impacts of LULC change on hydrology and were derived from Landsat TM images with a spatial resolution of 30 m using an interpretation method.The climate data, including daily values of precipitation and minimum-maximum temperature from January 1st, 1990, to December 31, 2000, were derived from 3 meteorological stations and 16 precipitation stations located in and around the Upper Fenhe River watershed (Table

Table 2 :
The warm-up, calibration, and validation period for the model.
performance of model calibration and validation.The Nash-Sutcliffe efficiency (NSE) is a normalized statistic that determines the relative magnitude of the residual variance compared to the measured data variance.NSE indicates how well the plot of observed versus simulated data fits the 1 : 1 line.Coefficient of determination (2

Table 3 :
Transition matrix of land uses from 1995 to 2010 (km 2 ).

Table 4 :
The results of parameter sensitivity analysis and calibration.

Table 5 :
Annual basin values of hydrologic features for the Upper Fenhe River watershed on the different land use from 1992 to 2000 (mm).

Table 6 :
LULC changes and average annual hydrologic features from 1992 to 2000 of three subbasins.