Impact Assessments of Rainfall – Runo ff Characteristics Response Based on Land Use Change via Hydrological Simulation

The hydrology response was studied considering the established fact of land use change in Dapoling basin. The whole period was divided into two (1965–1985 and 1986–2012) according to the major land use and land cover change in this region. Xinanjiang model was used to simulate discharge data in the two periods. The hydrologic response to the change could be evaluated by inspecting the response of model parameters and flood elements. The results show that the lag time varied, and the hydrologic elements including the mean runoff depth, flood peak and kurtosis coefficient varied with the rainfall depth. This result is significant for studying the response of runoff characteristic from land use and land cover change.


Introduction
Land use and land cover patterns change dramatically with the expansion of human activities, the impact of which on environment and sustainable development has attracted widespread attention [1][2][3][4].The mechanism of hydrologic process is complex, which relates not only to the spatial and temporal distribution of basin physiographic factors, but also to the basin's underlying surface conditions.Vörösmarty [2] showed that land use change and climate change are two factors that impact hydrology.On a long time scale, climate change impacts more than land use change.However, on a short time scale, land use change impacts more.The existing research [5][6][7] shows that the major land use variations that affect hydrology are afforestation and deforestation.Other significant factors influencing hydrology include agricultural development, wetland drainage, urbanization, land abandonment and changes due to forest fire.The immediate and direct impacts of these land use changes on water yield and water quality of basin runoff are multiple and different in temporal and spatial scale [5].
The hydrologic response to land use change reveals that land use change alters elements such as evapotranspiration [8,9], runoff depth [10,11], mean annual runoff [12][13][14], snowmelt [1,15], soil erosion [16,17] and the groundwater quality [18] to affect overall basin hydrology.It is well documented that land use change therefore can have a significant effect on river basin hydrology, which is mainly reflected on the change of the flood volume [1,9], peak flow discharge [10,13], flood frequency [19] and peak flow rate [20].De Roo et al. [21] explored the impact of land use change on flood in the Oder (The Czech Republic, Poland and Germany) and Meuse (France, Belgium, Germany, and The Netherlands) basins based on the distributed hydrological model LISFLOOD, which was developed by the flood group of the Natural Hazards Project of the Joint Research Centre of the European Commission.They concluded that the overall effects during the period experiencing land use change are: the flood risk of the basin tends to increase, soil storage capacity has decreased 12 mm before flood comes, the flood peak discharge has increased 0.2%, and water level has increased 1 cm.Mao et al. [22] studied the hydrologic response of land use change in the Great Lakes region of America using the Variable Infiltration Capacity (VIC) model.At regional scales, the simulated results due to land use change vary spatially and seasonally, which is largely related to vegetation conversion and the geographic location of each land use type.Compared to the condition before land use change, the deforestation and conversion of deciduous forest into woodland and row-crop agriculture in the research center leads to a reduced evapotranspiration and the increased total runoff; the conversion of majority evergreen into majority deciduous forest in the northern part leads to a reduced evapotranspiration and the increased total runoff; and the conversion of prairie grasslands into row agriculture in the southern and western part leads to an increased evapotranspiration and the reduced total runoff.In addition, the land degradation and vegetation cover changes affect soil quality [18].The authors [18] found that deforestation and agricultural expansion have significant impact on groundwater quality in Northern Kelantan, Malaysia.They tested the nitrate values and found that the concentrations show a significant positive trend in both the agricultural and residential wells, while forests and grasslands show an annual decrease.
Lag time is a fundamental descriptor of basin characteristics, reflecting the time difference between the rainfall and the resulting hydrograph [23].It is a comprehensive index related with some internal processes of the catchment and gives the information of the water storage, flow path, water resources quantity and concentration time.The calculation of lag time can reflect the information of the time of flood peak discharge, which is important for basin flood management.Lag time is the fingerprint of water that reflects the physiography and flow routing of basin, which changes conspicuously with the change of basin land surface [24].Previous research on lag time mainly considered the estimation of flood peak and the hydrograph shape [25,26], the comparison of the variation of slope flow and river flow [27], the impact of different urbanization degree on lag time [23,24,28], and the basis of choosing a flood event for rainfall-runoff model [29,30].
Recently, the studies of land use change in the Huai River basin mainly focused on runoff generation process [31], evapotranspiration [32], water yield [33,34], rainfall-runoff relationship [33], water quality [35], pollution characteristics [36] and sediment yield characteristics [37].Research on the impact of land use change on basin lag time in the Huai River basin is scarce.More understanding of the influence of land use change on lag time can provide guiding significance for flood forecasting and water resources management of the basin.Determining the lag time of a basin can provide significant guidance for forecasting the time of peak discharge for this basin and the downstream basins, which makes the search meaningful.
Based on previous studies and the available data of hydrology and land use, we chose Dapoling basin for this research.The land use change is an established fact of the selected basin.The mechanism of how land use change influences the hydrology has been widely discussed in other studies, thus the object of this research was to test the hydrology response based on the fact of land use change and provide guidance for flood forecasting.Dapoling basin is the head source of Huai River.Its hydrology condition is important for downstream basins.

Study Area
The authors of [38,39] used a detection method based on hydrology model to test the hydrology response to land use change in Xixian basin.Our study area was the head source sub-basin of Xixian.The sub-basin Dapoling is located in the upstream of the Huai River at the latitudes 113    The rainfall station of Dapoling basin is shown in Table 1 and Figure 1.According to the observed data in Table 1, the maximum water level stage is 104.86 m and the corresponding maximum discharge is 4200 m 3 /s at Dapoling station.The mean annual precipitation is 918 mm (derived from 1964 to 2005), half of which falls between June and September.The mean annual runoff depth is 375 mm.The rainfall station of Dapoling basin is shown in Table 1 and Figure 1.According to the observed data in Table 1, the maximum water level stage is 104.86 m and the corresponding maximum discharge is 4200 m 3 /s at Dapoling station.The mean annual precipitation is 918 mm (derived from 1964 to 2005), half of which falls between June and September.The mean annual runoff depth is 375 mm.

Land Use Data
The land use data are from Resource and Environment Data Cloud Platform of Chinese Academy of Sciences (Table 2 and Figure 2).The data show that the main land use type in the basin are farmland and woodland, accounting for 41.85% and 38.55% in the 1980s, 30.38% and 40.68% in the 1990s, and 41.81% and 38.06% in the 2000s, respectively.The principal variation is the cross change of farmland and paddy field.Compared with the 1990s, the paddy field reduced to 17.15% from 27.23%, and the farmland increased to 41.81% from 30.38%.Compared with the 1980s, the paddy field increased to 27.23% from 17.02%, and the farmland reduced to 30.38% from 41.85%.2).The data show that the main land use type in the basin are farmland and woodland, accounting for 41.85% and 38.55% in the 1980s, 30.38% and 40.68% in the 1990s, and 41.81% and 38.06% in the 2000s, respectively.The principal variation is the cross change of farmland and paddy field.Compared with the 1990s, the paddy field reduced to 17.15% from 27.23%, and the farmland increased to 41.81% from 30.38%.Compared with the 1980s, the paddy field increased to 27.23% from 17.02%, and the farmland reduced to 30.38% from 41.85%.

Xinanjiang Model
A physically-based hydrological model is an important tool for hydrology research and water resources management [40].Many hydrological models have been developed and applied to watersheds with various sizes.The Xinanjiang Model is widely used in China.It was developed by Zhao in 1973 (Hohai University, China) [41].The model is extensively used for runoff simulation in humid and semi-humid basins in China and has achieved high accuracy.The Xinanjiang Model is workable for variable input data from the minute scale to the daily scale and calculates the outflow at the corresponding time scale.Thus, it is suitable for calculating the flood discharge in small and medium watersheds with flood propagation time less than 24 h.The input data of Xinanjiang Model are mainly the rainfall and evaporation data.It does not consider Digital Elevation Model (DEM) data.Comparing to other distributed hydrological models, which need geographic data, Xinanjiang model only demands hydrology data, which makes the model easy to use.

Xinanjiang Model
A physically-based hydrological model is an important tool for hydrology research and water resources management [40].Many hydrological models have been developed and applied to watersheds with various sizes.The Xinanjiang Model is widely used in China.It was developed by Zhao in 1973 (Hohai University, China) [41].The model is extensively used for runoff simulation in humid and semi-humid basins in China and has achieved high accuracy.The Xinanjiang Model is workable for variable input data from the minute scale to the daily scale and calculates the outflow at the corresponding time scale.Thus, it is suitable for calculating the flood discharge in small and medium watersheds with flood propagation time less than 24 h.The input data of Xinanjiang Model are mainly the rainfall and evaporation data.It does not consider Digital Elevation Model (DEM) data.Comparing to other distributed hydrological models, which need geographic data, Xinanjiang model only demands hydrology data, which makes the model easy to use.
Considering the uneven distribution of rainfall and the underlying surface, the study basin is divided into several sub-units by the Thiessen polygon method in ArcGIS [42].The method assumes that the data of rainfall station located at the center of each Thiessen polygon could represent the areal mean rainfall in sub-units.
The evapotranspiration, the runoff production, the separation of runoff components and the flow concentration of hill slope are first simulated separately in each divided sub-unit successively, and then flood routing down the channels from the sub-unit outlets to the main basin outlet is obtained by applying the Muskingum method.The main characteristic of the model is the concept of runoff formation on repletion of storage, which means that runoff will not be produced until the soil moisture content of aeration zone reaches the field capacity, and thereafter the total runoff equals to the total rainfall.To describe the non-uniform distribution of tension water capacity in sub-basins, a tension water capacity curve is introduced (Figure 3): where f represents the area of runoff production; F represents the area of basin; W' represents the tension water capacity at a point; and WMM is the maximum value of W'.The areal mean tension water capacity, WM, is easy to obtain from Equation (2) by integration.B is the parameter shown in Table 3.
Considering the uneven distribution of rainfall and the underlying surface, the study basin is divided into several sub-units by the Thiessen polygon method in ArcGIS [42].The method assumes that the data of rainfall station located at the center of each Thiessen polygon could represent the areal mean rainfall in sub-units.
The evapotranspiration, the runoff production, the separation of runoff components and the flow concentration of hill slope are first simulated separately in each divided sub-unit successively, and then flood routing down the channels from the sub-unit outlets to the main basin outlet is obtained by applying the Muskingum method.
The main characteristic of the model is the concept of runoff formation on repletion of storage, which means that runoff will not be produced until the soil moisture content of aeration zone reaches the field capacity, and thereafter the total runoff equals to the total rainfall.To describe the nonuniform distribution of tension water capacity in sub-basins, a tension water capacity curve is introduced (Figure 3): where f represents the area of runoff production; F represents the area of basin; W' represents the tension water capacity at a point; and WMM is the maximum value of W'.The areal mean tension water capacity, WM, is easy to obtain from Equation ( 2) by integration.B is the parameter shown in Table 3.
. If Otherwise where P is the measured areal mean rainfall; EM is the measured pan evaporation; A is the ordinate value corresponding to WM; R is the runoff depth produced from permeable area;  is the initial tension water storage of sub-basin; and KC and WM are parameters shown in Table 3.
The Muskingum method is based on the principle of water balance and the relation of storage and outflow.It has two basic assumptions.One is that the water level is a straight line and flow varies linearly along the reach.The other is that there is a single relationship between the discharge Q at a section and the storage W in the river.
The water balance equation is as follows: where I is the inflow, Q is the outflow, W is the storage amount, and t is the computation time.Otherwise where P is the measured areal mean rainfall; EM is the measured pan evaporation; A is the ordinate value corresponding to WM; R is the runoff depth produced from permeable area; W 0 is the initial tension water storage of sub-basin; and KC and WM are parameters shown in Table 3.
The Muskingum method is based on the principle of water balance and the relation of storage and outflow.It has two basic assumptions.One is that the water level is a straight line and flow varies linearly along the reach.The other is that there is a single relationship between the discharge Q at a section and the storage W in the river.
The water balance equation is as follows: where I is the inflow, Q is the outflow, W is the storage amount, and t is the computation time.
The channel storage equation as follows: The flow proportion coefficient XE and the slope of the channel storage curve KE are used in the equation set.The flow routing equation can be obtained by solving the equation set, as below: Water 2019, 11, 866 where C 0 , C 1 , and C 2 are coefficients with the relationship: The structure of the Xinanjiang model is presented in Figure 4.The model parameters and the intermediate variables are displayed in Tables 3 and 4.

KC Calculation for evapotranspiration
The ratio of potential evapotranspiration to pan evaporation    [43].The 15 parameters can be divided into four categories: evapotranspiration parameters (KC, UM, LM, and C), runoff producing parameters (WM and B), water resources separating parameters (SM, EX, KG, and KI) and runoff concentration parameters (CS, CI, CG, KE, and XE).Each category is basically independent from the others.Nine of these parameters (KC, SM, KI, KG, CS, CI, CG, KE and XE) are sensitive.The remaining ones are insensitive.The values of insensitive parameters are given according to the basin information.The data series are separated into calibration and validation periods.The calibration process is for the sensitive ones and is conducted according to the order of parameter category shown in Table 3.
When calibrating the parameters, the model is divided into daily model and hourly model.The input data of daily model are in daily intervals.It mainly aims at water balance and calibrates the parameters KC, UM, LM, C, WM, and B. The input data of hourly model are in hourly intervals.It calibrates the remaining parameters.
The goodness-of-fit measures used in the hourly model are relative error and certainty coefficient.The relative error contained the runoff depth and flood peek discharge.
where Robs is observed runoff depth (mm), Rsim is calculated runoff depth (mm), Qm obs is observed peak discharge (m 3 /s), Qm sim is calculated peak discharge (m 3 /s), and ER and EQ are, respectively, related to runoff depth and peak discharge.The certainty coefficient as follows: where DC is the certainty coefficient, y c is the simulative value, y o is the observed value, y o is the average of the observed value, and n is the length of data series.

Daily Model
There is no evaporation station in Dapoling basin.We used the data from evaporation station named Xinyang of Xixian basin.The Dapoling basin is a sub-basin of Xixian basin.The distance of Dapoling to Xinyang is about 60 km.We used the observed daily precipitation data from 13 rainfall stations daily evapotranspiration data of Xinyang station between 1965 and 2012 for the input of daily model.The daily discharge data of Dapoling station were used for parameter calibration.

Hourly Model
We selected 37 flood events (Table 5) between 1965 and 2012 in the Dapoling basin.The small flood events between 1965 and 2012 were not considered for calibrating the flood forecast system used by the basin management agency.Therefore, we chose the flood events that were larger than 900 m 3 /s.To calculate the centroid of precipitation and runoff, we chose 28 single peak floods for lag time calculation.The observed precipitation data for model input were from the 13 rainfall stations listed in Table 1.The data were converted to average rainfall data by the Thiessen polygon method to represent the basin rainfall to be used conveniently in the calculation of lag time.The observed discharge data were from the Dapoling station.Initial precipitation and discharge data were recorded every half hour, every hour, every two hours, every four hours, or another time interval.The time interval depended on the measure level of the stations during 1965 and 2012.The rainfall and discharge data were calculated into the hourly interval data.The observed evaporation data were from Xinyang station outside the basin.Initial evaporation data were recorded every day.The data were calculated into the hourly interval by arithmetic mean method.The whole period was divided into two periods of 1965−1985 and 1986−2012 according to Qu et al. [38] and Zhang et al. [39].The comparison of runoff depth, flood peak, kurtosis coefficient, rainfall-runoff relationship and lag time was based on the two series.The lag time was calculated by the 28 single peak flood events.The others were calculated by all the flood events that participated in parameter calibration.

The Calculation of Lag Time
Four types of lag time are introduced in Figure 6 [23]: TLPC is the centroid lag-to-peak time from the centroid of precipitation to the peak discharge; TLC is the centroid lag time from the centroid of precipitation to the centroid of discharge; TLP is the lag-to-peak time from the beginning of precipitation to the peak discharge; and TLPP is the peak lag-to-peak time from the peak rainfall The comparison of runoff depth, flood peak, kurtosis coefficient, rainfall-runoff relationship and lag time was based on the two series.The lag time was calculated by the 28 single peak flood events.The others were calculated by all the flood events that participated in parameter calibration.

The Calculation of Lag Time
Four types of lag time are introduced in Figure 6 [23]: T LPC is the centroid lag-to-peak time from the centroid of precipitation to the peak discharge; T LC is the centroid lag time from the centroid of precipitation to the centroid of discharge; T LP is the lag-to-peak time from the beginning of precipitation to the peak discharge; and T LPP is the peak lag-to-peak time from the peak rainfall intensity to the peak discharge.The centroid lag time T LC equals to the mean flow concentration time of the basin that is usually meant.
The runoff depth, flood peak and kurtosis coefficient were calculated for all the selected floods.The runoff depth was determined by: R = 3.6QΔt/F (15) where R is runoff depth (mm), Q is mean discharge (m 3 /s), Δt is time interval (h), and F is area for the basin (km 2 ).The average value was determined by: The mean square error was determined by: where s is mean square error, xi is variable value, and ̅ is mean variable value.The kurtosis coefficient k was determined by:  =  ( − ̅ ) / ( − 1) The variables in Equation ( 18) have the same meaning as before.

The Calculation of Rainfall-Runoff Relationship
The rainfall-runoff correlogram reflects the quantitative relationship between one rainfall event and the total amount of corresponding flood.This research adopted the Antecedent Precipitation Index (Pa) to study the rainfall-runoff relationship in Dapoling basin.The linear regression equation was used to simulate the relationship between rainfall and runoff: where y is runoff depth (R/mm) and x is the sum of precipitation (P/mm) and its antecedent precipitation index (Pa/mm).The values of x and y were calculated for every event.
The slope of the line "a" reflects the runoff efficiency, and the x-intercept "b" reflects the rainfall threshold that produces runoff.The four types of lag time represent the time index of the runoff produced by the precipitation.The centroid lag time T LC equals the basin mean concentration time, therefore the information of concentration time can be obtained from the calculation of T LC .The other three lag times are all relative with the time of peak discharge.
The centroid of precipitation was determined by: where t wc is centroid of precipitation, w i is precipitation for period i, t i is time for period i, and n is the data size.The centroid of runoff was determined by: where t qc is centroid of runoff, Q i is mean runoff for period i.
The runoff depth, flood peak and kurtosis coefficient were calculated for all the selected floods.The runoff depth was determined by: R = 3.6Q∆t/F (15) where R is runoff depth (mm), Q is mean discharge (m 3 /s), ∆t is time interval (h), and F is area for the basin (km 2 ).The average value was determined by: average The mean square error was determined by: where s is mean square error, x i is variable value, and x is mean variable value.The kurtosis coefficient k was determined by: The variables in Equation ( 18) have the same meaning as before.

The Calculation of Rainfall-Runoff Relationship
The rainfall-runoff correlogram reflects the quantitative relationship between one rainfall event and the total amount of corresponding flood.This research adopted the Antecedent Precipitation Index (Pa) to study the rainfall-runoff relationship in Dapoling basin.The linear regression equation was used to simulate the relationship between rainfall and runoff: where y is runoff depth (R/mm) and x is the sum of precipitation (P/mm) and its antecedent precipitation index (Pa/mm).The values of x and y were calculated for every event.
The slope of the line "a" reflects the runoff efficiency, and the x-intercept "b" reflects the rainfall threshold that produces runoff.

Parameters Change of the Basin
The parameters before/after land use change are shown in Tables 6 and 7.The results in Table 7 show that the relative errors of runoff depth and peak discharge were all less than 0.25, and the certainty coefficients were all larger than 0.78.Twenty-eight of the 37 flood events qualified the inspection standard.The simulative result is good.Robs is observed runoff depth; Rsim is calculated runoff depth; Qm obs is observed peak discharge; Qm sim is calculated peak discharge, ER and EQ are related to runoff depth and peak discharge, respectively; DC is certainty coefficient.
The results in Table 6 show that the parameters of Xinanjiang Model varied largely after land use change, especially the sensitive ones.The ratio of potential evapotranspiration to pan evaporation KC enlarged from 0.77 to 0.95 after land use change, which means that the existing land use type promotes more evapotranspiration than runoff losses.The areal mean tension water capacity WM decreased from 150 to 145 mm.The three recession constants (surface water storage CS: 0.55-0.56;interflow storage CI: 0.85-0.9;and groundwater storage CG: 0.993-0.994)increased after land use change, which means that the runoff would recess slowly.The outflow coefficient of free water storage to the interflow KI (0.42-0.5) increased, and the outflow coefficient of free water storage to the groundwater flow KG (0.28-0.2) decreased after land use change.The land use data (Table 2 and Figure 2) show that, after 1985, the area of water and paddy field increased.The increase might have caused the above variation of runoff recession.

Runoff Depth
The calculated results of runoff depth before and after land use change in the basin are shown in Tables 8 and 9.The results in Table 8 indicate that the mean runoff depth after land use change (Series B) was slightly larger than that before land use change (Series A).The land use data in Table 2 express that the general variation trend of woodland and farmland decreased, while the urban land increased.This might have caused the slight increase of runoff depth.For the aspect of parameters of Xinanjiang Model in Table 3, the evaporation coefficient KC increased after land use change, indicating more consumption of water, while the areal mean tension water capacity WM reduced after land use change, indicating more production of water, which means that perhaps the WM influenced more.In addition, the mean square error after land use change was larger than that before land use change (before: 49.90; after: 50.04).The land use data in Table 2 show that, after 1985, the water and the paddy field increased and the farmland decreased, which may reduce the infiltration and produce more water.Other research [16] shows that the decreasing of farmland may cause runoff reducing.In Dapoling basin, the farmland decreased slightly (from 41.85% to 41.81%), which might have some impact on the variation of runoff.This means that, if the land use of an area tends to vary similarly with this condition, the runoff of a flood may increase.This has meaningful guidance to flood management.Among the range of rainfall exceeding 200 mm, five events were between 200 mm and 300 mm, one event was between 300 mm and 400 mm, and one event was greater than 400 mm.The same as below.
There were seven events with rainfall of less than 100 mm, 23 events with rainfall of between 100 mm and 200 mm, and seven events with rainfall of greater than 200 mm.Among the events with rainfall greater than 200 mm, five were between 200 mm and 300 mm, one was between 300 mm and 400 mm, and one was greater than 400 mm.Data in Table 9 indicate that the variation of runoff depth was influenced slightly by different size of rainfall events.The result for rainfall greater than 200 mm could not be obtained due to the lack of flood sample.This means that the result can only be implicated for the rainfall that smaller than 200 mm.

Flood Peak
The result of flood peak before and after land use change in Dapoling basin is shown in Tables 8  and 9. Table 8 indicates that the mean flood peak has increased due to land use change, which increased from 1821 m 3 /s before land use change (Series A) to 1916 m 3 /s after land use change (Series B).The maximum value of flood peak reduced after 1985 (from 4220 to 3740 m 3 /s).Land use data in Table 2 show that, after 1985, the water and paddy field increased, which means that the regulation and storage of water might have reduced the flood peak.
Data in Table 9 show the relationship between the variation of flood peak before/after land use change and the precipitation level.The mean flood peak after land use change was larger than that before land use change when rainfall was less than 200 mm.The variation of mean flood peak when rainfall was greater than 200 mm additionally could not be obtained obviously due to the lack of flood sample.However, Hood et al. [23] indicated that the variation of flood peak is related with the size of the rainfall event in low impact and traditional residential development areas [23].

Kurtosis Coefficient
The kurtosis coefficient reflects the shape of discharge hydrograph to a certain degree.The larger the kurtosis is, the more leptokurtic the hydrograph is, which has more values around the mean than a normal distribution.Leptokurtic hydrograph means that the flood rises and recesses rapidly.The result of kurtosis coefficient before and after land use change in Dapoling basin is shown in Tables 8 and 9.It can be concluded that the kurtosis coefficient after land use change (Series B) was larger than that before land use change (Series A).A larger kurtosis coefficient means more leptokurtic distribution of discharge hydrograph, which indicates that the flood after land use change rises and recesses more rapidly than that before land use change.This may be related with the increasing of urban area (from 0.63 to 0.88).8  and 9.It can be concluded that the kurtosis coefficient after land use change (Series B) was larger than that before land use change (Series A).A larger kurtosis coefficient means more leptokurtic distribution of discharge hydrograph, which indicates that the flood after land use change rises and recesses more rapidly than that before land use change.This may be related with the increasing of urban area (from 0.63 to 0.88).

Rainfall-Runoff Relationship
The rainfall-runoff correlograms before and after land use change are shown in Figure 7.The equation before land use change is as follows: The equation after land use change is as follows: y = 0.8033x − 58.623 (21) where y represents runoff depth (R/mm) of per event and x represents the sum of precipitation (P/mm) and its antecedent precipitation index (Pa/mm).
According to the result in Figure 7, the slope of the line after land use change was slightly larger than that before land use change, which means that the rainfall after land use change may produce more runoff than that before land use change.This is consistent with the change of runoff depth discussed above.The x-intercept after land use change was larger than that before land use change (Figure 7), which means that the rainfall after land use change produces runoff more difficult.

Mean Lag Time Change of the Basin
The results indicate that the mean basin lag time varied before/after 1985.The result of centroid lag-to-peak time TLPC is shown in Table 10.Compared with Series A, the mean centroid lag-to-peak The equation after land use change is as follows: y = 0.8033x − 58.623 (21) where y represents runoff depth (R/mm) of per event and x represents the sum of precipitation (P/mm) and its antecedent precipitation index (Pa/mm).
According to the result in Figure 7, the slope of the line after land use change was slightly larger than that before land use change, which means that the rainfall after land use change may produce more runoff than that before land use change.This is consistent with the change of runoff depth discussed above.The x-intercept after land use change was larger than that before land use change (Figure 7), which means that the rainfall after land use change produces runoff more difficult.

Mean Lag Time Change of the Basin
The results indicate that the mean basin lag time varied before/after 1985.The result of centroid lag-to-peak time T LPC is shown in Table 10.Compared with Series A, the mean centroid lag-to-peak time T LPC of Series B decreased from 11.92 h to 11.67 h.The mean square error has slightly increased from 9.97 to 10.21.The decreased T LPC means that the flood peak after land use change may come faster than that before land use change.The centroid lag-to-peak time T LPC in all selected floods is positive.Other studies found that the mean centroid lag-to-peak time T LPC may be negative [23].The negative values of the lag time are caused by the strong rainfall intensity occurring at the beginning of the precipitation.The intensity reaches its maximum soon after the precipitation.The long duration but low intensity precipitation process after the maximum leads to the centroid of precipitation towards right deviation, while the peak of the runoff appears soon after the peak rainfall intensity.Similarly, the calculated centroid lag time T LC also varied before/after 1985.The results are shown in Table 10.In Table 10, the average T LC slightly decreased from 25.28 h to 25.10 h.In addition, the variation range of centroid lag time T LC after land use change was smaller.Before land use change, the T LC varied from 6.6 h to 67.4 h, and its range was 60.8 h.The mean square error of T LC of Series A was 11.01.After land use change, the T LC varied from 6.9 h to 62.4 h, and its range was 55.5 h, with a mean square error of it was 11.73.T LC is the centroid lag time from the centroid of precipitation to the centroid of discharge, which means that T LC represents the mean concentration time of the basin.These variations of T LC indicate that the land use change resulted in the mean concentration time variation of the basin.The decreased average T LC means that the existing land use reduces the mean basin concentration time.The maximum value of T LC is smaller after land use change, while the minimum value is larger than that before 1985.This means that, after land use change, the distribution of the flood concentration time is more intensive.
The calculated results of lag-to-peak time T LP are shown in Table 10.The results indicate that the lag-to-peak time T LP after land use change was different from that before land use change.The average value under the condition after land use change was 1.33 times that before land use change.This means that the flood peak comes slower after land use change.In addition, the variation range of lag-to-peak time after land use change was larger than that before land use change.Before land use change, the T LP varied from 10 h to 61 h, the range was 51 h and the mean square error was 11.51.After land use change, the T LP varied from 10 h to 63 h, the range of which was 53 h and mean square error was 11.68.Lag-to-peak time T LP was positively correlated with precipitation duration (r: Series A: 0.834; Series B: 0.847) and runoff depth (r: Series A: 0.627; Series B: 0.605) before/after land use change.No relation was found among centroid lag-to-peak time T LPC , centroid lag time T LC or peak lag-to-peak time T LPP and precipitation duration or runoff depth (r < 0.35).
The calculated results of mean peak lag-to-peak time T LPP are shown in Table 10.The results show that the peak lag-to-peak time T LPP after land use change was slightly different from that before land use change.The average value after land use change was smaller than that before land use change (Series A: 11.11 h; Series B: 10.86 h), which means that the time of peak discharge will appear more quickly after land use change.In addition, the range of peak lag-to-peak time after land use change was similar to that before land use change.For the two series, the maximum value of T LPP was 28 h, and the minimum value was 4 h, but the mean square errors were different (Series A: 4.69; Series B: 4.86).The variation of mean square error indicates that the distribution range of the time from peak precipitation to peak discharge was slightly wider.
The mean basin lag times in different rainfall are shown in Table 11.The results in Table 11 show that the change of mean basin lag time varies with different rainfall depth.When the rainfall depth was less than 100 mm, the four kinds of lag times increased after land use change.When the rainfall was between 100 mm and 200 mm, the four kinds of lag times of Series B decreased compared with Series A. When the rainfall was greater than 200 mm, the T LC decreased, while the other three lag times seldom changed.Since the number of the events when rainfall was greater than 200 mm was few, the results of that condition need to be further studied.Synthesizing the results in Tables 10 and 11, different lag times were influenced by different size of rainfall events.The results in in Table 10 show that, for the all rainfall events, the values of centroid lag-to-peak time T LPC , centroid lag time T LC and peak lag-to-peak time T LPP decreased after land use change.The results in Table 11 show that, when the rainfall was between 100 and 200 mm, the above three lag times decreased after land use change, while, for other rainfall events, the three lag times increased or changed slightly.Therefore, the centroid lag-to-peak time T LPC , centroid lag time T LC and peak lag-to-peak time T LPP were influenced more by the rainfall events that were between 100 mm and 200 mm than other events.For the lag-to-peak time T LP , the mean value of all the events increased after land use change (Table 10).However, the results in Table 11 reveal that, when the rainfall was less than 100 mm, the value of T LP increased after land use change, which expressed that the T LP was influenced more by the rainfall events smaller than 100 mm.
This study compared the variation of four types of lag time before and after land use change in Dapoling basin by calculating the selected 28 flood events based on Xinanjiang Model.In the management of basin flood forecasting, different precipitation processes arrive at the time of peak discharge.Once the mean concentration time was determined, T LP , T LPP and T LPC could be reached from the time that the precipitation begins, the peak of precipitation appears and the precipitation ends in the whole process of precipitation.During the process of floods, the beginning time of precipitation is the most easily obtained.It is obvious that, if we want to obtain the time of peak discharge as early as possible, T LP is the most useful information, because we only need to know the beginning time of precipitation.The uncertainty of the future precipitation makes it hard to judge the peak of rainfall, thus the information of T LPP is difficult use.If the rainstorm center is in the upstream far away from the downstream, and the peak discharge does not appear in the downstream after the whole rainfall process, the time of peak discharge in the downstream can be obtained from T LPC .Every type of lag time has its characteristics, which means that it is necessary to choose appropriate types according to actual situation in basin flood management.The determination of average log time of a basin makes flood forecasting easier.When a flood event comes, the time of flood peak could be predicted.
The variation of the mean basin lag times for two periods was due to the land use change.Data in Table 2 indicate that land use types changed from the 1980s to the 1990s and form the 1990s to the 2000s.In addition, the fact that after 1985 the land use changed revealed by the variation of lag time distribution and the presence of outliers may relate with the land use data, which is a future aspect of this study.

Conclusions
Considering land use change, this study used the Xinanjiang Model to calibrate parameters before/after land use change and to compare the hydrology response to the change.The research showed that the parameters varied after the change.The change of parameters means that the existing land use promotes more evapotranspiration and consumes more water.The runoff discharge will recess more slowly.This may increase the difficulty of travel and rescue in the future.It is meaningful for the government to make contingency plans for the floods in this region.
Additionally, the variation results of the runoff depth, flood peak and kurtosis coefficient relationship showed that they have increased slightly after land use change for floods produced by rainfall smaller than 200 mm.The analysis of rainfall-runoff relationship revealed that the rainfall after land use change produced more runoff with less production threshold.Refering to the land use data, we found that the change of water, paddy field and farmland may be the main land uses that influence the change of these hydrology elements.
Due to the lack of huge floods, the variation of runoff depth, flood peak and kurtosis coefficient for floods produced by rainfall greater than 200 mm could not be well obtained.More flood events data should be collected to obtain a more useful conclusion.Dapoling basin is the head source of Huai River basin.There are some basins with similar environment in Huai River basin.The parameters of Dapoling can be the reference of other basins when calibrating hydrology model and the results obtained in Dapoling basin can be extended to other sub-basins in the Huai River basin, which can provide important references for soil and water conservation, water resources use and management of land use in the Huai River.
13 01"N to 32• 43 05"N, as shown in Figure1.The length of the main channel in Dapoling basin is about 73 km and the area is 1640 km2 , covered with various vegetation growing well.The Dapoling basin is located in the transition zone between subtropical zone and temperate zone, and is heavily affected by the monsoon climate in the flood season.The perennial average evaporation is about 1000 mm in the region and the annual average temperature is between 11 • C and 16 • C. The main soil types of the basin are light silty loam and sandy loam soil, and part is underlain by silt clay.Rice and wheat are the main crops grown in the area.Dapoling basin is the head source of the Huai River.Most of the basin is covered by Dabie Mountain.The rivers flow over a mountainous terrain with many tributaries and large slope.There are no reservoirs in the river.Water 2018, 10, x FOR PEER REVIEW 3 of 20 The sub-basin Dapoling is located in the upstream of the Huai River at the latitudes 113°16′04″E to 113°49′26″E, and longitudes 32°13′01″N to 32°43′05″N, as shown in Figure 1.The length of the main channel in Dapoling basin is about 73 km and the area is 1640 km 2 , covered with various vegetation growing well.The Dapoling basin is located in the transition zone between subtropical zone and temperate zone, and is heavily affected by the monsoon climate in the flood season.The perennial average evaporation is about 1000 mm in the region and the annual average temperature is between 11 °C and 16 °C.The main soil types of the basin are light silty loam and sandy loam soil, and part is underlain by silt clay.Rice and wheat are the main crops grown in the area.Dapoling basin is the head source of the Huai River.Most of the basin is covered by Dabie Mountain.The rivers flow over a mountainous terrain with many tributaries and large slope.There are no reservoirs in the river.

Figure 1 .
Figure 1.The study area of Dapoling basin.

Figure 1 .
Figure 1.The study area of Dapoling basin.

Figure 2 .
Figure 2. The land use type of the three periods.

Figure 2 .
Figure 2. The land use type of the three periods.

Figure 3 .
Figure 3.The tension water capacity curve in Xinanjiang Model.

Figure 3 .
Figure 3.The tension water capacity curve in Xinanjiang Model.

3. 3 .
Data Processing3.3.1.The Composing of Runoff SeriesThere are three test methods in the previous studies: the variation of the model errors of the whole period, the comparison of parameters in different periods, and the model simulation change compared with the reference period.One of the steps to test the change of data series is to run the model for some periods and compare the simulated runoff.In this study, the daily precipitation of rainfall stations and the daily evaporation were the input data for the model.The daily discharge was simulated by parameters of different periods.The authors found that the land use changed after 1985.Thus, the study period was divided into two parts of 1965−1985 and 1986−2012, according to the land use change tested in the previous study.We set two scenarios: Scenario A represented the condition before land use change, and Scenario B represented the condition after change.The hydrologic data were divided into two parts according to the two scenarios.Xinanjiang Model was used to calibrate parameters for the two parts separately and arrived at two sets of parameters.There were 18 flood events before 1985, of which 15 were for calibration and 3 for validation.There were 19 flood events after 1985, of which 16 were for calibration and 3 for validation.The set of parameters before land use change was named Parameter A, and the other was named Parameter B. Therefore, Parameter A represented the condition before land use change, and Parameter B represented after land use change.When simulating runoff, we used Parameter A to simulate the runoff of 1986 and 2012.The inputs were the precipitation and evapotranspiration data of 1986 and 2012.The output discharge and the observed discharge of 1965 and 1985 composed the whole period discharge and named Series A. Parameter B was used to simulate the runoff of 1965 and 1985.The inputs were the precipitation and evapotranspiration data of 1965 and 1985.The output discharge and the observed discharge of 1986 and 2012 composed the whole period discharge named Series B. The data processing is shown in Figure 5. Therefore, the two-discharge series were both based on the same precipitation and evapotranspiration data, thus Series A represented Scenario A, and Series B represented Scenario B. The difference of the two-runoff series was only caused by the land use change.Water 2018, 10, x FOR PEER REVIEW 10 of 20 model for some periods and compare the simulated runoff.In this study, the daily precipitation of rainfall stations and the daily evaporation were the input data for the model.The daily discharge was simulated by parameters of different periods.The authors found that the land use changed after 1985.Thus, the study period was divided into two parts of 1965−1985 and 1986−2012, according to the land use change tested in the previous study.We set two scenarios: Scenario A represented the condition before land use change (1965−1985), and Scenario B represented the condition after change (1986−2012).The hydrologic data were divided into two parts according to the two scenarios.Xinanjiang Model was used to calibrate parameters for the two parts separately and arrived at two sets of parameters.There were 18 flood events before 1985, of which 15 were for calibration and 3 for validation.There were 19 flood events after 1985, of which 16 were for calibration and 3 for validation.The set of parameters before land use change was named Parameter A, and the other was named Parameter B. Therefore, Parameter A represented the condition before land use change, and Parameter B represented after land use change.When simulating runoff, we used Parameter A to simulate the runoff of 1986 and 2012.The inputs were the precipitation and evapotranspiration data of 1986 and 2012.The output discharge and the observed discharge of 1965 and 1985 composed the whole period discharge and named Series A. Parameter B was used to simulate the runoff of 1965 and 1985.The inputs were the precipitation and evapotranspiration data of 1965 and 1985.The output discharge and the observed discharge of 1986 and 2012 composed the whole period discharge named Series B. The data processing is shown in Figure 5. Therefore, the two-discharge series were both based on the same precipitation and evapotranspiration data, thus Series A represented Scenario A, and Series B represented Scenario B. The difference of the two-runoff series was only caused by the land use change.

Figure 6 .
Figure 6.Definitions of terms used to describe hyetographs and response hydrographs.tw0, beginning of precipitation; twp, peak of precipitation; twc, centroid of precipitation; tpk, time of peak discharge; tqc, centroid of runoff; TLPC, centroid lag-to-peak; TLC, centroid lag; TLP, lag-to-peak; TLPP, peak lag-to-peak.The other variables are the same as above.

Figure 6 .
Figure6.Definitions of terms used to describe hyetographs and response hydrographs.t w0 , beginning of precipitation; t wp , peak of precipitation; t wc , centroid of precipitation; t pk , time of peak discharge; t qc , centroid of runoff; T LPC , centroid lag-to-peak; T LC , centroid lag; T LP , lag-to-peak; T LPP , peak lag-to-peak.The other variables are the same as above.

20 )
The rainfall-runoff correlograms before and after land use change are shown in Figure 7.The equation before land use change is as follows: y = 0.7544x − 51.314 (Water 2018, 10, x FOR PEER REVIEW 15 of 20 a normal distribution.Leptokurtic hydrograph means that the flood rises and recesses rapidly.The result of kurtosis coefficient before and after land use change in Dapoling basin is shown in Tables

Figure 7 .
Figure 7.The relationship between precipitation and runoff depth.

Figure 7 .
Figure 7.The relationship between precipitation and runoff depth.

Table 1 .
The data and stations of Dapoling basin.

Table 1 .
The data and stations of Dapoling basin.

Table 2 .
The land use data for three periods in the basin.The land use data are from Resource and Environment Data Cloud Platform of Chinese Academy of Sciences (Table2 and Figure

Table 2 .
The land use data for three periods in the basin.

Table 4 .
State variables in Xinanjiang model.
Zhao et al. introduced a detailed method of parameter calibration

Table 5 .
The flood events in Dapoling basin.
P refers to the observed rainfall, mm.R refers to the observed runoff depth, mm.Q m refers to the observed discharge, m 3 /s.The flood code with "*" was used for lag time calculation.

Table 6 .
Parameters of the two periods before/after land use change.

Table 7 .
The simulation results of Xinanjiang Model in Dapoling Basin.

Table 8 .
The result of hydrologic variables in Dapoling basin.
R refers to runoff depth (mm), Q m refers to flood peak (m 3 /s), s refers to mean square error, and the other variables are the same as above and below.

Table 9 .
The result of hydrologic variables of different rainfall in Dapoling basin.

Table 10 .
The average values of mean basin lag times in the basin.
T LPC means centroid lag-to-peak, T LC means centroid lag, T LP means lag-to-peak, T LPP means peak lag-to-peak, s means mean square error, the units in the table is h, the same as below.

Table 11 .
The lag times for different rainfall.