The particle tracking method (PTM) module was added into the open source Full Shallow-Water equations for Overland Flow in a two-dimensional (FullSWOF_2D) program, which has coupled rainfall–runoff and infiltration modules to determine the time of concentration (Tc) for impervious (Tci) and pervious (Tcp) surfaces. The updated program FullSWOF-PTM was tested using observed rainfall events with Nash–Sutcliffe efficiencies ranging from 0.60 to 0.95 (average of 0.75) for simulated runoff hydrographs. More than 400 impervious modeling cases with different surface slope (S0), roughness coefficient (n), length (L), and rainfall intensity (i) combinations were developed and simulated to obtain the Tci for developing the regression equation of Tci as a function of the four input parameters. More than 700 pervious modeling cases with different combinations of S0, n, L, i, and infiltration parameters including the saturated hydraulic conductivity, suction head, and moisture deficit were simulated to estimate the Tcp based on the travel time of 85% of particles arriving at the outlet and the ponding time. The regression equation of Tcp was developed as the sum of Tci and additional travel time as a function of infiltration parameters and i. The Tcp equation can be applied to wide ranges of input parameters in comparison to Akan’s equation.
An updated two-dimensional flow simulation program, FullSWOF-ZG, which fully (Full) solves shallow water (SW) equations for overland flow (OF) and includes submodules modeling infiltration by zones (Z) and flow interception by grate-inlet (G), was tested with 20 locally depressed curb inlets to validate the inlet efficiency (Eci), and with 80 undepressed curb inlets to validate the inlet lengths (LT) for 100% interception. Previous curb inlet equations were based on certain theoretical approximations and limited experimental data. In this study, 1000 road-curb inlet modeling cases from the combinations of 10 longitudinal slopes (S0, 0.1–1%), 10 cross slopes (Sx, 1.5–6%), and 10 upstream inflows (Qin, 6–24 L/s) were established and modeled to determine LT. The second 1000 modeling cases with the same 10 S0 and 10 Sx and 10 curb inlet lengths (Lci, 0.15–1.5 m) were established to determine Eci. The LT and Eci regression equations were developed as a function of input parameters (S0, Sx, and Qin) and Lci/LT with the multiple linear regression method, respectively. Newly developed regression equations were applied to 10,000 inlet design cases (10 S0, 10 Sx, 10 Qin, and 10 Lci combinations) and comprehensively compared with three equations in previous studies. The 100% intercepted gutter flow (Qg100) equations were derived, and over-prediction of Qg100 from previous methods was strongly correlated to smaller S0. Newly developed equations gave more accurate estimations of LT and Eci over a wide range of input parameters. These equations can be applied to designing urban drainage and road bioretention facilities, since they were developed using a large number of simulation runs with diverse input parameters, but previous methods often overpredict the gutter flow of total interception when the longitudinal slope S0 is small.
In the flat lowland agricultural areas of subtropical climate zones, the runoff process has a great influence on the regional water quantity and quality. In this study, field data about rainfall, evapotranspiration, soil moisture, groundwater table, and surface water dynamics were collected in two different experimental sites in the Taihu Basin, China. Results showed that densely distributed ditches contributed to shallow groundwater depths and persistent near-saturation soil. A correlation analysis was conducted to improve the understandings of runoff generation in humid lowland areas of the Taihu Basin. It was found that a Dunne overland flow was the dominant mechanism responsible for the rapid runoff generation. The total rainfall and runoff expressed a good linear relationship with an R2 of 0.95 in the Hongqiwei test site. The initial groundwater depth was considered as the indicator of the antecedent soil moisture estimation for the close relationship. The depression storage was suggested in a range from 4.72 to 8.03 mm for an estimation based on the water balance analysis for each rainfall event, which proves that the depression storage should not be neglected when calculating the runoff generation process in humid lowlands.
The two-dimensional overland flow simulation program, FullSWOF_2D, was revised to include submodules of determining infiltration by zones (Z) and grate-inlet (G) drainage from a 2D surface to a 1D pipe flow. The updated program, FullSWOF-ZG, was used to evaluate the performance of a road-bioretention strip (RBS) system and explore/understand key parameters of continuous RBS design. The program was validated using eight pervious surfaces under simulated rainfall events and tested with 20 experimental cases of a locally depressed curb inlet. The mean difference of simulated interception efficiencies (36.6%–86.0%) and observed interception efficiencies (34.8%–84.0%) of the curb inlet was 3.5%, which proves the program predicts the curb-inlet interception efficiency accurately. The 20 road-only and 20 RBS modeling cases were designed and modeled using the FullSWOF-ZG program. These case studies have different road lengths, curb inlet lengths, longitudinal slopes, cross slopes, bioretention-overflow inlet heights, and bioretention soil infiltration parameters. Only 34.6%–48.4% of the total runoff volume is intercepted by the RBS’s curb inlet under heavy rainfall (250 mm/h) and the remaining part of the runoff flows downstream along the road, which may cause local inundation and become a safety hazard. The curb inlet becomes the bottleneck of the RBS system that could impede the runoff flowing into the bioretention strip for detention and infiltration to improve the stormwater quality.
A distributed-framework hydrologic modeling system (DF-HMS) is a primary and significant component of a distributed-framework basin modeling system (DFBMS), which simulates the hydrological processes and responses after rainfall at the basin scale, especially for non-homogenous basins. The DFBMS consists of 11 hydrological feature units (HFUs) involving vertical and horizontal geographic areas in a basin. Appropriate hydrologic or hydraulic methods are adopted for different HFUs to simulate corresponding hydrological processes. The digital basin generation model is first developed to determine the essential information for hydrologic and hydraulic simulation. This paper mainly describes two significant HFUs contained in the DF-HMS for hydrologic modeling: Hilly sub-watershed and plain overland flow HFUs. A typical hilly area application case study in the Three Gorges area is introduced, which demonstrates DF-HMS’s good performance in comparison with the observed streamflow at catchment outlets.