Simulating saturation‐excess surface run‐off in a semi‐distributed hydrological model

The level of complexity, and the number of parameters, to include in a hydrological model is a relatively contentious issue in hydrological modelling. However, it can be argued that explicitly representing important run‐off generation processes can improve the practical value of a model's outputs. This paper explores the benefits of including a new function into an existing semi‐distributed hydrological model (the Pitman model) that is widely used in the sub‐Saharan Africa region. The new function was designed to represent saturation‐excess surface run‐off processes at subcatchment scales and was motivated by the evidence of dambo (low topography riparian areas) type features in many sub‐Saharan river basins. The results for uncertainty versions of the model, with and without the new function, were compared for 25 catchments, which were divided up into those where evidence of dambos exists and those where there is no such evidence. The results suggest that the new function certainly improves the model results for the catchments where dambos exist, but not in situations where saturation‐excess surface run‐off is not expected to occur. The overall conclusion is therefore that the addition of the new function is justified.     

Do not only connect: a model of infiltration‐excess overland flow based on simulation

The paper focusses on connectivity in the context of infiltration‐excess overland flow and its integrated response as slope‐base overland flow hydrographs. Overland flow is simulated on a sloping surface with some minor topographic expression and spatially differing infiltration rates. In each cell of a 128 × 128 grid, water from upslope is combined with incident rainfall to generate local overland flow, which is stochastically routed downslope, partitioning the flow between downslope neighbours. Simulations show the evolution of connectivity during simple storms. As a first approximation, total storm runoff is similar everywhere, discharge increasing proportionally with drainage area. Moderate differences in plan topography appear to have only a second‐order impact on hydrograph form and runoff amount. Total storm response is expressed as total runoff, runoff coefficient or total volume infiltrated; each plotted against total storm rainfall, and allowing variations in average gradient, overland flow roughness, infiltration rate and storm duration. A one‐parameter algebraic expression is proposed that fits simulation results for total runoff, has appropriate asymptotic behaviour and responds rationally to the variables tested. Slope length is seen to influence connectivity, expressed as a scale distance that increases with storm magnitude and can be explicitly incorporated into the expression to indicate runoff response to simple events as a function of storm size, storm duration, slope length and gradient. The model has also been applied to a 10‐year rainfall record, using both hourly and daily time steps, and the implications explored for coarser scale models. Initial trails incorporating erosion continuously update topography and suggest that successive storms produce an initial increase in erosion as rilling develops, while runoff totals are only slightly modified. Other factors not yet considered include the dynamics of soil crusting and vegetation growth. Copyright © 2014 John Wiley & Sons, Ltd.     

Run‐off processes from mountains to foothills: The role of soil stratigraphy and structure in influencing run‐off characteristics across high to low relief landscapes

The critical zone features that control run‐off generation, specifically at the regional watershed scale, are not well understood. Here, we addressed this knowledge gap by quantitatively and conceptually linking regional watershed‐scale run‐off regimes with critical zone structure and climate gradients across two physiographic provinces in the Southeastern United States. We characterized long‐term (~20 years) discharge and precipitation regimes for 73 watersheds with United States Geological Survey in‐stream gaging stations across the Appalachian Mountain and Piedmont physiographic provinces of North Carolina. Watersheds included in this analysis had <10% developed land and ranged in size from 14.1–4,390 km². Thirty‐four watersheds were located in the Piedmont physiographic province, which is typically classified as a low relief landscape with deep, highly weathered soils and regolith. Thirty‐nine watersheds were located in the Appalachian Mountain physiographic province, which is typically classified as a steeper landscape with highly weathered, but shallower soils and regolith. From the United States Geological Survey daily mean run‐off time series, we calculated annual and seasonal baseflow indices (BFI), minimum, mean, and maximum daily run‐off, and Pearson's correlation coefficients between precipitation and baseflow. Our results showed that Appalachian Mountain watersheds systematically had higher minimum daily flows and BFI values. Piedmont watersheds displayed much larger deviations from mean annual BFI in response to year‐to‐year variability in precipitation. A series of linear regression models between 21 landscape metrics and annual BFIs showed non‐linear and complex terrestrial–hydrological relationships across the two provinces. From these results, we discuss how distinct features of critical zone architecture, with specific focus on soil depth and stratigraphy, may be dominating the regulation of hydrological processes and run‐off regimes across these provinces.     

Estimating daily inflows of large lakes using a water‐balance‐based runoff coefficient scaling approach

Daily river inflow time series are highly valuable for water resources and water environment management of large lakes. However, the availability of continuous inflow data for large lakes is still relatively limited, especially for large lakes situated within humid plain regions with tens or even hundreds of tributaries. In this study, we choose the fifth largest freshwater Lake Chaohu in China as our study area to introduce a new approach to reconstruct historical daily inflows at ungauged subcatchments of large lakes. This approach makes use of water level, lake surface rainfall, evaporation from the lake, and catchment rainfall observations. Rainfall–runoff relationship at a reference catchment was analysed to select rainfall input and estimate run‐off coefficient firstly, and the run‐off coefficient was then transferred to ungauged subcatchments to initially estimate daily inflows. Run‐off coefficient was scaled to adjust daily inflows at ungauged subcatchments according to water balance of the lake. This approach was evaluated using sparsely measured inflows at eight subcatchments of Lake Chaohu and compared with the commonly used drainage area ratio method. Results suggest that the inflow time series reconstructed from this approach consistent well to corresponding observations, with mean R² and Nash–Sutcliffe efficiency values of 0.69 and 0.6, respectively. This approach outperforms drainage area ratio method in terms of mean R² and Nash–Sutcliffe efficiency values. Accuracy of this approach holds well when the number of water‐level station being used decreased from four to one.     

Hydrologic changes resulting from urban cover in seasonally snow‐covered catchments

There are few multibasin analyses of the effects of urban land cover on seasonal stream flow patterns within northern watersheds where winter snow cover is the norm. In this study, the effects of urban cover on stream flow were evaluated at nine catchments in southern Ontario, Canada, which vary greatly in urban impervious cover (1–84%) but cluster into two groups having ≥54% urban impervious area (‘urban’) and ≤11% impervious cover (‘rural’), respectively. Annual and seasonal run‐off totals (millimetres) were similar between the rural and urban groups and were relatively insensitive to urban cover. Instead, urban streams had significantly greater high flow frequency, flow variability and quickflow and lower baseflow compared with rural streams. Furthermore, differences in high flow frequency between urban and rural stream groups were largest in the summer and fall and less extreme in the winter and spring, perhaps because of the homogenizing effect of winter snow cover, frozen ground and spring melt on surface imperviousness. Although the clear clustering of streams into urban and rural groups precluded the identification of a threshold above which urban cover is the primary cause of flow differences, relatively high extreme flow frequency and flow variability in the two most urbanized of the rural streams (10–11% impervious) suggest that it may lie close to this range. Furthermore, whereas total run‐off volumes were not affected by urban cover, increases in stream flashiness and a greater frequency of high flow events particularly during the summer and fall may negatively impact stream biota and favour the transfer of surface‐deposited pollutants to urban streams. Copyright © 2014 John Wiley & Sons, Ltd.     

Effective run‐off flow length over biological soil crusts on silty loam soils in drylands

Surface hydrological behaviour is important in drylands because it affects the distribution of soil moisture and vegetation and the hydrological functioning of slopes and catchments. Microplot scale run‐off can be relatively easily measured, i.e. by rainfall simulations. However, slope or catchment run‐off cannot be deduced from microplots, requiring long‐time monitoring, because run‐off coefficients decrease with increasing drainage area. Therefore, to determine the slope length covered by run‐off (run‐off length) is crucial to connect scales. Biological soil crusts (BSCs) are good model systems, and their hydrology at slope scale is insufficiently known. This study provides run‐off lengths from BSCs, by field factorial experiments using rainfall simulation, including two BSC types, three rain types, three antecedent soil moistures and four plot lengths. Data were analysed by generalized linear modelling, including vascular plant cover as covariates. Results were the following: (i) the real contributing area is almost always much smaller than the topographical contributing area; (ii) the BSC type is key to controlling run‐off; run‐off length reached 3 m on cyanobacterial crust, but hardly over 1 m on lichen crust; this pattern remained through rain type or soil moisture; (iii) run‐off decreased with BSC development because soil sealing disappears; porosity, biomass and roughness increase and some changes occur in the uppermost soil layer; and (iv) run‐off flow increased with both rain type and soil moisture but run‐off coefficient only with soil moisture (as larger rains increased both run‐off and infiltration); vascular plant cover had a slight effect on run‐off because it was low and random. Copyright © 2014 John Wiley & Sons, Ltd.     

Experimental and analytical studies of solute transport during run‐off over vegetated surfaces

Solute transport in overland flow is considered as one of the main contributors to water pollution. Although many models of pollutant transport mechanism from soil to run‐off water have been proposed, the characteristics of solute transport accompanying the water run‐off over vegetated surface have not been well studied. In this study, a series of laboratory experiments were conducted to study the solute transport over vegetated surfaces. Based on the experimental results, an idea of the “stationary water layer” in run‐off was proposed. Applying the complete mixing theory in the stationary water layer, an analytical solute transport model was developed with the assumption that the upper run‐off completely mixes with the underlying water in the stationary water layer for each site. The results show that the predictions made by the present model are in good agreement with the measured experimental data. For the vegetated surfaces, the depth of stationary water layer is related to the rainfall intensity, bed slope, and vegetation density. The analytical solution shows that the maximum solute transport occurs at the time of concentration. This study advances our understanding of the mechanisms of solute transport over vegetated areas.     

Runoff‐generating processes in hydromorphic soils on a plot scale: free gravity‐driven versus pressure‐controlled flow

Draining soil water is an important runoff generator. This study aims to describe runoff‐generating processes on a plot scale (1 m²) in hydromorphic soils with different initial soil water contents. We irrigated 16 hydromorphic soils in the northern Pre‐Alps in Switzerland and recorded the variations in water content with time domain reflectometry (TDR) at five different depths per plot. Sprinkling was repeated three times at approximately 23‐h intervals and lasted for 1 h with a volume flux density of 70 mm h^?1. The comparison between the measured water content of the drainages with two physically based models revealed which of the flow processes dominated during water recessions. We distinguished between vertical drainage, lateral outflow and infiltration without drainage. Approximately 45% of all recorded time series of soil water content did not drain within approximately 20 h after the end of irrigation, about 25% drained laterally and 10% of the outflow was vertical. The drainage of the remaining 20% was the result of both lateral and vertical water flow (≈12%), or was not interpretable with the approaches applied (≈8%). Vertical flow was only observed in layers without any or with just a few hydromorphic features. Lateral draining horizons had approximately half the storage capacity and amplitude of water recession of those with vertical flow. Vertical flow was only observed in the upper soil horizons. Thus, vertical flow transmitted water to layers with lateral outflow and did not delay runoff by deep percolation. Increasing initial soil moisture increased the percentage of water content recordings according to a lateral outflow slightly, while vertical flow was less frequent. Copyright © 2010 John Wiley & Sons, Ltd.     

Modelling hydrological response to a fully‐monitored urban bioretention cell

Municipalities and agencies use green infrastructure to combat pollution and hydrological impacts (e.g., flooding) related to excess stormwater. Bioretention cells are one type of infiltration green infrastructure intervention that infiltrate and redistribute otherwise uncontrolled stormwater volume. However, the effects of these installations on the rest of the local water cycle is understudied; in particular, impacts on stormwater return flows and groundwater levels are not fully understood. In this study, full water cycle monitoring data were used to construct and calibrate a two‐dimensional Richards equation model (HYDRUS‐2D/3D) detailing hydrological implications of an unlined bioretention cell (Cleveland, Ohio) that accepts direct runoff from surrounding impervious surfaces. Using both preinstallation and postinstallation data, the model was used to (a) establish a mass balance to determine reduction in stormwater return flow, (b) evaluate green infrastructure effects on subsurface water dynamics, and (c) determine model sensitivity to measured soil properties. Comparisons of modelled versus observed data indicated that the model captured many hydrological aspects of the bioretention cell, including subsurface storage and transient groundwater mounding. Model outputs suggested that the bioretention cell reduced stormwater return flows into the local sewer collection system, though the extent of this benefit was attenuated during high inflow events that may have exhausted detention capacity. The model also demonstrated how, prior to bioretention cell installation, surface and subsurface hydrology were largely decoupled, whereas after installation, exfiltration from the bioretention cell activated a new groundwater dynamic. Still, the extent of groundwater mounding from the cell was limited in spatial extent and did not threaten other subsurface infrastructure. Finally, the sensitivity analysis demonstrated that the overall hydrological response was regulated by the hydraulics of the bioretention cell fill material, which controlled water entry into the system, and by the water retention parameters of the native soil, which controlled connectivity between the surface and groundwater.     

One‐ and two‐dimensional modelling of overland flow in semiarid shrubland,Jornada basin,New Mexico

Two distributed parameter models, a one‐dimensional (1D) model and a two‐dimensional (2D) model, are developed to simulate overland flow in two small semiarid shrubland watersheds in the Jornada basin, southern New Mexico. The models are event‐based and represent each watershed by an array of 1‐m² cells, in which the cell size is approximately equal to the average area of the shrubs. Each model uses only six parameters, for which values are obtained from field surveys and rainfall simulation experiments. In the 1D model, flow volumes through a fixed network are computed by a simple finite‐difference solution to the 1D kinematic wave equation. In the 2D model, flow directions and volumes are computed by a second‐order predictor–corrector finite‐difference solution to the 2D kinematic wave equation, in which flow routing is implicit and may vary in response to flow conditions. The models are compared in terms of the runoff hydrograph and the spatial distribution of runoff. The simulation results suggest that both the 1D and the 2D models have much to offer as tools for the large‐scale study of overland flow. Because it is based on a fixed flow network, the 1D model is better suited to the study of runoff due to individual rainfall events, whereas the 2D model may, with further development, be used to study both runoff and erosion during multiple rainfall events in which the dynamic nature of the terrain becomes an important consideration. Copyright © 2006 John Wiley & Sons, Ltd.     

Long‐term and seasonal hydrologic performance of an extensive green roof

Sustainable strategies such as green roofs have been implemented as stormwater management tools to mitigate disturbance of the hydrologic cycle resulting from urbanization. Green roofs, also referred to as vegetated roofs, can improve the urban landscape by reducing heat island effects, providing ecosystem services, and facilitating the retention and treatment of stormwater. Green roofs have received particular attention because they do not require acquisition and development of land and represent an application of biomimicry in design and construction. In this paper, we evaluate the effects of precipitation, evapotranspiration (ET), antecedent dry period (ADP), and seasonal variation on the run‐off quantity and distribution of an extensive, sedum covered, green roof on a commercial building in Syracuse, NY, USA. The green roof greatly facilitated retention of precipitation events without significant changes over the 4‐year study. The green roof retained on average 95.9 ± 3.6% (6.5 ± 5.6 mm) per rainfall event, with a range from 75% to 99.6% (33.2 to 3.3 mm). However, as precipitation quantity increased, the retention of water decreased. This high water retention capacity was the result of the combined effects of ET, stormwater storage (plants, growth media, and stormwater retention layer), and limited surface run‐off from the roof deck due to variation in the sloping of the green roof and the tapered insulation to the deck drains. The water retention capacity of the green roof did not change significantly between growing and nongrowing seasons. Slightly greater precipitation during the growing season coincided with increased ET. Average potential ET during the growing season was approximately 3 times greater than during the nongrowing season. The hydrologic performance of the green roof was not significantly impacted by an ADP greater than 2 days.     

Carbon loss by water erosion in drylands: implications from a study of vegetation change in the south‐west USA

Soil organic carbon (SOC) is an important component of the global carbon cycle yet is rarely quantified adequately in terms of its spatial variability resulting from losses of SOC due to erosion by water. Furthermore, in drylands, little is known about the effect of widespread vegetation change on changes in SOC stores and the potential for water erosion to redistribute SOC around the landscape especially during high‐magnitude run‐off events (flash floods). This study assesses the change in SOC stores across a shrub‐encroachment gradient in the Chihuahuan Desert of the south‐west USA. A robust estimate of SOC storage in surface soils is presented, indicating that more SOC is stored beneath vegetation than in bare soil areas. In addition, the change in SOC storage over a shrub‐encroachment gradient is shown to be nonlinear and highly variable within each vegetation type. Over the gradient of vegetation change, the heterogeneity of SOC increases, and newer carbon from C₃ plants becomes dominant. This increase in the heterogeneity of SOC is related to an increase in water erosion and SOC loss from inter‐shrub areas, which is self‐reinforcing. Shrub‐dominated drylands lose more than three times as much SOC as their grass counterparts. The implications of this study are twofold: (1) quantifying the effects of vegetation change on carbon loss via water erosion and the highly variable effects of land degradation on soil carbon stocks is critical. (2) If landscape‐scale understanding of carbon loss by water erosion in drylands is required, studies must characterize the heterogeneity of ecosystem structure and its effects on ecosystem function across ecotones subject to vegetation change. Copyright © 2013 John Wiley & Sons, Ltd.     

Modelling soil solute release and transport in run‐off on a loessial slope with and without surface stones

Stone covers on loessial slopes can increase the time of infiltration by slowing the velocity of the overland flow, which reduces the transport of solutes, but few mechanistic models have been tested under water‐scouring conditions. We carried out field experiments to test a previously proposed, physically based model of water and solute transport. The area of soil infiltration was calculated from the uncovered surface area, and Richards' equation and the kinematic wave equation were used to describe water infiltration and flow along slopes with stone covers. The transport of chemicals into the run‐off from the surface soil, presumably by diffusion, and their movement in the soil profile could be described by the convection–diffusion equations of the model. The simulated and measured data correlated well. The stones on the soil surface reduced the area available for infiltration but increased the Manning coefficient, eventually leading to increased water infiltration and decreased solute loss with run‐off. Our results indicated that the traditional model of water movement and solute migration could be used to simulate water transport and solute migration for stone‐covered soil on loessial slopes.     

Regression‐based approach to low flow prediction in the Maryland Piedmont region under joint climate and land use change

We examine the low flow records for six urbanized watersheds in the Maryland Piedmont region and develop regression equations to predict annual minimum low flow events. The effects of both future climate (based on precipitation and temperature projections from two climate models: Hadley and the Canadian Climate Centre (CCC)) and land use change are incorporated to illustrate possible future trends in low flows. A regression modelling approach is pursued to predict the minimum annual 7‐day low flow estimates for the proposed future scenarios. A regional regression model was calibrated with between 10 and 50 years of daily precipitation, daily average temperature, annual imperviousness, and the daily observed flow time‐series across six watersheds. Future simulations based on a 55 km² urbanizing watershed just north of Washington, DC, were performed. When land use and climate change were employed singly, the former predicted no trends in low flows and the latter predicted significant increasing trends under Hadley and no trends under CCC. When employed jointly, however, low flows were predicted to decrease significantly under CCC, whereas Hadley predicted no significant trends in low flows. Antecedent precipitation was the most influential predictor on low flows, followed by urbanization. Copyright © 2006 John Wiley & Sons, Ltd.     

An analytical study on three‐dimensional versus two‐dimensional water table‐induced flow patterns in a Tóthian basin

Although it has been increasingly acknowledged that groundwater flow pattern is complicated in the three‐dimensional (3‐D) domain, two‐dimensional (2‐D) water table‐induced flow models are still widely used to delineate basin‐scale groundwater circulation. However, the validity of 2‐D cross‐sectional flow field induced by water table has been seldom examined. Here, we derive the analytical solution of 3‐D water table‐induced hydraulic head in a Tóthian basin and then examine the validity of 2‐D cross‐sectional models by comparing the flow fields of selected cross sections calculated by the 2‐D cross‐sectional model with those by the 3‐D model, which represents the “true” cases. For cross sections in the recharge or discharge area of the 3‐D basin, even if head difference is not significant, the 2‐D cross‐sectional models result in flow patterns absolutely different from the true ones. For the cross section following the principal direction of groundwater flow, although 2‐D cross‐sectional models would overestimate the penetrating depth of local flow systems and underestimate the recharge/discharge flux, the flow pattern from the cross‐sectional model is similar to the true one and could be close enough to the true one by adjusting the decay exponent and anisotropy ratio of permeability. Consequently, to determine whether a 2‐D cross‐sectional model is applicable, a comparison of hydraulic head difference between 2‐D and 3‐D solutions is not enough. Instead, the similarity of flow pattern should be considered to determine whether a cross‐sectional model is applicable. This study improves understanding of groundwater flow induced by more natural water table undulations in the 3‐D domain and the limitations of 2‐D models accounting for cross‐sectional water table undulation only.     

Surface moisture area during rainfall–run‐off events to understand the hydrological dynamics of a basin in a plain region

The understanding of the hydrology of plain basins may be improved by the combined analysis of rainfall–run‐off records and remote sensed surface moisture data. Our work evaluates the surface moisture area (SMA) produced during rainfall–run‐off events in a plain watershed of the Argentine Pampas Region, and studies which hydrological variables are related to the generated SMA. The study area is located in the upper and middle basins of the Del Azul stream, characterized by the presence of small gently hilly areas surrounded by flat landscapes. Data from 9 rainfall–run‐off events were analysed. MODIS surface reflectance data were processed to calculate SMA subsequent to the peak discharge (post‐SMA), and previous to the rainfall events (prev‐SMA), to consider the antecedent wetness. Rainfall–run‐off data included total precipitation depth (P), maximum intensity of rainfall over 6 hr (I6max), surface run‐off registered between the beginning of the event and the day previous to the analysed MODIS scene (R), peak flow (Qp), and flood intensity (IF). In contrast with other works, post‐SMA showed a negative relationship with the R. Three groups of cases were identified: (a) Events of low I6max, high prev‐SMA, and low R were associated with slow and weakly channelized flow over plain areas, leading to saturated overland flow (SOF), with large SMA; (b) events of high I6max, low prev‐SMA, and medium to high R were rapidly transported along the gentle slopes of the basin, related to Hortonian overland flow (HOF) and low post‐SMA; and (c) events of medium to high I6max and prev‐SMA with medium R were related to heterogeneous input‐antecedent‐run‐off conditions combined: Local spatial conditions may have produced HOF or SOF, leading to an averaged response with medium SMA. The interactions between the geomorphology of the basin, the characteristics of the events, and the antecedent conditions may explain the obtained results. This analysis is relevant for the general knowledge of the hydrology of large plains, whose functioning studies are still in their early stages.     

Sensitivity of streamflow simulation in the Delaware River Basin to forecasted land‐cover change for 2030 and 2060

In order to simulate the potential effect of forecasted land‐cover change on streamflow and water availability, there has to be confidence that the hydrologic model used is sensitive to small changes in land cover (<10%) and that this land‐cover change exceeds the inherent uncertainty in forecasted conditions. To investigate this, a 26‐year streamflow record was simulated for 33 basins (54–928 km²) in the Delaware River Basin using three dates of land cover: the 2011 National Land‐Cover Dataset (Homer, Fry, & Barnes, 2012 ), 2030 land‐cover conditions representing median values from 101 equally‐likely forecasts, and 2060 land‐cover conditions corresponding to the same iterations used to represent 2030. Streamflow was simulated using a process‐based hydrologic model that includes both pervious and impervious methods as parameterized by three land‐cover‐based hydrologic response units (HRUs)—forested, agricultural, and developed land. Small, but significant differences in streamflow magnitude, variability, and seasonality were seen among the three time periods—2011, 2030, and 2060. Temporal differences were discernible from the range of conditions simulated with 101 equally likely forecasts for 2030. Development was co‐located with the most frequent landscape components, as characterized by topographic wetness index, resulting in a change in hydrology for each HRU, highlighting that knowing the location of disturbance is key to understanding potential streamflow changes. These results show that streamflow simulation using regional calibration that incorporates land‐cover‐based HRUs can be sensitive to relatively small changes in land‐cover and that temporal trends resulting from land‐cover change can be isolated in order to evaluate other changes that might affect water resources.     

Evaluation of distributed BMPs in an urban watershed—High resolution modeling for stormwater management

The urban environment modifies the hydrologic cycle resulting in increased runoff rates, volumes, and peak flows. Green infrastructure, which uses best management practices (BMPs), is a natural system approach used to mitigate the impacts of urbanization onto stormwater runoff. Patterns of stormwater runoff from urban environments are complex, and it is unclear how efficiently green infrastructure will improve the urban water cycle. These challenges arise from issues of scale, the merits of BMPs depend on changes to small‐scale hydrologic processes aggregated up from the neighborhood to the urban watershed. Here, we use a hyper‐resolution (1 m), physically based hydrologic model of the urban hydrologic cycle with explicit inclusion of the built environment. This model represents the changes to hydrology at the BMP scale (~1 m) and represents each individual BMP explicitly to represent response over the urban watershed. Our study varies both the percentage of BMP emplacement and their spatial location for storm events of increasing intensity in an urban watershed. We develop a metric of effectiveness that indicates a nonlinear relationship that is seen between percent BMP emplacement and storm intensity. Results indicate that BMP effectiveness varies with spatial location and that type and emplacement within the urban watershed may be more important than overall percent.     

A finite element model for simulating surface run‐off and unsaturated seepage flow in the shallow subsurface

This work introduces water–air two‐phase flow into integrated surface–subsurface flow by simulating rainfall infiltration and run‐off production on a soil slope with the finite element method. The numerical model is formulated by partial differential equations for hydrostatic shallow flow and water–air two‐phase flow in the shallow subsurface. Finite element computing formats and solution strategies are presented to obtain a numerical solution for the coupled model. An unsaturated seepage flow process is first simulated by water–air two‐phase flow under the atmospheric pressure boundary condition to obtain the rainfall infiltration rate. Then, the rainfall infiltration rate is used as an input parameter to solve the surface run‐off equations and determine the value of the surface run‐off depth. In the next iteration, the pressure boundary condition of unsaturated seepage flow is adjusted by the surface run‐off depth. The coupling process is achieved by updating the rainfall infiltration rate and surface run‐off depth sequentially until the convergence criteria are reached in a time step. A well‐conducted surface run‐off experiment and traditional surface–subsurface model are used to validate the new model. Comparisons with the traditional surface–subsurface model show that the initiation time of surface run‐off calculated by the proposed model is earlier and that the water depth is larger, thus providing values that are closer to the experimental results.