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.     

The role of overland flow and subsurface flow on the spatial distribution of soil moisture in the topsoil

Many investigations show relationships between topographical factors and the spatial distribution of soil moisture in catchments. However, few quantitative analyses have been carried out to elucidate the role of different hydrological processes in the spatial distribution of topsoil moisture in catchments. A spatially distributed rainfall—runoff model was used to investigate contributions of subsurface matric flow, macropore flow and surface runoff to the spatial distribution of soil moisture in a cultivated catchment. The model results show that lateral subsurface flow in the soil matrix or in macropores has a minor effect on the spatial distribution of soil moisture. Only when a perched groundwater table is maintained long enough, which is only possible if the subsurface is completely impermeable, may a spatial distribution in moisture content occur along the slope. Surface runoff, producing accumulations of soil moisture in flat flow paths of agricultural origin (field boundaries), was demonstrated to cause significant spatial variations in soil moisture within a short period after rainfall (<2 days). When significant amounts of surface runoff are produced, wetter moisture conditions will be generated at locations with larger upstream contributing areas. Copyright © 2001 John Wiley & Sons, Ltd.     

Investigating hydrologic connectivity and its association with threshold change in runoff response in a temperate forested watershed

Hydrological studies across varied climatic and physiographic regions have observed small changes in the ‘states of wetness’; based on average soil moisture, can lead to dramatic changes in the amount of water delivered to the stream channel. This non-linear behaviour of the storm response has been attributed to a critical switching in spatial organization of shallow soil moisture and hydrologic connectivity. However, much of the analysis of the role of soil moisture organization and connectivity has been performed in small rangeland catchments. Therefore, we examined the relationship between hydrologic connectivity and runoff response within a temperate forested watershed of moderate relief. We have undertaken spatial surveys of shallow soil moisture over a sequence of storms with varying antecedent moisture conditions. We analyse each survey for evidence of hydrologic connectivity and we monitor the storm response from the catchment outlet. Our results show evidence of a non-linear response in runoff generation over small changes in measures of antecedent moisture conditions; yet, unlike the previous studies of rangeland catchments, in this forested landscape we do not observe a significant change in geostatistical hydrologic connectivity with variations in antecedent moisture conditions. These results suggest that a priori spatial patterns in shallow soil moisture in forested terrains may not always be a good predictor of critical hydrologic connectivity that leads to threshold change in runoff generation, as has been the case in rangeland catchments. Copyright © 2007 John Wiley & Sons, Ltd.     

Controls of streamflow generation in small catchments across the snow–rain transition in the Southern Sierra Nevada,California

Processes controlling streamflow generation were determined using geochemical tracers for water years 2004–2007 at eight headwater catchments at the Kings River Experimental Watersheds in southern Sierra Nevada. Four catchments are snow‐dominated, and four receive a mix of rain and snow. Results of diagnostic tools of mixing models indicate that Ca²⁺, Mg²⁺, K⁺ and Cl^? behaved conservatively in the streamflow at all catchments, reflecting mixing of three endmembers. Using endmember mixing analysis, the endmembers were determined to be snowmelt runoff (including rain on snow), subsurface flow and fall storm runoff. In seven of the eight catchments, streamflow was dominated by subsurface flow, with an average relative contribution (% of streamflow discharge) greater than 60%. Snowmelt runoff contributed less than 40%, and fall storm runoff less than 7% on average. Streamflow peaked 2–4 weeks earlier at mixed rain–snow than snow‐dominated catchments, but relative endmember contributions were not significantly different between the two groups of catchments. Both soil water in the unsaturated zone and regional groundwater were not significant contributors to streamflow. The contributions of snowmelt runoff and subsurface flow, when expressed as discharge, were linearly correlated with streamflow discharge (R² of 0.85–0.99). These results suggest that subsurface flow is generated from the soil–bedrock interface through preferential pathways and is not very sensitive to snow–rain proportions. Thus, a declining of the snow–rain ratio under a warming climate should not systematically affect the processes controlling the streamflow generation at these catchments. Copyright © 2012 John Wiley & Sons, Ltd.     

Analysis of geomorphologic and hydrological characteristics in watershed saturated areas using topographic‐index threshold and geomorphology‐based runoff model

The objective of this paper is to investigate the variation of geomorphology and runoff characteristics in saturated areas under different partial contributing area (PCA) conditions. Geomorphologic information and hydrologic records from two mid‐size watersheds in northern Taiwan were selected for analysis. The PCA ratio in the watershed during a storm was assumed equal to the ratio of the surface‐flow volume to the direct runoff volume from measured hydrologic data. The extents of PCA regions were then determined by using a topographic‐index threshold. Consequently, the geomorphologic factors in saturated and unsaturated areas could be calculated using a digital elevation model, and these factors could then be linked to a geomorphology‐based IUH model for runoff simulation, which can consider both the surface‐ and subsurface‐flow processes in saturated and unsaturated areas, respectively. The results show that geomorphologic characteristics in the saturated areas vary significantly with different PCA ratios especially for higher order streams. A large PCA ratio results in a sharp hydrograph because the quick surface flow dominates the runoff process, whereas the hydrologic response in a low PCA case is dominated by the delayed subsurface flow. Copyright © 2007 John Wiley & Sons, Ltd.     

Thresholds for run‐off generation in a drained closed depression

Hydrological threshold behaviour has been observed across hillslopes and catchments with varying characteristics. Few studies, however, have evaluated rainfall–run‐off response in areas dominated by agricultural land use and artificial subsurface drainage. Hydrograph analysis was used to identify distinct hydrological events over a 9‐year period and examine rainfall characteristics, dynamic water storage, and surface and subsurface run‐off generation in a drained and farmed closed depression in north‐eastern Indiana, USA. Results showed that both surface flow and subsurface tile flow displayed a threshold relationship with the sum of rainfall amount and soil moisture deficit (SMD). Neither surface flow nor subsurface tile flow was observed unless rainfall amount exceeded the SMD. Timing of subsurface tile flow relative to soil moisture response on the shoulder slope of the depression indicated that the formation and drainage of perched water tables on depression hillslopes were likely the main mechanism that produced subsurface connectivity. Surface flow generation was delayed compared with subsurface tile flow during rainfall events due to differences in soil water storage along depression hillslopes and run‐off generation mechanisms. These findings highlight the substantial impact of subsurface tile drainage on the hydrology of closed depressions; the bottom of the depression, the wettest area prior to drainage installation, becomes the driest part of the depression after installation of subsurface drainage. Rapid connectivity of localized subsurface saturation zones during rainfall events is also greatly enhanced because of subsurface drainage. Thus, less fill is required to generate substantial spill. Understanding hydrologic processes in drained and farmed closed depressions is a critical first step in developing improved water and nutrient management strategies in this landscape.     

Hydrological functioning of mediterranean mountain basins in Vallcebre,Catalonia: Some challenges for hydrological modelling

The Vallcebre research catchments are located in the south-eastern Pyrenees, in an area of diverse land use and varying levels of degradation, including forested hillslopes, abandoned agricultural terraces and badland areas. Outside the badlands, the hydrological response is controlled by saturation mechanisms. Between September and June the spatial patterns of saturated areas and soil moisture are determined by subsurface flow, modified by the premature saturation of the inner parts of agricultural terraces, and the negative soil moisture anomalies induced by forest patches overgrown in grassland areas. During summer, this behaviour ceases because of soil moisture depletion and badland surfaces are the only hydrologically active areas, producing excess runoff in response to the intense rainstorms. During the beginning and the ending of the wet season, the saturation of active areas shows a hysteretic behaviour that breaks the linearity between the mean water reserve of the basin and the relative saturated area. © 1997 John Wiley & Sons, Ltd.     

Soil moisture and rainfall intensity thresholds for runoff generation in southwestern Wisconsin agricultural watersheds

The goal of this study was to improve understanding of the factors that influence runoff generation during non‐frozen ground periods in small agricultural watersheds in southwestern Wisconsin where the landscapes are controlled by dolostone bedrock in order to provide agricultural producers with a manure management tool. Six small watersheds (ranging from 6 to 17 ha) within two southwestern Wisconsin farm sites (Discovery Farms Program (DFP) and Pioneer Farm (PF)) were instrumented, and surface runoff was continuously monitored from 2004 to 2007. The soils in all watersheds were formed in deep (~1 m) loessial sites. A direct‐plant management strategy and corn‐soybean crop rotation were utilized within watersheds at DFP. A conventional tillage system (chisel plow in the fall followed by soil finisher in the spring) and a corn‐oat‐alfalfa crop rotation were utilized within watersheds at PF. At PF, the amount of precipitation leaving the landscape as surface runoff (1.8%) was two times greater compared to DFP (0.9%), indicating that the direct‐plant management system was better at retaining precipitation than the chisel plow/soil finisher system. Using breakpoint regression analysis, a non‐linear response in runoff generation with antecedent soil moisture (ASM) was observed with a threshold ASM of 0.39 cm³cm^?3 (approximately 80% of total porosity) for all six watersheds. Below this threshold, runoff coefficients were near zero. Above this threshold, runoff coefficients increased with ASM. A non‐linear response in runoff generation with maximum 30 min rainfall intensity (I30) was also observed, and threshold I30 values increased as ASM decreased and as crop cover increased. Copyright © 2012 John Wiley & Sons, Ltd.     

Field investigation and modelling of coupled stream discharge and shallow water‐table dynamics in a small Mediterranean catchment (Sardinia)

In the semi‐arid Mediterranean environment, the rainfall–runoff relationships are complex because of the markedly irregular patterns in rainfall, the seasonal mismatch between evaporation and rainfall, and the spatial heterogeneity in landscape properties. Watersheds often display considerable non‐linear threshold behavior, which still make runoff generation an open research question. Our objectives in this context were: to identify the primary processes of runoff generation in a small natural catchment; to test whether a physically based model, which takes into consideration only the primary processes, is able to predict spatially distributed water‐table and stream discharge dynamics; and to use the hydrological model to increase our understanding of runoff generation mechanisms. The observed seasonal dynamics of soil moisture, water‐table depth, and stream discharge indicated that Hortonian overland‐flow was negligible and the main mechanism of runoff generation was saturated subsurface‐flow. This gives rise to base‐flow, controls the formation of the saturated areas, and contributes to storm‐flow together with saturation overland‐flow. The distributed model, with a 1D scheme for the kinematic surface‐flow, a 2D sub‐horizontal scheme for the saturated subsurface‐flow, and ignoring the unsaturated flow, performed efficiently in years when runoff volume was high and medium, although there was a smoothing effect on the observed water‐table. In dry years, small errors greatly reduced the efficiency of the model. The hydrological model has allowed to relate the runoff generation mechanisms with the land‐use. The forested hillslopes, where the calibrated soil conductivity was high, were never saturated, except at the foot of the slopes, where exfiltration of saturated subsurface‐flow contributed to storm‐flow. Saturation overland‐flow was only found near the streams, except when there were storm‐flow peaks, when it also occurred on hillslopes used for pasture, where soil conductivity was low. The bedrock–soil percolation, simulated by a threshold mechanism, further increased the non‐linearity of the rainfall–runoff processes. Copyright © 2013 John Wiley & Sons, Ltd.     

A rainfall‐runoff model parameterized from GIS and runoff data

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

Estimating effective impervious area in urban watersheds using land cover,soil character and asymptotic curve number

Knowledge of the effective impervious area (EIA) or the degree to which impervious surfaces are hydraulically connected to the drainage system is useful for improving hydrological and environmental models and assessing the effectiveness of green stormwater infrastructure in urban watersheds. The goal of this research is to develop a method to estimate EIA fraction in urban watersheds using readily available data. Since EIA is dependent on rainfall–runoff response and cannot be solely determined based on the physical characteristics of a watershed, the EIA is linked with the asymptotic curve number (CN), a watershed index that represents runoff characteristics. In order for the method to be applicable to ungauged watersheds, the asymptotic CN is predicted using land cover and soil data from 35 urban catchments in Minnesota and Texas, USA. Similar data from 11 other urban catchments in Wisconsin and Texas, USA, are used to validate the results. A set of runoff depth versus EIA fraction curves is also developed to assess the impact of EIA reduction on discharge from an urban watershed in land-use planning studies.     

Appropriate vertical discretization of Richards' equation for two‐dimensional watershed‐scale modelling

A number of watershed‐scale hydrological models include Richards' equation (RE) solutions, but the literature is sparse on information as to the appropriate application of RE at the watershed scale. In most published applications of RE in distributed watershed‐scale hydrological modelling, coarse vertical resolutions are used to decrease the computational burden. Compared to point‐ or field‐scale studies, application at the watershed scale is complicated by diverse runoff production mechanisms, groundwater effects on runoff production, runon phenomena and heterogeneous watershed characteristics. An essential element of the numerical solution of RE is that the solution converges as the spatial resolution increases. Spatial convergence studies can be used to identify the proper resolution that accurately describes the solution with maximum computational efficiency, when using physically realistic parameter values. In this study, spatial convergence studies are conducted using the two‐dimensional, distributed‐parameter, gridded surface subsurface hydrological analysis (GSSHA) model, which solves RE to simulate vadose zone fluxes. Tests to determine if the required discretization is strongly a function of dominant runoff production mechanism are conducted using data from two very different watersheds, the Hortonian Goodwin Creek Experimental Watershed and the non‐Hortonian Muddy Brook watershed. Total infiltration, stream flow and evapotranspiration for the entire simulation period are used to compute comparison statistics. The influences of upper and lower boundary conditions on the solution accuracy are also explored. Results indicate that to simulate hydrological fluxes accurately at both watersheds small vertical cell sizes, of the order of 1 cm, are required near the soil surface, but not throughout the soil column. The appropriate choice of approximations for calculating the near soil‐surface unsaturated hydraulic conductivity can yield modest increases in the required cell size. Results for both watersheds are quite similar, even though the soils and runoff production mechanisms differ greatly between the two catchments. Copyright © 2003 John Wiley & Sons, Ltd.     

An integral-balance nonlinear model to simulate changes in soil moisture,groundwater and surface runoff dynamics at the hillslope scale

We present a system of ordinary differential equations (ODEs) capable of reproducing simultaneously the aggregated behavior of changes in water storage in the hillslope surface, the unsaturated and the saturated soil layers and the channel that drains the hillslope. The system of equations can be viewed as a two-state integral-balance model for soil moisture and groundwater dynamics. Development of the model was motivated by the need for landscape representation through hillslopes and channels organized following stream drainage network topology. Such a representation, with the basic discretization unit of a hillslope, allows ODEs-based simulation of the water transport in a basin. This, in turn, admits the use of highly efficient numerical solvers that enable space–time scaling studies. The goal of this paper is to investigate whether a nonlinear ODE system can effectively replicate observations of water storage in the unsaturated and saturated layers of the soil. Our first finding is that a previously proposed ODE hillslope model, based on readily available data, is capable of reproducing streamflow fluctuations but fails to reproduce the interactions between the surface and subsurface components at the hillslope scale. However, the more complex ODE model that we present in this paper achieves this goal. In our model, fluxes in the soil are described using a Taylor expansion of the underlying storage flux relationship. We tested the model using data collected in the Shale Hills watershed, a 7.9-ha forested site in central Pennsylvania, during an artificial drainage experiment in August 1974 where soil moisture in the unsaturated zone, groundwater dynamics and surface runoff were monitored. The ODE model can be used as an alternative to spatially explicit hillslope models, based on systems of partial differential equations, which require more computational power to resolve fluxes at the hillslope scale. Therefore, it is appropriate to be coupled to runoff routing models to investigate the effect of runoff and its uncertainty propagation across scales. However, this improved performance comes at the expense of introducing two additional parameters that have no obvious physical interpretation. We discuss the implications of this for hydrologic studies across scales.     

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.     

A simplified stochastic model for infiltration into a heterogeneous soil forced by random precipitation

A simplified stochastic infiltration model is presented, aimed at representing the rainfall-runoff transformation in the presence of heterogeneity in the soil and with precipitation as a random variable with complex temporal evolution. Such a model is based on a simple water mass balance of a surface soil layer, considered as a non-linear reservoir. The explicit inclusion of spatial heterogeneity allows the model to be used in sub-grid parametrizations at a variety of scales, from the distributed modelling of the hydrological response of small watersheds to the representation of surface mass fluxes in General Circulation Models. An approximate solution procedure is developed, which allows the estimation of statistical moments of the soil effective saturation and runoff inside discrete time steps where the hydraulic saturated conductivity and the rainfall intensity are taken as random variables with known probability density functions. As a first test of the proposed model, two different simulations, relative to two soils with different hydraulic conductivity distributions, are presented and discussed. A year long record of hourly averaged rainfall intensities, as measured by a tipping bucket gauge in Central Italy, is taken as the main input. The main finding is that the non-linear nature of the soil filter is such that, for random precipitation intensity, the coefficient of variation of the runoff is always higher than that of precipitation. Such a non-linear variability enhancement, due mainly to the threshold character of the soil mass balance equation, tends to be slightly dampened by the variability of the hydraulic saturated conductivity.     

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.     

On the non-uniqueness of the hydro-geomorphic responses in a zero-order catchment with respect to soil moisture

This study advances mechanistic interpretation of predictability challenges in hydro-geomorphology related to the role of soil moisture spatial variability. Using model formulations describing the physics of overland flow, variably saturated subsurface flow, and erosion and sediment transport, this study explores (1) why a basin with the same mean soil moisture can exhibit distinctly different spatial moisture distributions, (2) whether these varying distributions lead to non-unique hydro-geomorphic responses, and (3) what controls non-uniqueness in relation to the response type. Two sets of numerical experiments are carried out with two physically-based models, HYDRUS and tRIBS+VEGGIE+FEaST, and their outputs are analyzed with respect to pre-storm moisture state. The results demonstrate that distinct spatial moisture distributions for the same mean wetness arise because near-surface soil moisture dynamics exhibit different degrees of coupling with deeper-soil moisture and the process of subsurface drainage. The consequences of such variations are different depending on the type of hydrological response. Specifically, if the predominant runoff response is of infiltration excess type, the degree of non-uniqueness is related to the spatial distribution of near-surface moisture. If runoff is governed by subsurface stormflow, the extent of deep moisture contributing area and its “readiness to drain” determine the response characteristics. Because the processes of erosion and sediment transport superimpose additional controls over factors governing runoff generation and overland flow, non-uniqueness of the geomorphic response can be highly dampened or enhanced. The explanation is sediment composed by multi-size particles can alternate states of mobilization or surface shielding and the transient behavior is inherently intertwined with the availability of mobile particles. We conclude that complex nonlinear dynamics of hydro-geomorphic processes are inherent expressions of physical interactions. As complete knowledge of watershed properties, states, or forcings will always present the ultimate, if ever resolvable, challenge, deterministic predictability will remain handicapped. Coupling of uncertainty quantification methods and space-time physics-based approaches will need to evolve to facilitate mechanistic interpretations and informed practical applications.     

Controls on event runoff coefficients in the eastern Italian Alps

Analyses of event runoff coefficients provide essential insight on catchment response, particularly if a range of catchments and a range of events are compared by a single indicator. In this study we examine the effect of climate, geology, land use, flood types and initial soil moisture conditions on the distribution functions of the event runoff coefficients for a set of 14 mountainous catchments located in the eastern Italian Alps, ranging in size from 7.3 to 608.4 km². Runoff coefficients were computed from hourly precipitation, runoff data and estimates of snowmelt. A total of 535 events were analysed over the period 1989–2004. We classified each basin using a “permeability index” which was inferred from a geologic map and ranged from “low” to “high permeability”. A continuous soil moisture accounting model was applied to each catchment to classify ‘wet’ and ‘dry’ initial soil moisture conditions. The results indicate that the spatial distribution of runoff coefficients is highly correlated with mean annual precipitation, with the mean runoff coefficient increasing with mean annual precipitation. Geology, through the ‘permeability index’, is another important control on runoff coefficients for catchments with mean annual precipitation less than 1200 mm. Land use, as indexed by the SCS curve number, influences runoff coefficient distribution to a lesser degree. An analysis of the runoff coefficients by flood type indicates that runoff coefficients increase with event snowmelt. Results show that there exists an intermediate region of subsurface water storage capacity, as indexed by a flow–duration curve-based index, which maximises the impact of initial wetness conditions on the runoff coefficient. This means that the difference between runoff coefficients characterised by wet and dry initial conditions is negligible both for basins with very large storage capacity and for basins with small storage capacity. For basins with intermediate storage capacities, the impact of the initial wetness conditions may be relatively large.     

Towards the representation of groundwater in the Joint UK Land Environment Simulator

Groundwater is an important component of the hydrological cycle with significant interactions with soil hydrological processes. Recent studies have demonstrated that incorporating groundwater hydrology in land surface models (LSMs) considerably improves the prediction of the partitioning of water components (e.g., runoff and evapotranspiration) at the land surface. However, the Joint UK Land Environment Simulator (JULES), an LSM developed in the United Kingdom, does not yet have an explicit representation of groundwater. We propose an implementation of a simplified groundwater flow boundary parameterization (JULES-GFB), which replaces the original free drainage assumption in the default model (JULES-FD). We tested the two approaches under a controlled environment for various soil types using two synthetic experiments: (1) single-column and (2) tilted-V catchment, using a three-dimensional (3-D) hydrological model (ParFlow) as a benchmark for JULES’ performance. In addition, we applied our new JULES-GFB model to a regional domain in the UK, where groundwater is the key element for runoff generation. In the single-column infiltration experiment, JULES-GFB showed improved soil moisture dynamics in comparison with JULES-FD, for almost all soil types (except coarse soils) under a variety of initial water table depths. In the tilted-V catchment experiment, JULES-GFB successfully represented the dynamics and the magnitude of saturated and unsaturated storage against the benchmark. The lateral water flow produced by JULES-GFB was about 50% of what was produced by the benchmark, while JULES-FD completely ignores this process. In the regional domain application, the Kling-Gupta efficiency (KGE) for the total runoff simulation showed an average improvement from 0.25 for JULES-FD to 0.75 for JULES-GFB. The mean bias of actual evapotranspiration relative to the Global Land Evaporation Amsterdam Model (GLEAM) product was improved from −0.22 to −0.01 mm day⁻¹. Our new JULES-GFB implementation provides an opportunity to better understand the interactions between the subsurface and land surface processes that are dominated by groundwater hydrology.