Next Article in Journal
Optimal Spatial Design of Capacity and Quantity of Rainwater Harvesting Systems for Urban Flood Mitigation
Next Article in Special Issue
Spatio-Temporal Patterns of Water Table and Vegetation Status of a Deserted Area
Previous Article in Journal
Yaya-Beas, R.-E., et al. Helminth Egg Removal Capacity of UASB Reactors under Subtropical Conditions. Water 2015, 7, 2402–2421
Previous Article in Special Issue
Vapor Flow Resistance of Dry Soil Layer to Soil Water Evaporation in Arid Environment: An Overview
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Case Report

Flood-Runoff in Semi-Arid and Sub-Humid Regions, a Case Study: A Simulation of Jianghe Watershed in Northern China

1
College of Water Resources Science and Engineering, Taiyuan University of Technology, Taiyuan 030024, China
2
Institute of Water Resources and Environmental Geology, Taiyuan University of Technology, Taiyuan 030024, China
3
Institute of Hydraulic Engineering and Water Resources Management, RWTH Aachen University, Aachen D 52056, Germany
4
School of Computing, Science & Engineering, University of Salford, Manchester M5 4WT, UK
*
Author to whom correspondence should be addressed.
Water 2015, 7(9), 5155-5172; https://doi.org/10.3390/w7095155
Submission received: 6 July 2015 / Revised: 5 September 2015 / Accepted: 14 September 2015 / Published: 22 September 2015

Abstract

:
This paper presents a modeling application of surface runoff using the Hydrologic Modelling System (HEC-HMS). A case study was carried out for the Jianghe watershed, a typical semi-arid and sub-humid geo-climatic region in northern China. Two modeling schemes using different descriptive sub-mechanism models provided by HEC-HMS for runoff volume, direct runoff and routing (channel flow) were investigated. The modeling results were compared with historical observation data. This work shows that HEC-HMS can be a suitable modeling tool for specific situations in China. With the appropriate selection of the sub-mechanism models, HEC-HMS can be applied to various situations, including the typical semi-arid and sub-humid conditions in northern China.

1. Introduction

Surface runoff or flood runoff is the flow of water over the ground surface. This happens when soil is fully saturated or heavy rain falls over a short period of time such that it cannot be completely absorbed by the ground surface. Surface runoff plays an important role in the water cycle and resource distribution in geological ecosystems. It also has a direct effect on soil erosion and consequently, flooding. There are many factors that influence surface runoff, including rainfall dynamics, ground geo-topography, types of soil and vegetation, land use and development, etc. [1,2].
In China, the modeling of surface runoff has been studied since the early 1960s. Specific models for different climatic and environmental conditions have been proposed [3], such as the Xinanjiang-model used for the humid regions [4,5,6,7] and the Shanbei-model used for the arid regions in China [3,8]. However, these two models are not able to represent all of the typical geo-climatic characteristics of the vast and diverse territory of China. For example, in northern China, where the climate presents both semi-arid and sub-humid characteristics, the rainfall and runoff processes are much more complicated than those that occur in solely humid or arid regions.
The Hydrologic Modeling System (HEC-HMS) is an integrated modeling and simulation tool for all hydrologic processes of dendritic watershed systems. It consists of different sub-mechanism models for rainfall loss, direct runoff and routing, i.e., the component processes involved in the transfer of rainfall to runoff. This paper reports a case study of the application of the HEC-HMS in the flood modeling of a typical semi-arid and sub-humid region, the Jianghe watershed, in northern China. It aims to investigate the suitability of the sub-mechanism models provided by HEC-HMS and to identify an effective modeling strategy using the appropriate sub-models for flooding processes under these climatic and environmental conditions.

2. Study Area and Data

The Jianghe watershed is located in the midwest of Changzhi city, Shanxi Province, China (Figure 1). It comprises a catchment area of about 270 Km2, with an elevation range from 977 to 1544 m. The altitude is high in the north, west and south, but low in the east and the center. The watershed is located in a typical semi-arid and sub-humid climatic region in northern China, where the annual average precipitation is approximately 564 mm, with a recorded maximum of 872.2 mm in 1971 and a minimum of 299.1 mm in 1983. The precipitation in a single year varies substantially but typically, intensive rainfalls occur in summer (June to September). The annual average water surface evaporation is about 1775.8 mm and the annual average temperature is about 9.1 °C.
The watershed was discretized into 12 sub-basins as shown in Figure 2. Digital Elevation Model (DEM) data with a 30-m resolution, obtained from the Computer Network Information Center, Chinese Academy of Sciences, were employed to characterize the stream network and the topographic attributes (e.g., area, length and slope) of each sub-basin based on a 1:50,000 topographical map. The characterized topographical data are listed in Table 1.
The land use data were obtained from the United States Geological Survey. The main type of land in the Jianghe watershed is grassland of low to medium coverage, comprising about half of the total area. The other land uses are forest, shrub, open forest and flood plain. The land use data are listed in Table 2. The soil type data were obtained from the China Soil Science Database. The major type in the region is cinnamon soil belonging to the C soil type according to the Soil Conservation Service (SCS) classification.
Figure 1. The location and area of Jianghe watershed.
Figure 1. The location and area of Jianghe watershed.
Water 07 05155 g001
Table 1. The Topographical Parameters of the Sub-basins.
Table 1. The Topographical Parameters of the Sub-basins.
Sub-basins No.Land PlaneStream/Channel
Area (km2)Length (m)Slope (m/m)Length (m)Slope (m/m)
117.152870.36456180.022
212.336000.26837300.016
335.757130.38353790.017
412.639490.39737290.038
535.561260.26753740.008
613.143580.34720090.057
720.243450.34828910.045
88.223350.28623150.013
932.449830.35265800.022
1055.957560.36356410.014
117.731800.42418780.045
1219.342140.36448890.010
Table 2. Land Use Types of Jianghe Watershed.
Table 2. Land Use Types of Jianghe Watershed.
No.Land UseGridArea
km2Percentage %
1Forest land11,14310.143.76
2Shrub land34,65531.5511.69
3Open forestland47,19242.9615.91
4Low coverage grassland69,04662.8623.28
5Middle coverage grassland76,94270.0525.94
6Flood plain57,58752.4319.42
Total296,565270100
Figure 2. The hydrological model implemented in HEC-HMS.
Figure 2. The hydrological model implemented in HEC-HMS.
Water 07 05155 g002
The precipitation and discharge data were obtained from the Hydrology and Water Resource Survey Administration of Changzhi, Shanxi Province, China. There are five observation stations evenly distributed in the catchment area as shown in Figure 1. They are Lizhuang, Zhongcun, Xishangzhuang and Baquan and Beizhangdian; Beizhangdian also occurs within the stream gauge station. The precipitation data from the five stations were weighted at the center of each sub-basin using an inverse-distance-squared method [9]. The historical precipitation and discharge data show that rainfall and flooding in the Jianghe catchment area present three main characteristics, i.e.:
  • Flooding formed by localized rainstorms: the runoff was mainly generated by infiltration-excess when the flooding had a rapid rise and recession with a single peak value but the total flood volume was small.
  • Flooding formed by uniform rainfall across the whole area: the runoff was dominated by saturation-excess when the flooding changed little over the period of time and the total flood volume was large.
  • Flooding formed by mixed rainfalls: the runoff was due to both infiltration- and saturation-excess when the flooding possessed the two previous characteristics.

3. Modeling Schemes and Methodology

HEC-HMS has previously been used to simulate the rainfall-runoff processes in the arid, semi-arid, humid and sub-humid geo-climatic environments in many countries [10,11,12,13]. It has optional models for the sub-mechanisms. The following models are used to estimate the runoff volume: the Initial and Constant-rate model, the Deficit and Constant-rate model, the Soil Conservation Service curve number (SCS-CN) model, and the Green and Ampt model [9,14,15,16]. For the direct runoff, there is the SCS Unit Hydrography (SCS-UH) model, the Snyder’s UH model, the Clark’s UH model and the Kinematic-wave model [9,14,15]. For the routing (channel flow), there is the Lag model, the Muskingum model, the Modified Kinematic-wave model and the Muskingum Cunge model [9].
Previous studies in China have shown that the runoff volume models such as the Initial and Constant-Rate and the Soil Conservation Service curve number (SCS-CN) models, and the direct runoff models such as the Kinematic wave and the SCS Unit Hydrograph (SCS-UH) models have been used successfully to model flooding [17,18,19,20]. The Muskingum model is the conventional model for channel flow in arid, semi-arid, humid and sub-humid regions.
According to the summarized characteristics of the rainfall and flood, the geo-topography, the soil types and the land use, two schemes shown in Table 3 were planned for the investigation. Given that the flood durations are relatively short and the scale of the individual sub-basins to be studied is relatively small, the influence of the sub-ground or base flow was neglected in the study. Jianghe watershed has been under observation for more than 50 years. A previous study [21] found that the condition of the river channels was stable with no significant change. Based on this, in the current study, the Muskingum model was calibrated using Scheme I only. The calibrated parameters were directly employed in Scheme II later. The details of the models employed in the two schemes are described in the following sections.
Table 3. Planned modeling analysis.
Table 3. Planned modeling analysis.
Modeling SchemeRunoff-Volume ModelDirect-Runoff ModelRouting Model
ІInitial and constant-rateKinematic waveMuskingum
IISoil Conservation Service curve number (SCS-CN)SCS Unit Hydrograph (SCS-UH)Muskingum

3.1. Scheme І

3.1.1. The Runoff-Volume Model: Initial and Constant-Rate

The runoff-volume model of the initial and constant-rate has been successfully applied in a fair number of previous studies [22,23,24]. The basic concept of the model is that the maximum potential infiltration rate, fc, is constant throughout a flood event. It includes an initial loss, Ia, to reflect the interception and surface retention. The rainfall-excess or precipitation-excess, pet, is defined by Equation (1):
p e t = { 0      i f   p i < I a p t f c i f   p i < I a   a n d   p t > f c 0    i f   p i > I a   a n d   p t > f c
where, p t is the average precipitation during a time interval Δt, pi is the observation data at time ti, and p i is the total precipitation observed in the flooding period. In this study, the initial value of f c was selected from a range of 1.3–3.9 mm/h, while I a was within a 10%–20% range of the total precipitation according to [25]. It was also assumed that there was no impervious area in terms of the land uses discussed previously.

3.1.2. The Direct-Runoff Model: Kinematic Wave

The Kinematic-wave model is for the surface flow. This model uses Equation (2) as the flow governing Equation:
{ h t + q = i e   q = 1 n S 0 1 / 2 h 5 / 3
where h is the water height or depth above the surface, t is the time, q is the discharge per unit width, i e is the net rainfall, n is the Manning’s roughness coefficient, and S0 is the land slope.
To estimate direct-runoff using the Kinematic-wave model, the catchment area needs to be discretized into series of elementary components including the planes of surface flow, sub-channels of collection, channels of collection and the main channels. In this study, each sub-basin was characterized into one surface flow plane and one main channel (trapezoidal area). The corresponding parameters are the length, slope and roughness (N) for planes, and the length, slope, bottom width (W), side slope (K) and Manning’s roughness coefficient (n) for channels. The length and slope data are listed in Table 1. The values of the parameters W and K were estimated through calibration. The values of the parameters N and n were estimated following [25] and [26] within the range of 0.025–0.05.

3.1.3. The Routing Model: Muskingum

The Muskingum model is used to describe the outflows from the reaches of flows using the following Equation:
O 2 = C 0 I 2 + C 1 I 1 + C 2 O 1
{ C 0 = 0.5 t K x K K x + 0.5 t C 1 = 0.5 t + K x K K x + 0.5 t C 2 = K K x 0.5 t K K x + 0.5 t
C 0 + C 1 + C 2 = 1
where, I1 and I2 are the inflows to a flow reach at the start and the end of a period of time, respectively, O1 and O2 are the outflows from the reach at the same moment, ∆t is the length of the period of time, K is the travel time throughout the reach area, and x is the Muskingum weighting factor. In this study, the initial value of K was evaluated using the following definition:
K = L V w
where L is the length of flow reaches, Vw is the velocity of the flood waves, which is about 1.33–1.67 times the normal flow rate [25]. The initial value of the factor x is in the range of 0–0.5 according to statistical analysis.

3.2. Scheme II

3.2.1. The Runoff-Volume Model: SCS Curve Number

The SCS-CN model assumes that the accumulated rainfall-excess depends upon the cumulative precipitation, soil type, land use and the previous moisture conditions [27]. This model defines the accumulated rainfall-excess using the following equation:
p e = ( p I a ) 2 p I a + S
where pe is the accumulated rainfall-excess, p is the accumulated rainfall depth at a particular moment, Ia is the initial loss, and S is the potential maximum retention (water infiltrated into the ground). An empirical linear relationship between Ia and S is defined as Ia = 0.2S. The maximum retention, S, is estimated using the following Equation:
S = 25400 C N 254
where CN (called the SCS curve number) is used to represent the combined effects of the primary characteristics of the catchment area, including soil type, land use and the previous moisture condition. It takes a value in the range of 30–98.

3.2.2. The Direct-Runoff Model: SCS Unit Hydrograph

The SCS Unit Hydrograph (UH) is a parametric model based on the average UH derived from gauged rainfall and runoff data of a large number of small agricultural watersheds throughout the US.
A dimensionless UH discharge, Ut, which is defined as a ratio of the momentary discharge to the single peak discharge Up, is related to another dimensionless time parameter, t, which is a fraction of the momentary time and the time of the peak discharge, Tp. The Up and Tp are defined by the two equations below:
U p = C A T p
T p = Δ t 2 +   t l a g
where A is the watershed catchment area, C is a conversion constant, Δt is the excess precipitation duration (computational interval in HEC-HMS), tlag is the basin lag, i.e., the time difference between the center of the mass of rainfall excess and the peak of the UH.

4. Calibration and Modeling

In general, flood analysis tended to use large events. However, considering that the scale of the events should reflect the range of the applicability and performance of a model [28], this study selected thirteen flood events recorded in the period from 1980 to 1998 for the modeling test. The thirteen events had scales from small to large with flood peak discharges in a range of 28.5–461 m3/s. The precipitation data from the five stations were weighted at the center of each sub-basin in terms of an inverse-distance-squared method [9]. Thereafter, the precipitation data were interpolated into a sequential data series with 15-minute time intervals. Finally, the data series was used for the Time-series Database of the HEC-HMS.
Seven flood events in the period from 1980 to 1990 were selected for model calibration to determine the characteristic parameters. Thereafter, the determined parameters were used to predict the other six flood events recorded from 1990 to 1998. The characteristic parameters are the fc for the Initial and Constant-rate model, CN for the SCS-CN model, N, W, K and n for the Kinematic-wave model, the lag time, tlag, for the SCS-UH model, and the travel time K and weighting factor x for the Muskingum model.
The flood-volume, the flood-peak discharge, the peak-time and the Nash efficiency coefficient are the four objective variables assessed in the calibration. The flood volume relative error, REv, the flood peak discharge relative error, REp, and the time difference of the peak appearance, ∆T, are defined by Equations (9)–(11), respectively, while the Nash efficiency coefficient, DC, is defined by Equation (12).
R E v = Q s Q 0 Q 0 × 100 %
R E p = q s q 0 q 0 × 100 %
T = T s T 0
D C = 1 i = 1 n ( q s ( i ) q 0 ( i ) ) 2 i = 1 n ( q 0 ( i ) q 0 , m e a n ) 2 × 100 %
where Qs is the modeled flood volume, Q0 is the observed flood volume, qs is the modeled value of peak discharge and q0 is the observed peak discharge, Ts is the modeled time of the peak appearance and T0 is the observed time of the peak appearance, i indicates the sequence of intervals, and q0, mean is the average flood volume observed in an event. The best result is that the REv and REp have the lowest absolute values, while ∆T is small and DC is close to 1.
Trial-and-error and the objective function are two conventional methods adopted for calibration. The trial-and-error is a manual control method, which has an advantage to be able to follow the physical meaning of the parameters defined in the models, but it is time-consuming and labor-intensive. The objective function is an automatic control method, which optimizes the characteristic parameters for a targeted result in terms of a predefined algorithm. The algorithm may generate unrealistic values of the parameters, which are beyond the range of physical meaning. In this work, a combined strategy was adopted for model calibration. At first, it used the objective function to obtain a set of the optimized values of the parameters. When the obtained values of the parameters did not meet the physical meaning, they were further revised through the trial-and-error method. The physical meaning of the hydrological parameters, such as fc, K and x, is based on the sub-basin scale, while the physical meaning of the geographic parameters, such as N, W, K and n, is based on the land use grid data scale.
The HEC-HMS provides two algorithms and seven objective functions for optimization, in which the Nelder and Mead algorithm and the Peak-weighted Root Mean Square Error objective function were employed in this study due to their simplicity and good performance [29]. The results of the parameter calibration are listed in Table 4, Table 5 and Table 6.
Table 4. The Calibration Result of the Run-off Volume and the Direct-runoff, Scheme І.
Table 4. The Calibration Result of the Run-off Volume and the Direct-runoff, Scheme І.
Sub-Basin No.Initial and Constant-Rate ModelKinematic Wave Model
fc (mm/h)NW (m)K (m/m)n
11.30.451.50.120.031
21.60.191.70.130.035
30.80.124.50.210.028
41.70.581.50.130.043
51.60.196.30.250.037
61.30.253.90.270.032
71.50.217.10.320.027
81.30.1612.50.430.032
90.70.272.40.150.026
101.90.244.20.210.031
111.50.731.30.150.026
121.70.225.30.180.035
Table 5. The Calibration Result of the Routing, Schemes I and II.
Table 5. The Calibration Result of the Routing, Schemes I and II.
Muskingum ModelReach 1Reach 2Reach 3Reach 4Reach 5Reach 6Reach 7Reach 8
K(h)1.82.30.81.22.01.40.60.5
x0.400.250.350.400.300.350.300.30
Table 6. The Calibration Result of the Run-off Volume and the Direct-runoff, Scheme II.
Table 6. The Calibration Result of the Run-off Volume and the Direct-runoff, Scheme II.
Sub-Basins No.SCS-CN ModelSCS UH Model
CNLag Time (min)
17530
28320
38832
46522
58335
66825
78525
89013
97828
106532
115518
127723

5. Results and Discussion

Table 7 and Table 8 list the modeling results of the flood volume, peak discharge, their relative errors with respect to the observation data, the difference between the modeled and the observed time of peak appearance, and the Nash efficiency coefficients for the two modeling schemes. Figure 3, Figure 4, Figure 5 and Figure 6 show the modeling and observed variation in the flow rate over time.
Table 7. The Modeling Results of Scheme I.
Table 7. The Modeling Results of Scheme I.
Recorded EventsQsQ0REvqsq0REp∆TDC
104 m3104 m3%m3/sm3/s%min
Events for Calibration
10 August 198058.267.2−13.429.732.4−8.3240.86
27 July 198271.787.9−18.5104.3123−15.2220.87
15 July 198582.294.3−12.851.459.1−13.1480.74
23 July 198579.992.8−13.998.9103−3.9250.92
14 July 198760.674.8−18.951.857.2−9.3150.91
18 July 1988579.5933.2−37.9333.7461−27.6450.84
13 June 199060.867.19.633.6366.5540.79
Events for Prediction
11 July 1990128.7118.88.383.698.8−15.4300.77
16 August 1991177.5169.74.6248.1300−17.3260.89
3 July 199343.1947.5−9.345.750−8.5340.93
28 August 199551.961.5−15.636.138.9−7.2230.79
12 July 199648.856.1−12.826.128.5−8.5280.93
08 July 1998332.7437.2−23.942.949.8−13.8350.87
Table 8. The Modeling Results of Scheme II.
Table 8. The Modeling Results of Scheme II.
Recorded EventsQsQ0REvqsq0REp∆TDC
104 m3104 m3%m3/sm3/s%min
Events for Calibration
10 August 198072.967.28.627.732.4−14.3−730.81
27 July 198281.387.9−7.5118.6123−3.5180.86
15 July 1985102.794.38.960.959.13.2−150.71
23 July 1985100.492.88.293.4103−9.3−250.81
14 July 198767.574.8−9.746.757.2−18.2−650.76
18 July 19881152.4933.223.5547.246118.7−350.71
13 June 199073.467.19.437.1363.2−150.88
Events for Prediction
11 July 1990129.5118. 98.9101.698.82.8−250.93
16 August 1991150169.7−11.6290.4300−3.2−120.78
3 July 199351.447.58.252.5504.3−150.94
28 August 199556.161.6−8.936.838.9−5.5−480.71
12 July 199661.956.110.630.628.57.3−360.77
8 July 1998529.9437.221.252.449.85.3−960.84
In the modeling, the parameters fc and N showed a high sensitivity in Scheme I, while the CN showed a high sensitivity in Scheme II. Meanwhile, the Muskingum parameter K showed a high sensitivity in both schemes. These results demonstrate that the soil types, land use, basin topography and the river length play a major role in the flood run-off processes.

5.1. The Runoff-Volume Models

Schemes I and II employ the Initial and Constant-rate model and the SCS-CN model, respectively, to simulate runoff generation. According to the result of the flood volume relative error, REv, it can be seen that the error of Scheme I is relatively small for the four flood events, 13 June 1990, 11 July 1990, 16 August 1991, and 3 July 1993, with absolute values less than 10%, but the error increases significantly for the events of 27 July 1982, 14 July 1987, 28 July 1995, and 08 July 1998, with absolute values greater than 15%. However, the error of Scheme II has a small variation and a relatively small average absolute value for most of these events.
This can be explained by the following: due to the semi-arid and sub-humid attributes, the flood runoff is dominated by infiltration-excess under the condition of short but intensive rainfall, or by combined infiltration- and saturation-excess under the condition of long-lasting rainfall of various intensities. In Figure 3, Figure 4, Figure 5 and Figure 6, it can be seen that the floods on 13 June 1990, 11 July 1990 and 3 July 1993 were caused by rainfall of short duration, while the flood on 16 August 1991 was caused by rainfall with a strong intensity. The relative error of the flood volume of these events in Scheme I is relatively small and shows that the Initial and Constant-rate model is effective for modeling the infiltration-excess runoff. However, the increase in the error on 27 July 1982, 14 July 1987, 28 July 1995 and 8 July 1998 indicates that there was additional saturation-excess runoff, a fact that can been seen in Figure 3, Figure 4, Figure 5 and Figure 6, which shows that the four events had a long, continuous duration. Although the SCS-CN model in Scheme II demonstrated a better performance for runoff volume prediction, it needs to be pointed out that the SCS-CN is an empirical model, which does not consider extreme situations.
Figure 3. Flow Rate vs Time, Scheme I Calibration.
Figure 3. Flow Rate vs Time, Scheme I Calibration.
Water 07 05155 g003

5.2. The Direct-Runoff and the Routing Models

Schemes I and II employ the Kinematic-wave model and the SCS-UH model, respectively, to simulate direct runoff. However, both schemes use the Muskingum model to simulate the routing (channel flow). For the floods on 15 July 1985, 13 June 1990, 11 July 1990, 16 August 1991 and 3 July 1993, which were short in duration and strong in intensity, as shown in Figure 3, Figure 4, Figure 5 and Figure 6, the difference between the predicted and observed (ΔT) peak time was within the range of 26–54 min in Scheme I compared with 15–25 min in Scheme II. For the floods on 10 August 1980, 14 July 1987, 28 August 1995 and 8 July 1998, which had long durations, the difference was within the range of 15–35 min using Scheme I and 48–96 min using Scheme II. The results show that the SCS-UH model of Scheme II is suitable for the situations of infiltration-excess, while the Kinematic-wave model of Scheme I is suitable for the situations of infiltration- and saturation-excess for semi-arid and sub-humid environmental conditions.
Figure 4. Flow Rate vs Time, Scheme I Prediction.
Figure 4. Flow Rate vs Time, Scheme I Prediction.
Water 07 05155 g004
Figure 5. Flow rate vs time, scheme II calibration.
Figure 5. Flow rate vs time, scheme II calibration.
Water 07 05155 g005
Figure 6. Flow Rate vs Time, Scheme II Prediction.
Figure 6. Flow Rate vs Time, Scheme II Prediction.
Water 07 05155 g006

5.3. The Overall Comparison of the Two Modeling Schemes

From the modeling results shown in Table 7 and Table 8, it can be seen that 84.6% of the flood volume predictions of Scheme I had relative errors with absolute values less than 20%, while 92.3% of the predictions of Scheme II had relative errors less than 20%. For the flood peak discharge, 84.6% of the predictions of Scheme I had relative errors with absolute values less than 20%, while 100% of the predictions of Scheme II had relative errors with absolute values less than 20%. Meanwhile, the DCs were higher than 0.7 for all prediction cases. The results show that both schemes are viable to describe the flood-runoff processes in the Jianghe area. Overall, the average absolute values of the four performance indicator parameters in Scheme I are REv = 15.34, REp = 11.89, ΔT = 31.46 and DC = 0.85, compared with REv = 9.67, REp = 7.6, ΔT = 36.77 and DC = 0.8 in Scheme II. Overall, Scheme II performed better than Scheme I in the case study of the Janghe watershed.

6. Conclusions

This paper presents a case study of the modeling of the rainfall-runoff processes in a typical semi-arid and sub-humid region in northern China. Two schemes using the modeling tool HEC-HMS were investigated for the suitability and accuracy of the sub-component models provided. The proposed hydrological model is based on the hydrological characteristics, geo-topography, soil type and land use in the area. The models and predictions were validated and compared using observation data. The following four conclusions can be drawn from the study:
  • The HEC-HMS is a suitable tool for modeling the rainfall-runoff processes in the Jianghe watershed, which is located in a typical region of semi-arid and sub-humid climate in northern China.
  • The SCS-CN model performed better than the initial and constant-rate model in the estimation of runoff generation.
  • The kinematic wave model demonstrated a relatively good result for the prediction of the flood peak time for situations of long-lasting rainfall duration with varied intensity, while the SCS-UH model demonstrated a relatively good result for situations with short rainfall duration but high intensity.
  • An appropriate combination of the volume generation model and the direct runoff model is capable of producing reliable modeling estimates using the HEC-HMS for various applications. In general, Scheme II can be used for semi-arid and sub-humid environmental conditions.

Acknowledgments

The research is funded by Shanxi Province Science Foundation (No. 2010011030-1), China.

Author Contributions

This work was conducted in collaboration between all authors. Hua Jin led the project and defined the research theme. Hua Jin and Rui Liang designed methods and modelling, conducted the modelling investigations. Hua Jin, Rui Liang and Yu Wang analyzed the data, interpreted the results and wrote the paper. Prasad Tumula contributed to discussion, analysis, interpretation. All authors have contributed to the revision and approved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Beven, K.; Robert, E. Horton’s perceptual model of infiltration processes. Hydrol. Process. 2004, 18, 3447–3460. [Google Scholar] [CrossRef]
  2. McMahon, T.A.; Finlayson, B.L.; Haines, A.T.; Srikanthan, R. Global Runoff: Continental Comparisons of Annual Flows and Peak Discharges; Catena Verlag: Cremlingen-Destedt, Germany, 1992. [Google Scholar]
  3. Zhao, R.J. Watershed Hydrological Model-Xinanjiang Model and Shanbei Model; Water and Power Press: Beijing, China, 1984. (In Chinese) [Google Scholar]
  4. Song, X.M.; Kong, F.Z.; Zhan, C.S.; Han, J.W. Hybrid optimization rainfall-runoff simulation based on Xinanjiang model and artificial neural network. J. Hydrol. Eng. 2012, 17, 1033–1041. [Google Scholar] [CrossRef]
  5. Zhu, Q.A.; Zhang, W.C. The applicability study of Xinanjiang model on simulation of rainfall-runoff and flooding hydrographs in Jiangkou basin. J. Water Resour. Water Eng. 2004, 15, 19–23. [Google Scholar]
  6. Cheng, C.T.; Ou, C.P.; Chau, K.W. Combining a fuzzy optimal model with a genetic algorithm to solve multi-objective rainfall-runoff model calibration. J. Hydrol. 2002, 268, 72–86. [Google Scholar] [CrossRef]
  7. Zhao, R.J. The Xinanjiang model applied in China. J. Hydrol. 1992, 135, 371–381. (In Chinese) [Google Scholar]
  8. Dai, R.X.; Li, L.; Liu, X.M.; Wang, W.; Wu, W. Application of a hydrology model based on DEM in Ju River Watershed. China Rural. Water Hydropower 2007, 3, 1–3. (In Chinese) [Google Scholar]
  9. United States Army Corps of Engineers (USACE). Hydrologic Modeling System HEC-HMS; United States Army Corps of Engineers: Washington, DC, USA, 2000. [Google Scholar]
  10. Unucka, J.; Adamec, M. Modeling of the land cover impact on the rainfall-runoff relations in the Olse catchment. J. Hydrol. Hydromech. 2008, 56, 257–271. [Google Scholar]
  11. Cydzik, K.; Hogue, T.S. Modeling postfire response and recovery using the hydrologic engineering center hydrologic modeling system (HEC-HMS). J. Am. Water Resour. Assoc. (JAWRA) 2009, 45, 702–714. [Google Scholar] [CrossRef]
  12. Sherif, M.; Mohamed, M.; Shetty, A.; Almulla, M. Rainfall-runoff modeling of three wadis in the Northern area of UAE. J. Hydrol. Eng. 2011, 16, 10–20. [Google Scholar] [CrossRef]
  13. Meenu, R.; Rehana, S.; Mujumdar, P.P. Assessment of hydrologic impacts of climate change in Tunga–Bhadra river basin, India with HEC-HMS and SDSM. Hydrol. Process. 2013, 27, 1572–1589. [Google Scholar] [CrossRef]
  14. De Silva, M.M.G.T.; Weerakoon, S.B.; Herath, S. Modeling of event and continuous flow hydrographs with HEC-HMS: Case study in the Kelani River Basin, Sri Lanka. J. Hydrol. Eng. 2014, 19, 800–806. [Google Scholar] [CrossRef]
  15. Kaffas, K.; Ssanthou, V. Application of a continuous rainfall-runoff model to the basin of Kosynthos River using the hydrologic software HEC-HMS. Glob. NEST J. 2014, 16, 188–203. [Google Scholar]
  16. Kowalik, T.; Walega, A. Estimation of CN parameter for small agricultural watersheds using asymptotic functions. Water 2015, 7, 939–955. [Google Scholar] [CrossRef]
  17. Wang, Z. Distributed Hydrological Model Research in Upper Fenhe Watershed. Master’s Thesis, Taiyuan University of Technology, Taiyuan, China, May 2003. [Google Scholar]
  18. Zhang, J.J.; Na, L.; Zhang, B. Applicability of the distributed hydrological model of HEC-HMS in a small watershed of the Loess Plateau area. J. Beijing For. Univ. 2009, 31, 52–57. (In Chinese) [Google Scholar]
  19. Li, Z.B.; Lu, K.X. Approximate analytical solution of rainfall runoff process on permeable slope. J. Hydraul. Eng. 2009, 6, 8–21. (In Chinese) [Google Scholar]
  20. Chen, F.; Lin, F.; Chen, X.W. Application of HEC-HMS distributed hydrological model to the rainflood simulation in Jinjiang river basin. J. Huaqiao Univ. Nat. Sci. 2012, 33, 325–329. (In Chinese) [Google Scholar]
  21. Ren, H.L. Study on the Hydrological Model of the Semi-Arid and Semi-humid Region—An Example of the Beizhangdian Watershed. Master’s Thesis, Taiyuan University of Technology, Taiyuan, China, May 2006. [Google Scholar]
  22. Halwatura, D.; Najim, M.M.M. Application of the HEC-HMS model for runoff simulation in a tropical catchment. Environ. Model. Softw. 2013, 46, 155–162. [Google Scholar] [CrossRef]
  23. Du, J.K.; Li, Q.; Rui, H.Y.; Zuo, T.H.; Zheng, D.P.; Xu, Y.P.; Xu, C.Y. Assessing the effects of urbanization on annual runoff and flood events using an integrated hydrological modeling system for Qinhuai River basin. China J. Hydrol. 2012, 464, 127–139. [Google Scholar] [CrossRef]
  24. Ali, M.; Khan, S.J.; Aslam, I.; Khan, Z. Simulation of the impacts of land-use change on surface runoff of Lai Nullah Basin in Islamabad, Pakistan. Landsc. Urban Plan. 2011, 102, 271–279. [Google Scholar] [CrossRef]
  25. United States Army Corps of Engineers (USACE). Engineering and Design Flood-Runoff Analysis; United States Army Corps of Engineers: Washington, DC, USA, 1994. [Google Scholar]
  26. Wu, C.G. Hydraulics; Higher Education Press: Beijing, China, 2008. (In Chinese) [Google Scholar]
  27. United States Department of Agriculture. Estimation of direct runoff from storm rainfall. In National Engineering Handbook; United States Department of Agriculture: Washington, DC, USA, 2004. [Google Scholar]
  28. Reed, S.; Koren, V.; Smith, M.; Zhang, Z.Y.; Moreda, F.; Seo, D.J.; Participants, D. Overall distributed model intercomparison project results. J. Hydrol. 2004, 298, 27–60. [Google Scholar] [CrossRef]
  29. Deng, X.; Dong, X.H.; Bo, H.J. Research on influence of objective function on HEC-HMS model parameter calibration. Water Resour. Power 2010, 28, 17–19. (In Chinese) [Google Scholar]

Share and Cite

MDPI and ACS Style

Jin, H.; Liang, R.; Wang, Y.; Tumula, P. Flood-Runoff in Semi-Arid and Sub-Humid Regions, a Case Study: A Simulation of Jianghe Watershed in Northern China. Water 2015, 7, 5155-5172. https://doi.org/10.3390/w7095155

AMA Style

Jin H, Liang R, Wang Y, Tumula P. Flood-Runoff in Semi-Arid and Sub-Humid Regions, a Case Study: A Simulation of Jianghe Watershed in Northern China. Water. 2015; 7(9):5155-5172. https://doi.org/10.3390/w7095155

Chicago/Turabian Style

Jin, Hua, Rui Liang, Yu Wang, and Prasad Tumula. 2015. "Flood-Runoff in Semi-Arid and Sub-Humid Regions, a Case Study: A Simulation of Jianghe Watershed in Northern China" Water 7, no. 9: 5155-5172. https://doi.org/10.3390/w7095155

APA Style

Jin, H., Liang, R., Wang, Y., & Tumula, P. (2015). Flood-Runoff in Semi-Arid and Sub-Humid Regions, a Case Study: A Simulation of Jianghe Watershed in Northern China. Water, 7(9), 5155-5172. https://doi.org/10.3390/w7095155

Article Metrics

Back to TopTop