Understanding plot‐scale hydrology of Lesser Himalayan watershed—A field study and HYDRUS‐2D modelling approach

Soil moisture dynamics have a significant effect on overland flow generation. Catchment aspect is one of the major controlling factors of overland flow and soil moisture behaviour. A few experimental studies have been carried out in the uneven topography of the Himalayas. This study presents plot‐scale experiments using portable rainfall simulator at an altitude of 1,230 m above mean sea level and modelling of overland flow using observed datasets. Two plots were selected in 2 different aspects of Aglar watershed of Lesser Himalaya; the agro‐forested (AF) plot was positioned at the north aspect whereas the degraded (DE) plot was located at the south aspect of the hillslope. HS flumes and rain gauges were installed to measure the runoff at the outlet of the plot and the rainfall depth during rainfall simulation experiments. Moreover, 10 soil moisture sensors were installed at upslope and downslope locations of both the plots at 5, 15, 25, 35, and 45 cm depth from ground level to capture the soil moisture dynamics. The tests were conducted at intensities of 79.8 and 75 mm/hr in AF plot and 82.2 and 72 mm/hr in the DE plot during Test 1 and Test 2, respectively. The observed data indicate the presence of reinfiltration process only in the AF plot. The high water holding capacity and the presence of reinfiltration process results in less runoff volume in the AF plot compared with the DE plot. The Hortonian overland flow mechanism was found to be the dominant overland flow mechanism as only a few layers of top soil get saturated during all of the rainfall–runoff experiments. The runoff, rainfall, and soil moisture data were subsequently used to calibrate the parameters of HYDRUS‐2D overland flow module to simulate the runoff hydrograph and soil moisture. The components of hydrograph were evaluated in terms of peak discharge, runoff volume and time of concentration, the results were found to be within the satisfactory range. The goodness of fit of simulated hydrographs were more than 0.85 and 0.95 for AF and DE plot, respectively. The model produced satisfactory simulation results of soil moisture for all of the rainfall–runoff experiments. The HYDRUS‐2D overland flow module was found promising to simulate the runoff hydrograph and soil moisture in plot‐scale research.     

A simple model of hillslope response for overland flow generation

This paper deals with the derivation of the hydrological response of a hillslope on the assumption of quick runoff by surface runoff generation. By using the simple non‐linear storage based model, first proposed by Horton, an analytical solution of the overland flow equations over a plane hillslope was derived. This solution establishes a generalization for different flow regimes of Horton's original solution, which is valid for the transitional flow regime only. The solution proposed was compared successfully with that of Horton and, for the turbulent flow regime, to the one derived from kinematic wave theory. This solution can be applied easily to both stationary and non‐stationary rainfall excess events. An analytical solution for the instantaneous response function (IRF) was also derived. Finally, simple expressions to compute peak and time to peak of IRF are proposed. Copyright © 2001 John Wiley & Sons, Ltd.     

Simulating hydrological response for the R‐5 catchment: comparison of two models and the impact of the roads

In this paper the performance of two hydrological‐response models is evaluated and compared based upon simulations for a single rainfall–runoff event. The two models are QPBRRM, a relatively simple model of Horton overland flow, and InHM, a comprehensive physics‐based model of each of the known streamflow generation mechanisms. The rainfall–runoff event focused upon in this study is from the small rangeland catchment in Oklahoma known as R‐5. When calibrated, both QPBRRM and InHM are shown to effectively simulate the R‐5 event. The calibration procedures used in this study for QPBRRM and InHM were quite different. The calibration of QPBRRM was a curve fitting exercise, whereas the calibration of InHM was based upon an internally valid estimate of the continuous head field. In this study QPBRRM did not perform well outside of the calibrated range. The impact of the roads cutting across the R‐5 catchment is simulated with InHM and discussed for the first time in the study reported here. The relative merits of QPBRRM and InHM are each discussed. Copyright © 2002 John Wiley & Sons, Ltd.     

A new surface runoff parameterization with subgrid-scale soil heterogeneity for land surface models

Soil heterogeneity plays an important role in determining surface runoff generation mechanisms. At the spatial scales represented by land surface models used in regional climate model and/or global general circulation models (GCMs) for numerical weather prediction and climate studies, both infiltration excess (Horton) and saturation excess (Dunne) runoff may be present within a studied area or a model grid cell. Proper modeling of surface runoff is essential to a reasonable representation of feedbacks in the land–atmosphere system. In this paper, a new surface runoff parameterization that dynamically represents both Horton and Dunne runoff generation mechanisms within a model grid cell is presented. The new parameterization takes into account of effects of soil heterogeneity on Horton and Dunne runoff. A series of numerical experiments are conducted to study the effects of soil heterogeneity on Horton and Dunne runoff and on soil moisture storage under different soil and precipitation conditions. The new parameterization is implemented into the current version of the hydrologically based variable infiltration capacity (VIC) land surface model and tested over three watersheds in Pennsylvania. Results show that the new parameterization plays a very important role in partitioning the water budget between surface runoff and soil moisture in the atmosphere–land coupling system. Significant underestimation of the surface runoff and overestimation of subsurface runoff and soil moisture could be resulted if the Horton runoff mechanism were not taken into account. Also, the results show that the Horton runoff mechanism should be considered within the context of subgrid-scale spatial variability of soil properties and precipitation. An assumption of time-invariant spatial distribution of potential infiltration rate may result in large errors in surface runoff and soil moisture. In addition, the total surface runoff from the new parameterization is less sensitive to the choice of the soil moisture shape parameter of the distribution.     

Use of the Connectivity of Runoff Model (CRUM) to investigate the influence of storm characteristics on runoff generation and connectivity in semi‐arid areas

Much attention has been given to the surface controls on the generation and transmission of runoff in semi‐arid areas. However, the surface controls form only one part of the system; hence, it is important to consider the effect that the characteristics of the storm event have on the generation of runoff and the transmission of flow across the slope. The impact of storm characteristics has been investigated using the Connectivity of Runoff Model (CRUM). This is a distributed, dynamic hydrology model that considers the hydrological processes relevant to semi‐arid environments at the temporal scale of a single storm event. The key storm characteristics that have been investigated are the storm duration, rainfall intensity, rainfall variability and temporal structure. This has been achieved through the use of a series of defined storm hydrographs and stochastic rainfall. Results show that the temporal fragmentation of high‐intensity rainfall is important for determining the travel distances of overland flow and, hence, the amount of runoff that leaves the slope as discharge. If the high‐intensity rainfall is fragmented, then the runoff infiltrates a short distance downslope. Longer periods of high‐intensity rainfall allow the runoff to travel further and, hence, become discharge. Therefore, storms with similar amounts of high‐intensity rainfall can produce very different amounts of discharge depending on the storm characteristics. The response of the hydrological system to changes in the rainfall characteristics can be explained using a four‐stage model of the runoff generation process. These stages are: (1) all water infiltrating, (2) the surface depression store filling or emptying without runoff occurring, (3) the generation and transmission of runoff and (4) the transmission of runoff without new runoff being generated. The storm event will move the system between the four stages and the nature of the rainfall required to move between the stages is determined by the surface characteristics. This research shows the importance of the variable‐intensity rainfall when modelling semi‐arid runoff generation. The amount of discharge may be greater or less than the amount that would have been produced if constant rainfall intensity is used in the model. Copyright © 2006 John Wiley & Sons, Ltd.     

Measurements of the initiation of post‐wildfire runoff during rainstorms using in situ overland flow detectors

Overland flow detectors (OFDs) were deployed in 2012 on a hillslope burned by the 2010 Fourmile Canyon fire near Boulder, Colorado, USA. These detectors were simple, electrical resistor‐type instruments that output a voltage (0–2·5 V) and were designed to measure and record the time of runoff initiation, a signal proportional to water depth, and the runoff hydrograph during natural convective rainstorms. Initiation of runoff was found to be spatially complex and began at different times in different locations on the hillslope. Runoff started first at upstream detectors 56% of the time, at the mid‐stream detectors 6%, and at the downstream detectors 38% of the time. Initiation of post‐wildfire runoff depended on the time‐to‐ponding, travel time between points, and the time to fill surface depression storage. These times ranged from 0·5–54, 0·4–1·1, and 0·2–14 minutes, respectively, indicating the importance of the ponding process in controlling the initiation of runoff at this site. Time‐to‐ponding was modeled as a function of the rainfall acceleration (i.e. the rate of change of rainfall intensity) and either the cumulative rainfall at the start of runoff or the soil–water deficit. Measurements made by the OFDs provided physical insight into the spatial and temporal initiation of post‐wildfire runoff during unsteady flow in response to time varying natural rainfall. They also provided data that can be telemetered and used to determine critical input parameters for hydrologic rainfall–runoff models. Copyright © 2015 John Wiley & Sons, Ltd.     

Error analysis for the evaluation of model performance: rainfall–runoff event summary variables

This paper provides a procedure for the evaluation of model performance for rainfall–runoff event summary variables, such as total discharge or peak runoff. The procedure is based on the analysis of model errors, defined as the differences between observed values and values predicted by a simulation model. Model errors can (i) indicate whether and where the model can be improved, (ii) be used to measure the performance of a model, and (iii) be used to compare model simulations. In this paper, both statistical and graphical methods are used to characterize model errors. We explore model recalibration by relating model errors to the model predictions, and to external, independent variables. The R‐5 catchment data sets that we used in this study include summary variables for 72 rainfall–runoff events. The simulations used in this study were previously conducted with the quasi‐physically based rainfall–runoff model QPBRRM for 11 different characterizations of the R‐5 catchment, each with increasing information or a refined spatial discretization of the overland flow planes. This paper is about proposing model diagnostics and not about procedures for using diagnostics for model modification. Copyright © 2007 John Wiley & Sons, Ltd.     

Seasonal rainfall–runoff relationships in a lowland forested watershed in the southeastern USA

Hydrological processes of lowland watersheds of the southern USA are not well understood compared to a hilly landscape due to their unique topography, soil compositions, and climate. This study describes the seasonal relationships between rainfall patterns and runoff (sum of storm flow and base flow) using 13 years (1964–1976) of rainfall and stream flow data for a low‐gradient, third‐order forested watershed. It was hypothesized that runoff–rainfall ratios (R/P) are smaller during the dry periods (summer and fall) and greater during the wet periods (winter and spring). We found a large seasonal variability in event R/P potentially due to differences in forest evapotranspiration that affected seasonal soil moisture conditions. Linear regression analysis results revealed a significant relationship between rainfall and runoff for wet (r² = 0·68; p < 0·01) and dry (r² = 0·19; p = 0·02) periods. Rainfall‐runoff relationships based on a 5‐day antecedent precipitation index (API) showed significant (r² = 0·39; p < 0·01) correspondence for wet but not (r² = 0·02; p = 0·56) for dry conditions. The same was true for rainfall‐runoff relationships based on 30‐day API (r² = 0·39; p < 0·01 for wet and r² = 0·00; p = 0·79 for dry). Stepwise regression analyses suggested that runoff was controlled mainly by rainfall amount and initial soil moisture conditions as represented by the initial flow rate of a storm event. Mean event R/P were higher for the wet period (R/P = 0·33), and the wet antecedent soil moisture condition based on 5‐day (R/P = 0·25) and 30‐day (R/P = 0·26) prior API than those for the dry period conditions. This study suggests that soil water status, i.e. antecedent soil moisture and groundwater table level, is important besides the rainfall to seasonal runoff generation in the coastal plain region with shallow soil argillic horizons. Copyright © 2011 John Wiley & Sons, Ltd.     

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

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

Potential of time‐lapse photography for identifying saturation area dynamics on agricultural hillslopes

Mapping saturation areas during rainfall events is important for understanding the dynamics of overland flow. In this study, we evaluate the potential of high temporal resolution time‐lapse photography for mapping the dynamics of saturation areas (i.e., areas where water is visually ponding on the surface) on the hillslope scale during natural rainfall. We take 1 image per minute over a 100 × 15 m² depression area on an agricultural field in the Hydrological Open Air Laboratory, Austria. The images are georectified and classified by an automated procedure, using grey intensity as a threshold to identify saturation area. The optimum threshold T is obtained by comparing saturation areas from the automated analysis with the manual analysis of 149 images. T is found to be highly correlated with an image brightness characteristic defined as the greyscale image histogram mode M (Pearson correlation r = 0.91). We estimate T as T = M + C where C is a calibration parameter assumed to be constant during each event. The automated procedure estimates the total saturation area close to the manual analysis with mean normalized root mean square error of 9% and 21% if C is calibrated for each event and taken constant for all events, respectively. The spatial patterns of saturation are estimated with a geometric mean accuracy index of 94% as compared to the manual analysis of the same photos. The patterns are tested against field observations for one date as a preliminary demonstration, which yields a root mean square error of the shortest distance between the measured boundary points and the automatically classified boundary as 23 cm. The usefulness of the patterns is illustrated by exploring run‐off generation processes of an example event. Overall, the proposed classification method based on grey intensity is found to process images with highly varying brightnesses well. It is more efficient than the manual tracing for a large number of images, which allows the exploration of surface flow processes at high temporal resolution.     

An approach to analysing plot scale infiltration and runoff responses to rainfall of fluctuating intensity

Simulated rainfall of fluctuating intensity was applied to runoff plots on bare dryland soils in order to explore a new method for analysing the non‐steady‐state responses of infiltration and overland flow. The rainfall events all averaged 10 mm/h but included intensity bursts of up to 70 mm/h and lasting 5–15 min, as well as periods of low intensity and intermittency of up to 25 min. Results were compared with traditional steady‐state estimates of infiltrability made under simulated rainfall sustained at a fixed intensity of 10 mm/h. Mean event infiltration rate averaged 13.6% higher under fluctuating intensities, while runoff ratios averaged only 63% of those seen under constant intensity. In order to understand the changing soil infiltrability, up to three affine Horton infiltration equations were fitted to segments of each experiment. All equations had the same final infiltrability f_c, but adjusted values for coefficients f₀ (initial infiltrability) and K_f (exponential decay constant) were fitted for periods of rainfall that followed significant hiatuses in rainfall, during which subsurface redistribution allowed near‐surface soil suction to recover. According to the fitted Horton equations, soil infiltrability recovered by up 10–24 mm/h during intra‐event rainfall hiatuses of 15 to 20‐min duration, contributing to higher overall event infiltration rates and to reduced runoff ratios. The recovery of infiltrability also reduced the size of runoff peaks following periods of low intensity rainfall, compared with the predictions based on single Horton infiltration equations, and in some cases, no runoff at all was recorded from late intensity peaks. The principal finding of this study is that, using a set of affine equations, the intra‐event time variation of soil infiltrability can be tracked through multiple intensity bursts and hiatuses, despite the lack of steady‐state conditions. Copyright © 2016 John Wiley & Sons, Ltd.     

Overland runoff erosion dynamics on steep slopes with forages under field simulated rainfall and inflow

Although numerous studies have acknowledged that vegetation can reduce erosion, few process-based studies have examined how vegetation cover affect runoff hydraulics and erosion processes. We present field observations of overland flow hydraulics using rainfall simulations in a typical semiarid area in China. Field plots (5 × 2 m²) were constructed on a loess hillslope (25°), including bare soil plot as control and three plots with planted forage species as treatments—Astragalus adsurgens, Medicago sativa and Cosmos bipinnatus. Both simulated rainfall and simulated rainfall + inflow were applied. Forages reduced soil loss by 55–85% and decreased overland flow rate by 12–37%. Forages significantly increased flow hydraulic resistance expressed by Darcy–Weisbach friction factor by 188–202% and expressed by Manning's friction factor by 66–75%; and decreased overland flow velocity by 28–30%. The upslope inflow significantly increased overland flow velocity by 67% and stream power by 449%, resulting in increased sediment yield rate by 108%. Erosion rate exhibited a significant linear relationship with stream power. M. sativa exhibited the best in reducing soil loss which probably resulted from its role in reducing stream power. Forages on the downslope performed better at reducing sediment yield than upslope due to decreased rill formation and stream power. The findings contribute to an improved understanding of using vegetation to control water and soil loss and land degradation in semiarid environments.     

Nutrient export via overland flow from a cultivated field of an Ultisol in southern China

Understanding the influence of complex interactions among hydrological factors, soil characteristics and biogeochemical functions on nutrient dynamics in overland flow is important for efficiently managing agricultural nonpoint pollution. Experiments were conducted to assess nutrient export from Ultisol soils in the Sunjia catchment, Jiangxi province, southern China, between 2003 and 2005. Four plots were divided into two groups: two peanut plots and two agroforestry (peanut intercropped with citrus) plots. During the study period, we collected water samples for chemical analyses after each rainfall event that generated overland flow to assess nutrient export dynamics. The concentrations of potassium (K) and nitrate‐N (NO₃^––N) in overland flow were higher during the wetting season (winter and early spring). This reflects the solubility of K and NO₃^––N, the accumulation of NO₃^––N during the dry season and an increase in desorption processes and mixing with pre‐event water caused by prolonged contact with soil in areas with long‐duration, low‐intensity rainfall. In contrast, concentrations of total nitrogen (TN) and total phosphorus (TP) were higher during the wet season (late March to early July) and during the dry season (mid‐July to the end of September or early October). This was due to the interaction between specific hydrological regimes, the properties of the Ultisol and particulate transport processes. Variations in nutrient concentrations during storm events further identified that event water was the dominant source of total nitrogen and total phosphorus, and pre‐event water was the dominant source of NO₃^––N. In addition, the results obtained for the different land uses suggest that agroforestry practices reduce nutrient loss via overland flow. Copyright © 2012 John Wiley & Sons, Ltd.     

A hydrochemical approach to quantify the role of return flow in a surface flow‐dominated catchment

Stormflow generation in headwater catchments dominated by subsurface flow has been studied extensively, yet catchments dominated by surface flow have received less attention. We addressed this by testing whether stormflow chemistry is controlled by either (a) the event‐water signature of overland flow, or (b) the pre‐event water signature of return flow. We used a high‐resolution hydrochemical data set of stormflow and end‐members of multiple storms in an end‐member mixing analysis to determine the number of end‐members needed to explain stormflow, characterize and identify potential end‐members, calculate their contributions to stormflow, and develop a conceptual model of stormflow. The arrangement and relative positioning of end‐members in stormflow mixing space suggest that saturation excess overland flow (26–48%) and return flow from two different subsurface storage pools (17–53%) are both similarly important for stormflow. These results suggest that pipes and fractures are important flow paths to rapidly release stored water and highlight the value of within‐event resolution hydrochemical data to assess the full range and dynamics of flow paths.     

Spatial and temporal variability of the hydrological response in a small Mediterranean research catchment (Vallcebre,Eastern Pyrenees)

This paper analyses the spatial and temporal variability of the hydrological response in a small Mediterranean catchment (Cal Rodó). The first part of the analysis focuses on the rainfall–runoff relationship at seasonal and monthly scale, using an 8‐year data set. Then, using storm‐flow volume and coefficient, the temporal variability of the rainfall–runoff relationship and its relationship with several hydrological variables are analysed at the event scale from hydrographs observed over a 3‐year period. Finally, the spatial non‐linearity of the hydrological response is examined by comparing the Cal Rodó hydrological response with the Can Vila sub‐catchment response at the event scale. Results show that, on a seasonal and monthly scale, there is no simple relationship between rainfall and runoff depths, and that evapotranspiration is a factor that introduced some non‐linearity in the rainfall–runoff relationship. The analysis of monthly values also reveals the existence of a threshold in the relationship between rainfall and runoff depths, denoting a more contrasted hydrological response than the one usually observed in humid catchments. At the event scale, the storm‐flow coefficient has a clear seasonal pattern with an alternance between a wet period, when the catchment is hydrologically responsive, and a dry summer period, when the catchment is much less reactive to any rainfall. The relationship between the storm‐flow coefficient and rainfall depth, rainfall maximum intensity and base‐flow shows that observed correlations are the same as those observed for humid conditions, even if correlation coefficients are notably lower. Comparison with the Can Vila sub‐catchment highlights the spatial heterogeneity of the rainfall‐runoff relationship at the small catchment scale. Although interpretation in terms of runoff processes remains delicate, heterogeneities between the two catchments seem to be related to changes in the ratio between infiltration excess and saturation processes in runoff formation. Copyright © 2007 John Wiley & Sons, Ltd.     

A simple two‐parameter model for scaling hillslope surface runoff

The objective of this research was to develop and parameterise a physically justified yet low‐parameter model to quantify observed changes in surface runoff ratios with hillslope length. The approach starts with the assumption that a unit of rainfall‐excess runoff generated at a point is a fraction β of precipitation P (m) which travels some variable distance down a slope before reinfiltrating, depending on the local rainfall, climate, soils, etc. If this random distance travelled Y is represented by a distribution, then a survival function will describe the probability of this unit of runoff travelling further than some distance x (m). The total amount of per unit width runoff Q (m²) flowing across the lower boundary of a slope of length λ (m) may be considered the sum of all the proportions of the units of rainfall excess runoff integrated from the lower boundary x = 0 to the upper boundary x = λ of the slope. Using these assumptions we derive a model Q(λ) = βPμλ/(μ + λ), (μ > 0, 0 ≤ β ≤ 1, λ ≥ 0) that describes the change in surface runoff with slope length, where μ (m) is the mean of the random variable Y. Dividing both sides of this equation by Pλ yields a simple two‐parameter equation for the dimensionless hillslope runoff ratio Q_h(λ) = βμ/(μ + λ). The model was parameterised with new rainfall and runoff data collected from three replicates of bounded 2 m wide plots of four different lengths (0.5, 1.0, 2.0 and 4.0 m) for 2 years from a forested SE Australian site, and with 32 slope length–runoff data sets from 12 other published studies undertaken between 1934 and 2010. Using the parameterised model resulted in a Nash and Sutcliffe statistic between observed and predicted runoff ratio (for all data sets combined) of 0.93, compared with –2.1 when the runoff ratio was fixed at the value measured from the shortest plot. Copyright © 2014 John Wiley & Sons, Ltd.     

Further testing of the Integrated Hydrology Model (InHM): event‐based simulations for a small rangeland catchment located near Chickasha,Oklahoma

In the paper that is the foundation for this study, VanderKwaak and Loague (2001. Water Resources Research 37 : 999–1013) reported a demonstration of a fully coupled comprehensive physics‐based hydrologic‐response model, InHM (Integrated Hydrology Model), for two rainfall‐runoff events from the small rangeland catchment known as R‐5. The InHM simulations reported herein address (in three phases) limitations in the VanderKwaak and Loague (2001. Water Resources Research 37 : 999–1013) simulations. In Phase I, a new finite‐element mesh was selected to represent R‐5. In Phase II, with the new mesh in place, evaporation was considered for the R‐5 events. In Phase III, with the new mesh in place and evaporation considered, the geology of R‐5 was approximated. Each phase, compared with the results reported by VanderKwaak and Loague (2001. Water Resources Research 37 : 999–1013), shows a change in the simulated near‐surface response. The performance of InHM for 15 R‐5 events is also reported herein. The results from two stages of model calibration are presented. The uncertainty in initial soil‐water content estimates for event‐based simulation is shown to be a major limitation for physics‐based models. The performance of InHM, relative to past event‐based simulation efforts with a quasi‐physically based rainfall‐runoff model, is better for both peak stormflow and the time to peak stormflow, but worse for stormflow depth. The InHM simulations reported here set the stage for continuous simulation of near‐surface response for the R‐5 catchment with InHM. Copyright © 2004 John Wiley & Sons, Ltd.     

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

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

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

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

Mechanisms and pathways of lateral flow on aspen‐forested,Luvisolic soils,Western Boreal Plains,Alberta, Canada

Rainfall simulation experiments by Redding and Devito ( 2008 , Hydrological Processes 23: 4287–4300) on two adjacent plots of contrasting antecedent soil moisture storage on an aspen‐forested hillslope on the Boreal Plain showed that lateral flow generation occurred only once large soil storage capacity was saturated combined with a minimum event precipitation of 15–20 mm. This paper extends the results of Redding and Devito ( 2008 , Hydrological Processes 23: 4287–4300) with detailed analysis of pore pressure, soil moisture and tracer data from the rainfall simulation experiments, which is used to identify lateral flow generation mechanisms and flow pathways. Lateral flow was not generated until soils were wet into the fine textured C horizon. Lateral flow occurred dominantly through the clay‐rich Bt horizon by way of root channels. Lateral flow during the largest event was dominated by event water, and precipitation intensity was critical in lateral flow generation. Lateral flow was initiated as preferential flow near the soil surface into root channels, followed by development of a perched water table at depth, which also interacted with preferential flow pathways to move water laterally by the transmissivity feedback mechanism. The results indicate that lateral flow generated by rainfall on these hillslopes is uncommon because of the generally high available soil moisture storage capacity and the low probability of rainfall events of sufficient magnitude and intensity. Copyright © 2010 John Wiley & Sons, Ltd.