Relationships between rainfall intensity, duration and suspended particle washoff from an urban road surface

A basic understanding of the relationships between rainfall intensity, duration of rainfall and the amount of suspended particles in stormwater runoff generated from road surfaces has been gained mainly from past washoff experiments using rainfall simulators. Simulated rainfall was generally applied at constant intensities, whereas rainfall temporal patterns during actual storms are typically highly variable. This paper discusses a rationale for the application of the constant-intensity washoff concepts to actual storm event runoff. The rationale is tested using suspended particle load data collected at a road site located in Toowoomba, Australia. Agreement between the washoff concepts and measured data is most consistent for intermediate-duration storms (duration <5 h and >1 h). Particle loads resulting from these storm events increase linearly with average rainfall intensity. Above a threshold intensity, there is evidence to suggest a constant or plateau particle load is reached. The inclusion of a peak discharge factor (maximum 6 min rainfall intensity) enhances the ability to predict particle loads.

Indian Ocean sea surface salinity variations in a coupled model

The variability of the sea surface salinity (SSS) in the Indian Ocean is studied using a 100-year control simulation of the Community Climate System Model (CCSM 2.0). The monsoon-driven seasonal SSS pattern in the Indian Ocean, marked by low salinity in the east and high salinity in the west, is captured by the model. The model overestimates runoff int the Bay of Bengal due to higher rainfall over the Himalayan-Tibetan regions which drain into the Bay of Bengal through Ganga-Brahmaputra rivers. The outflow of low-salinity water from the Bay of Bengal is to strong in the model. Consequently, the model Indian Ocean SSS is about 1 less than that seen in the climatology. The seasonal Indian Ocean salt balance obtained from the model is consistent with the analysis from climatological data sets. During summer, the large freshwater input into the Bay of Bengal and its redistribution decide the spatial pattern of salinity tendency. During winter, horizontal advection is the dominant contributor to the tendency term. The interannual variability of the SSS in the Indian Ocean is about five times larger than that in coupled model simulations of the North Atlantic Ocean. Regions of large interannual standard deviations are located near river mouths in the Bay of Bengal and in the eastern equatorial Indian Ocean. Both freshwater input into the ocean and advection of this anomalous flux are responsible for the generation of these anomalies. The model simulates 20 significant Indian Ocean Dipole (IOD) events and during IOD years large salinity anomalies appear in the equatorial Indian Ocean. The anomalies exist as two zonal bands: negative salinity anomalies to the north of the equator and positive to the south. The SSS anomalies for the years in which IOD is not present and for ENSO years are much weaker than during IOD years. Significant interannual SSS anomalies appear in the Indian Ocean only during IOD years.

Atrazine degradation and transport in runoff on a Black Vertosol

In Australia communities are concerned about atrazine being detected in drinking water supplies. It is important to understand mechanisms by which atrazine is transported from paddocks to waterways if we are to reduce movement of agricultural chemicals from the site of application. Two paddocks cropped with grain sorghum on a Black Vertosol were monitored for atrazine, potassium chloride (KCl) extractable atrazine, desethylatrazine (DEA), and desisopropylatrazine (DIA) at 4 soil depths (0-0.05, 0.05-0.10, 0.10-0.20, and 0.20-0.30 m) and in runoff water and runoff sediment. Atrazine + DEA + DIA (total atrazine) had a half-life in soil of 16-20 days, more rapid dissipation than in many earlier reports. Atrazine extracted in dilute potassium chloride, considered available for weed control, was initially 34% of the total and had a half-life of 15-20 days until day 30, after which it dissipated rapidly with a half life of 6 days. We conclude that, in this region, atrazine may not pose a risk for groundwater contamination, as only 0.5% of applied atrazine moved deeper than 0.20 m into the soil, where it dissipated rapidly. In runoff (including suspended sediment) atrazine concentrations were greatest during the first runoff event (57 days after application) (85 μg/L) and declined with time. After 160 days, the total atrazine lost in runoff was 0.4% of the initial application. The total atrazine concentration in runoff was strongly related to the total concentration in soil, as expected. Even after 98% of the KCl-extractable atrazine had dissipated (and no longer provided weed control), runoff concentrations still exceeded the human health guideline value of 40 μg/L. For total atrazine in soil (0-0.05 m), the range for coefficient of soil sorption (Kd) was 1.9-28.4 mL/g and for soil organic carbon sorption (KOC) was 100-2184 mL/g, increasing with time of contact with the soil and rapid dissipation of the more soluble, available phase. Partition coefficients in runoff for total atrazine were initially 3, increasing to 32 and 51 with time, values for DEA being half these. To minimise atrazine losses, cultural practices that maximise rain infiltration, and thereby minimise runoff, and minimise concentrations in the soil surface should be adopted.

Impact of river runoff into the ocean on Indian summer monsoon

Rivers of the world discharge about 36000 km 3 of freshwater into the ocean every year. To investigate the impact of river discharge on climate, we have carried out two 100 year simulations using the Community Climate System Model (CCSM3), one including the river runoff into the ocean and the other excluding it. When the river discharge is shut off, global average sea surface temperature (SST) rises by about 0.5 degrees C and the Indian Summer Monsoon Rainfall (ISMR) increases by about 10% of the seasonal total with large increase in the eastern Bay of Bengal and along the west coast of India. In addition, the frequency of occurrence of La Nina-like cooling events in the equatorial Pacific increases and the correlation between ISMR and Pacific SST anomalies become stronger. The teleconnection between the SST anomalies in the Pacific and monsoon is effected via upper tropospheric meridional temperature gradient and the North African-Asian Jet axis.

Modeling of Hillslope Runoff and Soil Erosion at Rainfall Events Using Cellular Automata Approach

National Science Fund for Distinguished Young Scholars of China [40225004]; National Natural Science Foundation of China [40471048]

An Integrated Approach to Quantifying Groundwater and Surface Water Contributions of Stream Flow

The management of water resources in Ireland prior to the Water Framework Directive (WFD) has focussed on surface water and groundwater as separate entities. A critical element to the successful implementation of the
WFD is to improve our understanding of the interaction between the two and flow mechanisms by which groundwaters discharge to surface waters. An improved understanding of the contribution of groundwater to surface water is required for the classification of groundwater body status and the determination of groundwater quality thresholds. The results of the study will also have a wider application to many areas of the WFD.
A subcommittee of the WFD Groundwater Working Group (GWWG) has been formed to develop a methodology to estimate the groundwater contribution to Irish Rivers. The group has selected a number of analytical techniques to quantify components of stream flow in an Irish context (Master Recession Curve, Unit Hydrograph, Flood Studies Report methodologies and
hydrogeological analytical modelling). The components of stream flow that can be identified include deep groundwater, intermediate and overland. These analyses have been tested on seven pilot catchments that have a variety of hydrogeological settings and have been used to inform and constrain a mathematical model. The mathematical model used was the NAM (NedbØr-AfstrØmnings-Model) rainfall-runoff model which is a module of DHIs MIKE 11 modelling suite. The results from these pilot catchments have been used to develop a decision model based on catchment descriptors from GIS datasets for the selection of NAM parameters. The datasets used include the mapping of aquifers, vulnerability and subsoils, soils, the Digital Terrain Model, CORINE and lakes. The national coverage of the GIS datasets has allowed the extrapolation of the mathematical model to regional catchments across Ireland.

Assessing the combined effect of land management and wildfire on runoff and soil erosion in north central Portugal

Strong and sometimes extreme responses in runoff and soil erosion following wildfires have been reported worldwide. However, in the case of North-Central Portugal, little research had been carried out regarding the hydrologic and erosive impacts of several land management activities in recently burnt areas (such as ground preparation, post-fire logging or post-fire mitigation treatments). This study aims to assess post-fire runoff and soil erosion response on Eucalypt and Maritime pine plantations during the first, second and third years following wildfires. The effect of several pre-fire ground preparation operations (ploughed down-slope, contour ploughed and inclined terraces), post-fire logging activities (on both the eucalypt and pine plantations), as well as the application of hydromulch (a post-fire emergency treatment) on overland flow and soil erosion were compared to burnt but undisturbed and untreated areas. The intensive monitoring of runoff, soil erosion and selected soil properties served to determine the main factors involved in post-fire runoff and soil erosion and their spatial and temporal variation. Soil water repellency deserved special attention, due to its supposed important role for overland flow generation. Repeated rainfall simulation experiments (RSE’s), micro-scale runoff plots and bounded sediment fences were carried out and/or installed immediately after the wildfire on seven burnt slopes. Micro-scale runoff plots results under natural rainfall conditions were also compared to the RSE’s results, which was useful for assessing the representativeness of the data obtained with artificial rainfall. The results showed comparable runoff coefficient (20-60%) but lower sediment losses (125-1000 g m-2) than prior studies in Portugal, but especially outside Portugal. Lower sediment losses were related with the historic intensive land use in the area. In evaluating these losses, however, the shallowness and stoniness of the soils, as well as the high organic matter fraction of the eroded sediments (50%) must not be overlooked. Sediment limited erosion was measured in all the ploughed sites, probably due to the time since ploughing (several years). The disturbance of the soil surface cover due to post-fire logging and wood extraction substantially increased sediment losses at both the pine and eucalypt sites. Hydromulch effectiveness in reducing the runoff (70%) and sediment losses (83%) was attributed to the protective high coverage provided by hydromulch. The hydromulch significantly affected the soil cover and other soil properties and these changes also reduced the soil erosion risk. The rainfall amount was the main factor explaining the variance in runoff. However, a shift from rainfall amount to rainfall intensity was detected when either the surface cover or the infiltration capacity (hydrophilic conditions) increased. Sediment losses were controlled by rainfall intensity and surface cover. The role of soil water repellency on runoff generation was not consistent; the overall repellency levels alone were not enough to assess its hydrological impact. Soil water repellency explained runoff generation in the specific-sites model better than in the overall model. Additionally, soil moisture content was a better predictor for soil water repellency than antecedent rainfall. The natural rainfall results confirmed that RSE’s were able to capture the specific sediment losses and its organic matter content as well as the differences between the ploughed and unploughed sites. Repeated RSE’s also captured the seasonal variations in runoff and sediment losses attributed to soil water repellency. These results have implications for post-fire soil erosion modelling and soil conservation practices in the region, or areas with the same land use, climate and soil characteristics. The measured sediment loss, as well as the increasing frequency of ploughing in recently burnt and unburnt eucalypt stands, suggests ploughing is not an effective as a soil conservation measure. Logging activities with less impact are recommended in order to maintain the forest litter protecting the soil surface. Due to its high effectiveness in reducing runoff and soil erosion, hydromulch is recommended for highly sensitive and vulnerable areas.

Simulating current global river runoff with a global hydrological model: model revisions, validation, and sensitivity analysis

Global hydrological models (GHMs) model the land surface hydrologic dynamics of continental-scale river basins. Here we describe one such GHM, the Macro-scale - Probability-Distributed Moisture model.09 (Mac-PDM.09). The model has undergone a number of revisions since it was last applied in the hydrological literature. This paper serves to provide a detailed description of the latest version of the model. The main revisions include the following: (1) the ability for the model to be run for n repetitions, which provides more robust estimates of extreme hydrological behaviour, (2) the ability of the model to use a gridded field of coefficient of variation (CV) of daily rainfall for the stochastic disaggregation of monthly precipitation to daily precipitation, and (3) the model can now be forced with daily input climate data as well as monthly input climate data. We demonstrate the effects that each of these three revisions has on simulated runoff relative to before the revisions were applied. Importantly, we show that when Mac-PDM.09 is forced with monthly input data, it results in a negative runoff bias relative to when daily forcings are applied, for regions of the globe where the day-to-day variability in relative humidity is high. The runoff bias can be up to - 80% for a small selection of catchments but the absolute magnitude of the bias may be small. As such, we recommend future applications of Mac-PDM.09 that use monthly climate forcings acknowledge the bias as a limitation of the model. The performance of Mac-PDM.09 is evaluated by validating simulated runoff against observed runoff for 50 catchments. We also present a sensitivity analysis that demonstrates that simulated runoff is considerably more sensitive to method of PE calculation than to perturbations in soil moisture and field capacity parameters.

Global runoff routing with the hydrological component of the ECMWF NWP system

A global river routing scheme coupled to the ECMWF land surface model is implemented and tested within the framework of the Global Soil Wetness Project II, to evaluate the feasibility of modelling global river runoff at a daily time scale. The exercise is designed to provide benchmark river runoff predictions needed to verify the land surface model. Ten years of daily runoff produced by the HTESSEL land surface scheme is input into the TRIP2 river routing scheme in order to generate daily river runoff. These are then compared to river runoff observations from the Global Runoff Data Centre (GRDC) in order to evaluate the potential and the limitations. A notable source of inaccuracy is bias between observed and modelled discharges which is not primarily due to the modelling system but instead of to the forcing and quality of observations and seems uncorrelated to the river catchment size. A global sensitivity analysis and Generalised Likelihood Uncertainty Estimation (GLUE) uncertainty analysis are applied to the global routing model. The ground water delay parameter is identified as being the most sensitive calibration parameter. Significant uncertainties are found in results, and those due to parameterisation of the routing model are quantified. The difficulty involved in parameterising global river discharge models is discussed. Detailed river runoff simulations are shown for the river Danube, which match well observed river runoff in upstream river transects. Results show that although there are errors in runoff predictions, model results are encouraging and certainly indicative of useful runoff predictions, particularly for the purpose of verifying the land surface scheme hydrologicly. Potential of this modelling system on future applications such as river runoff forecasting and climate impact studies is highlighted. Copyright © 2009 Royal Meteorological Society.

PERSiST: a flexible rainfall-runoff modelling toolkit for use with the INCA family of models

Runoff generation processes and pathways vary widely between catchments. Credible simulations of solute and pollutant transport in surface waters are dependent on models which facilitate appropriate, catchment-specific representations of perceptual models of the runoff generation process. Here, we present a flexible, semi-distributed landscape-scale rainfall-runoff modelling toolkit suitable for simulating a broad range of user-specified perceptual models of runoff generation and stream flow occurring in different climatic regions and landscape types. PERSiST (the Precipitation, Evapotranspiration and Runoff Simulator for Solute Transport) is designed for simulating present-day hydrology; projecting possible future effects of climate or land use change on runoff and catchment water storage; and generating hydrologic inputs for the Integrated Catchments (INCA) family of models. PERSiST has limited data requirements and is calibrated using observed time series of precipitation, air temperature and runoff at one or more points in a river network. Here, we apply PERSiST to the river Thames in the UK and describe a Monte Carlo tool for model calibration, sensitivity and uncertainty analysis

Calibration of structure in a distributed forecasting model for a semiarid flash flood: dynamic surface storage and channel roughness

Flash floods pose a significant danger for life and property. Unfortunately, in arid and semiarid environment the runoff generation shows a complex non-linear behavior with a strong spatial and temporal non-uniformity. As a result, the predictions made by physically-based simulations in semiarid areas are subject to great uncertainty, and a failure in the predictive behavior of existing models is common. Thus better descriptions of physical processes at the watershed scale need to be incorporated into the hydrological model structures. For example, terrain relief has been systematically considered static in flood modelling at the watershed scale. Here, we show that the integrated effect of small distributed relief variations originated through concurrent hydrological processes within a storm event was significant on the watershed scale hydrograph. We model these observations by introducing dynamic formulations of two relief-related parameters at diverse scales: maximum depression storage, and roughness coefficient in channels. In the final (a posteriori) model structure these parameters are allowed to be both time-constant or time-varying. The case under study is a convective storm in a semiarid Mediterranean watershed with ephemeral channels and high agricultural pressures (the Rambla del Albujón watershed; 556 km 2 ), which showed a complex multi-peak response. First, to obtain quasi-sensible simulations in the (a priori) model with time-constant relief-related parameters, a spatially distributed parameterization was strictly required. Second, a generalized likelihood uncertainty estimation (GLUE) inference applied to the improved model structure, and conditioned to observed nested hydrographs, showed that accounting for dynamic relief-related parameters led to improved simulations. The discussion is finally broadened by considering the use of the calibrated model both to analyze the sensitivity of the watershed to storm motion and to attempt the flood forecasting of a stratiform event with highly different behavior.

Imbalanced land-surface water budgets in a numerical weather prediction system

There has been a significant increase in the skill and resolution of numerical weather prediction models (NWPs) in recent decades, extending the time scales of useful weather predictions. The land-surface models (LSMs) of NWPs are often employed in hydrological applications, which raises the question of how hydrologically representative LSMs really are. In this paper, precipitation (P), evaporation (E) and runoff (R) from the European Centre for Medium-Range Weather Forecasts (ECMWF) global models were evaluated against observational products. The forecasts differ substantially from observed data for key hydrological variables. In addition, imbalanced surface water budgets, mostly caused by data assimilation, were found on both global (P-E) and basin scales (P-E-R), with the latter being more important. Modeled surface fluxes should be used with care in hydrological applications and further improvement in LSMs in terms of process descriptions, resolution and estimation of uncertainties is needed to accurately describe the land-surface water budgets.

Effect of tillage on fractal indices describing soil surface microrelief of a Brazilian Alfisol

Soil surface roughness is known to influence water infiltration, runoff and erosion. Soil surface roughness changes with management and weather and its mathematical description still remains an important issue. The main objective of this study was to investigate the effect of tillage on the two fractal indices, fractal dimension, D, and crossover length, 1, currently used in characterizing soil surface microrelief. The statistical index random roughness, RR, was also assessed. Field experiments were done on an Alfisol located at Rio Grande do Sul State (Brazil). Two tillage treatments (conventional versus direct drilling) were tested. The soil surface microrelief was assessed by point elevation measurements in 16 plots for each treatment. The sampling scheme was a square grid with 20 x 20 mm between point spacing and the plot size was 280 x 280 mm, so that each data set consisted of 225 individual elevation points. All indices were calculated after trend removal, both by slope correction, i.e., oriented microrelief, and by slope plus tillage marks correction, i.e., random microrelief. The implemented algorithm for estimating D and 1 consisted in evaluating the roughness around the local root mean square deviation (RMS) of the point elevation values. Irrespective of tillage treatment and detrending procedure, fractal behavior extended only over a bounded range of scales, from 40 to 100 mm, due to the experimental setup. In these conditions, assessing fractal indices was not always straightforward. The statistical index RR and the fractal index I were significantly different between tillage treatments for oriented and random surface conditions. D values of random soil surfaces were not affected by tillage treatment, whereas D values of oriented microrelief were significantly lower in the direct drilled plots. Removal of tillage marks trend resulted in a significant increase in D values. Within each tillage treatment, 1 and D were significantly correlated. (c) 2006 Elsevier B.V. All rights reserved.

Experimental and numerical analyses about the efficiency of flow through devices for the sediment controll in urban runoff

As land is developed, the impervious surfaces that are created increase the amount of runoff during rainfall events, disrupting the natural hydrologic cycle, with an increment in volume of runoff and in pollutant loadings. Pollutants deposited or derived from an activity on the land surface will likely end up in stormwater runoff in some concentration, such as nutrients, sediment, heavy metals, hydrocarbons, gasoline additives, pathogens, deicers, herbicides and pesticides. Several of these pollutants are particulate-bound, so it appears clear that sediment removal can provide significant water-quality improvements and it appears to be important the knowledge of the ability of stromwater treatment devices to retain particulate matter. For this reason three different units which remove sediments have been tested through laboratory. In particular a roadside gully pot has been tested under steady hydraulic conditions, varying the characteristics of the influent solids (diameter, particle size distribution and specific gravity). The efficiency in terms of particles retained has been evaluated as a function of influent flow rate and particles characteristics; results have been compared to efficiency evaluated applying an overflow rate model. Furthermore the role of particles settling velocity in efficiency determination has been investigated. After the experimental runs on the gully pot, a standard full-scale model of an hydrodynamic separator (HS) has been tested under unsteady influent flow rate condition, and constant solid concentration at the input. The results presented in this study illustrate that particle separation efficiency of the unit is predominately influenced by operating flow rate, which strongly affects the particles and hydraulic residence time of the system. The efficiency data have been compared to results obtained from a modified overflow rate model; moreover the residence time distribution has been experimentally determined through tracer analyses for several steady flow rates. Finally three testing experiments have been performed for two different configurations of a full-scale model of a clarifier (linear and crenulated) under unsteady influent flow rate condition, and constant solid concentration at the input. The results illustrate that particle separation efficiency of the unit is predominately influenced by the configuration of the unit itself. Turbidity measures have been used to compare turbidity with the suspended sediments concentration, in order to find a correlation between these two values, which can allow to have a measure of the sediments concentration simply installing a turbidity probe.

A high-resolution map of direct and indirect connectivity of erosion risk areas to surface waters in Switzerland: A risk assessment tool for planning and policy-making

Off-site effects of soil erosion are becoming increasingly important, particularly the pollution of surface waters. In order to develop environmentally efficient and cost effective mitigation options it is essential to identify areas that bear both a high erosion risk and high connectivity to surface waters. This paper introduces a simple risk assessment tool that allows the delineation of potential critical source areas (CSA) of sediment input into surface waters concerning the agricultural areas of Switzerland. The basis are the erosion risk map with a 2 m resolution (ERM2) and the drainage network, which is extended by drained roads, farm tracks, and slope depressions. The probability of hydrological and sedimentological connectivity is assessed by combining soil erosion risk and extended drainage network with flow distance calculation. A GIS-environment with multiple-flow accumulation algorithms is used for routing runoff generation and flow pathways. The result is a high resolution connectivity map of the agricultural area of Switzerland (888,050 ha). Fifty-five percent of the computed agricultural area is potentially connected with surface waters, 45% is not connected. Surprisingly, the larger part of 34% (62% of the connected area) is indirectly connected with surface waters through drained roads, and only 21% are directly connected. The reason is the topographic complexity and patchiness of the landscape due to a dense road and drainage network. A total of 24% of the connected area and 13% of the computed agricultural area, respectively, are rated with a high connectivity probability. On these CSA an adapted land use is recommended, supported by vegetated buffer strips preventing sediment load. Even areas that are far away from open water bodies can be indirectly connected and need to be included in planning of mitigation measures. Thus, the connectivity map presented is an important decision-making tool for policy-makers and extension services. The map is published on the web and thus available for application.