Runoff generation mechanisms in pastures of the Sand Mountain region of Alabama—a field investigation

Excessive application of poultry litter to pastures in the Sand Mountain region of north Alabama has resulted in phosphorus (P) contamination of surface water bodies and buildup of P in soils of this region. Since surface runoff is recognized as the primary mechanism of P transport, understanding surface runoff generation mechanisms are crucial for alleviating water quality problems in this region. Identification of surface runoff generation mechanisms is also important for delineation of hydrologically active areas (HAAs). Therefore, the specific objective of this study was to identify surface runoff generation mechanisms (infiltration excess versus saturation excess) using distributed surface and subsurface sensors and rain gauge. Results from three rainfall events (2·13–3·43 cm) of differing characteristics, and sensor data at four locations with differing soil hydraulic properties along the hillslope showed that the main surface runoff generation mechanism in this region is infiltration excess. Because of this, rainfall intensity and soil hydraulic conductivity were found to play dominant roles in surface runoff generation in this region. Further, only short periods of a few rainfall events during which the rainfall intensity is high produce surface runoff. This study indicates that perhaps subsurface flows and transport of P in subsurface flows need to be quantified to reduce P contamination of surface water bodies in this region. Current studies at this location are identifying spatial and temporal distribution of HAAs, quantifying rainfall characteristics that generate runoff, and estimating runoff volume that results from connected HAAs. Copyright © 2008 John Wiley & Sons, Ltd.     

Dynamics of stormflow generation—A hillslope‐scale field study in east‐central Pennsylvania,USA

A 40 m × 20 m mowed, grass hillslope adjacent to a headwater stream within a 26‐ha watershed in east‐central Pennsylvania, USA, was instrumented to identify and map the extent and dynamics of surface saturation (areas with the water table at the surface) and surface runoff source areas. Rainfall, stream flow and surface runoff from the hillslope were recorded at 5‐min intervals from 11 August to 22 November 1998, and 13 April to 12 November 1999. The dynamics of the water table (0 to 45 cm depth from the soil surface) and the occurrence of surface runoff source areas across the hillslope were recorded using specially designed subsurface saturation and surface runoff sensors, respectively. Detailed data analyses for two rainfall events that occurred in August (57·7 mm in 150 min) and September (83·6 mm in 1265 min) 1999, illustrated the spatial and temporal dynamics of surface saturation and surface runoff source areas. Temporal data analyses showed the necessity to measure the hillslope dynamics at time intervals comparable to that of rainfall measurements. Both infiltration excess surface runoff (runoff caused when rainfall intensity exceeds soil infiltration capacity) and saturation excess surface runoff (runoff caused when soil moisture storage capacity is exceeded) source areas were recorded during these rainfall events. The August rainfall event was primarily an infiltration excess surface runoff event, whereas the September rainfall event produced both infiltration excess and saturation excess surface runoff. Occurrence and disappearance of infiltration excess surface runoff source areas during the rainfall events appeared scattered across the hillslope. Analysis of surface saturation and surface runoff data showed that not all surface saturation areas produced surface runoff that reached the stream. Emergence of subsurface flow to the surface during the post‐rainfall periods appeared to be a major flow process dominating the hillslope after the August rainfall event. Copyright © 2002 John Wiley & Sons, Ltd.     

Effects of hillslope topography on hydrological responses in a weathered granite mountain,Japan: comparison of the runoff response between the valley‐head and the side slope

To evaluate the effects of hillslope topography on storm runoff in a weathered granite mountain, discharge rate, soil pore water pressures, and water chemistry were observed on two types of hillslope: a valley‐head (a concave hillslope) and a side slope (a planar hillslope). Hydrological responses on the valley‐head and side slope reflected their respective topographic characteristics and varied with the rainfall magnitude. During small rainfall events (<35 mm), runoff from the side slope occurred rapidly relative to the valley‐head. The valley‐head showed little response in storm runoff. As rainfall amounts increased (35–60 mm), the valley‐head yielded a higher flow relative to the side slope. For large rainfall events (>60 mm), runoff from both hillslopes increased with rainfall, although that from the valley‐head was larger than that from the side slope. The differences in the runoff responses were caused by differences in the roles of lower‐slope soils and the convergence of the hillslope. During small rainfall events, the side slope could store little water; in contrast, all rainwater could be stored in the soils at the valley‐head hollow. As the amount of rainfall increased, the subsurface saturated area of the valley‐head extended from the bottom to the upper portion of the slope, with the contributions of transient groundwater via lateral preferential flowpaths due to the high concentration of subsurface water. Conversely, saturated subsurface flow did not contribute to runoff responses, and the subsurface saturated area at the side slope did not extend to the upper slope for the same storm size. During large rainfall events, expansion of the subsurface saturated area was observed in both hillslopes. Thus, differences in the concentration of subsurface water, reflecting hillslope topography, may create differences in the extension of the subsurface saturated area, as well as variability in runoff responses. Copyright © 2007 John Wiley & Sons, Ltd.     

Effects of differential hillslope‐scale water retention characteristics on rainfall–runoff response at the Landscape Evolution Observatory

Hillslopes turn precipitation into runoff and thus exert important controls on various Earth system processes. It remains difficult to collect reliable data necessary for understanding and modeling these Earth system processes in real catchments. To overcome this problem, controlled experiments are being conducted at the Landscape Evolution Observatory at Biosphere 2, The University of Arizona. Previous experiments have revealed differences in hydrological response between 2 landscapes within Landscape Evolution Observatory, even though both landscapes were designed to be identical. In an attempt to discover where the observed differences stem from, we use a fully 3‐dimensional hydrological model (CATchment HYdrology) to show the effect of soil water retention characteristics and saturated hydraulic conductivity on the hydrological response of these 2 hillslopes. We also show that soil water retention characteristics can be derived at hillslope scale from experimental observations of soil moisture and matric potential. It is found that differences in soil packing between the 2 landscapes may be responsible for the observed differences in hydrological response. This modeling study also suggests that soil water retention characteristics and saturated hydraulic conductivity have a profound effect on rainfall–runoff processes at hillslope scale and that parametrization of a single hillslope may be a promising step in modeling rainfall–runoff response in real catchments.     

Upscaling of saturated conductivity for Hortonian runoff modelling

Stochastic and deterministic upscaling techniques are developed that upscale saturated conductivity at the support of 0.04 m² to representative actual infiltration (I_b) for support units (blocks) of 10¹–10⁴ m², as a function of steady state rainfall and runon to the block, under Hortonian runoff (infiltration excess overland flow). Parameters in the upscaling techniques represent the surface runoff flow pattern and the spatial probability distribution of saturated conductivity within the 10¹–10⁴ m² block. The stochastic upscaling technique represents the spatial process of infiltration and runoff using a simple process-imitating model, estimating I_b using Monte Carlo simulation. The deterministic upscaling technique aggregates these processes by a deterministic function relating rainfall and runon to I_b. The stochastic upscaling technique is shown to be capable to upscale saturated conductivity derived from ring infiltrometers to I_b values of plots (1 m²) corresponding to measured I_b values using rainfall simulators. It is shown that both upscaling techniques can be used to estimate I_b for each time step and each block in transient rainfall–runoff models, giving better estimates of cumulative runoff from a hillslope and a small catchment than model runs that do not use upscaling techniques.     

Spatial–temporal variability and hydrologic connectivity of runoff generation areas in a North Alabama pasture—implications for phosphorus transport

This study delineated spatially and temporally variable runoff generation areas in the Sand Mountain region pasture of North Alabama under natural rainfall conditions, and demonstrated that hydrologic connectivity is important for generating hillslope response when infiltration‐excess (IE) runoff mechanism dominates. Data from six rainfall events (13·7–32·3 mm) on an intensively instrumented pasture hillslope (0·12 ha) were analysed. Analysis of data from surface runoff sensors, tipping bucket rain gauge and HS‐flume demonstrated spatial and temporal variability in runoff generation areas. Results showed that the maximum runoff generation area, which contributed to runoff at the outlet of the hillslope, varied between 67 and 100%. Furthermore, because IE was the main runoff generation mechanism on the hillslope, the data showed that as the rainfall intensity changed during a rainfall event, the runoff generation areas expanded or contracted. During rainfall events with high‐intensity short‐ to medium‐duration, 4–8% of total rainfall was converted to runoff at the outlet. Rainfall events with medium‐ to low‐intensity, medium‐duration were found less likely to generate runoff at the outlet. In situ soil hydraulic conductivity (k) was measured across the hillslope, which confirmed its effect on hydrologic connectivity of runoff generation areas. Combined surface runoff sensor and k‐interpolated data clearly showed that during a rainfall event, lower k areas generate runoff first, and then, depending on rainfall intensity, runoff at the outlet is generated by hydrologically connected areas. It was concluded that in IE‐runoff‐dominated areas, rainfall intensity and k can explain hydrologic response. The study demonstrated that only connected areas of low k values generate surface runoff during high‐intensity rainfall events. Identification of these areas would serve as an important foundation for controlling nonpoint source pollutant transport, especially phosphorus. The best management practices can be developed and implemented to reduce transport of phosphorus from these hydrologically connected areas. Copyright © 2009 John Wiley & Sons, Ltd.     

Hillslopes in humid-tropical climates aren't always wet: Implications for hydrologic response and landslide initiation in Puerto Rico

The devastating impacts of the widespread flooding and landsliding in Puerto Rico following the September 2017 landfall of Hurricane Maria highlight the increasingly extreme atmospheric disturbances and enhanced hazard potential in mountainous humid-tropical climate zones. Long-standing conceptual models for hydrologically driven hazards in Puerto Rico posit that hillslope soils remain wet throughout the year, and therefore, that antecedent soil wetness imposes a negligible effect on hazard potential. Our post-Maria in situ hillslope hydrologic observations, however, indicate that while some slopes remain wet throughout the year, others exhibit appreciable seasonal and intra-storm subsurface drainage. Therefore, we evaluated the performance of hydro-meteorological (soil wetness and rainfall) versus intensity-duration (rainfall only) hillslope hydrologic response thresholds that identify the onset of positive pore-water pressure, a predisposing factor for widespread slope instability in this region. Our analyses also consider the role of soil-water storage and infiltration rates on runoff generation, which are relevant factors for flooding hazards. We found that the hydro-meteorological thresholds outperformed intensity-duration thresholds for a seasonally wet, coarse-grained soil, although they did not outperform intensity-duration thresholds for a perennially wet, fine-grained soil. These end-member soils types may also produce radically different stormflow responses, with subsurface flow being more common for the coarse-grained soils underlain by intrusive rocks versus infiltration excess and/or saturation excess for the fine-grained soils underlain by volcaniclastic rocks. We conclude that variability in soil-hydraulic properties, as opposed to climate zone, is the dominant factor that controls runoff generation mechanisms and modulates the relative importance of antecedent soil wetness for our hillslope hydrologic response thresholds.     

Flood forecasting and infiltration modeling/Prévision de crue et modélisation de l’infiltration

Abstract

Abstract The role of accuracy in the representation of infiltration on the effectiveness of real-time flood forecasting models was investigated. A simple semi-distributed model of conceptual type with adaptive estimate of hydraulic characteristics included in the infiltration component was selected. Infiltration was described by a very accurate approach recently formulated for complex rainfall patterns, or alternatively through a simpler formulation known as an extension of the classical time compression approximation. The results indicated that, for situations involving a significant rainfall variability in space, the inaccuracy in the representation of infiltration cannot be corrected by the adaptive component of the rainfall–runoff model. A preliminary analysis of the role of an approximation of saturated hydraulic conductivity to be used in each homogeneous area of the semi-distributed model used both in non-adaptive version and in real-time is also presented.     

A simplified model for estimating field‐scale surface runoff hydrographs

The problem of obtaining field‐scale surface response to rainfall events is complicated by the spatial variability of infiltration characteristics of the soil and rainfall. In this paper, we develop and test a simplified model for generating surface runoff over fields with spatial variation in both rainfall rate and saturated hydraulic conductivities. The model is able to represent the effects of local variation in infiltration, as well as the run‐on effect that controls infiltration of excess water from saturated upstream areas. The effective rainfall excess is routed to the slope outlet using a simplified solution of the kinematic wave approximation. Model results are compared to averaged hydrographs from numerically‐intensive Monte–Carlo simulations for observed and design rainfall events and soil patterns that are typical of Central Italy. The simplified model is found to yield satisfactory results at a relatively small computational expense. A proposal to include a simple channel routing scheme is also presented as a prelude to extend this conceptualization to watershed scales. Copyright © 2006 John Wiley & Sons, Ltd.     

Delineating runoff processes and critical runoff source areas in a pasture hillslope of the Ozark Highlands

The identification of runoff contributing areas would provide the ideal focal points for water quality monitoring and Best Management Practice (BMP) implementation. The objective of this study was to use a field‐scale approach to delineate critical runoff source areas and to determine the runoff mechanisms in a pasture hillslope of the Ozark Highlands in the USA. Three adjacent hillslope plots located at the Savoy Experimental Watershed, north‐west Arkansas, were bermed to isolate runoff. Each plot was equipped with paired subsurface saturation and surface runoff sensors, shallow groundwater wells, H‐flumes and rain gauges to quantify runoff mechanisms and rainfall characteristics at continuous 5‐minute intervals. The spatial extent of runoff source areas was determined by incorporating sensor data into a geographic information‐based system and performing geostatistical computations (inverse distance weighting method). Results indicate that both infiltration excess runoff and saturation excess runoff mechanisms occur to varying extents (0–58% for infiltration excess and 0–26% for saturation excess) across the plots. Rainfall events that occurred 1–5 January 2005 are used to illustrate the spatial and temporal dynamics of the critical runoff source areas. The methodology presented can serve as a framework upon which critical runoff source areas can be identified and managed for water quality protection in other watersheds. Copyright © 2008 John Wiley & Sons, Ltd.     

Capabilities and limitations of detailed hillslope hydrological modelling

Hillslope hydrological modelling is considered to be of great importance for the understanding and quantification of hydrological processes in hilly or mountainous landscapes. In recent years a few comprehensive hydrological models have been developed at the hillslope scale which have resulted in an advanced representation of hillslope hydrological processes (including their interactions), and in some operational applications, such as in runoff and erosion studies at the field scale or lateral flow simulation in environmental and geotechnical engineering. An overview of the objectives of hillslope hydrological modelling is given, followed by a brief introduction of an exemplary comprehensive hillslope model, which stimulates a series of hydrological processes such as interception, evapotranspiration, infiltration into the soil matrix and into macropores, lateral and vertical subsurface soil water flow both in the matrix and preferential flow paths, surface runoff and channel discharge. Several examples of this model are presented and discussed in order to determine the model's capabilities and limitations. Finally, conclusions about the limitations of detailed hillslope modelling are drawn and an outlook on the future prospects of hydrological models on the hillslope scale is given.The model presented performed reasonable calculations of Hortonian surface runoff and subsequent erosion processes, given detailed information of initial soil water content and soil hydraulic conditions. The vertical and lateral soil moisture dynamics were also represented quite well. However, the given examples of model applications show that quite detailed climatic and soil data are required to obtain satisfactory results. The limitations of detailed hillslope hydrological modelling arise from different points: difficulties in the representations of certain processes (e.g. surface crusting, unsaturated–saturated soil moisture flow, macropore flow), problems of small‐scale variability, a general scarcity of detailed soil data, incomplete process parametrization and problems with the interdependent linkage of several hillslopes and channel–hillslope interactions. Copyright © 1999 John Wiley & Sons, Ltd.     

Macroporosity and infiltration in blanket peat: the implications of tension disc infiltrometer measurements

Little is known about the processes of infiltration and water movement in the upper layers of blanket peat. A tension infiltrometer was used to measure hydraulic conductivity in a blanket peat in the North Pennines, England. Measurements were taken from the surface down to 20 cm in depth for peat under four different vegetation covers. It was found that macropore flow is a significant pathway for water in the upper layers of this soil type. It was also found that peat depth and surface vegetation cover were associated with macroporosity and saturated hydraulic conductivity. The proportion of macropore flow was found to be greater at 5 cm depth than at 0, 10 and 20 cm depth. Peat beneath a Sphagnum cover tends to be more permeable and a greater proportion of macropore flow can occur beneath this vegetation type. Functional macroporosity and matrix flow in the near‐surface layers of bare peat appear to have been affected by weathering processes. Comparision of results with rainfall records demonstrates that infiltration‐excess overland flow is unlikely to be a common runoff‐generating mechanism on blanket peat; rather, a saturation‐excess mechanism combined with percolation‐excess above much less permeable layers dominates the runoff response. Copyright © 2001 John Wiley & Sons, Ltd.     

Development and parameterization of an infiltration model accounting for water depth and rainfall intensity

Experimental work has clearly shown that the effective hydraulic conductivity (K_e) or effective infiltration rate (f_e) on the local scale of a plot cannot be considered as constant but are dependent on water depth and rainfall intensity because non‐random microtopography‐related variations in hydraulic conductivity occur. Rainfall–runoff models generally do not account for this: models assume that excess water is uniformly spread over the soil surface and within‐plot variations are neglected. In the present study, we propose a model that is based on the concepts of microtopography‐related water depth‐dependent infiltration and partial contributing area. Expressions for the plot scale K_e and f_e were developed that depend on rainfall intensity and runon from upslope (and thus on water depth). To calibrate and validate the model, steady state infiltration experiments were conducted on maize fields on silt loam soils in Belgium, with different stages and combinations of rainfall intensity and inflow, simulating rainfall and runon. Water depth–discharge and depth–inundation relationships were established and used to estimate the effect of inundation on K_e. Although inflow‐only experiments were found to be unsuitable for calibration, the model was successfully calibrated and validated with the rainfall simulation data and combined rainfall–runon data (R²: 0.43–0.91). Calibrated and validated with steady state infiltration experiments, the model was combined with the Green–Ampt infiltration equation and can be applied within a two‐dimensional distributed rainfall–runoff model. The effect of water depth–dependency and rainfall intensity on infiltration was illustrated for a hillslope. Copyright © 2012 John Wiley & Sons, Ltd.     

The role of run-on for overland flow and the characteristics of runoff generation in the Loess Plateau,China

Abstract

The Loess Plateau in China is overlain by deep and loose soil. As in other semi-arid regions, convective precipitation produces storms, typically of short duration, relatively high intensity and limited areal extent. Infiltration excess (Hortonian mechanism) of precipitation is conventionally assumed to be more prominent than saturation excess (Dunne mechanism) for storm runoff generation. This assumption is true at a point during the storm. However, the runoff generation mechanism is altered when the runoff is conditioned by a lateral redistribution movement of water, i.e. run-on, as the spatial scale increases. In the Loess Plateau, the effects of run-on may be significant, because of the deep and loose surface soil layer. In this study, the role of run-on for overland flow in the Upper Wei River basin, located in the Loess Plateau, is evaluated by means of a simple numerical model at the hillslope scale. The results show that almost all the Hortonian overland flow infiltrates into the soil along the flat hillslope and dry gully before it reaches the river channel. Most of the runoff is generated from the saturated soil near the river channel and from the subsurface. The run-on process takes much longer than the infiltration, facilitating rainfall–runoff modelling at a daily time step. A hydrological model is employed to investigate the characteristics of runoff generation in the Upper Wei River basin. The analysis shows that the subsurface flow contribution to total streamflow is more than 53% from October to March, while the overland flow contribution exceeds 72% from April to September.

Editor D. Koutsoyiannis; Associate editor Dawen Yang

Citation Liu, D.F., Tian, F.Q., Hu, H.C., and Hu, H.P., 2012. The role of run-on for overland flow and the characteristics of runoff generation in the Loess Plateau, China. Hydrological Sciences Journal, 57 (6), 1107–1117.     

Rainfall–runoff responses on Arctic hillslopes underlain by continuous permafrost,North Slope,Alaska, USA

The Arctic hydrologic cycle is intensifying, as evidenced by increased rates of precipitation, evapotranspiration, and riverine discharge. However, the controls on water fluxes from terrestrial to aquatic systems in upland Arctic landscapes are poorly understood. Upland landscapes account for one third of the Arctic land surface and are often drained by zero‐order geomorphic flowpath features called water tracks. Previous work in the region attributed rapid runoff response at larger stream orders to water tracks, but models suggest water tracks are hydrologically disconnected from the surrounding hillslope. To better understand the role of water tracks in upland landscapes, we investigated the surface and subsurface hydrologic responses of 6 water tracks and their hillslope watersheds to natural patterns of rainfall, soil thaw, and drainage. Between storms, both water track discharge and the water table in the hillslope watersheds exhibited diel fluctuations that, when lagged by 5 hr, were temporally correlated with peak evapotranspiration rate. Water track soils remained saturated for more of the summer season than soils in their surrounding hillslope watersheds. When rainfall occurred, the subsurface response was nearly instantaneous, but the water tracks took significantly longer than the hillslopes to respond to rainfall, and longer than the responses previously observed in nearby larger order Arctic streams. There was also evidence for antecedent soil water storage conditions controlling the magnitude of runoff response. Based on these observations, we used a broken stick model to test the hypothesis that runoff production in response to individual storms was primarily controlled by rainfall amount and antecedent water storage conditions near the water track outlet. We found that the relative importance of the two factors varied by site, and that water tracks with similar watershed geometries and at similar landscape positions had similar rainfall–runoff model relationships. Thus, the response of terrestrial water fluxes in the upland Arctic to climate change depends on the non‐linear interactions between rainfall patterns and subsurface water storage capacity on hillslopes. Predicting these interactions across the landscape remains an important challenge.     

Difference infiltrometer: a method to measure temporally variable infiltration rates during rainstorms

We developed a difference infiltrometer to measure time series of non‐steady infiltration rates during rainstorms at the point scale. The infiltrometer uses two, tipping bucket rain gages. One gage measures rainfall onto, and the other measures runoff from, a small circular plot about 0.5‐m in diameter. The small size allows the infiltration rate to be computed as the difference of the cumulative rainfall and cumulative runoff without having to route water through a large plot. Difference infiltrometers were deployed in an area burned by the 2010 Fourmile Canyon Fire near Boulder, Colorado, USA, and data were collected during the summer of 2011. The difference infiltrometer demonstrated the capability to capture different magnitudes of infiltration rates and temporal variability associated with convective (high intensity, short duration) and cyclonic (low intensity, long duration) rainstorms. Data from the difference infiltrometer were used to estimate saturated hydraulic conductivity of soil affected by the heat from a wildfire. The difference infiltrometer is portable and can be deployed in rugged, steep terrain and does not require the transport of water, as many rainfall simulators require, because it uses natural rainfall. It can be used to assess infiltration models, determine runoff coefficients, identify rainfall depth or rainfall intensity thresholds to initiate runoff, estimate parameters for infiltration models, and compare remediation treatments on disturbed landscapes. The difference infiltrometer can be linked with other types of soil monitoring equipment in long‐term studies for detecting temporal and spatial variability at multiple time scales and in nested designs where it can be linked to hillslope and basin‐scale runoff responses. Published 2012. This article is a U.S. Government work and is in the public domain in the USA.     

An analysis of the impact of spatial variability in rainfall on runoff and sediment predictions from a distributed model

This paper investigates the effect of introducing spatially varying rainfall fields to a hydrological model simulating runoff and erosion. Pairs of model simulations were run using either spatially uniform (i.e. spatially averaged) or spatially varying rainfall fields on a 500‐m grid. The hydrological model used was a simplified version of Thales which enabled runoff generation processes to be isolated from hillslope averaging processes. Both saturation excess and infiltration excess generation mechanisms were considered, as simplifications of actual hillslope processes. A 5‐year average recurrence interval synthetic rainfall event typical of temperate climates (Melbourne, Australia) was used. The erosion model was based on the WEPP interrill equation, modified to allow nonlinear terms relating the erosion rate to rainfall or runoff‐squared. The model results were extracted at different scales to investigate whether the effects of spatially varying rainfall were scale dependent. A series of statistical metrics were developed to assess the variability due to introducing the spatially varying rainfall field. At the catchment (approximately 150 km²) scale, it was found that particularly for saturation excess runoff, model predictions of runoff were insensitive to the spatial resolution of the rainfall data. Generally, erosion processes at smaller sub‐catchment scales, particularly when the sediment generation equation had non linearity, were more sensitive to spatial rainfall variability. Introducing runon infiltration reduced the total runoff and sediment yield at all scales, and this process was also most sensitive to the rainfall resolution. Copyright © 2011 John Wiley & Sons, Ltd.     

Infiltration and runoff generation processes in fire‐affected soils

Post‐wildfire runoff was investigated by combining field measurements and modelling of infiltration into fire‐affected soils to predict time‐to‐start of runoff and peak runoff rate at the plot scale (1 m²). Time series of soil‐water content, rainfall and runoff were measured on a hillslope burned by the 2010 Fourmile Canyon Fire west of Boulder, Colorado during cyclonic and convective rainstorms in the spring and summer of 2011. Some of the field measurements and measured soil physical properties were used to calibrate a one‐dimensional post‐wildfire numerical model, which was then used as a ‘virtual instrument’ to provide estimates of the saturated hydraulic conductivity and high‐resolution (1 mm) estimates of the soil‐water profile and water fluxes within the unsaturated zone. Field and model estimates of the wetting‐front depth indicated that post‐wildfire infiltration was on average confined to shallow depths less than 30 mm. Model estimates of the effective saturated hydraulic conductivity, K_s, near the soil surface ranged from 0.1 to 5.2 mm h^?1. Because of the relatively small values of K_s, the time‐to‐start of runoff (measured from the start of rainfall), t_p, was found to depend only on the initial soil‐water saturation deficit (predicted by the model) and a measured characteristic of the rainfall profile (referred to as the average rainfall acceleration, equal to the initial rate of change in rainfall intensity). An analytical model was developed from the combined results and explained 92–97% of the variance of t_p, and the numerical infiltration model explained 74–91% of the variance of the peak runoff rates. These results are from one burned site, but they strongly suggest that t_p in fire‐affected soils (which often have low values of K_s) is probably controlled more by the storm profile and the initial soil‐water saturation deficit than by soil hydraulic properties. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

Assessing the impact of the hydraulic properties of a crusted soil on overland flow modelling at the field scale

Soil surface crusts are widely reported to favour Hortonian runoff, but are not explicitly represented in most rainfall‐runoff models. The aim of this paper is to assess the impact of soil surface crusts on infiltration and runoff modelling at two spatial scales, i.e. the local scale and the plot scale. At the local scale, two separate single ring infiltration experiments are undertaken. The first is performed on the undisturbed soil, whereas the second is done after removal of the soil surface crust. The HYDRUS 2D two‐dimensional vertical infiltration model is then used in an inverse modelling approach, first to estimate the soil hydraulic properties of the crust and the subsoil, and then the effective hydraulic properties of the soil represented as a single uniform layer. The results show that the crust hydraulic conductivity is 10 times lower than that of the subsoil, thus illustrating the limiting role the crust has on infiltration. Moving up to the plot scale, a rainfall‐runoff model coupling the Richards equation to a transfer function is used to simulate Hortonian overland flow hydrographs. The previously calculated hydraulic properties are used, and a comparison is undertaken between a single‐layer and a double‐layer representation of the crusted soil. The results of the rainfall‐runoff model show that the soil hydraulic properties calculated at the local scale give acceptable results when used to model runoff at the plot scale directly, without any numerical calibration. Also, at the plot scale, no clear improvement of the results can be seen when using a double‐layer representation of the soil in comparison with a single homogeneous layer. This is due to the hydrological characteristics of Hortonian runoff, which is triggered by a rainfall intensity exceeding the saturated hydraulic conductivity of the soil surface. Consequently, the rainfall‐runoff model is more sensitive to rainfall than to the subsoil's hydrodynamic properties. Therefore, the use of a double‐layer soil model to represent runoff on a crusted soil does not seem necessary, as the increase of precision in the soil discretization is not justified by a better performance of the model. Copyright © 2005 John Wiley & Sons, Ltd.     

Near-surface water fluxes and their controls in a sloping heterogeneously layered volcanic soil beneath a supra-wet tropical montane cloud forest (NW Costa Rica)

Tropical montane cloud forests (TMCF) receive additional (‘occult’) inputs of water from fog and wind-driven rain. Together with the concomitant reduction in evaporative losses, this typically leads to high soil moisture levels (often approaching saturation) that are likely to promote rapid subsurface flow via macropores. Although TMCF make up an estimated 6.6% of all remaining montane tropical forest and occur mostly in steep headwater areas that are protected in the expectation of reduced downstream flooding, TMCF hillslope hydrological functioning has rarely been studied. To better understand the hydrological response of a supra-wet TMCF (net precipitation up to 6535 mm y⁻¹) on heterogeneously layered volcanic ash soils (Andosols), we examined temporal and spatial soil moisture dynamics and their contribution to shallow subsurface runoff and stormflow for a year (1 July 2003–30 June 2004) in a small headwater catchment on the Atlantic (windward) slope near Monteverde, NW Costa Rica. Particular attention was paid to the partitioning of water fluxes into lateral subsurface flow and vertical percolation. The presence of a gravelly layer (C-horizon) at ~25 cm depth of very high hydraulic conductivity (geometric mean: 502 mm h⁻¹) intercalated between two layers of much lower conductivity (7.5 and 15.7 mm h⁻¹ above and below, respectively), controlled both surface infiltration and delayed vertical water movement deeper into the soil profile. Soil water fluxes during rainfall were dominated by rapid lateral flow in the gravelly layer, particularly at high soil moisture levels. In turn, this lateral subsurface flow controlled the magnitude and timing of stormflow from the catchment. Stormflow amount increased rapidly once topsoil moisture content exceeded a threshold value of ~0.58 cm³ cm⁻³. Responses were not affected appreciably by rainfall intensity because soil hydraulic conductivities across the profile largely exceeded prevailing rainfall intensities.