Spatio‐temporal variability of piezometric response on two steep alpine hillslopes

Information on the main drivers of subsurface flow generation on hillslopes of alpine headwater catchments is still missing. Therefore, the dominant factors controlling the water table response to precipitation at the hillslope scale in the alpine Bridge Creek Catchment, Northern Italy, were investigated. Two steep hillslopes of similar size, soil properties and vegetation cover but contrasting topography were instrumented with 24 piezometric wells. Sixty‐three (63) rainfall‐runoff events were selected over three years in the snow‐free months to analyse the influence of rainfall depth, antecedent moisture conditions, hillslope topographic characteristics and soil depth on shallow water table dynamics. Piezometric response, expressed as percentage of well activation and water peak magnitude, was strongly correlated with soil moisture status, as described by an index combining antecedent soil moisture and rainfall depth. Hillslope topography was found to be a dominant control only for the convex‐divergent hillslope and during wet conditions. Timing of water table response depended primarily on soil depth and topographic position, with piezometric peak response occurring later and showing a greater temporal variability at the hillslope bottom, characterized by thicker soil. The relationship between mean hillslope water table level and standard deviation for all wells reflected the timing of the water table response at the different locations along the hillslopes. The outcomes of this research contribute to a better understanding of the controls on piezometric response at the hillslope scale in steep terrain and its role on the hydrological functioning of the study catchment and of other sites with similar physiographic characteristics. Copyright © 2014 John Wiley & Sons, Ltd.     

Gypsies in the palace: experimentalist's view on the use of 3‐D physics‐based simulation of hillslope hydrological response

As a fundamental unit of the landscape, hillslopes are studied for their retention and release of water and nutrients across a wide range of ecosystems. The understanding of these near‐surface processes is relevant to issues of runoff generation, groundwater–surface water interactions, catchment export of nutrients, dissolved organic carbon, contaminants (e.g. mercury) and ultimately surface water health. We develop a 3‐D physics‐based representation of the Panola Mountain Research Watershed experimental hillslope using the TOUGH2 sub‐surface flow and transport simulator. A recent investigation of sub‐surface flow within this experimental hillslope has generated important knowledge of threshold rainfall‐runoff response and its relation to patterns of transient water table development. This work has identified components of the 3‐D sub‐surface, such as bedrock topography, that contribute to changing connectivity in saturated zones and the generation of sub‐surface stormflow. Here, we test the ability of a 3‐D hillslope model (both calibrated and uncalibrated) to simulate forested hillslope rainfall‐runoff response and internal transient sub‐surface stormflow dynamics. We also provide a transparent illustration of physics‐based model development, issues of parameterization, examples of model rejection and usefulness of data types (e.g. runoff, mean soil moisture and transient water table depth) to the model enterprise. Our simulations show the inability of an uncalibrated model based on laboratory and field characterization of soil properties and topography to successfully simulate the integrated hydrological response or the distributed water table within the soil profile. Although not an uncommon result, the failure of the field‐based characterized model to represent system behaviour is an important challenge that continues to vex scientists at many scales. We focus our attention particularly on examining the influence of bedrock permeability, soil anisotropy and drainable porosity on the development of patterns of transient groundwater and sub‐surface flow. Internal dynamics of transient water table development prove to be essential in determining appropriate model parameterization. Copyright © 2010 John Wiley & Sons, Ltd.     

Hillslope hydrology and wetland response of two small zero‐order boreal catchments on the Precambrian Shield

Two Precambrian Shield zero‐order catchments were monitored from January 2003 to July 2004 to characterize their hydrological and biogeochemical characteristics prior to a forest management experiment. Hydrometric observations were used to examine temporal trends in hillslope‐wetland connectivity and the hillslope runoff processes that control wetland event response. The hillslope groundwater flux from the longer transect (E1) was continuous throughout the study period. Groundwater fluxes from a shorter and steeper hillslope (E0) were intermittent during the study period. Large depression storage elements (termed micro‐basins) located on the upper hillslope of the E1 catchment appeared to be at least partly responsible for the observed rapid wetland runoff responses. These micro‐basins were hydrologically connected to a downslope wetland by a subsurface channel of glacial cobbles that functioned as a macropore channel during episodic runoff events. The runoff response from the hilltop micro‐basins is controlled by antecedent water table position and water is quickly piped to the wetland fringe through the cobble channel during high water table conditions. During periods of low water table position, seepage along the bedrock–soil interface from the hilltop micro‐basin and other hillslopes maintained hillslope–wetland connectivity. The micro‐basins create a dynamic variable source‐area runoff system where the contributing area expands downslope during episodic runoff events. The micro‐basins occupied 30% of the E1 catchment and are a common feature on the Precambrian Shield. Copyright © 2007 John Wiley & Sons, Ltd.     

Physically‐based modelling of double‐peak discharge responses at Slapton Wood catchment

Heavy winter rainfall produces double‐peak hydrographs at the Slapton Wood catchment, Devon, UK. The first peak is saturation‐excess overland flow in the hillslope hollows and the second (i.e. the delayed peak) is subsurface stormflow. The physically‐based spatially‐distributed model SHETRAN is used to try to improve the understanding of the processes that cause the double peaks. A three‐stage (multi‐scale) approach to calibration is used: (1) water balance validation for vertical one‐dimensional flow at arable, grassland and woodland plots; (2) two‐dimensional flow for cross‐sections cutting across the stream valley; and (3) three‐dimensional flow in the full catchment. The main data are for rainfall, stream discharge, evaporation, soil water potential and phreatic surface level. At each scale there was successful comparison with measured responses, using as far as possible parameter values from measurements. There was some calibration but all calibrated values at one scale were used at a larger scale. A large proportion of the subsurface runoff enters the stream from three dry valleys (hillslope hollows), and previous studies have suggested convergence of the water in the three large hollows as being the major mechanism for the production of the delayed peaks. The SHETRAN modelling suggests that the hillslopes that drain directly into the stream are also involved in producing the delayed discharges. The model shows how in the summer most of the catchment is hydraulically disconnected from the stream. In the autumn the catchment eventually ‘wets up’ and shallow subsurface flows are produced, with water deflected laterally along the soil‐bedrock interface producing the delayed peak in the stream hydrograph. Copyright © 2007 John Wiley & Sons, Ltd.     

A low‐dimensional hillslope‐based catchment model for layered groundwater flow

Despite the strong interaction between surface and subsurface waters, groundwater flow representation is often oversimplified in hydrological models. For instance, the interplay between local or shallow aquifers and deeper regional‐scale aquifers is typically neglected. In this work, a novel hillslope‐based catchment model for the simulation of combined shallow and deep groundwater flow is presented. The model consists of the hillslope‐storage Boussinesq (hsB) model representing shallow groundwater flow and an analytic element (AE) model representing deep regional groundwater flow. The component models are iteratively coupled via a leakage term based on Darcy's law, representing delayed recharge to the regional aquifer through a low conductivity layer. Simulations on synthetic single hillslopes and on a two‐hillslope open‐book catchment are presented, and the results are compared against a benchmark three‐dimensional Richards equation model. The impact of hydraulic conductivity, hillslope plan geometry (uniform, convergent, divergent), and hillslope inclination (0.2%, 5%, and 30%) under drainage and recharge conditions are examined. On the single hillslopes, good matches for heads, hydrographs, and exchange fluxes are generally obtained, with the most significant differences in outflows and heads observed for the 30% slope and for hillslopes with convergent geometry. On the open‐book catchment, cumulative outflows are overestimated by 1–4%. Heads in the confined and unconfined aquifers are adequately reproduced throughout the catchment, whereas exchange fluxes are found to be very sensitive to the hillslope drainable porosity. The new model is highly efficient computationally compared to the benchmark model. The coupled hsB/AE model represents an alternative to commonly used groundwater flow representations in hydrological models, of particular appeal when surface–subsurface exchanges, local aquifer–regional aquifer interactions, and low flows play a key role in a watershed's dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

On the interrelations between topography,soil depth,soil moisture,transpiration rates and species distribution at the hillslope scale

Relations between the spatial patterns of soil moisture, soil depth, and transpiration and their influence on the hillslope water balance are not well understood. When determining a water balance for a hillslope, small scale variations in soil depth are often ignored. In this study we found that these variations in soil depth can lead to distinct patterns in transpiration rates across a hillslope. We measured soil moisture content at 0.05 and 0.10 m depth intervals between the soil surface and the soil–bedrock boundary on 64 locations across the trenched hillslope in the Panola Mountain Research Watershed, Georgia, USA. We related these soil moisture data to transpiration rates measured in 14 trees across the hillslope using 28 constant heat sapflow sensors. Results showed a lack of spatial structure in soil moisture across the hillslope and with depth when the hillslope was in either the wet or the dry state. However, during the short transition period between the wet and dry state, soil moisture did become spatially organized with depth and across the hillslope. Variations in soil depth and thus total soil water stored in the soil profile at the end of the wet season caused differences in soil moisture content and transpiration rates between upslope and midslope sections at the end of the summer. In the upslope section, which has shallower soils, transpiration became limited by soil moisture while in the midslope section with deeper soils, transpiration was not limited by soil moisture. These spatial differences in soil depth, total water available at the end of the wet season and soil moisture content during the summer appear responsible for the observed spatial differences in basal area and species distribution between the upslope and midslope sections of the hillslope.     

A hillslope-based hydrological model using catchment area and width functions

Abstract

A flow-interval hillslope discretization scheme is proposed for catchment hydrological modelling. By this scheme, a two-dimensional catchment is simplified into a one-dimensional cascade of flow intervals linked by the main stream. Each flow interval comprises a set of parallel hillslopes. The hillslope is the fundamental computational unit in the hydrological model providing lateral inflow to the main stream. The size of hillslope is determined by the catchment area and width functions. Catchment runoff is the total of hillslope responses through the river routing. Tests in four Japanese catchments showed that the model performed well on simulating the overall water balance, general flow pattern, and daily and hourly hydrographs of a whole catchment, as well as simultaneous simulation in different subcatchments. Characteristics of catchment hydrological responses and model applicability are discussed.     

SULPHATE DYNAMICS IN RELATION TO GROUNDWATER–SURFACE WATER INTERACTIONS IN HEADWATER WETLANDS OF THE SOUTHERN CANADIAN SHIELD

The spatial and temporal distribution of sulphate (SO₄) concentrations in peat pore water and the outlet streams of two forested swamps was related to variations in the magnitude of upland runoff, wetland water levels and flow path. The swamps were located in headwater catchments with contrasting till depths typical of the southern Canadian Shield. Inputs of SO₄ from shallow hillslope tills and streams showed little seasonal variation in either source or concentration in both swamps. Sulphate dynamics at the outlet stream reflected hydrological and biogeochemical processes within the valley wetlands, which in turn were partly controlled by catchment hydrogeology. During high runoff, maximum water table elevations and peak surface flow in the swamps resulted in upland inputs largely bypassing anoxic peat. Consequently, SO₄ concentrations of 8–10 mg/l at the swamp outlets were similar to stream and groundwater inputs. During periods of low flow, concentrations of SO₄ at the swamp outlets declined to less than 3 mg/l. At this time lower water table elevations resulted in increased interaction of input water with anoxic peats, and therefore, SO₄ reduction. Contrasts in till depth and the nature of groundwater flow between catchments resulted in differences in SO₄ dynamics between years and swamps. In dry summers the absence of groundwater inputs to the swamp in the catchment with thin till resulted in a large water table drawdown and re-oxidation of accumulated S, which contributed to maximum SO₄ concentrations (up to 35 mg/l) during storm runoff. Continuous groundwater input to the swamp in the catchment with deeper till was critical to maintaining saturated surfaces and efficient SO₄ retention during both dry and wet summers. A conceptual model of wetland SO₄ retention and export, based on catchment hydrogeology, is developed to generalize the SO₄ dynamics of valley bottom wetlands at the landscape scale. © 1997 by John Wiley & Sons, Ltd.     

Sub‐grid scale parameterization of hillslope runoff and erosion processes for catchment‐scale models of semi‐arid landscapes

The processes of hillslope runoff and erosion are typically represented at coarse spatial resolution in catchment‐scale models due to computational limitations. Such representation typically fails to incorporate the important effects of topographic heterogeneity on runoff generation, overland flow, and soil erosion. These limitations currently undermine the application of distributed catchment models to understand the importance of thresholds and connectivity on hillslope and catchment‐scale runoff and erosion, particularly in semi‐arid environments. This paper presents a method for incorporating high‐resolution topographic data to improve sub‐grid scale parameterization of hillslope overland flow and erosion models. Results derived from simulations conducted using a kinematic wave overland flow model at 0.5 m spatial resolution are used to parameterize the depth–discharge relationship in the overland flow model when applied at 16 m resolution. The high‐resolution simulations are also used to derive a more realistic parameterization of excess flow shear stress for use in the 16 m resolution erosion model. Incorporating the sub‐grid scale parameterization in the coarse‐resolution model (16 m) leads to improved predictions of overland flow and erosion when evaluated using results derived from high‐resolution (0.5 m) model simulations. The improvement in performance is observed for a range of event magnitudes and is most notable for erosion estimates due to the non‐linear dependency between the rates of erosion and overland flow. Copyright © 2013 John Wiley & Sons, Ltd.     

Water storage dynamics in a till hillslope: the foundation for modeling flows and turnover times

Studies on hydrology, biogeochemistry, or mineral weathering often rely on assumptions about flow paths, water storage dynamics, and transit times. Testing these assumptions requires detailed hydrometric data that are usually unavailable at the catchment scale. Hillslope studies provide an alternative for obtaining a better understanding, but even on such well‐defined and delimited scales, it is rare to have a comprehensive set of hydrometric observations from the water divide down to the stream that can constrain efforts to quantify water storage, movement, and turnover time. Here, we quantified water storage with daily resolution in a hillslope during the course of almost an entire year using hydrological measurements at the study site and an extended version of the vertical equilibrium model. We used an exponential function to simulate the relationship between hillslope discharge and water table; this was used to derive transmissivity profiles along the hillslope and map mean pore water velocities in the saturated zone. Based on the transmissivity profiles, the soil layer transmitting 99% of lateral flow to the stream had a depth that ranged from 8.9 m at the water divide to under 1 m closer to the stream. During the study period, the total storage of this layer varied from 1189 to 1485 mm, resulting in a turnover time of 2172 days. From the pore water velocities, we mapped the time it would take a water particle situated at any point of the saturated zone anywhere along the hillslope to exit as runoff. Our calculations point to the strengths as well as limitations of simple hydrometric data for inferring hydrological properties and water travel times in the subsurface.     

Identifying and linking source areas of flow and P transport in dairy‐grazed headwater catchments,North Island,New Zealand

We examined the applicability of the critical‐source area (CSA) concept to the dairy‐grazed 192‐ha Upper Toenepi catchment and its 8·7‐ha Kiwitahi sub‐catchment, New Zealand. We evaluated if phosphorus (P) transport from land into stream is dominated by saturation‐excess (SE) and infiltration‐excess (IE) runoff during stormflow and by sub‐surface (<1·5 m depth) flows during baseflow. We measured stream flow and shallow groundwater levels, collected monthly stream, tile drain (TDA) and groundwater samples, and flow‐proportional stream samples from the Kiwitahi sub‐catchment, and determined their dissolved reactive phosphorus (DRP) and total phosphorus (TP) concentrations. In the Kiwitahi sub‐catchment, during storm events, IE contributions were significant. Contributions from SE appeared significant in the Upper Toenepi catchment. However, in both catchments, sub‐surface contributions dominated stormflow and baseflow periods. Absence of water table at the surface and the water table gradient towards the stream indicated that P transport during events was not limited to surface runoff. The dynamics of the groundwater table and the occurrence of SE areas were influenced by proximity to the stream and hillslope positions. Baseflow accounted for 42% of the annual flow in the Kiwitahi sub‐catchment, and contributed 37 and 52% to the DRP and TP loads, respectively. The P transport during baseflow appeared equally important as P losses from CSAs during stormflow. The close resemblance in P levels between groundwater and stream samples during baseflow demonstrates the importance of shallow groundwater for stream flow. In the Upper Toenepi catchment, contributions from effluent ponds (EFFs) dominated P loads. Management strategies should focus on controlling P release from EFFs, and on decreasing Olsen P concentrations in soil to minimize leaching of P via sub‐surface flow to streams. Research is needed to quantify the role of sub‐surface flow as well as to expand management strategies to minimize P transfers during stormflow and baseflow conditions. Copyright © 2010 John Wiley & Sons, Ltd.     

The roles of channels and hillslopes in rainfall/run‐off lag times during intense storms in a steep catchment

The response time (lag time) between rainfall input and run‐off output in headwater catchments is a key parameter for flood prediction. Lag times are expected to be controlled by run‐off processes, both on hillslopes and in channels. To demonstrate these effects on peak lag times within a 4.5‐km² catchment, we measured stream water levels at up to 16 channel locations at 1‐min intervals and compared the lag times with topographic indices describing the length and gradient of the hillslope and channel flow path. We captured storm events with a total precipitation of 38–198 mm and maximum hourly precipitation intensity of 9–90 mm/hr. There were positive relationships between lag time and flow path length as well as the ratio of the flow path length and the square root of the gradient of channels for the most intense storms, demonstrating that channel flow paths generally defined the variation in lag times. Topographic analysis showed that hillslope flow path lengths were similar among locations, whereas channel flow path length increased almost one order of magnitude with a 100‐fold increase in catchment area. Thus, the relative importance of hillslope flow path decreased with increasing catchment area. Our results indicate that the variation in lag times is small when hillslopes are sufficiently wet; thus, catchment‐scale variation in lag times can be explained almost entirely by channel processes. Detailed topographic channel information can improve prediction of flood peak timing, whereas hillslopes can be treated as homogeneous during large flood events.     

Precipitation‐runoff processes in Shimen hillslope micro‐catchment of Taihang Mountain,north China

Precipitation runoff is a critical hillslope hydrological process for downslope streamflow and piedmont/floodplain recharge. Shimen hillslope micro‐catchment is strategically located in the central foothill region of Taihang Mountains, where runoff is crucial for water availability in the piedmont corridors and floodplains of north China. This study analyzes precipitation‐runoff processes in the Shimen hillslope micro‐catchment for 2006–2008 using locally designed runoff collection systems. The study shows that slope length is a critical factor, next only to precipitation, in terms of runoff yield. Regression analysis also shows that runoff is related positively to precipitation, and negatively to slope length. Soil mantle in the study area is generally thin and is therefore not as critical a runoff factor as slope length. The study shows a significant difference between overland and subsurface runoff. However, that between the 0–10 and 10–20 cm subsurfaces is insignificant. Runoff hardly occurs under light rains (<10 mm), but is clearly noticeable under moderate‐to‐rainstorm events. In the hillslope catchment, vertical infiltration (accounting for 42–84% of the precipitation) dominates runoff processes in subsurface soils and weathered granite gneiss bedrock. A weak lateral flow (at even the soil/bedrock interface) and the generally small runoff suggest strong infiltration loss via deep percolation. This is critical for groundwater recharge in the downslope piedmont corridors and floodplains. This may enhance water availability, ease water shortage, avert further environmental degradation, and reduce the risk of drought/flood in the event of extreme weather conditions in the catchment and the wider north China Plain. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling of sediment dynamics in a laboratory‐scale experimental catchment

The variability of hillslope form and function is examined experimentally using a simple model catchment in which most landscape development parameters are either known or controlled. It is demonstrated that there is considerable variability in sediment output from similar catchments, subjected to the same hydrological processes, and for which the initial hillslope profiles are the same. The results demonstrate that, in the case of catchments with a linear initial hillslope profile, the sediment output is initially high but reduces through time, whereas for a concave initial profile the sediment output was smaller and relatively constant. Concave hillslope profiles also displayed reduced sediment output when compared with linear slopes with the same overall slope. Using this experimental model catchment data, the SIBERIA landscape evolution model was tested for its ability to predict temporal sediment transport. When calibrated for the rainfall and erodible material, SIBERIA is able to simulate mean temporal sediment output for the experimental catchment over a range of hillslope profiles and rainfall intensities. SIBERIA is also able to match the hillslope profile of the experimental catchments. The results of the study provide confidence in the ability of SIBERIA to predict temporal sediment output. The experimental and modelling data also demonstrate that, even with all geomorphic and hydrological variables being known and/or controlled, there is still a need for long‐term stream gauging to obtain reliable assessments of field catchment hydrology and sediment transport. Copyright © 2005 John Wiley & Sons, Ltd.     

Lumped hysteretic model for subsurface stormflow developed using downward approach

The terms ‘downward’ and ‘upward’ (synonymous with ‘top‐down’ and ‘bottom‐up’ respectively) are sometimes used when describing methods for developing hydrological models. A downward approach is used here to develop a lumped catchment‐scale model for subsurface stormflow at the 0·94 km² Slapton Wood catchment. During the development, as few assumptions as possible are made about the behaviour of subsurface stormflow at the catchment scale, and no assumptions are made about its behaviour at smaller scales. (In an upward approach, in contrast, the modelling would be based on assumptions about, and data for, the behaviour at smaller scales, such as the hillslope, plot, and point scales.) The model has a single store with a relatively simple relationship between discharge and storage, based on equations describing hysteretic patterns seen in a graph of discharge against storage. Double‐peaked hydrographs have been observed at the catchment outlet. Rainfall on the channel and infiltration‐excess and saturation‐excess runoff give a rapid response, and shallow subsurface stormflow gives a delayed response. Hydrographs are successfully simulated for the large delayed responses observed in 1971–1980 and 1989–1991, then a lumped model for the rapid response is coupled to the lumped hysteretic model and some double‐peaked hydrographs simulated. A physical interpretation is developed for the lumped hysteretic model, making use of information on patterns of perched saturation observed in 1982 on a hillslope at the Slapton Wood catchment. Downward and upward approaches are complementary, and the most robust way to develop and improve lumped catchment models is to iterate between downward and upward steps. Possible next steps are described. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

Seasonal and interannual variations of nitrate and chloride in stream waters related to spatial and temporal patterns of groundwater concentrations in agricultural catchments

Nitrate concentrations in streamwater of agricultural catchments often exhibit interannual variations, which are supposed to result from land‐use changes, as well as seasonal variations mainly explained by the effect of hydrological and biogeochemical cycles. In catchments on impervious bedrock, seasonal variations of nitrate concentrations in streamwater are usually characterized by higher nitrate concentrations in winter than in summer. However, intermediate or inverse cycles with higher concentrations in summer are sometimes observed. An experimental study was carried out to assess the mechanisms that determine the seasonal cycles of streamwater nitrate concentrations in intensive agricultural catchments. Temporal and spatial patterns of groundwater concentrations were investigated in two adjacent catchments located in south‐western Brittany (France), characterized by different seasonal variations of streamwater nitrate concentrations. Wells were drilled across the hillslope at depths ranging from 1·5 to 20 m. Dynamics of the water table were monitored and the groundwater nitrate and chloride concentrations were measured weekly over 2 years. Results highlighted that groundwater was partitioned into downslope domains, where denitrification induced lower nitrate concentrations than into mid‐slope and upslope domains. For one catchment, high subsurface flow with high nitrate concentrations during high water periods and active denitrification during low water periods explained the higher streamwater nitrate concentrations in winter than in summer. For the other catchment, the high contribution of groundwater with high nitrate concentrations smoothed or inverted this trend. Increasing bromide/chloride ratio and nitrate concentrations with depth argued for an effect of past agricultural pressure on this catchment. The relative contribution of flows in time and correlatively the spatial origin of waters, function of the depth and the location on the hillslope, and their chemical characteristics control seasonal cycles of streamwater nitrate concentrations and can influence their interannual trends. Copyright © 2004 John Wiley & Sons, Ltd.     

Infiltration into frozen soil: From core‐scale dynamics to hillslope‐scale connectivity

Infiltration into frozen soil is a key hydrological process in cold regions. Although the mechanisms behind point‐scale infiltration into frozen soil are relatively well understood, questions remain about upscaling point‐scale results to estimate hillslope‐scale run‐off generation. Here, we tackle this question by combining laboratory, field, and modelling experiments. Six large (0.30‐m diameter by 0.35‐m deep) soil cores were extracted from an experimental hillslope on the Canadian Prairies. In the laboratory, we measured run‐off and infiltration rates of the cores for two antecedent moisture conditions under snowmelt rates and diurnal freeze–thaw conditions observed on the same hillslope. We combined the infiltration data with spatially variable data from the hillslope, to parameterise a surface run‐off redistribution model. We used the model to determine how spatial patterns of soil water content, snowpack water equivalent (SWE), and snowmelt rates affect the spatial variability of infiltration and hydrological connectivity over frozen soil. Our experiments showed that antecedent moisture conditions of the frozen soil affected infiltration rates by limiting the initial soil storage capacity and infiltration front penetration depth. However, shallow depths of infiltration and refreezing created saturated conditions at the surface for dry and wet antecedent conditions, resulting in similar final infiltration rates (0.3 mm hr^?1). On the hillslope‐scale, the spatial variability of snowmelt rates controlled the development of hydrological connectivity during the 2014 spring melt, whereas SWE and antecedent soil moisture were unimportant. Geostatistical analysis showed that this was because SWE variability and antecedent moisture variability occurred at distances shorter than that of topographic variability, whereas melt variability occurred at distances longer than that of topographic variability. The importance of spatial controls will shift for differing locations and winter conditions. Overall, our results suggest that run‐off connectivity is determined by (a) a pre‐fill phase, during which a thin surface soil layer wets up, refreezes, and saturates, before infiltration excess run‐off is generated and (b) a subsequent fill‐and‐spill phase on the surface that drives hillslope‐scale run‐off.