A numerical investigation of bedrock groundwater recharge and exfiltration on soil mantled hillslopes

Here we use Richards Equation models of variably saturated soil and bedrock groundwater flow to investigate first-order patterns of the coupling between soil and bedrock flow systems. We utilize a Monte Carlo sensitivity analysis to identify important hillslope parameters controlling bedrock recharge and then model the transient response of bedrock and soil flow to seasonal precipitation. Our results suggest that hillslopes can be divided into three conceptual zones of groundwater interaction, (a) the zone of lateral unsaturated soil moisture accumulation (upper portion of hillslope), (b) the zone of soil saturation and bedrock recharge (middle of hillslope) and (c) the zone of saturated-soil lateral flow and bedrock groundwater exfiltration (bottom of hillslope). Zones of groundwater interaction expand upslope during periods of precipitation and drain downslope during dry periods. The amount of water partitioned to the bedrock groundwater system a can be predicted by the ratio of bedrock to soil saturated hydraulic conductivity across a variety of hillslope configurations. Our modelled processes are qualitatively consistent with observations of shallow subsurface saturation and groundwater fluctuation on hillslopes studied in our two experimental watersheds and support a conceptual model of tightly coupled shallow and deep subsurface circulation where groundwater recharge and discharge continuously stores and releases water from longer residence time storage.     

Snowmelt‐driven macropore flow and soil saturation in a semiarid forest

Lateral subsurface flow is generally assumed to occur as a result of the development of a saturated zone above a low‐permeability interface such as at the soil–bedrock contact, and it is often augmented by macropore flow. Our objective was to evaluate the development of lateral subsurface flow and soil saturation at a semiarid ponderosa pine forest in New Mexico with respect to the conceptual model of saturation building above the soil–bedrock contact. At this site, we have long‐term observations of the water budget components, including lateral flow. A 1·5 m deep by 7 m long trench was constructed to observe lateral subsurface flow and development of saturation directly. Our observations are based on flow resulting from a melting snowdrift. The edge of the drift was about 7 m upslope from the trench. Lateral subsurface flow only occurred from root macropores in the Bt soil horizon. Saturation developed and grew outward from flowing root macropores, rather than growing upward from the soil–bedrock interface. This macropore‐centred saturation resulted in a highly heterogeneous distribution of water content until enough macropores began flowing and individual macropore saturated zones grew large enough to coalesce and saturate large volumes of the soil. Our observations are based on one snowmelt event and a relatively short hillslope flow path, and thus do not represent a full range of hydrologic conditions. Nevertheless, the observed behaviour did not conform to the traditional model of soil–bedrock control of saturation and lateral flow. Copyright © 2004 John Wiley & Sons, Ltd.     

PRINCIPES NOUVEAUX POUR DISTINGUER CE QUI EST Dû A L'ÉCOULEMENT SOUTERRAIN SUR L'HYDROGRAPHE DES FLEUVES

Abstract

Throughflow has been measured from three soil horizons on a 12 slope with impermeable, bedrock. Storm flow comes from the 10–45 cm horizon and is controlled by the upslope extent of saturated conditions. Base flow comes from the 45–75 cm horizon and is supplied by slow unsaturated flow from the whole soil mass to a small constant zone of saturation.

Differences between input and output stream hydrographs over 270 metres of channel are attributed to throughflow and correlate well with measured values providing a basis for separating throughflow components from the stream hydrograph. Observed stream flows contain no overland flow or ground water flow components. The main basin flood peak is not generated within this control section of channel but is produced in the headwater zone (0.1 km²) by the faster runoff characteristics of the soils in that area and by topographic factors which lead to rapid channel extension.     

Mechanisms responsible for streamflow generation on a small,salt-affected and deeply weathered hillslope

This paper considers the contributions of overland flow, throughflow and deep seepage to the generation of streamflow in a salt-affected, deeply weathered landscape. Runoff mechanisms on a small hillslope in south-western Australia were dependent on the extent and development of variable source areas. In winter, streamflow generation was controlled by returnflow, saturation overland flow and throughflow. In summer, post-ponding, infiltration-excess and saturation overland flow dominated. The extent of the variable source area and the magnitude of streamflow were due to antecedent soil moisture, rainfall and slope morphology. Concave hillslope sections accumulated soil moisture due to both saturated and unsaturated lateral flow processes. Throughflow provided the mechanism and vehicle for solute movement from the groundwater discharge area to the stream. However, discharge from the deep aquifer was the primary mechanism responsible for soil salinity and maintaining the core of the variable source area. Estimates of throughflow which only take account of soil-water movement and disregard returnflow, will underestimate the magnitude of throughflow.     

A low‐dimensional physically based model of hydrologic control of shallow landsliding on complex hillslopes

Hillslopes have complex three‐dimensional shapes that are characterized by their plan shape, profile curvature of surface and bedrock, and soil depth. To investigate the stability of complex hillslopes (with different slope curvatures and plan shapes), we combine the hillslope‐storage Boussinesq (HSB) model with the infinite slope stability method. The HSB model is based on the continuity and Darcy equations expressed in terms of storage along the hillslope. Solutions of the HSB equation account explicitly for plan shape by introducing the hillslope width function and for profile curvature through the bedrock slope angle and the hillslope soil depth function. The presented model is composed of three parts: a topography model conceptualizing three‐dimensional soil mantled landscapes, a dynamic hydrology model for shallow subsurface flow and water table depth (HSB model) and an infinite slope stability method based on the Mohr–Coulomb failure law. The resulting hillslope‐storage Boussinesq stability model (HSB‐SM) is able to simulate rain‐induced shallow landsliding on hillslopes with non‐constant bedrock slope and non‐parallel plan shape. We apply the model to nine characteristic hillslope types with three different profile curvatures (concave, straight, convex) and three different plan shapes (convergent, parallel, divergent). In the presented model, the unsaturated storage has been calculated based on the unit head gradient assumption. To relax this assumption and to investigate the effect of neglecting the variations of unsaturated storage on the assessment of slope stability in the transient case, we also combine a coupled model of saturated and unsaturated storage and the infinite slope stability method. The results show that the variations of the unsaturated zone storage do not play a critical role in hillslope stability. Therefore, it can be concluded that the presented dynamic slope stability model (HSB‐SM) can be used safely for slope stability analysis on complex hillslopes. Our results show that after a certain period of rainfall the convergent hillslopes with concave and straight profiles become unstable more quickly than others, whilst divergent convex hillslopes remain stable (even after intense rainfall). In addition, the relation between subsurface flow and hillslope stability has been investigated. Our analyses show that the minimum safety factor (FS) occurs when the rate of subsurface flow is a maximum. In fact, by increasing the subsurface flow, stability decreases for all hillslope shapes. Copyright © 2008 John Wiley & Sons, Ltd.     

Hillslope hydrographs by the finite element method

Simulation models may be used to explore the implications of making specific assumptions about the nature of a real world system, and then to make predictions of the behaviour of that system under a set of naturally occurring conditions. It is important that understanding generated by the former should be gained before predictive use of the system model. This paper describes and uses a finite-element model of transient, partially saturated water flow within a hillslope soil mantle overlying an impermeable bedrock, to make an investigation into the effects of parameter variations and initial conditions on the hillslope hydrograph. The results clearly demonstrate that the response of the hillslope system to rainfall is highly non-linear and that the initial conditions, particularly in the unsaturated zone, are of paramount importance in governing the timing and magnitude of the hydrograph peak. Hillslope convergence appears as the dominant topographic parameter but the non-linearity of the response and the complex interdependence between the soil and topographic parameters restrict the possibility of further definite conclusions about the relative sensitivity of the simulated hillslope hydrograph to changes in these parameters.     

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.     

A low‐dimensional model of bedrock weathering and lateral flow coevolution in hillslopes: 2. Controls on weathering and permeability profiles,drainage hydraulics,and solute export pathways

The advance of a chemical weathering front into the bedrock of a hillslope is often limited by the rate weathering products that can be carried away, maintaining chemical disequilibrium. If the weathering front is within the saturated zone, groundwater flow downslope may affect the rate of transport and weathering—however, weathering also modifies the rock permeability and the subsurface potential gradient that drives lateral groundwater flow. This feedback may help explain why there tends to be neither “runaway weathering” to great depth nor exposed bedrock covering much of the earth and may provide a mechanism for weathering front advance to keep pace with incision of adjacent streams into bedrock. This is the second of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory. Here, we show how a simplified kinetic model of 1‐D rock weathering can be extended to consider lateral flow in a 2‐D hillslope. Exact and approximate analytical solutions for the location and thickness of weathering within the hillslope are obtained for a number of cases. A location for the weathering front can be found such that lateral flow is able to export weathering products at the rate required to keep pace with stream incision at steady state. Three pathways of solute export are identified: “diffusing up,” where solutes diffuse up and away from the weathering front into the laterally flowing aquifer; “draining down,” where solutes are advected primarily downward into the unweathered bedrock; and “draining along,” where solutes travel laterally within the weathering zone. For each pathway, a different subsurface topography and overall relief of unweathered bedrock within the hillslope is needed to remove solutes at steady state. The relief each pathway requires depends on the rate of stream incision raised to a different power, such that at a given incision rate, one pathway requires minimal relief and, therefore, likely determines the steady‐state hillslope profile.     

Physical controls on septic leachate movement in the vadose zone at the hillslope scale,Putnam County,New York,USA

The fate and transport of contaminants in the vicinity of septic fields remains poorly understood in many hydrogeomorphological environments. We report hydrometric data from an intensive hillslope‐scale experiment conducted between 29 August and 11 November 1998 at a residential leach field in New York State. The objective of our study was to characterize water flux within the vadose zone, understand the physical controls on the flux, and predict how this ultimately will affect subsurface water quality. Soil‐water flux was calculated using matric potential measurements from a network of 25 tensiometer nests, each nest consisting of three tensiometers installed to depths of 10, 50 and 130 cm. Unsaturated hydraulic conductivity curves were derived at each depth from field‐determined time‐domain reflectometry–tensiometry moisture‐release curves and borehole permeametry measurements. Flownets indicated that a strong upward flux of soil water occurred between rainstorms. Following the onset of (typically convective) rainfall, low near‐surface matric potentials were rapidly converted to near‐saturated and saturated conditions, promoting steep vertical gradients through the near‐surface horizons of the hillslope. Lateral hydraulic gradients were typically 10 times smaller than the vertical gradients. Resultant flow vectors showed that the flux was predominantly vertical through the vadose zone, and that the flux response to precipitation was short‐lived. The flux response was controlled primarily by the shape of the unsaturated hydraulic conductivity curves, which indicated a rapid loss of conductivity below saturation. Thus, soil water had a very high residence time in the vadose zone. The absence of rapid wetting at 130 cm and the delayed and small phreatic zone response to rainfall indicated that water movement through macropores did not occur on this hillslope. These results are consistent with a Cl tracing experiment, which demonstrated that the tracer was retained in the vadose zone for several months after injection to the system. Copyright © 2002 John Wiley & Sons, Ltd.     

The distribution of solute processes on an acid hillslope and the delivery of solutes to a stream: II. Exchangeable Al3+

In order to identify the distribution of aluminium (Al) within an acid hillslope and its release to a stream, the spatial distribution of acid ammonium oxalate extractable Al (Al_o) and exchangeable Al³⁺ have been investigated on a podzolized hillslope in Bicknoller Combe, Somerset, UK. The eluviated Al from topsoils is mainly deposited in the lower soil horizons forming podzolic B horizons, but some Al flows downslope carried by lateral throughflow. Al oxides may provide the main source of exchangeable Al³⁺ on the study slope due to high soil acidity. Examination of the spatial distribution of exchangeable Al³⁺ suggests that the slope hollow, where active convergent throughflow occurs, and the saturation wedge at the base of the slope are the main delivery routes of dissolved Al³⁺ to the stream. Divalent base cations (Ca²⁺ and Mg²⁺), supplied from atmospheric input and organic decomposition and carried by throughflow, exchange Al³⁺ via cation exchange reactions under high water content. Laterally illuviated Al oxides in the lower hollow adjacent to the saturation wedge probably provide a pool for continuous delivery of Al either as soluble or complexed forms to the stream via the saturated wedge. Copyright © 1999 John Wiley & Sons, Ltd.     

Lateral water flux in the unsaturated zone: A mechanism for the formation of spatial soil heterogeneity in a headwater catchment

Measurements of soil water potential and water table fluctuations suggest that morphologically distinct soils in a headwater catchment at the Hubbard Brook Experimental Forest in New Hampshire formed as a result of variations in saturated and unsaturated hydrologic fluxes in the mineral soil. Previous work showed that each group of these soils had distinct water table fluctuations in response to precipitation; however, observed variations in soil morphology also occurred above the maximum height of observed saturation. Variations in unsaturated fluxes have been hypothesized to explain differences in soil horizon thickness and presence/absence of specific horizons but have not been explicitly investigated. We examined tensiometer and shallow groundwater well records to identify differences in unsaturated water fluxes among podzols that show distinct morphological and chemical differences. The lack of vertical hydraulic gradients at the study sites suggests that lateral unsaturated flow occurs in several of the soil units. We propose that the variations in soil horizon thickness and presence/absence observed at the site are due in part to slope‐parallel water flux in the unsaturated portion of the solum. In addition, unsaturated flow may be involved in the translocation of spodic material that primes those areas to contribute water with distinct chemistry to the stream network and represents a potential source/sink of organometallic compounds in the landscape.     

Seasonal soil water dynamics in the jarrah forest,Western Australia. II: Results from a site with fine-textured soil profiles

Seasonal soil water dynamics were measured at a fine-textured, upslope site within the jarrah forest of southwest Western Australia and compared to the results from a coarse-textured hillslope transect. Gravity drainage dominated during winter and early spring. This reversed in early summer and an upward potential gradient was observed to 7 m depth. A shallow ephemeral saturation zone was observed above a clay pan at 1.5 m depth. This saturation zone persisted through late winter and early spring, contrasting with the short-lived saturation in the duricrust on the hillslope transect. The annual maximum to minimum unsaturated soil water storage was about 530 mm, 50 mm greater than the hillslope transect and higher than most values reported elsewhere in Australia. Significant soil water content changes following winter rain were generally restricted to 6 m but at one site occurred to 9 m. These depths were significantly less than the coarser-textured hillslope transect. Soil water drying rates averaged 5 mm day^?1 during extended dry periods compared to 3.5 mm day^?1 on the hillslope transect. The drying rate occurred uniformly through the profile until late summer when a significant decrease in the upper 3 m was observed.     

Parametric study of a physically-based,plot-scale hillslope hydrological model through virtual experiments

Abstract

A physically-based hillslope hydrological model with shallow overland flow and rapid subsurface stormflow components was developed and calibrated using field experiments conducted on a preferential path nested hillslope in northeast India. Virtual experiments were carried out to perform sensitivity analysis of the model using the automated parameter estimation (PEST) algorithm. Different physical parameters of the model were varied to study the resulting effects on overland flow and subsurface stormflow responses from the theoretical hillslopes. It was observed that topographical shapes had significant effects on overland flow hydrographs. The slope profiles, surface storage, relief, rainfall intensity and infiltration rates primarily controlled the overland flow response of the hillslopes. Prompt subsurface stormflow responses were mainly dominated by lateral preferential flow, as soil matrix flow rates were very slow. Rainfall intensity and soil macropore structures were the most influential parameters on subsurface stormflow. The number of connected soil macropores was a more sensitive parameter than the size of macropores. In hillslopes with highly active vertical and lateral preferential pathways, saturation excess overland flow was not evident. However, saturation excess overland flow was generated if the lateral macropores were disconnected. Under such conditions, rainfall intensity, duration and preferential flow rate governed the process of saturation excess overland flow generation from hillslopes.

Editor D. Koutsoyiannis; Associate editor C. Perrin     

Characterization of subsurface fluxes at the plot scale during flash floods in the Valescure catchment,France

This study focuses on a 10-m² plot within a granitic hillslope in Cevennes mountainous area in France, in order to study infiltration and subsurface hydrological processes during heavy rainfalls and flash floods. The monitoring device included water content at several depths (0–70 cm for the shallow soil water; 0–10 m for the deep water) during both intense artificial and natural rainfall events, chemical and physical tracers, time-lapse electrical resistivity tomography. During the most intense events, the infiltrated water was estimated to be some hundreds of millimetres, which largely exceeds the topsoil capacity (≤40 cm deep in most of the cases). The weathered/fractured rock area below the soil clearly has an active role in the water storage and sub-surface flow dynamics. Vertical flow was dominant in the first 0–10 m, and lateral flow was effective at 8–10 m depth, at the top of the saturated area. The speed of the vertical flow was estimated between 1 and 10 m/hr, whereas it was estimated between 0.1 and 1 m/hr for the lateral flow. The interpretation of the experiments might lead to a local pattern of the 2D-hydrological processes and profile properties, which could be generic for most of the mountainous catchments under Mediterranean climate. It suggests that fast triggering of floods at the catchment scale cannot be explained by a mass transfer within the hillslope, but should be due to a pressure wave propagation through the bedrock fractures, which allows exfiltration of the water downstream the hillslope.     

Snowmelt runoff and soil moisture dynamics on steep subalpine hillslopes

Snowmelt water supplies streamflow and growing season soil moisture in mountain regions, yet pathways of snowmelt water and their effects on moisture patterns are still largely unknown. This study examined how flow processes during snowmelt runoff affected spatial patterns of soil moisture on two steep sub‐alpine hillslope transects in Rocky Mountain National Park, CO, USA. The transects have northeast‐facing and east‐facing aspects, and both extend from high‐elevation bedrock outcrops down to streams in valley bottoms. Spatial patterns of both snow depth and near‐surface soil moisture were surveyed along these transects in the snowmelt and summer seasons of 2008–2010. To link these patterns to flow processes, soil moisture was measured continuously on both transects and compared with the timing of discharge in nearby streams. Results indicate that both slopes generated shallow lateral subsurface flow during snowmelt through near‐surface soil, colluvium and bedrock fractures. On the northeast‐facing transect, this shallow subsurface flow emerged through mid‐slope seepage zones, in some cases producing saturation overland flow, whereas the east‐facing slope had no seepage zones or overland flow. At the hillslope scale, earlier snowmelt timing on the east‐facing slope led to drier average soil moisture conditions than on the northeast‐facing slope, but within hillslopes, snow patterns had little relation to soil moisture patterns except in areas with persistent snow drifts. Results suggest that lateral flow and exfiltration processes are key controls on soil moisture spatial patterns in this steep sub‐alpine location. Copyright © 2014 John Wiley & Sons, Ltd.     

The soil-moisture zone in a physically-based hydrologic model

The present paper describes an approach to modelling the unsaturated soil-moisture zone in the framework of an integrated physically-based hydrologic response model. It is supposed that the subsurface flow regime may be viewed as two separate entities — a saturated flow system which may be modelled by standard two-dimensional regional techniques, and a single overlying unsaturated zone in which the flow is essentially vertical. Coupling takes place via the definition of saturation at the lower boundary of the unsaturated zone, and via a conservative water balance. Attention is focused on the computational procedure for the unsaturated zone as a self-contained module. The major difficulties are the definition of the interface between the saturated and unsaturated zones, the nonlinear character of the equation used to describe unsaturated flow, the inclusion of realistic atmospheric boundary conditions, and, the interaction between water uptake by plants and available soil-moisture. Each of these points is discussed, in turn, with the emphasis on mathematically formulating the problem in such a way that the most important physical features are reproduced with a minimal amount of computational effort. The text concludes with a few illustrative examples.     

Runoff generation processes estimated from hydrological observations on a steep forested hillslope with a thin soil layer

In order to understand runoff generation processes on a forested hillslope involving large heterogeneities, this study monitored runoff from a steep hillslope with a thin soil layer as well as matric potential within it and analyzed their responses to storm rainfall. A comparison of storm runoff responses from the study slope with those from two adjacent catchments, one of which includes it, showed that physical properties of the slope reflected the runoff characteristics: although no responses occurred in very dry conditions because of the absence of wet zones near the stream, the area contributing to storm runoff more rapidly extended to the whole slope due to its topographic properties. They also caused its steep hydrographs produced in the wettest condition where almost all the rainfall contributed to storm runoff. In this wettest condition, tensiometric responses near bedrock showed the vertical quick propagation of the rainfall pulse, and a good agreement of storm hydrograph simulated through a kinematic wave runoff model suggested that runoff from the slope was produced by a lateral flow on the bedrock receiving the quick propagation. In a transition process from dry to the wettest conditions, the development of the lateral flow producing smaller responses at the downslope end was estimated from decreasing of matric potential near bedrock from high negative to low values with increasing cumulative rainfall.     

A quasi three-dimensional model for flow and transport in unsaturated and saturated zones: 1. Implementation of the quasi two-dimensional case

A quasi three-dimensional (QUASI 3-D) model is presented for simulating the subsurface water flow and solute transport in the unsaturated and in the saturated zones of soil. The model is based on the assumptions of vertical flow in the unsaturated zone and essentially horizontal groundwater flow. The 1-D Richards equation for the unsaturated zone is coupled at the phreatic surface with the 2-D flow equation for the saturated zone. The latter was obtained by averaging 3-D flow equation in the saturated zone over the aquifer thickness. Unlike the Boussinesq equation for a leaky-phreatic aquifer, the developed model does not contain a storage term with specific yield and a source term for natural replenishment. Instead it includes a water flux term at the phreatic surface through which the Richards equation is linked with the groundwater flow equation. The vertical water flux in the saturated zone is evaluated on the basis of the fluid mass balance equation while the horizontal fluxes, in that equation, are prescribed by Darcy law. A 3-D transport equation is used to simulate the solute migration. A numerical algorithm to solve the problem for the general quasi 3-D case was developed. The developed methodology was exemplified for the quasi 2-D cross-sectional case (QUASI2D). Simulations for three synthetic problems demonstrate good agreement between the results obtained by QUASI2D and two fully 2-D flow and transport codes (SUTRA and 2DSOIL). Yet, simulations with the QUASI2D code were several times faster than those by the SUTRA and the 2DSOIL codes.     

Mixing of event and pre‐event water in a shallow Entisol in sloping farmland based on isotopic and hydrometric measurements,SW China

Water percolation and flow processes in subsurface geologic media play an important role in determining the water source for plants and the transport of contaminants or nutrients, which is essential for water resource management and the development of measures for pollution mitigation. During June 2013, the dynamics of the rainwater, soil water, subsurface flows and groundwater in a shallow Entisol on sloping farmland were monitored using a hydrometric and isotopic approach. The results showed that effective mixing of rainwater and soil water occurred in hours. The rebound phenomenon of δD profiles in soils showed that most isotope‐depleted rainwater largely bypassed the soil matrix when the water saturation in the soil was high. Preferential‐flow, which was the dominant water movement pattern in the vadose zone, occurred through the whole soil profile, and infrequent piston‐flow was mainly found at 20–40 cm in depth. The interflow in the soil layer, composed of 75.2% rainwater, was only generated when the soil profile had been saturated. Underflow in the fractured mudrock was the dominant flow type in this hillslope, and outflow was dominated by base flow (groundwater flow) with a mean contribution of 76.7%. The generation mechanism of underflow was groundwater ridging, which was superimposed upon preferential‐flow composed mainly of rainwater. The quick mixing process of rainwater and soil water and the rapid movement of the mixture through preferential channels in the study soil, which shows a typical bimodal pore size distribution, can explain the prompt release of pre‐event water in subsurface flow. Water sources of subsurface flows at peak discharge could be affected by the antecedent soil water content, rain characteristics and antecedent groundwater levels. Copyright © 2016 John Wiley & Sons, Ltd.