An Investigation into Sub-Basin Rainfall Losses in Different Underlying Surface Conditions Using HEC-HMS : A Case Study of a Loess Hilly Region in Gedong Basin in the Western Shanxi Province of China

Basins located in loess hilly–gully regions often suffer flood disasters during the flood season. Meanwhile, the underlying surface of the region can increase the rainfall losses, thereby reducing the flood volume. Therefore, the prediction of rainfall losses on the underlying surface is necessary for scientifically and reasonably forecasting the flood volume. The relationship between the rainfall losses and underlying characteristics was investigated and a method for predicting the rainfall losses using HEC-HMS was presented in this paper with a case study in the Gedong basin, a typical loess hilly region of western Shanxi Province in northern China. Results showed that HEC-HMS could be applied to loess hilly–gully regions. The loss computation results suggested that the losses of sub-basins varied with the density of rainfall. The analysis of influences of rainfall losses, including forestland percentage and slope, indicated that the former had a positive impact, while the latter had a negative influence. The impact of forestland percentage is larger than that of slope. Furthermore, with the increase of forestland percentage, its correlation with rainfall losses was enhanced, and the correlation coefficient ranged between 0.64 and 0.84 from the 1970s to the 2010s.


Introduction
Rainfall losses, which occur during flood events, play a crucial role in real-time flood forecasting and flood estimation [1,2].The losses are generated through canopy interception, infiltration, depression storage, evaporation, and evapotranspiration [3].Rainfall losses by canopy interception are a significant part of hydrological losses from forested ecosystems.Rainfall interception losses mainly rely on the rainfall characteristics, forest structure, and climatic changes governing the rates of evaporation and evapotranspiration during and after rainfall events [4,5].Infiltration is a major process in flood generation.Infiltration rate is a function of antecedent soil moisture content that decides the magnitude of flood peak and volume [6].In semi-arid agricultural areas, evaporation and evapotranspiration from the soil have a significant impact on rainfall losses [7].However, the underlying surface is a paramount factor affecting rainfall losses during flood events.The effect of the underlying surface on rainfall losses is different because of its vegetation, slope, and tillage measures.It is considered that vegetation might prevent rainfall from reaching the surface due to canopy interception and infiltration; the surface slope is large, the flow speed is fast, and its retention time is short, leading to a small infiltration amount.Tillage measures (terrace) change the slope micro-terrain and increase the roughness so that soil moisture infiltration performance changes, thus affecting rainfall losses.In short, the underlying surface can affect the rainfall losses by interception and infiltration in the flood process [8].
In recent years, many studies have been conducted to examine the relationship between rainfall losses and land use type with different models.The Gash model [9] was applied to the basin covered by sugarcane and riparian forest in order to assess the influence of rapid sugarcane canopy changes on rainfall interception losses and explain the main factors determining the rainfall interception losses [10].The Liu model [11] was used in an experimental multispecies (Acacia mangium, Gliricidia sepium, Guazuma ulmifolia, Ochroma pyramidale, and Pachira quinata) tree plantation in Soberania National Park to predict rainfall interception losses and estimate the water storage capacities of tree boles; the results revealed that the interspecific differences between observed and simulated cumulative interception loss were significant, with Acacia mangium intercepting more rainfall than other species.The Liu model was most sensitive to variations of evaporation rate [12].In a mixed evergreen and deciduous broadleaved forest, the rainfall interception loss predictions of the revised Gash model were evaluated and the magnitude of gross precipitation and its distribution into interception losses, through fall and stemflow were quantified [13].In a word, many investigations concerning rainfall losses focused on vegetation interception, while infiltration research was relatively scarce.Gedong basin is a typical loess hilly-gully region; a number of studies in loess hilly-gully regions focus on soil erosion and soil loss.Feng et al. [14] compared the soil and water conservation performances under three single land use types (cultivated land, CL; switchgrass, SG; and abandoned land, AL) and two composite land use types (CL-SG and CL-AL).The results indicated a general trend in the number of runoff and soil loss events for the five land-use types: CL = CL-SG > CL-AL > SG > AL, and the vegetation coverage was the primary factor controlling soil erosion.Zhang et al. [15] analyzed the influences of the thickness of an aeolian sand layer overlying a loess slope on runoff and sediment production processes by eight simulated rainstorms in the Wind-Water Erosion Crisscross Region of the northern Loess Plateau.Zhu et al. [16] combined three sets of plot data (short slope plots, long slope plots, and soil conservation plots) to evaluate the effectiveness of different conservation measures in reducing runoff and soil loss.Yan et al. [17] explored the effect of watershed management practices on the relationship between runoff and sediment by analyses in the hilly-gully regions of the Loess Plateau; the results suggested that a combination of hillslope and gully erosion control practices effectively reduced sediment delivery and erosion.Li et al. [18] determined the effect of different land use types (artificial forestland, native grassland, and artificial grassland) on soil organic carbon and total soil nitrogen by experimental and statistical analysis; the results showed that land use types had great influence on soil organic carbon and total soil nitrogen, and artificial grassland was the optimal choice to mitigate soil carbon and nitrogen loss in the loess hilly-gully region.
In summary, a majority of previous studies in the loess hilly-gully region concentrated on soil loss and erosion; only a few studies have analyzed rainfall losses with a focus on vegetation interception.In the Gedong basin, most rainfall losses are generated through interception and infiltration.Hence the sub-basins' rainfall losses (interception and infiltration) under different underlying surface conditions are determined in this study by the HEC-HMS model (Hydrologic Engineering Center-Hydrological Model System, U.S. Army Corps of Engineers-Hydrologic Engineering Center, Washington, DC, USA).
The HEC-HMS model has different methods to calculate the rainfall losses such as Soil and Conservation Service, Green and Ampt, Initial-Constant, Deficit-Constant, Exponential and Soil Moisture Accounting [6].The Soil and Conservation Service (SCS) Curve Number (CN) loss model estimates precipitation excess as a function of cumulative precipitation, land use, soil cover, and antecedent soil moisture [6].The CN can reflect the underlying condition of the basin.Considering CN can reflect the underlying characteristics of the basin, the SCS-CN loss method is applied to this study.
The objective of this study is to analyze the relationship between the rainfall losses and underlying characteristics in the hilly-gully loess region.The underlying characteristics include vegetation, soil, slope, and so on.In this study, the SCS-CN loss model is used to calculate the rainfall losses of sub-basins in different underlying surface conditions, and then the relationship between the rainfall losses and underlying characteristics (the percentage of forestland area and surface slope in sub-basins) is analyzed on the basis of rainfall losses calculation results and existing experimental results in the loess hilly-gully region.

Meteorology and Hydrology
Gedong basin, with an area of 724.6 km 2 , is located in the western part of Shanxi Province, China (Figure 1).The geographical coordinates are east longitude 111 • 03 -111 • 35 , latitude 37 • 38 -38 • 11 , about 30 km from east to west, 30 km from north to south.The basin is a typical hilly and gully region of the Loess Plateau.The area of soil erosion is 524 km 2 , and the area of soil and water conservation is 621 km 2 .The climate of the basin is arid and semi-arid warm temperature, the average annual temperature is 8.7 • C, and the average annual rainfall is 517 mm.Rainfall in the basin is rainstorm type, which is characterized by a small area of heavy rain, large intensity, short duration, and the largest rainstorms concentrated in July or August.According to Gedong hydrological station statistics, the average annual runoff is 109 million m 3 .

Soil and Land Use
The soil types are relatively simple; most of the study areas are covered with clay loam, sandy loam, and sandy clay.Land use data of the Gedong basin are extracted from remote sensing images, mainly of forest land, grassland, construction land, and cultivated land (the terraces are located on cultivated land); the situation of land use is shown in Figure 2. Through years of comprehensive

Soil and Land Use
The soil types are relatively simple; most of the study areas are covered with clay loam, sandy loam, and sandy clay.Land use data of the Gedong basin are extracted from remote sensing images, mainly of forest land, grassland, construction land, and cultivated land (the terraces are located on cultivated land); the situation of land use is shown in Figure 2. Through years of comprehensive management of soil and water conservation, the land use efficiency of Gedong basin is significantly changed.From the remote sensing image data, the area of the terrace is 98.48 km 2 , accounting for 13.6% of the basin area.Land use in different stages is shown in Table 1.After 2010, the area of cultivated land accounts for 24.7% of the total area; forestland is the main land use type, accounting for 60.89% of the total area.Compared with the 1980s, the forest land increased by 17.58%, and the cultivated land decreased by 11.58%.From the 1980s to 2010s, part of the cultivated land was turned into forest.

Surface Slope Characteristics
Slope has a great impact on rainfall losses, runoff, infiltration, soil erosion, and so on.The slope map of Gedong basin is shown in Figure 3.

Surface Slope Characteristics
Slope has a great impact on rainfall losses, runoff, infiltration, soil erosion, and so on.The slope map of Gedong basin is shown in Figure 3.

Surface Slope Characteristics
Slope has a great impact on rainfall losses, runoff, infiltration, soil erosion, and so on.The slope map of Gedong basin is shown in Figure 3.The minimum slope of the Gedong basin is 0 • and the maximum slope of the Gedong basin is 54.76 • .The average slope of the 11 sub-basins is 13.53 • , the maximum is 14.97 • , and the minimum is 12.44 • .Figure 4 shows the scatter plot of average slope with the change in basin area.It can be seen that with the increase in basin area, the average slope first increases, then gradually stabilizes; the relationship can be expressed as follows: Y average slope = 1.6815lnX + 6.6423R 2 = 0.6877. ( Water 2017, 9, 870 6 of 28 The minimum slope of the Gedong basin is 0° and the maximum slope of the Gedong basin is 54.76°.The average slope of the 11 sub-basins is 13.53°, the maximum is 14.97°, and the minimum is 12.44°.Figure 4 shows the scatter plot of average slope with the change in basin area.It can be seen that with the increase in basin area, the average slope first increases, then gradually stabilizes; the relationship can be expressed as follows: Yaverage slope = 1.6815lnX + 6.6423R 2 = 0.6877. ( According to the "TD/T1014-2007 second national land survey technical regulations" provisions, the slope can be categorized into five levels; the break points are 2, 6, 15, and 25° [30].According to this regulation, combined with the actual situation of the basin, the slope is divided into 0-6°, 6-15°, 15-25°, and 25-55°.Figure 5 reflects the percentage of the area under different slopes, varying with the basin area.It can be seen in the figure that a slope area of 0-6° is about 15% of the total basin area; 25-55° is about 6% of the total basin area; 6-15° and 15-25° area percentages vary from 30% to 50%.When the basin area is less than 75 km 2 , the 6-15° slope area percentage is higher than that for 15-25°; however, when the basin area is greater than 75 km 2 , the 6-15° slope area percentage is lower than that of 15-25°.It can be seen in the figure that a slope area of 0-6 • is about 15% of the total basin area; 25-55 • is about 6% of the total basin area; 6-15 • and 15-25 • area percentages vary from 30% to 50%.When the basin area is less than 75 km 2 , the 6-15 • slope area percentage is higher than that for 15-25 • ; however, when the basin area is greater than 75 km 2 , the 6-15 • slope area percentage is lower than that of 15-25 • .
Combined with changes in land use, the average slope of cultivated land is 8.93 • in the 2010s, and 15.2 • in 1980s, and most of the cultivated land with slope > 15 • is converted into forest.Returning cultivated land into forest has a great impact on cultivated land with slope > 15 • [31,32].

Data Collection
The HEC-HMS model was developed by the U.S. Army Corps of Engineers-Hydrologic Engineering Center to simulate rainfall and runoff.It is used in this study to compute the rainfall losses in the sub-basins.Model is composed of basin models, meteorological models, control specifications, and input data components.These components are used to simulate the hydrologic process in a watershed.Basin characteristics such as study area, streams parameters, digital elevation, slope, land use, soil types, and curve number are extracted from ArcGIS [33].
The input parameters used in the model are listed in Table 2.A simulation calculates the rainfall-runoff response in the basin based on input data from the meteorological model.Meteorological data mainly include precipitation and discharge [34].Precipitation data is obtained from the rainfall stations installed in the study area.In the model, these data are converted into effective precipitation data using the Thiessen polygon method by default (Figure 6).However, compared with spatial interpolation methods, the Thiessen polygon method has some limitations in calculating rainfall [34,35].The information of temporal and spatial precipitation distribution is significant when selecting data for the model calibration and validation.Discharge data are provided by Gedong hydrological station, which is the only hydrological station at the outlet of the basin.

3
Routing Muskingum K (travel time), X (weighting factor) A simulation calculates the rainfall-runoff response in the basin based on input data from the meteorological model.Meteorological data mainly include precipitation and discharge [34].Precipitation data is obtained from the rainfall stations installed in the study area.In the model, these data are converted into effective precipitation data using the Thiessen polygon method by default (Figure 6).However, compared with spatial interpolation methods, the Thiessen polygon method has some limitations in calculating rainfall [34,35].The information of temporal and spatial precipitation distribution is significant when selecting data for the model calibration and validation.Discharge data are provided by Gedong hydrological station, which is the only hydrological station at the outlet of the basin.A total of 18 flood events from 1970 to 2012 were simulated, and a parameter sensitivity analysis was carried out to provide the basis for the next analysis of rainfall loss in each sub-basin under In the control specifications, the time period and time step of the simulation run are defined (15 min) [33].Input data components contain time-series data, paired data, and grid data; these sections are used as boundary conditions or parameters in the model.

The Simulation of the HEC-HMS Model
The simulation consists of three steps.The first step is watershed delineation.In this step, the basic data such as precipitation, discharge, geological, soil data, land use, and topographic data are entered into the model.After importing the data, the model is calibrated to define the parameters and then these parameters are validated.Model calibration and validation are important to determine the factors that influence the characteristics of the study area.After the above steps are completed, the simulation results are obtained [36,37].

Rainfall Losses Calculation
In this study, the SCS Curve Number method is applied to calculate the losses in the catchment.This method describes the loss of precipitation by vegetation interception, depression storage, evaporation, evapotranspiration, and infiltration, and indicates whether runoff is generated.The SCS method calculates the precipitation excess as a function of cumulative rainfall, soil cover, land use, and antecedent moisture content [38,39].The equation is expressed as follows: where Q is the cumulative rainfall excess, P is the total rainfall, I a represents the initial abstraction, and S denotes the maximum potential storage.The formula is established when P > I a .
The initial abstraction can be computed with: The maximum potential storage has a relationship with the curve number, which can be computed as follows: CN (curve number) can be estimated as a function of land use, soil type, and antecedent moisture content in a basin; it comprehensively reflects the characteristics of the basin before precipitation.The CN layer is determined in ArcGIS by utilizing the soil type and land use map.In order to develop the CN layer, the soil map is divided into different hydrological groups such as A, B, C, and D [3,39].It is found that the CN varies from 57 to 88.The final CN for the sub-catchments is computed using the following formula: where Ai represents the drainage area of sub-catchment i and CN i is the CN of sub-catchment i.

Direct Runoff Calculation
In this study, the direct runoff in the catchment is calculated by utilizing the SCS Unit Hydrograph.This describes the process of converting the excess rainfall in the catchment into surface runoff.With the progress of precipitation, excess rainfall forms surface runoff along the surface soil into the river, then reaches the sub-basin exports, and ultimately forms runoff.The lag time for SCS unit hydrograph for each sub-catchment is determined as follows [3,40].
The relationship between the peak flow and the peak arrival time of the runoff curve is expressed as follows: where U p is the peak flow (m 3 /s), C is the conversion coefficient, A is the drainage area (km 2 ), and T p is the peak appearance time (h).
The relationship between the time of the arrival peak and the duration of the unit rainfall is where ∆t is the duration of the unit rainfall (the time interval of the simulation process in the HEC-HMS model) and t lag is lag time of the peak.

Flow Routing
The movement of a flood wave through a river reach is described by the Muskingum method.The key to using the Muskingum method is to determine K and X [20,21,41,42].K and X are calculated as follows: Water 2017, 9, 870 where K is the spread time of the flood wave through the river reach (s), ∆L is the length of routing (m), V w is the flood wave celerity (m/s), β is the exponent, V av is the average velocity (m/s), n is Manning's roughness, R is the hydraulic radius (m), and S is the bed slope of the channel (m/m) [43,44].
The value of β depends on river characteristics such as channel geometry, slope, and roughness.The rivers in Gedong basin are mostly wide, rectangular channels, so the value of 5/3 for β is obtained from Manning's equation for a wide rectangular channel where the hydraulic radius is equal to the average depth [43,44].
where, X is the weighting factor, Q 0 is the reference flow from the inflow hydrograph (m 3 /s), Q min is the minimum flow (m 3 /s), Q p is the peak flow (m 3 /s), and P is the bottom width of flow area or average width (m).

Model Construction and Simulation
The Gedong basin is divided into 11 sub-basins in order to better represent the spatial variation of parameters.The hydrological model requires that each sub-basin has at least one rainfall node that could represent the sub-basin.In this study, the centroid of the sub-basin was selected as the rainfall node of each sub-basin.The convergence lines of the two sub-basins formed the river channel, until the basin exit section.A generalized model of Gedong basin is shown in Figure 7.

Sensitivity Analysis
In this study, 18 flood events in different stages were applied to the simulation.A sensitivity analysis was used to examine the relative changes in the model outputs with respect to the change in the model input parameters [

Sensitivity Analysis
In this study, 18 flood events in different stages were applied to the simulation.A sensitivity analysis was used to examine the relative changes in the model outputs with respect to the change in the model input parameters [41].Theoretically, if x 1 , x 2 , . . ., x n were model input variables and y = y(x 1 , x 2 , . . ., x n ) was the model output, then the relative sensitivity (flexibility) of y with respect to the ith variable at (x 1 , . . ., x i , . . ., x n ) was equal to e(y) = ∂y ∂x i x 1 , . . ., x i , . . ., x n x i y x 1 , . . ., x i , . . ., If the absolute value of e was equal to or greater than 1, the model input parameter was flexible.Otherwise, the model input parameter was weakly flexible or inflexible [45,46].Sensitivity analysis was performed in two stages.The impact of CN (curve number), I a (initial abstraction), RC (attenuation coefficient), R (peak ratio), K (travel time), X (weighting factor of flow), and t lag (lag time) on the peak discharge (P) and the impact of CN, I a , RC, R, K, X, and t lag on flood volume (V) were assessed in the first and second stages, respectively [46].The results of the sensitivity analysis are shown in Table 3.The sensitivity analysis illustrated that CN, I a , t lag , and K were sensitive to the peak discharge and flood volume in a loess hilly region such as Gedong basin.
The parameters that were consistent with the flood events of the 1970s and 1980s observed in the hydrograph are shown in Table 4.The observed and simulated hydrographs of the 1970s and 1980s are shown in Figure 8. Overall, the figures described the shape and trend of the hydrographs of eight flood events as being similar, except 1970 and 1977.The total volume was slightly overestimated based on the observed hydrograph in 1977.The peak was slightly underestimated based on the observed hydrograph in 1982, 1986, and 1987, while a relatively perfect match was obtained in 1979.The parameters that were consistent with the flood events of the 1990s in the observed hydrograph are shown in Table 5.The observed and simulated hydrographs of 1990s are shown in Figure 9.
On the whole, the figures described the shape and trend of the hydrographs of four flood events as being similar.However, the peak of all events was slightly underestimated based on the observed hydrograph.
The parameters that were consistent with the flood events of the 2000s in the observed hydrograph are shown in Table 6.The observed and simulated hydrographs of 1990s are shown in Figure 9. On the whole, the figures described the shape and trend of the hydrographs of four flood events as being similar.However, the peak of all events was slightly underestimated based on the observed hydrograph.
The parameters that were consistent with the flood events of the 2000s in the observed hydrograph are shown in Table 6.The observed and simulated hydrographs of the 2000s are shown in Figure 10.The observed and simulated hydrographs of the 2000s are shown in Figure 10.The parameters that were consistent with the flood events of the 2010s in the observed In contrast to the hydrographs of the 1970s, 1980s, and 1990s, similar results were obtained in the 2000s.The shape and trend of the hydrographs of four flood events were similar.The total volume was slightly overestimated based on the observed hydrographs of 2002 and 2003.The peak was slightly overestimated based on the observed hydrographs in 2005 and 2007.
The parameters that were consistent with the flood events of the 2010s in the observed hydrograph are shown in Table 7.The observed and simulated hydrographs of the 2010s are shown in Figure 11.It can be seen that the fitness between the simulated and observed hydrographs was less precise than other periods, especially in 2012.This is because, after 2010, the forestland area changed greatly, resulting in great changes of parameters and their underlying characteristics.Moreover, in the SCS model, soil initial moisture was classified into dry, moderate, and wet; such a simple and coarse classification might also lead to increased error during the flood simulation process [6].
Model simulated results and evaluation criteria [47] for 18 flood events are summarized in Table 8.It can be seen from Table 8 that the peak flood absolute error between the observed and simulated peak discharge ranged from 0.56% to 18.49%, and the peak current difference ranged from 0 to 1 h.Moreover, the Nash coefficient was higher than 0.7, and the maximum value was 0.81.In summary, the HEC-HMS model was suitable for flood simulation of a loess hilly region in Gedong basin.

Calculation of Losses in Sub-Basins
Losses in the study area were calculated in the sub-basin for 1972, 1987, 1991, 1996, 2006, and 2012.They represented the underlying characteristics of different ages.The percentage of forestland in each sub-basin and sub-basin area is shown in Table 9.According to the "Standard for hydrological information and hydrological forecasting (GB/T 22482-2008)," in this work we proposed four model performance classes as a guidance on reference Nash coefficient range, which were denoted as Unsatisfactory (Nash coefficient < 0.5), Acceptable (0.5 ≤ Nash coefficient < 0.7), Good (0.7 ≤ Nash coefficient < 0.9), or Very good (Nash Water 2017, 9, 870 19 of 27 coefficient ≥ 0.9) [48].The Nash coefficients of 19 July 1972, 30 June 1987, 15 September 1991, 9 August 1996, 14 August 2006, and 30 July 2012 used for the calculation of rainfall losses, as shown in Table 10.Table 10 illustrates that the probability of the model fit being considered Unsatisfactory, Acceptable, Good, and Very good was 0%, 16.67%, 83.33%, and 0%, respectively.So the calculation of rainfall losses was feasible.
The results of losses calculation are shown in Tables 11-16.Table 11 demonstrates that the total precipitation volume generated in the whole basin was 590.47 mm and the loss volume was 474.63 mm, with an average percentage of 80.38% per sub-basin.The losses of W630, W520, W460, and W830 were below the average value.Table 12 showed that the total precipitation volume generated in the whole basin was 487.79 mm and the loss volume was 395.64 mm, with an average percentage of 81.11% per sub-basin.The losses of W630, W520, W460, and W830 were below average.Table 13 shows that the total precipitation volume generated in the whole basin was 623.08 mm and the loss volume was 516.81 mm, with an average percentage of 82.94% per sub-basin.The losses of W630, W520, W460, and W830 were below average value.
Table 14 shows that the total precipitation volume generated in the whole basin was 706.35 mm and the loss volume was 588.74 mm, with an average percentage of 83.35% per sub-basin.The losses of W630, W520, W460, and W830 were below average value.
Table 15 illustrates that the total precipitation volume generated in the whole basin was 499.77 mm and the loss volume was 439.7 mm, with an average percentage of 87.98% per sub-basin.The losses of W740, W630, W520, W460, and W830 were below average value.
Table 16 shows that the total precipitation volume generated in the whole basin was 680.47 mm and the loss volume was 612.87 mm, with an average percentage of 90.07%per sub-basin.The losses of W740, W630, W520, W460, and W830 were below average value.
The spatial distribution of the sub-basin rainfall losses is shown in Figure 12.
Water 2017, 9, 870 22 of 28 Table 15 illustrates that the total precipitation volume generated in the whole basin was 499.77 mm and the loss volume was 439.7 mm, with an average percentage of 87.98% per sub-basin.The losses of W740, W630, W520, W460, and W830 were below average value.Table 16 shows that the total precipitation volume generated in the whole basin was 680.47 mm and the loss volume was 612.87 mm, with an average percentage of 90.07%per sub-basin.The losses of W740, W630, W520, W460, and W830 were below average value.
The spatial distribution of the sub-basin rainfall losses is shown in Figure 12.
(a) (b) On the whole, with the increase of forestland percentage, the rainfall losses in the watershed increased from 1972 to 2012.In 1972, the variation range of the rainfall losses of sub-basins was 70.61% to 90.96% and the average percentage of the rainfall losses of sub-basins was 80.38%; in 1987, the variation range was 71.50% to 90.98% and the average percentage of the rainfall losses of sub-basins was 81.11%; in 1991, the variation range was 72.25% to 91.56% and the average percentage of the rainfall losses of sub-basins was 82.94%; in 1996, the variation range was 72.29% to 92.11% and the On the whole, with the increase of forestland percentage, the rainfall losses in the watershed increased from 1972 to 2012.In 1972, the variation range of the rainfall losses of sub-basins was 70.61% to 90.96% and the average percentage of the rainfall losses of sub-basins was 80.38%; in 1987, the variation range was 71.50% to 90.98% and the average percentage of the rainfall losses of sub-basins was 81.11%; in 1991, the variation range was 72.25% to 91.56% and the average percentage of the rainfall losses of sub-basins was 82.94%; in 1996, the variation range was 72.29% to 92.11% and the average percentage of the rainfall losses of sub-basins was 83.35%; in 2006, the variation range was 74.45% to 98.02% and the average percentage of the rainfall losses of sub-basins was 87.98%; in 2012, the variation range was 74.56% to 98.69% and the average percentage of the rainfall losses of sub-basins was 90.07%.

The Relationship between Rainfall Losses and Forestland Percentage and Slope
According to the experiment in Wangjiagou (a semi-arid hilly loess region of Shanxi Province of China), with an increase of slope angles, the runoff per unit area slightly increased on a short slope (7 m long), but decreased after reaching a maximum at 15 • and then decreased with slope angle on a long slope (20 m long), which might be related to the complicated effect of several factors (e.g., rainfall conditions, rill development) on soil infiltrability [16].In general, with the increase in slope, the rainfall losses gradually decreased.Moreover, for different rainfall periods, about 15 • below the slope, slope had a greater impact on infiltration, while a slope greater than 15 • had less influence on infiltration [49].The results revealed that the rainfall losses' decline was largest at 15 • , and the effect of slope on rainfall losses was complex.Thus, rainfall losses were influenced by forestland percentage and slope in Gedong basin.Multiple regression analysis was used to analyze the effects of forestland percentage and slope on rainfall losses in different sub-basins.The results suggested that the effect of forestland on rainfall losses was greater than that of slope, and rainfall losses increased as the forestland percentage increased and slope decreased.The regression equation was as follows: When the slope was [0, 15 where y represents rainfall losses, x 1 represents the surface slope, and x 2 represents forestland.
As shown in regression Equations ( 14) and ( 15), the result of multiple regression analysis was similar to the experimental result in loess hilly regions [49].

The Impact of Forestland Percentage on Rainfall Losses
In the Wangjiagou basin (a typical loess hilly region), Li Gang et al. [49] obtained the infiltration characteristics of different land types through a rainfall infiltration experiment.The steady infiltration rate of forestland was 0.96-0.99mm/min, and the steady infiltration rate of cultivated land was 0.39-0.83mm/min (the higher the slope of cultivated land, the lower the infiltration rate).This indicated that the infiltration rate of forestland was higher than that of cultivated land.So forestland was the main factor influencing rainfall losses in loess hilly regions.
According to the slope analysis, the areas of W740, W630, W520, W460, and W830 were all less than 75 km 2 , so the slope of these sub-basins was mainly 6-15 • .In these sub-basins, rainfall losses should be higher than in other sub-basins.However, the losses of W740, W630, W520, W460, and W830 were below the average value in Tables 11-16, which was caused by the percentage of forestland in these sub-basins being relatively lower than in others.
The correlation between rainfall losses and percentage of forestland in sub-basins was analyzed under the same rainfall level and different underlying conditions.According to the level of rainfall (30)(31)(32)(33)(34)(35)(36)(37)(38)(39)(40)  13 mm more than 9 August 1996, respectively.The results were consistent with the view that the rainfall losses had a positive correlation with forestland percentage.In contrast, the flood volume had a negative correlation with forestland percentage.The information of losses per sub-basin was similar to the information on flood contributing areas.The results of this information were helpful to determine the water availability of different areas.The correlation between rainfall losses and percentage of forestland in sub-basins is shown in Figure 13.From the 1970s to the 2010s, with the increase of forestland, the correlation coefficient between the percentage of rainfall losses and the percentage of forestland gradually increased.The correlation coefficient of 1972 was lower (0.64).In 2012, the correlation coefficient reached 0.84, which indicates that the rainfall losses of the basin had a positive correlation with the percentage of forestland.

Conclusions
This paper presents a case study of the modeling of the rainfall losses and flood processes in a typical loess hilly region in northern China.The SCS-CN model is used to compute the per sub-basin rainfall losses of Gedong basin in HEC-HMS.The model and flood simulation results are validated From the 1970s to the 2010s, with the increase of forestland, the correlation coefficient between the percentage of rainfall losses and the percentage of forestland gradually increased.The correlation coefficient of 1972 was lower (0.64).In 2012, the correlation coefficient reached 0.84, which indicates that the rainfall losses of the basin had a positive correlation with the percentage of forestland.

Conclusions
This paper presents a case study of the modeling of the rainfall losses and flood processes in a typical loess hilly region in northern China.The SCS-CN model is used to compute the per sub-basin rainfall losses of Gedong basin in HEC-HMS.The model and flood simulation results are validated and compared using observation data from different stages.Losses in the study area are calculated per sub-basin for different stages.The following conclusions can be drawn from the study: (1) From the 1970s to the 2010s, when the rainfall was 40-50 mm, the losses of sub-basins were in the range of 71.50%-98.02%;when the rainfall was 50-60 mm, the losses of sub-basins were in the range of 70.61%-91.56%;and when the rainfall was 60-70 mm, the losses of sub-basins were in the range of 72.29%-98.96%.(2) In the loess hilly region, rainfall losses are influenced by forestland percentage and slope.Multiple regression analysis suggests that the rainfall losses increase as the forestland area increases and the slope decreases, and the effect of forestland percentage on rainfall losses is greater than that of slope.From 1972 to 2012, with the increase of forestland percentage, the average percentage of rainfall losses in the Gedong basin increased from 80.38% to 90.07%.(3) Under the same rainfall level, the rainfall losses are positively correlated with forestland percentage.With the increase of forestland, the correlation coefficient between the percentage of rainfall losses and the percentage of forestland gradually increases.From the 1970s to the 2010s, the correlation coefficient ranged between 0.64 and 0.84.(4) Sensitivity analysis indicates that curve number (CN), initial abstraction (I a ), lag time (t lag ), and travel time (K) are the sensitive parameters determining rainfall losses in the basin.
The HEC-HMS is appropriate for modeling the flood processes in this typical loess hilly region in northern China.However, the simulation accuracy of the model may be affected by the limitations of the Thiessen polygon method in calculating rainfall.In the future, other spatial interpolation methods such as the Kriging method can be used to calculate the rainfall for improving the simulation accuracy of the model.

Figure 1 .
Figure 1.The location and area of Gedong basin.

Figure 1 .
Figure 1.The location and area of Gedong basin.

Figure 2 .
Figure 2. The land use types in Gedong basin.

Figure 2 .
Figure 2. The land use types in Gedong basin.

Figure 2 .
Figure 2. The land use types in Gedong basin.

Figure 3 .
Figure 3.The slope map of Gedong basin.

Figure 3 .
Figure 3.The slope map of Gedong basin.

Figure 4 .
Figure 4. Average slope varying with basin area.Figure 4. Average slope varying with basin area.

Figure 4 .
Figure 4. Average slope varying with basin area.Figure 4. Average slope varying with basin area.

Figure 4 .
Figure 4. Average slope varying with basin area.

Figure 5 .
Figure 5. Percentage of area under different slopes varying with basin area.

Figure 5 .
Figure 5. Percentage of area under different slopes varying with basin area.

Figure 6 .
Figure 6.Runoff and rainfall gauging station as well as Thiessen polygon in Gedong basin.

Figure 6 .
Figure 6.Runoff and rainfall gauging station as well as Thiessen polygon in Gedong basin.
41].Theoretically, if x1, x2,…, xn were model input variables and y = y(x1, x2,…, xn) was the model output, then the relative sensitivity (flexibility) of y with respect to the ith

Figure 8 .
Figure 8. Observed and simulated hydrographs of the 1970s and 1980s.

Figure 8 .
Figure 8. Observed and simulated hydrographs of the 1970s and 1980s.

Figure 9 .
Figure 9. Observed and simulated hydrographs of the 1990s.

Figure 10 .
Figure 10.Observed and simulated hydrographs of the 2000s.

Figure 11 .
Figure 11.Observed and simulated hydrographs of the 2010s.

Figure 11 .
Figure 11.Observed and simulated hydrographs of the 2010s.

Figure 12 .
Figure 12.Sub-basin rainfall losses of the study area.

Figure 13 .
Figure 13.The correlation between the percentage of rainfall loss and the percentage of forestland.

Figure 13 .
Figure 13.The correlation between the percentage of rainfall loss and the percentage of forestland.

Table 1 .
Land use information in different stages.

Table 3 .
Flexibility coefficients of the model.

Table 6 .
Calibrated parameters of the 2000s.

Table 6 .
Calibrated parameters of the 2000s.

Table 7 .
Calibrated parameters of the 2010s.

Table 8 .
Simulated results and evaluation criteria.

Table 8 .
Simulated results and evaluation criteria.

Table 9 .
Sub-basin area and forestland percentage.

Table 10 .
The Nash coefficients of flood events.