Simulation and Prediction of Climate Variability and Assessment of the Response of Water Resources in a Typical Watershed in China

The assessment of water resource responses to climate change is required in water resource planning and management, protecting environmental quality, and managing watersheds. This study modeled surface runoff and baseflow responses to variations in precipitation (0%, ±10%, and ±20%) and temperature (0 ◦C, ±1 ◦C, and ±2 ◦C) in 25 types of scenarios in the Lanhe Watershed (1140 km2), which possesses the typical hydrological and meteorological characteristics of the Loess Plateau in China. The study is based on the Soil and Water Assessment Tool (SWAT), which was calibrated and validated using the coefficient of determination (R2), Nash-Suttcliffe (Ens), and Percent bias (PBIAS), using the observed streamflow of Shangjingyou Station, a unique gauging station in the study area. The model was calibrated with daily streamflow, from 1967 to 1996, and then validated from 1997 to 2011. R2, Ens, and PBIAS were 0.95 and 0.84, 0.78 and 0.72, and 0.6% and −9.1% in annual and monthly calibration periods, 0.90 and 0.78, 0.74 and 0.67, and 22.1% and 18.8% in annual and monthly validation periods, and the overall performance ratings was “satisfactory”. The assessment indicates that surface runoff is likely to be more affected than baseflow when altering temperatures and precipitation, and the noticeable changes of surface and baseflow are from June to September and October to November, respectively. Results also indicate that surface runoff and baseflow are very sensitive to the projected reduction in temperature, rather than to an increase of temperature, while precipitation is a constant. In turn, when the temperature is a constant, the surface runoff is sensitive to the projected increase in precipitation and the baseflow is sensitive to the decrease in precipitation.


Introduction
The impacts of climate change on hydrology cycling and water resource quantity and timing have attracted more and more researchers to the fields of hydrology, geohydrology, ecology, et al. [1][2][3].Anticipating changes in water resources related to hydrological processes is particularly important for regions where water resources are scarce owing to population growth and economic prosperity [4][5][6][7][8], especially in regions of the upper and middle reaches of the Yellow River, where the shortage of water resources has resulted in a series of ecological and environmental issues, as well as the problem of social economy [9][10][11][12][13][14].
Most studies indicate that hydrological processes may be sensitive to changes in precipitation and temperature [15][16][17][18], which are the main climate variables.Changes of temperature have a strong impact on evapotranspiration, soil moisture, and the thermal circulation of the atmosphere, and, thus, on surface runoff and baseflow, as well as on crop production.The variability of precipitation directly changes the water resources of surface runoff, baseflow, and interflow, which then cause floods and droughts, which result from issues of social economy and the ecosystem and affect the sustainable use of water resources [19][20][21][22].Thus, temperature and precipitation are usually chosen as the main influencing factors in order to design climate variability scenarios [23,24].The climate change scenarios in most research are often projected using global or regional climate model experiments, of which a major limitation is large-scale resolution [25].In this study, based on the annual and monthly mean values (P for precipitation, T for temperature), the arbitrary values of precipitation (P, P ± 10%, and P ± 20%) and temperature (T, T ± 1 • C, and T ± 2 • C) were selected to simulate climate variability in order to bracket the range of possible precipitations and temperatures.
The main target of this paper is study of the water resources of the Lanhe Watershed under climate variability using the Soil and Water Assessment Tool (SWAT) model, which performs well all around the world, in different climatic situations, for solving hydrological environmental problems [26][27][28].The Lanhe Watershed is a typical watershed of the Loess Plateau, where the climate presents semi-arid characteristics, with most precipitation falling from June to September [29].Precipitation stations and the stream gauge station in the study area can offer long-term continuous data sets, from 1967 to 2011, including daily, monthly, and annual precipitation and discharge data as a basis.In this paper, the first aim is to calibrate and validate the SWAT model using the observed discharge from gauging stations.The second aim is to apply the model by inspecting the sensitivity of surface runoff and baseflow occurring with the change of temperature and precipitation, while noting the better performance of the model.

Study Area
The Lanhe Watershed, located in the upper and middle reaches of the Yellow River, is a representative region of the Loess Plateau.The Lanhe River is a secondary tributary of the Yellow River, and lies in the middle of Shanxi Province, China (Figure 1).The study area is about 1140 km 2 , the total river length is 57.6 km, and the mean of the river slope is 4.05‰.It has semi-arid and arid climate characteristics, with less precipitation in winter and spring, and more rain in summer and autumn, the precipitation of which accounts for more than 60% of the annual precipitation [30,31].The weather is affected by local and regional climatic influences, with northwestern winds in late autumn, winter, and spring, and southeastern winds in summer.The annual precipitation is 477.3 mm, the annual temperature is 8.8 • C, and the annual evapotranspiration is 1310.7 mm.In the Lanhe Watershed, about 80% of the annual precipitation falls from June to September (the flood season), and is usually concentrated in several rain events during July and August [32].Streamflow in the flood season accounts for 60%-80% of the annual total streamflow, and runoff generation is dominated by surface runoff, produced by the infiltration-excess mechanism; baseflow is commonly from 0.1 m 3 /s to 0.3 m 3 /s [33].The location of the Lanhe Watershed is shown in Figure 1.

SWAT (Soil and Water Assessment Tool) Model
SWAT was developed by the United States Department of Agricultural Research Service and Texas A&M University laboratories [34].It is a physically based, spatially distributed, continuous simulation watershed model, and it has been utilized to simulate the hydrological processes of interception, infiltration, surface runoff, evapotranspiration, percolation, lateral flow, and groundwater recharge, based on the water balance equation at a daily time step [35].The SWAT model also provides several software tools, which include the pre-and post-processing software, automatic and manual calibration procedures, and an ArcGIS SWAT interface, as well as abundant theoretical and technical documentation, which consist of theoretical documentation, a user's manual and guide, a developer's manual, etc. [36][37][38].Most of the software and documentation are conveniently and freely obtainable from [39].Thus, it has been widely used around the world for hydrological processes and soil movement, the assessment of water resources and water quality, and in the movement and transformation of nutrients (such as nitrogen, phosphorus) and pesticides [40,41].In recent years, the model has been increasingly applied in evaluating climate impacts on water resources at the basin scale, for small to large catchments [42][43][44][45][46][47].
In the SWAT model, both the Soil Conservation Service (SCS) curve number and the Green and Ampt infiltration method are developed to estimate surface runoff, and the three potential evapotranspiration (PET) methods, the Penman-Monteith method, the Priestley-Taylor method, and the Hargreaves method, are used to simulate PET [48,49].In addition, the SWAT model subdivides groundwater systems as shallow or deep aquifers.The shallow aquifer is an unconfined aquifer, returning groundwater to the streams, and a deep aquifer is a confined aquifer, contributing flow to the outside of the subbasin [35].The SCS curve number and the Penman-Monteith method were selected in this paper to calculate surface runoff and PET.

SWAT (Soil and Water Assessment Tool) Model
SWAT was developed by the United States Department of Agricultural Research Service and Texas A&M University laboratories [34].It is a physically based, spatially distributed, continuous simulation watershed model, and it has been utilized to simulate the hydrological processes of interception, infiltration, surface runoff, evapotranspiration, percolation, lateral flow, and groundwater recharge, based on the water balance equation at a daily time step [35].The SWAT model also provides several software tools, which include the pre-and post-processing software, automatic and manual calibration procedures, and an ArcGIS SWAT interface, as well as abundant theoretical and technical documentation, which consist of theoretical documentation, a user's manual and guide, a developer's manual, etc. [36][37][38].Most of the software and documentation are conveniently and freely obtainable from [39].Thus, it has been widely used around the world for hydrological processes and soil movement, the assessment of water resources and water quality, and in the movement and transformation of nutrients (such as nitrogen, phosphorus) and pesticides [40,41].In recent years, the model has been increasingly applied in evaluating climate impacts on water resources at the basin scale, for small to large catchments [42][43][44][45][46][47].
In the SWAT model, both the Soil Conservation Service (SCS) curve number and the Green and Ampt infiltration method are developed to estimate surface runoff, and the three potential evapotranspiration (PET) methods, the Penman-Monteith method, the Priestley-Taylor method, and the Hargreaves method, are used to simulate PET [48,49].In addition, the SWAT model subdivides groundwater systems as shallow or deep aquifers.The shallow aquifer is an unconfined aquifer, returning groundwater to the streams, and a deep aquifer is a confined aquifer, contributing flow to the outside of the subbasin [35].The SCS curve number and the Penman-Monteith method were selected in this paper to calculate surface runoff and PET.

Data Collection and Analysis
Digital Elevation Model (DEM) data using a 30 m resolution were downloaded from the Computer Network Information Center of the Chinese Academy of Sciences.It was used to delineate the subwatershed and stream network, as well as to derive the sub-watershed parameters of the area, length, width, mean slope, and average elevation.
For both land use and soil data, the resolution was 1000 m, and was acquired from the Environmental and Ecological Data Center of the Chinese Academy of Sciences.The main land use types were evergreen coniferous forest, deciduous broad-leaved forest, shrubwood, grassland, and farmland; the main soil types were calcareous cambisol, eutric cambisol, calcareous fluvisol, haplustalf, calcareous alfisol, and calcareous rhogosol.
The meteorological data of the Lanhe Watershed were obtained from the China Meteorological Data Sharing Service System.The daily climate data included maximum and minimum temperatures, wind speed, and relative humidity, as well as daily sunshine duration, from 1967 to 2011.The mean annual values of temperature, evapotranspiration, wind speed, relative humidity, and sunshine duration were 8.8 • C, 1310.7 mm, 2.1 m/s, 59%, and 2542.3 h, respectively.
The precipitation and streamflow data were obtained from the Hydrology and Water Resource Survey Administration of Shanxi, China, and contained daily precipitation and runoff.There are seven precipitation stations (Pingshang, 21350; Puming, 21250; Yanjiagou, 21300; Baijiazhuang, 21450; Louzi, 21500; Caocheng, 21600; and Shangjingyou, 21650), which are evenly distributed in the study area, and one stream gauge station (Shangjingyou), which is the outlet of the Lanhe catchment.The longest time span of meteorological and hydrological data (45-year period, from 1967 to 2011) was used for the research.The average annual precipitation was 477.3 mm, with the maximum and minimum rainfalls being 816.9 mm and 266.5 mm; the mean annual flow was estimated to be about 43.41 mm 3 /year, 38.1 mm for runoff depth, and 0.07 for the mean runoff coefficient; the maximum and minimum temperatures were 10.4 • C and 7.5 • C, respectively.

Characterization and Delineation of Hydrological Response Units
In this study, model setup and parameterization were done using the ArcSWAT interface.First, the Lanhe Watershed was divided into 25 subbasins, decided by a threshold area of 1000 ha (Figure 2) with DEM data of 30 m in resolution.Then, all subbasins were fully divided into 81 hydrologic response units (HRUs), based on a combination of land use, soil class (type), and land slope, representing over 15% threshold of land use, soil, and slope, respectively.

Model Calibration and Validation
The determination of the most sensitive parameters is the key, and first step, for model calibration and validation at the watershed scale [50].The sensitivity analysis can classify the sensitive parameters, from most sensitive to least sensitive for the input parameters.Within the model, a number of parameters are specified in order to identify the spatial units of the subbasins, and the hydrologic response units (HRUs) may be assessed from available data.Latin Hypercube One-factor-At-a-Time (LH-OAT) was used for a preliminary sensitivity analysis of the parameters [51].The Shuffle Complex Evolution Method (SCEM-UA), proposed by Song et al. [52], was then employed for a formal calibration using the streamflow data of Shangjingyou Gauge Station, from 1967 to 2011.The coefficient of determination (R 2 ), Nash-Suttcliffe (Ens), and Percent bias (PBIAS) were used to assess model performance.The coefficient of determination (R 2 ), Nash-Suttcliffe (Ens), and Percent bias (PBIAS), defined by Equations ( 1)-( 3), were done using a statistics performance for the model [53,54].R 2 commonly describes the degree of agreement between simulated values and measured values.A value of 1.0 indicates a perfect prediction, and a value of 0.0 suggests an unsatisfactory prediction.In general, when a value is greater than 0.5 it is considered to be acceptable [55].Ens indicates how well the plot of observed versus simulated data fits the 1:1 line.A value of 0.75 < Ens ≤ 1.00 is very good, 0.65 < Ens ≤ 0.75 is good, 0.5 < Ens ≤ 0.65 is satisfactory, and Ens ≤ 0.5 is unsatisfactory [56].
where: n is the total number of observations.Qi obs is the observed data; Qi avg is the average of the observed data.Pi sim is the simulated data; Pi avg is the average of the simulated data.

Scenarios of Climate Change
Precipitation and temperature variables significantly affect the hydrological process; thus, in this study, they were selected in order to design climate change scenarios with the arbitrary scenarios method, based on the tendencies of temperature and precipitation in the study area.

Design of the Climate Scenarios
Historical temperature and precipitation data (from 1967 to 2011) were collected in order to investigate climatic change and the hydrometeorological characteristics of the study area for the design of climate change scenarios.The Mann-Kendall test (M-K test) assessed the trends and abrupt changes in a time series [59], and linear trend analyses were used for temperature analyses [60]; and Ens indicates how well the plot of observed versus simulated data fits the 1:1 line.A value of 0.75 < Ens ≤ 1.00 is very good, 0.65 < Ens ≤ 0.75 is good, 0.5 < Ens ≤ 0.65 is satisfactory, and Ens ≤ 0.5 is unsatisfactory [56].
where: n is the total number of observations.Q i obs is the observed data; Q i avg is the average of the observed data.P i sim is the simulated data; P i avg is the average of the simulated data.

Scenarios of Climate Change
Precipitation and temperature variables significantly affect the hydrological process; thus, in this study, they were selected in order to design climate change scenarios with the arbitrary scenarios method, based on the tendencies of temperature and precipitation in the study area.

Design of the Climate Scenarios
Historical temperature and precipitation data (from 1967 to 2011) were collected in order to investigate climatic change and the hydrometeorological characteristics of the study area for the design of climate change scenarios.The Mann-Kendall test (M-K test) assessed the trends and abrupt changes in a time series [59], and linear trend analyses were used for temperature analyses [60]; and the five-year moving average and accumulating anomalies were applied to the precipitation analysis.
UF and UB are standardized statistical series which can be written as Equations ( 4) and ( 5) in the M-K test [61].When the UF curve and the UB curve intersect in the confidence zone, the cross point indicates that it is an abrupt point [61].
Figures 3 and 4 show temperature tendency and change, and Figure 5 demonstrates precipitation tendency and change.According to Figures 3 and 4, the mean annual temperature presented a wave of increasing tendency, which experienced a sudden change, which was the intersection point between UF and UB curves in 1990.The annual average temperature, maximum temperature, and minimum temperature were 8.8 • C, 10.4 • C, and 7.5 • C, and the rate was 0.29 • C every 10 years.From Figure 5, the precipitation for a five-year moving average anomaly showed a decreasing tendency.In the periods of 1967-1972, 1979-1984, and 1993-2006, precipitation presented a decreasing tendency.However, in 1973-1978, 1985-1992, and 2007-2011, the precipitation indicated an increasing tendency.The annual average, maximum, and minimum precipitations were 477.3 mm, 816.9 mm (1967), and 266.5 mm (1972).
Water 2016, 8, 490 6 of 15 UF and UB are standardized statistical series which can be written as Equations ( 4) and ( 5) in the M-K test [61].When the UF curve and the UB curve intersect in the confidence zone, the cross point indicates that it is an abrupt point [61].
Figures 3 and 4 show temperature tendency and change, and Figure 5 demonstrates precipitation tendency and change.According to Figure 3 and Figure 4, the mean annual temperature presented a wave of increasing tendency, which experienced a sudden change, which was the intersection point between UF and UB curves in 1990.The annual average temperature, maximum temperature, and minimum temperature were 8.8 °C, 10.4 °C, and 7.5 °C, and the rate was 0.29 °C every 10 years.From Figure 5, the precipitation for a five-year moving average anomaly showed a decreasing tendency.In the periods of 1967-1972, 1979-1984, and 1993-2006, precipitation presented a decreasing tendency.However, in 1973-1978, 1985-1992, and 2007-2011, the precipitation indicated an increasing tendency.The annual average, maximum, and minimum precipitations were 477.3 mm, 816.9 mm (1967), and 266.5 mm (1972).From previous investigations, average annual temperatures will increase from 1.7 °C to 2.3 °C, and the average annual precipitation will increase from 5% to 23%, according to a report on comprehensive environmental evolution assessment in the west of China.Meanwhile, considering the hydrometeorological characteristics of the study area, generally, the annual precipitation varied from −20% to +20%, and the annual mean temperature varied from −2 °C to 2 °C.From previous investigations, average annual temperatures will increase from 1.7 • C to 2.3 • C, and the average annual precipitation will increase from 5% to 23%, according to a report on comprehensive environmental evolution assessment in the west of China.Meanwhile, considering the  From previous investigations, average annual temperatures will increase from 1.7 °C to 2.3 °C, and the average annual precipitation will increase from 5% to 23%, according to a report on comprehensive environmental evolution assessment in the west of China.Meanwhile, considering the hydrometeorological characteristics of the study area, generally, the annual precipitation varied from −20% to +20%, and the annual mean temperature varied from −2 °C to 2 °C.As a result, climate scenarios were formulated by changing the temperatures and precipitations of the annual mean.Temperature was constructed by altering the annual average temperature (8.8 °C) in steps of 1 °C for ±2 °C (such as +10.8 °C, +9.8 °C, 8.8 °C, 7.8 °C, and 6.8 °C), and the annual mean precipitation (477.3 mm) changing by steps of 10% for ±20% (such as 572.76 mm, 525.03 mm, 477.3 mm, 429.57mm, and 381.84 mm).Twenty-five scenarios of climate change were formulated using the arbitrary scenario method for SWAT model simulation (see Table 1).The base scenario is S33, P being the annual average precipitation and T being the annual average temperature.

Methods of Sensitivity Analysis
For assessing climate change impacts on water resources, the rate of change of hydrological factors η ， , such as surface runoff and baseflow were defined using Equation (6).As a result, climate scenarios were formulated by changing the temperatures and precipitations of the annual mean.Temperature was constructed by altering the annual average temperature (8.8 • C) in steps of 1  1).The base scenario is S33, P being the annual average precipitation and T being the annual average temperature.

Methods of Sensitivity Analysis
For assessing climate change impacts on water resources, the rate of change of hydrological factors η ∆P,∆T , such as surface runoff and baseflow were defined using Equation (6).
η ∆P,∆T = W P+∆P,T+∆T − W P,T W P,T × 100 W P,T is the value of the hydrological factors obtained for the surface runoff and baseflow using simulations of the SWAT model in the base scenarios.W P+∆P,T+∆T is the simulation result of surface runoff and baseflow, changing the annual average precipitation and temperature, by steps of ∆P (10%) for precipitation, or ∆T (1 • C) for temperature.

SWAT Model Calibration and Validation Results
The results of model calibration and validation are shown in Table 2 and Figures 6 and 7. From the calibration and validation results of annual and monthly streamflows, R 2 values are 0.95 and 0.84, and 0.78 and 0.72, and all results greater than 0.5 are considered to be acceptable; Ens values are 0.78 and 0.72, and 0.74 and 0.67, and the simulated results of the ratings are "very good" or "good"; PBIAS values are 0.6% and −9.1%, and 22.1% and 18.8%, with a performance rating of "very good" for the calibration period, and a rating of "satisfactory" for the validation period.According to the evaluation criteria, the performance rating of the overall model simulated was "satisfactory".Figures 6 and 7 show that the measured and simulated annual and monthly streamflow were closely followed most of time.To sum up, the statistic values proved that the model performed well, not only on an annual scale but also on a monthly scale in the Lanhe Watershed, and that it can be used to predict water resources (surface runoff and baseflow).evaluation criteria, the performance rating of the overall model simulated was "satisfactory".Figures 6 and 7 show that the measured and simulated annual and monthly streamflow were closely followed most of time.To sum up, the statistic values proved that the model performed well, not only on an annual scale but also on a monthly scale in the Lanhe Watershed, and that it can be used to predict water resources (surface runoff and baseflow).From the calibration and validation results, the optimal parameters of calibration and validation are listed in Table 3. From the calibration and validation results, the optimal parameters of calibration and validation are listed in Table

Sensitivity Analysis of Surface Runoff
Table 4 shows the annual surface runoff obtained by the SWAT model with respect to the different temperatures and precipitations in each scenario.Figure 8 shows the rate of change relative to the base scenario (S33), of which the value is 43,411 × 10 3 m 3 .From Table 4 and Figure 8, both the precipitation and temperature variability changed the surface runoff, and the contribution of precipitation variability was far greater than that of temperature change.Comparing the surface flow of S33 with other scenarios, the smallest relative change of the surface flow decreased by 40.8% (17,729 × 10 3 m 3 ) when the precipitation decreased by 20% and the temperature increased by 2 • C (S51 scenario).The smallest relative change in the surface flow increased by 78.9% (34,268.7 × 10 3 m 3 ) when the precipitation increased by 20% and the temperature decreased by 2 • C (S15).
surface runoff, and the contribution of precipitation variability was far greater than that of temperature change.Comparing the surface flow of S33 with other scenarios, the smallest relative change of the surface flow decreased by 40.8% (17,729 × 10 3 m 3 ) when the precipitation decreased by 20% and the temperature increased by 2 °C (S51 scenario).The smallest relative change in the surface flow increased by 78.9% (34,268.7 × 10 3 m 3 ) when the precipitation increased by 20% and the temperature decreased by 2 °C (S15).With respect to the same precipitation, the relative changes of surface runoff increased in the range of 1.5% (645 × 10 3 m 3 , S11-S21) to 15.1% (6565.8× 10 3 m 3 , S15-S25), while the temperature decreased by 1 °C, and the relative changes in surface runoff decreased in the range of 1.6% (708.4 × With respect to the same precipitation, the relative changes of surface runoff increased in the range of 1.5% (645 × 10 3 m 3 , S11-S21) to 15.1% (6565.8× 10 3 m 3 , S15-S25), while the temperature decreased by 1 • C, and the relative changes in surface runoff decreased in the range of 1.6% (708.4 × 10 3 m 3 , S42-S52) to 9% (3900.7 × 10 3 m 3 , S45-S55) while the temperature increased by 1 • C, compared to the S33 scenario.The results show that the contribution of the reduction in temperature was greater than that of the increase in temperature for surface runoff.
With respect to situations where the same temperature is used, the relative changes of surface runoff increased in the range of 20.3% (8791.9× 10 3 m 3 , S54-S53) to 34.3% (14,910.2× 10 3 m 3 , S15-S14), while the precipitation change increased by 10%, and the relative changes of surface runoff decreased in the range of 14.3% (6205.2× 10 3 m 3 , S42-S41) to 26.0% (11,293.9 × 10 3 m 3 , S13-S12) while the precipitation decreased by 10%, relative to the S33 scenarios.The performance indicated that the contribution of increasing the relative precipitation variation rates was greater than that of the decrease for surface runoff.
In order to learn more about the changes in surface runoff and baseflow within a year, it is necessary to analyze periodic features of surface runoff and baseflow.In this paper, and based on flood and irrigation water management, the hydrological year was divided into four periods: the first period (the pre-crop period) from March to May, the second period (the flood season) from June to September, the third period (the post-flood season) from October to November, and the fourth period (the ice period) from December to February.Monthly surface flow changes were simulated under four scenarios: S31, (1 − 20%)P, T; S33, P, T; S51, (1 − 20%)P, T + 2 • C; and S53, P, T + 2 • C.
Figure 9 depicts that the surface runoff of S33 is greater than that of the S31, S51, and S53 scenarios, with a significant change in the average monthly surface water being in the flood season, with August showing the largest decrease, and a slight change in the ice period.The changes were in accordance with the monthly changes in precipitation.It should also be noted that the surface runoff reduction was greater under a 20% decrease in precipitation rather than a 2 • C increase in temperature for the flood season, and there were slight variations relative to other periods.
Figure 9 depicts that the surface runoff of S33 is greater than that of the S31, S51, and S53 scenarios, with a significant change in the average monthly surface water being in the flood season, with August showing the largest decrease, and a slight change in the ice period.The changes were in accordance with the monthly changes in precipitation.It should also be noted that the surface runoff reduction was greater under a 20% decrease in precipitation rather than a 2 °C increase in temperature for the flood season, and there were slight variations relative to other periods.

Sensitivity Analysis of Baseflow
Table 5 presents the annual baseflow obtained by the SWAT model, in accordance with each scenario.The baseflow of S33 was 6463.8 × 10 3 m 3 (runoff depth 5.67 mm), occurring with the annual average precipitation and temperature.

Sensitivity Analysis of Baseflow
Table 5 presents the annual baseflow obtained by the SWAT model, in accordance with each scenario.The baseflow of S33 was 6463.8 × 10 3 m 3 (runoff depth 5.67 mm), occurring with the annual average precipitation and temperature.Table 5 and Figure 10 show the differences and rate of change compared with the base S33 scenario.In contrast to the S33 scenario, the minimum relative change of the baseflow decreased by 26.8% (4731.5 × 10 3 m 3 , S51), while the precipitation decreased by 20% and the temperature increased by 2 • C. The maximum relative change of the baseflow increased by 30% (8402.9× 10 3 m 3 , S15), while the precipitation increased by 20% and the temperature decreased by 2   10 show the differences and rate of change compared with the base S33 scenario.In contrast to the S33 scenario, the minimum relative change of the baseflow decreased by 26.8% (4731.5 × 10 3 m 3 , S51), while the precipitation decreased by 20% and the temperature increased by 2 °C.The maximum relative change of the baseflow increased by 30% (8402.9× 10 3 m 3 , S15), while the precipitation increased by 20% and the temperature decreased by 2 °C, respectively.For the same precipitation, when the temperature decreases by 1 °C, the baseflow will increase from 362 × 10 3 m 3 to 465.4 × 10 3 m 3 , by 5.6% to 7.2%; when the temperature increases by 1 °C, the baseflow will decrease from 271.5 × 10 3 m 3 to 381.4 × 10 3 m 3 , by 4.2% to 5.9%.
For the same temperature, the relative changes in baseflow caused by decreasing precipitation are greater than those caused by an increasing precipitation.For instance, when the precipitation For the same precipitation, when the temperature decreases by 1 • C, the baseflow will increase from 362 × 10 3 m 3 to 465.4 × 10 3 m 3 , by 5.6% to 7.2%; when the temperature increases by 1 • C, the baseflow will decrease from 271.5 × 10 3 m 3 to 381.4 × 10 3 m 3 , by 4.2% to 5.9%.
For the same temperature, the relative changes in baseflow caused by decreasing precipitation are greater than those caused by an increasing precipitation.For instance, when the precipitation change decreases by 10%, the baseflow will decrease from 530 × 10 3 m 3 to 652.8 × 10 3 m 3 , by 8.2% to 10.1%; when the precipitation change increases by 10%, the baseflow will increase from 491.3 × 10 3 m 3 to 562.3 × 10 3 m 3 , by 7.6% to 8.7%.
Results also show that climate change has more of a contribution to the surface runoff relative to the baseflow.
The monthly baseflow changes are shown in Figure 11.The overall results for each operative scenario (S31, S33, S51, S53) showed that the baseflow yield changes implied significant changes in the post-flood season (October and November), with October showing the largest difference, and slight differences for the ice period (from December to February) and the pre-crop period (from March to May).The baseflow changes appeared to be a time-lag phenomenon with respect to the monthly changes in precipitation.Furthermore, it should also be noted that the contribution of precipitation variability was far greater than that of temperature change.For the same precipitation, when the temperature decreases by 1 °C, the baseflow will increase from 362 × 10 3 m 3 to 465.4 × 10 3 m 3 , by 5.6% to 7.2%; when the temperature increases by 1 °C, the baseflow will decrease from 271.5 × 10 3 m 3 to 381.4 × 10 3 m 3 , by 4.2% to 5.9%.
For the same temperature, the relative changes in baseflow caused by decreasing precipitation are greater than those caused by an increasing precipitation.For instance, when the precipitation change decreases by 10%, the baseflow will decrease from 530 × 10 3 m 3 to 652.8 × 10 3 m 3 , by 8.2% to 10.1%; when the precipitation change increases by 10%, the baseflow will increase from 491.3 × 10 3 m 3 to 562.3 × 10 3 m 3 , by 7.6% to 8.7%.
Results also show that climate change has more of a contribution to the surface runoff relative to the baseflow.
The monthly baseflow changes are shown in Figure 11.The overall results for each operative scenario (S31, S33, S51, S53) showed that the baseflow yield changes implied significant changes in the post-flood season (October and November), with October showing the largest difference, and slight differences for the ice period (from December to February) and the pre-crop period (from March to May).The baseflow changes appeared to be a time-lag phenomenon with respect to the monthly changes in precipitation.Furthermore, it should also be noted that the contribution of precipitation variability was far greater than that of temperature change.

Conclusions
This study used the SWAT model, which was calibrated and validated with R 2 , Ens, and PBIAS, and presented the impacts of climate variability on surface runoff and baseflow in the Lanhe Watershed, which has similar hydrological and meteorological characteristics to those of the Loess Plateau.In the calibration period , the annual and monthly R 2 , Ens, and PBIAS were 0.95 and 0.84, 0.78 and 0.72, and 0.6% and −9.1%; in the validation period (1997-2011), the annual and monthly R 2 , Ens, and PBIAS were 0.90 and 0.78, 0.74 and 0.67, and 22.1% and 18.8%.The performance rating of the calibrated and validated SWAT model was better for the prediction variability of surface runoff and baseflow with variability in precipitation and temperature.
The simulated results indicated that the surface runoff was likely to be more affected than baseflow by changes in precipitation and temperature.Additionally, there were distinct variations in seasonal and monthly surface runoff and baseflows.Surface runoff had noticeable changes during the flood season, from June to September, and baseflow had significant variability in the post-flood season, from October to November.Furthermore, the surface runoff and the baseflow had slight variations relative to other periods.
When the precipitation was a constant, the decrease in temperature was more significant than an increase of the surface runoff, and the same results were found for baseflow.Nevertheless, when the temperature was a constant, an increase in the rate of precipitation was more significant than a decrease of the surface runoff; the decrease of precipitation was more remarkable than the increase in baseflow.

Figure 3 .
Figure 3.The M-K test for average annual temperature in the Lanhe Watershed.

Figure 3 .
Figure 3.The M-K test for average annual temperature in the Lanhe Watershed.
the study area, generally, the annual precipitation varied from −20% to +20%, and the annual mean temperature varied from −2 • C to 2 • C.

Figure 3 .
Figure 3.The M-K test for average annual temperature in the Lanhe Watershed.

Figure 4 .
Figure 4. Average annual temperature and tendencies in the Lanhe Watershed.

Figure 5 .
Figure 5.The annual average precipitation departure and accumulated anomalies of the Lanhe Watershed.

Figure 6 .
Figure 6.Observed and simulated annual streamflows at Shangjingyou Station.Figure 6. Observed and simulated annual streamflows at Shangjingyou Station.

Figure 6 .
Figure 6.Observed and simulated annual streamflows at Shangjingyou Station.Figure 6. Observed and simulated annual streamflows at Shangjingyou Station.

Figure 7 .
Figure 7. Observed and simulated monthly average streamflows at Shangjingyou Station.

Figure 7 .
Figure 7. Observed and simulated monthly average streamflows at Shangjingyou Station.

Figure 8 .
Figure 8. Rate of annual surface runoff change from base and other scenarios.

Figure 8 .
Figure 8. Rate of annual surface runoff change from base and other scenarios.

Figure 10 .
Figure 10.Rate of annual baseflow change of the base and other scenarios.

Figure 10 .
Figure 10.Rate of annual baseflow change of the base and other scenarios.

Figure 10 .
Figure 10.Rate of annual baseflow change of the base and other scenarios.

Table 1 .
Climate scenarios designed for the Lanhe Watershed.

Table 1 .
Climate scenarios designed for the Lanhe Watershed.

Table 2 .
Performance statistics in the calibration and validation periods.

Table 2 .
Performance statistics in the calibration and validation periods.

Table 3 .
Results of parameter calibration.

Table 4 .
Annual surface runoff simulated in different scenarios.

Table 5 .
Annual baseflow simulated in different scenarios.

Table 5 .
Annual baseflow simulated in different scenarios.

Table 5 and
Figure