A two-layer model for studying 2D dissolved pollutant runoff over impermeable surfaces

Dissolved pollutants in stormwater are a main contributor to water pollution in urban environments. However, many existing transport models are semi-empirical and only consider one-dimensional flows, which limit their predictive capacity. Combining the shallow water and the advection–diffusion equations, a two-dimensional physically based model is developed for dissolved pollutant transport by adopting the concept of a ‘control layer’. A series of laboratory experiments has been conducted to validate the proposed model, taking into account the effects of buildings and intermittent rainfalls. The predictions are found to be in good agreement with experimental observations, which supports the assumption that the depth of the control layer is constant. Based on the validated model, a parametric study is conducted, focusing on the characteristics of the pollutant distribution and transport rate over the depth. The hyetograph, including the intensity, duration and intermittency, of rainfall event has a significant influence on the pollutant transport rates. The depth of the control layer, rainfall intensity, surface roughness and area length are dominant factors that affect the dissolved pollutant transport. Finally, several perspectives of the new pollutant transport model are discussed. This study contributes to an in-depth understanding of the dissolved pollutant transport processes on impermeable surfaces and urban stormwater management.     

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.     

Kinematic wave solutions for pollutant transport by runoff over an impervious plane,with instantaneous or finite‐period mixing

Kinematic wave solutions are derived for transport of a conservative non‐point‐source pollutant during a rainfall‐runoff event over an impervious plane. Rainfall is assumed to be steady, uniform and finite in duration. Prior to the start of rainfall, the pollutant is distributed uniformly over the plane. When rainfall occurs, the pollutant is washed off in a particular manner and the mixing of pollutant in the runoff water occurs either instantaneously or in a finite period of time under the assumption that the pollutant is soluble and is mixed completely in the runoff water. Copyright © 2002 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.     

Kinematic wave solutions for pollutant transport over an infiltrating plane with finite‐period mixing and mixing zone

Kinematic wave solutions are derived for transport of a conservative non‐point‐source pollutant during a rainfall‐runoff event over an infiltrating plane for two cases: (i) finite‐period mixing and (ii) soil‐mixing zone. Rainfall is assumed to be steady, uniform and finite in duration, and it is assumed to have zero concentration of pollutants. Infiltration is assumed constant in time and space. Prior to the start of rainfall, the pollutant is distributed uniformly over the plane. In the first case, when rainfall occurs, the mixing of pollutant in the runoff water occurs in a finite period of time. In the second case, the chemical concentration is assumed to be a linearly decreasing function of rainfall intensity and overland flow. The solute concentration and discharge are found to depend on the flow characteristics as well as the solute concentration parameters. The characteristics of solute concentration and discharge graphs seem to be similar to those reported in the literature and observed in laboratory experiments. Copyright © 2002 John Wiley & Sons, Ltd.     

Water repellency as run‐off and soil detachment controlling factor in a dry‐Mediterranean hillslope (South of Spain)

This study examines the effect of water repellency on controlling temporal variability of runoff generation mechanisms and soil detachment on metamorphic derived soils under dry‐Mediterranean climate. The research is carried out in an unburnt Mediterranean hillslope in souther Spain characterized by a patchy vegetation pattern and shallow soils. The Water Drop Penetration Time test (WDPT) is applied to measure water repellency at the end of summer (Sep‐2008), mid autumn (Nov‐2008) and mid winter (Feb‐2009). Rainfall simulations were used to obtain runoff generation and soil detachment in the same periods of time. The main shrub specie is Cistus monspeliensis which leaves a load of litter during the summer due to the lack of water. This great amount of organic material is accumulated under the shrubs triggering an extreme water repellence (WDPT > 6,000 s) that limits infiltration processes. This process is enforced due to the low soil water content at the end of dry season. Certain water repellency (WDPT > 1,500 s) is also observed on bare soil as consequence of their sandier texture and the accumulation of annual plants which die at the end of the wet season. Soil moisture increases during the autumn and water repellency disappears in both shrub and bare soil at the middle of the wet season (WDPT < 5 s). The main consequence is that the temporal trend of water repellency controls the mechanism and frequency of runoff generation and, hence, soil detachment. At the end of the summer, Hortonian mechanisms predominates when water repellency is extreme, even in soils under Cistus monspeliensis where runoff generation can reach higher peaks of overland flow and sediment concentration. Conversely, only the saturation of soil could generate runoff during the wet season being this quite less frequent in bare soil and absent in shrub. Copyright © 2010 John Wiley & Sons, Ltd.     

Comparing theoretically supported rainfall‐runoff erosivity factors at the Sparacia (South Italy) experimental site

Interpreting rainfall‐runoff erosivity by a process‐oriented scheme allows to conjugate the physical approach to soil loss estimate with the empirical one. Including the effect of runoff in the model permits to distinguish between detachment and transport in the soil erosion process. In this paper, at first, a general definition of the rainfall‐runoff erosivity factor REF_e including the power of both event runoff coefficient Q_R and event rainfall erosivity index EI₃₀ of the Universal Soil Loss Equation (USLE) is proposed. The REF_e factor is applicable to all USLE‐based models (USLE, Modified USLE [USLE‐M] and Modified USLE‐M [USLE‐MM]) and it allows to distinguish between purely empirical models (e.g., Modified USLE‐M [USLE‐MM]) and those supported by applying theoretical dimensional analysis and self‐similarity to Wischmeier and Smith scheme. This last model category includes USLE, USLE‐M, and a new model, named USLE‐M based (USLE‐MB), that uses a rainfall‐runoff erosivity factor in which a power of runoff coefficient multiplies EI₃₀. Using the database of Sparacia experimental site, the USLE‐MB is parameterized and a comparison with soil loss data is carried out. The developed analysis shows that USLE‐MB (characterized by a Nash–Sutcliffe Efficiency Index NSEI equal to 0.73 and a root mean square error RMSE = 11.7 Mg ha^?1) has very similar soil loss estimate performances as compared with the USLE‐M (NSEI = 0.72 and RMSE = 12.0 Mg ha^?1). However, the USLE‐MB yields a maximum discrepancy factor between predicted and measured soil loss values (176) that is much lower than that of USLE‐M (291). In conclusion, the USLE‐MB should be preferred in the context of theoretically supported USLE type models.     

A rainfall‐runoff model parameterized from GIS and runoff data

A guiding principle in hydrological modelling should be to keep the number of calibration parameters to a minimum. A reduced number of parameters to be calibrated, while maintaining the accuracy and detail required by modern hydrological models, will reduce parameter and model structure uncertainty and improve model diagnostics. In this study, the dynamics of runoff are derived from the distribution of distances from points in the catchments to the nearest stream. This distribution is unique for each catchment and can be determined from a geographical information system. The distribution of distances, will, when a celerity of (subsurface) flow is introduced, provide a distribution of travel times, or a unit hydrograph (UH). For spatially varying levels of saturation deficit, we have different celerities and, hence, different UHs. Runoff is derived from the superposition of the different UHs. This study shows how celerities can be estimated if we assume that recession events represent the combined UHs for different levels of saturation deficit. A new soil moisture routine which estimates saturated and unsaturated volumes of subsurface water and with only one parameter to calibrate is included in the new model. The performance of the new model is compared with that of the Swedish HBV model and is found to perform equally well for eight Norwegian catchments although the number of parameters to be calibrated in the module concerning soil moisture and runoff dynamics is reduced from seven in the HBV model to one in the new model. It is also shown that the new model has a more realistic representation of the subsurface hydrology. Copyright © 2013 John Wiley & Sons, Ltd.     

Influence of long‐period filter cut‐off on elastic spectral displacements

The effect of the long‐period filter cut‐off, T_c, on elastic spectral displacements is investigated using a strong ground‐motion database from Europe and the Middle East. The relation between the filter and oscillator responses is considered to observe the influence of T_c for both analogue and digital records, and the variations with site classification, magnitude, filter order and viscous damping. Robust statistics are derived using the re‐processed European data to generalize the effects of the long‐period filter cut‐off on maximum oscillator deformation demands as a function of these seismological and structural features. Statistics with a 95% confidence interval are derived to suggest usable period ranges for spectral displacement computations as a function of T_c. The results indicate that the maximum period at which spectral displacements can be confidently calculated depend strongly on the site class, magnitude and filter order. The period range where reliable long‐period information can be extracted from digital accelerograms is twice that of analogue records. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

Rainfall‐runoff analysis based on competing linear impulse responses: decomposition of rainfall‐runoff processes

Many recent studies have successfully used neural networks for non‐linear rainfall‐runoff modelling. Due to fundamental limitation of linear structures, approaches employing linear models have been generally considered inferior to the neural network approaches in this area. However, the authors believe that with an appropriate extension, the concept of linear impulse responses can be a viable tool since it enables one to understand underlying dynamics of rainfall‐runoff processes. In this paper, the use of competing impulse responses for rainfall‐runoff analysis is proposed. The proposed method is based on the switch over of competing linear impulse‐responses, each of which satisfies the constraints of non‐negativity and uni‐modality. The computational analyses performed for the rainfall‐runoff data in the Seolma‐Chun experimental basin, Korea showed that the proposed method can yield promising results. Considering the basin characteristics as well as the results from this study, it may be concluded that three impulse responses are enough for rainfall‐runoff analysis. Copyright © 2007 John Wiley & Sons, Ltd.     

An analytical approach to ascertain saturation‐excess versus infiltration‐excess overland flow in urban and reference landscapes

Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration‐excess (where precipitation rate exceeds infiltration capacity) or saturation‐excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one‐dimensional model that distinguishes between infiltration‐excess overland flow (IEOF) and saturation‐excess overland flow (SEOF) using Green–Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low‐ and high‐intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth‐weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high‐intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.     

Effects of conditional parameterization on performance of rainfall‐runoff model regarding hydrologic non‐stationarity

Better parameterization of a hydrological model can lead to improved streamflow prediction. This is particularly important for seasonal streamflow forecasting with the use of hydrological modelling. Considering the possible effects of hydrologic non‐stationarity, this paper examined ten parameterization schemes at 12 catchments located in three different climatic zones in east Australia. These schemes are grouped into four categories according to the period when the data are used for model calibration, i.e. calibration using data: (1) from a fixed period in the historical records; (2) from different lengths of historical records prior to prediction year; (3) from different climatic analogue years in the past; and (4) data from the individual months. Parameterization schemes were evaluated according to model efficiency in both the calibration and verification period. The results show that the calibration skill changes with the different historic periods when data are used at all catchments. Comparison of model performance between the calibration schemes indicates that it is worth calibrating the model with the use of data from each individual month for the purpose of seasonal streamflow forecasting. For the catchments in the winter‐dominant rainfall region of south‐east Australia, a more significant shift in rainfall‐runoff relationships at different periods was found. For those catchments, model calibration with the use of 20 years of data prior to the prediction year leads to a more consistent performance. Copyright © 2011 John Wiley & Sons, Ltd.     

Automatic calibration of HEC‐HMS using single‐objective and multi‐objective PSO algorithms

This study presents single‐objective and multi‐objective particle swarm optimization (PSO) algorithms for automatic calibration of Hydrologic Engineering Center‐ Hydrologic Modeling Systems rainfall‐runoff model of Tamar Sub‐basin of Gorganroud River Basin in north of Iran. Three flood events were used for calibration and one for verification. Four performance criteria (objective functions) were considered in multi‐objective calibration where different combinations of objective functions were examined. For comparison purposes, a fuzzy set‐based approach was used to determine the best compromise solutions from the Pareto fronts obtained by multi‐objective PSO. The candidate parameter sets determined from different single‐objective and multi‐objective calibration scenarios were tested against the fourth event in the verification stage, where the initial abstraction parameters were recalibrated. A step‐by‐step screening procedure was used in this stage while evaluating and comparing the candidate parameter sets, which resulted in a few promising sets that performed well with respect to at least three of four performance criteria. The promising sets were all from the multi‐objective calibration scenarios which revealed the outperformance of the multi‐objective calibration on the single‐objective one. However, the results indicated that an increase of the number of objective functions did not necessarily lead to a better performance as the results of bi‐objective function calibration with a proper combination of objective functions performed as satisfactorily as those of triple‐objective function calibration. This is important because handling multi‐objective optimization with an increased number of objective functions is challenging especially from a computational point of view. Copyright © 2012 John Wiley & Sons, Ltd.     

The effect of estimating areal rainfall using self‐similarity topography method on the simulation accuracy of runoff prediction

Rainfall input for hydrologic modelling is assumed uniformly distributed over the entire catchment. This can lead to significant errors. Investigations of areal rainfall in mountain areas are typically limited by a lack of adequate meteorological and hydrogeological records. This study focuses on areal rainfall in mountain areas within the Kaidu River Basin, China, with the aim of analyzing the influence of areal rainfall on the simulation accuracy of runoff prediction. We conducted a simulation using MIKE 11/NAM rainfall‐runoff model over 92 days of the rain season and compared the simulation error in different methods. On the basis of properties of self‐similarity degree (SSD) in analyzing the detailed characteristics of terrain, areal rainfall was calculated to model the runoff. The results of the model simulations are generally consistent with observed data, indicating that the self‐similarity topography method is able to reflect the spatial change of rainfall. This indicates that the proposed methodology is applicable for the management of water resources in mountain area. The modelling and self‐similarity topography method study allowed quantification of the spatial rainfall and provided an insight into their implications in hydrological forecasting. Copyright © 2011 John Wiley & Sons, Ltd.     

Projecting climate change impacts on stream flow regimes with tracer‐aided runoff models ‐ preliminary assessment of heterogeneity at the mesoscale

The northern mid‐high latitudes form a region that is sensitive to climate change, and many areas already have seen – or are projected to see – marked changes in hydroclimatic drivers on catchment hydrological function. In this paper, we use tracer‐aided conceptual runoff models to investigate such impacts in a mesoscale (749 km²) catchment in northern Scotland. The catchment encompasses both sub‐arctic montane sub‐catchments with high precipitation and significant snow influence and drier, warmer lowland sub‐catchments. We used downscaled HadCM3 General Circulation Model outputs through the UKCP09 stochastic weather generator to project the future climate. This was based on synthetic precipitation and temperature time series generated from three climate change scenarios under low, medium and high greenhouse gas emissions. Within an uncertainty framework, we examined the impact of climate change at the monthly, seasonal and annual scales and projected impacts on flow regimes in upland and lowland sub‐catchments using hydrological models with appropriate process conceptualization for each landscape unit. The results reveal landscape‐specific sensitivity to climate change. In the uplands, higher temperatures result in diminishing snow influence which increases winter flows, with a concomitant decline in spring flows as melt reduces. In the lowlands, increases in air temperatures and re‐distribution of precipitation towards autumn and winter lead to strongly reduced summer flows despite increasing annual precipitation. The integration at the catchment outlet moderates these seasonal extremes expected in the headwaters. This highlights the intimate connection between hydrological dynamics and catchment characteristics which reflect landscape evolution. It also indicates that spatial variability of changes in climatic forcing combined with differential landscape sensitivity in large heterogeneous catchments can lead to higher resilience of the integrated runoff response. Copyright © 2012 John Wiley & Sons, Ltd.     

Hydrological processes of a closed catchment‐lake system in a semi‐arid Mediterranean environment

Data collected in 4 years of field observations were used in conjunction with continuous simulation models to study, at the small‐basin scale, the water balance of a closed catchment‐lake system in a semi‐arid Mediterranean environment. The open water evaporation was computed with the Penman equation, using the data set collected in the middle of the lake. The surface runoff was partly measured at the main tributary and partly simulated using a distributed, catchment, hydrological model, calibrated with the observed discharge. The simplified structure of the developed modelling mainly concerns soil moisture dynamics and bedrock hydraulics, whereas the flow components are physically based. The calibration produced high efficiency coefficients and showed that surface runoff is greatly affected by soil water percolation into fractured bedrock. The bedrock reduces the storm‐flow peaks and the interflow and has important multi‐year effects on the annual runoff coefficients. The net subsurface outflow from the lake was calculated as the residual of the lake water balance. It was almost constant in the dry seasons and increased in the wet seasons, because of the moistening of the unsaturated soil. During the years of observation, rainfall 30% higher than average caused abundant runoff and a continuous rise in the lake water levels. The analysis allows to predict that, in years with lower than the average rainfall, runoff will be drastically reduced and will not be able to compensate for negative balance between precipitation and lake evaporation. Such highly unsteady situations, with great fluctuations in lake levels, are typical of closed catchment‐lake systems in the semi‐arid Mediterranean environment. Copyright © 2012 John Wiley & Sons, Ltd.     

Vegetation controls on small‐scale runoff and erosion dynamics in a degrading dryland environment

This paper investigates the controls of vegetation on runoff and erosion dynamics in the dryland environment of Jornada, New Mexico, USA. As the American southwest has seen significant shifts in the dominant vegetation species in the past 150 years, an understanding of the vegetation effects on hydrological and erosional processes is vital for understanding and managing environmental change. Small‐scale rainfall simulations were carried out to identify the hydrological and erosional processes resulting from the grassland and shrubland vegetation species. Results obtained using tree‐regression analysis suggested that the primary vegetation control on runoff and erosion is the shrub type and canopy density, which directly affects the local microtopographic gradient of mounds beneath the shrubs. Significant interactions and feedbacks were found to occur among the local mound gradient, crust cover, soil aggregate stability and antecedent soil moisture between the different vegetation species for both the runoff and erosion responses. Although some of the shrub species were found to produce higher sediment yields than the grass species, the distinguishing feature of the grassland was the significantly higher enrichment in the fine sediment fraction compared to all other surface cover types. This enrichment in fines has important implications for nutrient movement in such environments. Copyright © 2009 John Wiley & Sons, Ltd.     

Soil detachment and transport on field‐ and laboratory‐scale interrill areas: erosion processes and the size‐selectivity of eroded sediment

Field‐ and laboratory‐scale rainfall simulation experiments were carried out in an investigation of the temporal variability of erosion processes on interrill areas, and the effects of such variation upon sediment size characteristics. Poorly aggregated sandy soils from the semi‐arid environment of Senegal, West Africa, were used on both a 40 m² field plot and a 0·25 m² laboratory plot; rainfall intensity for all experiments was 70 mm h^?1 with a duration of 1 to 2 hours. Time‐series measurements were made of the quantity and the size distribution of eroded material: these permitted an estimate of the changing temporal balance between the main erosion processes (splash and wash). Results from both spatial scales showed a similar temporal pattern of runoff generation and sediment concentration. For both spatial scales, the dominant erosional process was detachment by raindrops; this resulted in a dynamic evolution of the soil surface under raindrop impact, with the rapid formation of a sieving crust followed by an erosion crust. However, a clear difference was observed between the two scales regarding the size of particles detached by both splash and wash. While all measured values were lower than the mean weight diameter (MWD) value of the original soil (mean 0·32 mm), demonstrating the size‐selective nature of wash and splash processes, the MWD values of washed and splashed particles at the field scale ranged from 0·08 to 0·16 mm and from 0·12 to 0·30 mm respectively, whereas the MWD values of washed and splashed particles at the laboratory scale ranged from 0·13 to 0·29 mm and from 0·21 to 0·32 mm respectively. Thus only at the field scale were the soil particles detached by splash notably coarser than those transported by wash. This suggests a transport‐limited erosion process at the field scale. Differences were also observed between the dynamics of the soil loss by wash at the two scales, since results showed wider scatter in the field compared to the laboratory experiments. This scatter is probably related to the change in soil surface characteristics due to the size‐selectivity of the erosion processes at this spatial scale. Copyright © 2006 John Wiley & Sons, Ltd.