Macropores and preferential flow—a love‐hate relationship

Preferential flow is of high relevance for runoff generation, transport of chemicals and nutrients, and the transit time distribution of water in the soil or watershed. However, preferential flow effects are generally ignored in lumped hydrological models. And even most physically‐based models ignore macropores and preferential flow features at the soil and hillslope scale. Keith Beven was never satisfied with this situation and he tried again and again to convince the scientific community to focus their research on the complex topic of macropore and preferential flow. Although he recognized how difficult it is to correctly include preferential flow in hydrological models, he made substantial progress defining and describing macropore flow and showing its relevance, developing models to simulate preferential flow, and in particular, the interaction between macropores and the soil matrix. In this short commentary, I reflect on these achievements and outline a vision for research in preferential flow experiments and modeling.     

Integrating observations and models to determine the effect of seasonally frozen ground on hydrologic partitioning in alpine hillslopes in the Colorado Rocky Mountains,USA

This study integrated spatially distributed field observations and soil thermal models to constrain the impact of frozen ground on snowmelt partitioning and streamflow generation in an alpine catchment within the Niwot Ridge Long-Term Ecological Research site, Colorado, USA. The study area was comprised of two contrasting hillslopes with notable differences in topography, snow depth and plant community composition. Time-lapse electrical resistivity surveys and soil thermal models enabled extension of discrete soil moisture and temperature measurements to incorporate landscape variability at scales and depths not possible with point measurements alone. Specifically, heterogenous snowpack thickness (~0–4 m) and soil volumetric water content between hillslopes (~0.1–0.45) strongly influenced the depths of seasonal frost, and the antecedent soil moisture available to form pore ice prior to freezing. Variable frost depths and antecedent soil moisture conditions were expected to create a patchwork of differing snowmelt infiltration rates and flowpaths. However, spikes in soil temperature and volumetric water content, as well as decreases in subsurface electrical resistivity revealed snowmelt infiltration across both hillslopes that coincided with initial decreases in snow water equivalent and early increases in streamflow. Soil temperature, soil moisture and electrical resistivity data from both wet and dry hillslopes showed that initial increases in streamflow occurred prior to deep soil water flux. Temporal lags between snowmelt infiltration and deeper percolation suggested that the lateral movement of water through the unsaturated zone was an important driver of early streamflow generation. These findings provide the type of process-based information needed to bridge gaps in scale and populate physically based cryohydrologic models to investigate subsurface hydrology and biogeochemical transport in soils that freeze seasonally.     

The relevance of preferential flow in catchment scale simulations: Calibrating a 3D dual-permeability model using DREAM

The occurrence of preferential flow in the subsurface has often been shown in field experiments. However, preferential flow is rarely included in models simulating the hydrological response at the catchment scale. If it is considered, preferential flow parameters are typically determined at the plot scale and then transferred to larger-scale simulations. Here, we successfully used the optimization algorithm DiffeRential Evolution Adaptive Metropolis (DREAM) to calibrate a 3D physics-based dual-permeability model directly at the catchment scale. In order to keep computational costs of the optimization routine at a reasonable level, we limited the number of parameters to be calibrated to the ones that had been shown before to be most influential for the simulation of discharge. We also calibrated parameters of the matrix domain and the macropore domain with a fixed parameter ratio between soil layers instead of calibrating every layer separately. These ratios reflected observed depth profiles of soil hydraulic properties at our study site. The dual-permeability parameter sets identified during calibration were able to simulate observed discharge time series satisfactorily but did not outperform a calibrated single-domain reference model scenario. Saturated hydraulic conductivities of the macropore domain were calibrated such that they became very similar to matrix saturated hydraulic conductivities, thereby effectively removing the effect of macropores. This suggests that the incorporation of vertical preferential flow as represented by the dual-permeability approach was not relevant for reproducing the hydrometric response reasonably well in the studied catchment. We also tested the scale-invariance of the calibrated dual-permeability parameter sets by using the parameter sets performing best at catchment scale to simulate plot-scale bromide depth profiles obtained from tracer irrigation experiments. This parameter transfer proved to be not successful, indicating that soil hydraulic parameters are scale-variant, independent of the direction of parameter transfer.     

Directions of preferential flow in a hillslope soil, 1. Quasi‐steady flow

Preferred infiltration is mainly perceived as vertically down whereas subsurface storm flow is thought to occur parallel to slopes. The transition from vertical to lateral flow in a layered hillslope soil is the focus of the contribution. Transient flow is assumed to move as a wetting front. Three time‐domain reflectometry (TDR) wave‐guides, each 0·15 m long, were mounted in the shape of a truncated tetrahedron with its peak pointing down. Each wave‐guide focuses the front velocity along its axis. The three front‐velocity vectors are decomposed into their x, y and z components, which are then assembled to the resultant velocity vector. The volume density flux of preferred flow is the product of the front velocity and the mobile water content. The latter is the amplitude of transient soil moisture measured with each wave‐guide. The resultant vector of the volume flux density is computed similarly to the velocity vector. The experimental approach allows for the rapid assessment of transient flows without relying on the variation of water potentials. The experiments indicate that the directions of the resultant vectors of velocity and volume flux density can be estimated if the moisture variations of the three TDR wave‐guides are strongly correlated during the passing of the wetting front. Copyright © 2004 John Wiley & Sons, Ltd.     

Effects of antecedent moisture and macroporosity on infiltration and water flow in frozen soil

Infiltration into frozen soil plays an important role in soil freeze–thaw and snowmelt-driven hydrological processes. To better understand the complex thermal energy and water transport mechanisms involved, the influence of antecedent moisture content and macroporosity on infiltration into frozen soil was investigated. Ponded infiltration experiments on frozen macroporous and non-macroporous soil columns revealed that dry macroporous soil produced infiltration rates reaching 10³ to 10⁴ mm day⁻¹, two to three orders of magnitude larger than dry non-macroporous soil. Results suggest that rapid infiltration and drainage were a result of preferential flow through initially air-filled macropores. Using recorded flow rates and measured macropore characteristics, calculations indicated that a combination of both saturated flow and unsaturated film flow likely occurred within macropores. Under wet conditions, regardless of the presence of macropores, infiltration was restricted by the slow thawing rate of pore ice, producing infiltration rates of 2.8 to 5.0 mm day⁻¹. Reduced preferential flow under wet conditions was attributed to a combination of soil swelling, due to smectite-rich clay (that reduced macropore volume), and pore ice blockage within macropores. In comparison, dry soil column experiments demonstrated that macropores provided conduits for water and thermal energy to bypass the frozen matrix during infiltration, reducing thaw rates compared with non-macroporous soils. Overall, results showed the dominant control of antecedent moisture content on the initiation, timing, and magnitude of infiltration and flow in frozen macroporous soils, as well as the important role of macropore connectivity. The study provides an important data set that can aid the development of hydrological models that consider the interacting effects of soil freeze–thaw and preferential flow on snowmelt partitioning in cold regions.     

Quantifying annual groundwater recharge and storage in the central Sierra Nevada using naturally occurring 35S

Identifying aquifer vulnerability to climate change is of vital importance in the Sierra Nevada and other snow‐dominated basins where groundwater systems are essential to water supply and ecosystem health. Quantifying the component of new (current year's) snowmelt in groundwater and surface water is useful in evaluating aquifer vulnerability because significant annual recharge may indicate that streamflow will respond rapidly to annual variability in precipitation, followed by more gradual decreases in recharge as recharge declines over decades. Hydrologic models and field‐based studies have indicated that young (<1 year) water is an important component of streamflow. The goal of this study was to utilize the short‐lived, naturally occurring cosmogenic isotope sulfur‐35 (³⁵S) to quantify new snowmelt contribution to groundwater and surface waters in Sagehen Creek Basin (SCB) and Martis Valley Groundwater Basin (MVGB) located within the Tertiary volcanics of the central Sierra Nevada, CA. Activities of ³⁵S were measured in dissolved sulfate (³⁵SO₄^2?) in SCB and MVGB snowpack, groundwater, springs, and streamflow. The percent of new snowmelt (PNS) in SCB streamflow ranged from 0.2 ± 6.6% during baseflow conditions to 14.0 ± 3.4% during high‐flow periods of snowmelt. Similar to SCB, the PNS in MVGB groundwater and streamflow was typically <30% with the largest fractions occurring in late spring or early summer following peak streamflow. The consistently low PNS suggests that a significant fraction of annual snowmelt in SCB and MVGB recharges groundwater, and groundwater contributions to streamflow in these systems have the potential to mitigate climate change impacts on runoff.     

Identifying the flow pattern and natural recharge at a strategic groundwater resource in the Dornogobi Province,Mongolia

The strategic project of economic development in the Dornogobi Province in Mongolia is dependent on water supply. Thus a comprehensive hydrogeological characterization was focused on the Upper Cretaceous multi-aquifer system north of Sainshand city. A conceptual model was developed to discover the groundwater flow pattern essential to correct the setting of the numerical model of groundwater flow created using MODFLOW to assess the natural recharge of the aquifer. The conceptualization was based on geological and hydrogeological characterization. However, the evaluation of hydrochemistry proved to be the key factor revealing the principal feature of the groundwater flow pattern, which is the presence of preferential flow zones. These zones allow for intensive transfer of relatively fresh Na(Mg,Ca)?HCO₃-dominated groundwater into discharge areas, where it leaks into the Quaternary aquifer. The numerical model suggested an enormous natural recharge of 22 100 m³/d, originating in 64% of the preferential flow zones.     

High-frequency monitoring of the occurrence of preferential flow on hillslopes and its relationship with rainfall features,soil moisture and landscape

ABSTRACT

High-frequency monitoring was conducted to quantify the frequency and controlling factors of preferential flow (PF) in a monsoon-influenced sub-humid mountainous catchment (6.48 km²) of Northern China. Rainfall was measured using nine bucket raingauges. Soil moisture probes were set up at 12 sites to observe the PF. Overall, 129 rainfall events were identified during the years 2014–2016. The average PF occurrence was 41%, which increased to 71% during heavy rainfall events (>20 mm) revealing a strong influence of the amount and intensity of rainfall. The study also revealed that the PF increased with antecedent soil moisture. Soil moisture was much higher on flat sites compared to sloping sites, providing evidence that the topography has a strong influence on rainfall infiltration and runoff which, subsequently, influence soil moisture variation and the occurrence of PF. Our findings provide valuable insights into the hydrological processes for studies in regions with similar environmental conditions.     

Subsurface storm flow formation at different hillslopes and implications for the ‘old water paradox’

Although many studies over the past several decades have documented the importance of subsurface stormflow (SSF) in hillslopes, its formation is still not well understood. Therefore, we studied SSF formation in the vadose soil zone at four different hillslopes during controlled sprinkling experiments and natural rainfall events. Event and pre‐event water fractions were determined using artificially traced sprinkling water and ²²²Rn as natural tracer. SSF formation and the fraction of pre‐event water varied substantially at different hillslopes. Both intensity of SSF and fraction of pre‐event water depended on whether SSF in preferential flow paths was fed directly from precipitation or was fed indirectly from saturated parts of the soil. Soil water was rapidly mobilized from saturated patches in the soil matrix and was subsequently released into larger pores, where it mixed with event water. Substantial amounts of pre‐event water, therefore, were contained in fast flow components like subsurface storm flow and also in overland flow. This finding has consequences for commonly used hydrograph separation methods and might explain part of the ‘old water paradox’. Copyright © 2007 John Wiley & Sons, Ltd.     

Evaporation from seasonally frozen bare and vegetated ground at various groundwater table depths in the Ordos Basin,Northwest China

In cold climates, the process of freezing–thawing significantly affects the ground surface heat balance and water balance. To better understand the mechanism of evaporation from seasonally frozen soils, we performed field experiments at different water table depths on vegetated and bare ground in a semiarid region in China. Soil moisture and temperature, air temperature, precipitation, and water table depths were measured over a 5‐month period (November 1, 2016, to March 14, 2017). The evaporation, which was calculated by a mass balance method, was high in the periods of thawing and low in the periods of freezing. Increased water table depth in the freezing period led to high soil moisture in the upper soil layer, whereas lower initial groundwater levels during freezing–thawing decreased the cumulative evaporation. The extent of evaporation from the bare ground was the same in summer as in winter. These results indicate that a noteworthy amount of evaporation from the bare ground is present during freezing–thawing. Finally, the roots of Salix psammophila could increase the soil temperature. This study presents an insight into the joint effects of soil moisture, temperature, ground vegetation, and water table depths on the evaporation from seasonally frozen soils. Furthermore, it also has important implications for water management in seasonally frozen areas.     

Differences in gravity‐dominated unsaturated flow during autumn rains and snowmelt

Five tracer experiments have been performed in a coarse‐textured soil near the new main airport at Gardermoen, Norway. In two lysimeter walls, 30 and 40 measuring points form the basis for spatial moment calculations. Although experiments were performed under different meteorological conditions (autumn and snowmelt) and at two different sites, the ratios of centres of vertical mass over cumulative infiltration were of the same order of magnitude, indicating a gravity‐dominated flow. Two‐dimensional transport simulations with SUTRA (Voss, 1984), with a priori estimated input parameters and random fields of soil hydraulic properties revealed a relatively good agreement with the experimental results. Three possible sources of heterogeneity affecting the vertical displacement of solute during snowmelt were identified: variability of soil physical properties, soil surface elevations and variability of ground frost during the melting period. To obtain accurate predictions, soil heterogeneity was the most important factor to characterize for the coarse‐textured soil under consideration. Copyright © 1999 John Wiley & Sons, Ltd.     

A physically based distributed subsurface–surface flow dynamics model for forested mountainous catchments

This study was designed to develop a physically based hydrological model to describe the hydrological processes within forested mountainous river basins. The model describes the relationships between hydrological fluxes and catchment characteristics that are influenced by topography and land cover. Hydrological processes representative of temperate basins in steep terrain that are incorporated in the model include intercepted rainfall, evaporation, transpiration, infiltration into macropores, partitioning between preferential flow and soil matrix flow, percolation, capillary rise, surface flow (saturation‐excess and return flow), subsurface flow (preferential subsurface flow and baseflow) and spatial water‐table dynamics. The soil–vegetation–atmosphere transfer scheme used was the single‐layer Penman–Monteith model, although a two‐layer model was also provided. The catchment characteristics include topography (elevation, topographic indices), slope and contributing area, where a digital elevation model provided flow direction on the steepest gradient flow path. The hydrological fluxes and catchment characteristics are modelled based on the variable source‐area concept, which defines the dynamics of the watershed response. Flow generated on land for each sub‐basin is routed to the river channel by a kinematic wave model. In the river channel, the combined flows from sub‐basins are routed by the Muskingum–Cunge model to the river outlet; these comprise inputs to the river downstream. The model was applied to the Hikimi river basin in Japan. Spatial decadal values of the normalized difference vegetation index and leaf area index were used for the yearly simulations. Results were satisfactory, as indicated by model efficiency criteria, and analysis showed that the rainfall input is not representative of the orographic lifting induced rainfall in the mountainous Hikimi river basin. Also, a simple representation of the effects of preferential flow within the soil matrix flow has a slight significance for soil moisture status, but is insignificant for river flow estimations. Copyright © 2005 John Wiley & Sons, Ltd.     

Land cover effects on soil infiltration capacity measured using plot scale rainfall simulation in steep tropical lowlands of Central Panama

Understanding land use/land cover (LULC) effects on tropical soil infiltration is crucial for maximizing watershed scale hydro-ecosystem services and informing land managers. This paper reports results from a multiyear investigation of LULC effects on soil bulk infiltration in steep, humid tropical, and lowland catchments. A rainfall simulator applied water at measured rates on 2 × 6 m plots producing infiltration through structured, granulated, and macroporous Ferralsols in Panama's central lowlands. Time-lapse electrical resistivity tomography (ERT) helped to visualize infiltration depth and bulk velocity. A space-for-time substitution methodology allowed a land-use history investigation by considering the following: (a) a continuously heavy-grazed cattle pasture, (b) a rotationally grazed traditional cattle pasture, (c) a 4-year-old (y.o.) silvopastoral system with nonnative improved pasture grasses and managed intensive rotational grazing, (d) a 7 y.o. teak (Tectona grandis) plantation, (e) an approximately 10 y.o. secondary succession forest, (f) a 12 y.o. coffee plantation (Coffea canephora), (g) an approximately 30 y.o. secondary succession forest, and (h) a >100 y.o. secondary succession forest. Within a land cover, unique plot sites totalled two at (a), (c), (d), (e), and (g); three at (b); and one at (f) and (h). Our observations confirmed measured infiltration scale dependency by comparing our 12 m² plot-scale measurements against 8.9 cm diameter core-scale measurements collected by others from nearby sites. Preferential flow pathways (PFPs) significantly increased soil infiltration capacity, particularly in forests greater than or equal to 10 y.o. Time-lapse ERT observations revealed shallower rapid bulk infiltration and increased rapid lateral subsurface flow in pasture land covers when compared with forest land covers and highlighted how much subsurface flow pathways can vary within the Ferralsol soil class. Results suggest that LULC effects on PFPs are the dominant mechanism by which LULC affects throughfall partitioning, runoff generation, and flow pathways.     

Subsurface lateral flow generation in aspen and conifer‐dominated hillslopes of a first order catchment in northern Utah

Mountain headwater catchments in the semi‐arid Intermountain West are important sources of surface water because these high elevations receive more precipitation than neighboring lowlands. This study examined subsurface runoff in two hillslopes, one aspen dominated, the other conifer dominated, adjacent to a first order stream in snow‐driven northern Utah. Snow accumulation, soil moisture, trenchflow and streamflow were examined in hillslopes and their adjacent stream. Snow water equivalents (SWEs) were greater under aspen stands compared to conifer, the difference increasing with higher annual precipitation. Semi‐variograms of shallow spatial soil moisture patterns and transects of continuous soil moisture showed no increase in soil moisture downslope, suggesting the absence of subsurface flow in shallow (～12 cm) soil layers of either vegetation type. However, a clear threshold relationship between soil moisture and streamflow indicated hillslope–stream connectivity, deeper within the soil profile. Subsurface flow was detected at ～50 cm depth, which was sustained for longer in the conifer hillslope. Soil profiles under the two vegetation types varied, with deep aspen soils having greater water storage capacity than shallow rocky conifer soils. Though SWEs were less under the conifers, the soil profile had less water storage capacity and produced more subsurface lateral flow during the spring snowmelt. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

Soil moisture states,lateral flow,and streamflow generation in a semi‐arid,snowmelt‐driven catchment

Hydraulic connectivity on hillslopes and the existence of preferred soil moisture states in a catchment have important controls on runoff generation. In this study we investigate the relationships between soil moisture patterns, lateral hillslope flow, and streamflow generation in a semi‐arid, snowmelt‐driven catchment. We identify five soil moisture conditions that occur during a year and present a conceptual model based on field studies and computer simulations of how streamflow is generated with respect to the soil moisture conditions. The five soil moisture conditions are (1) a summer dry period, (2) a transitional fall wetting period, (3) a winter wet, low‐flux period, (4) a spring wet, high‐flux period, and (5) a transitional late‐spring drying period. Transitions between the periods are driven by changes in the water balance between rain, snow, snowmelt and evapotranspiration. Low rates of water input to the soil during the winter allow dry soil regions to persist at the soil–bedrock interface, which act as barriers to lateral flow. Once the dry‐soil flow barriers are wetted, whole‐slope hydraulic connectivity is established, lateral flow can occur, and upland soils are in direct connection with the near‐stream soil moisture. This whole‐slope connectivity can alter near‐stream hydraulics and modify the delivery of water, pressure, and solutes to the stream. Copyright © 2005 John Wiley & Sons, Ltd.     

Spatio‐temporal controls of snowmelt and runoff generation during rain‐on‐snow events in a mid‐latitude mountain catchment

A network of 30 standalone snow monitoring stations was used to investigate the snow cover distribution, snowmelt dynamics, and runoff generation during two rain‐on‐snow (ROS) events in a 40 km² montane catchment in the Black Forest region of southwestern Germany. A multiple linear regression analysis using elevation, aspect, and land cover as predictors for the snow water equivalent (SWE) distribution within the catchment was applied on an hourly basis for two significant ROS flood events that occurred in December 2012. The available snowmelt water, liquid precipitation, as well as the total retention storage of the snow cover were considered in order to estimate the amount of water potentially available for the runoff generation. The study provides a spatially and temporally distributed picture of how the two observed ROS floods developed in the catchment. It became evident that the retention capacity of the snow cover is a crucial mechanism during ROS. It took several hours before water was released from the snowpack during the first ROS event, while retention storage was exceeded within 1 h from the start of the second event. Elevation was the most important terrain feature. South‐facing terrain contributed more water for runoff than north‐facing slopes, and only slightly more runoff was generated at open compared to forested areas. The results highlight the importance of snowmelt together with liquid precipitation for the generation of flood runoff during ROS and the large temporal and spatial variability of the relevant processes. Copyright © 2015 John Wiley & Sons, Ltd.     

On the controls of preferential flow in soils of different hillslope position and lithological origin

Soils derived from different lithologies and their controls on preferential flow remain underexplored in forested landscapes. In the same lithology, the propensity for preferential flow occurrence at different hillslope positions also remains largely elusive. By utilizing a soil moisture response time method, we compared preferential flow occurrence between a shale site (Shale Hills, silt loam soils) and a sandstone site (Garner Run, sandy loam soils) at four hillslope positions: ridge-top, North- and South-facing mid-slopes and toe slope, for over 2 years. The catchments are neighbouring and covered by temperate forest. For the four hillslope positions, Shale Hills had higher preferential flow frequencies compared to Garner Run. Between these two catchments, the South-facing mid-slope sites showed the highest contrasts in preferential flow frequency (33.5% of events at Shale Hills vs. 8.8% at Garner Run) while the ridge-top sites showed the lowest contrasts (18.7 vs. 13.2%). Additionally, over the unfrozen period, for seven out of eight monitoring sites, drier antecedent conditions tended to be more favourable for preferential flows to occur, with significant (p < .01) relationships at two sites. Except for the South-facing mid-slope sites, both Shale Hills and Garner Run had two preferential flow pathways. The characteristic preferential flow pathways at Shale Hills were the B_w and C horizons, and for Garner Run, preferential flow moved from the E/AE horizon to the B_w horizon. This study shows that shale-derived soils tended to have higher preferential flow occurrence than sandstone soils, but hillslope positions exhibit different levels of contrasts. More effort should be paid to study the impact of lithology on preferential flows in the context of land surface modelling and biogeochemical reactions to improve ecosystem services of headwater catchments.     

Tracer‐based evidence of heterogeneity in subsurface flow and storage within a boreal hillslope

Runoff from boreal hillslopes is often affected by distinct soil boundaries, including the frozen boundary and the organic‐mineral boundary (OMB), where highly porous and hydraulically conductive organic material overlies fine‐grained mineral soils. Viewed from the surface, ground cover appears as a patchwork on sub‐meter scales, with thick, moss mats interspersed with lichen‐covered, silty soils with gravel inclusions. We conducted a decameter‐scale subsurface tracer test on a boreal forest hillslope in interior Alaska to quantify locations and mechanisms of transport and storage in these soils, focusing on the OMB. A sodium bromide tracer was added as a slug addition to a pit and sampled at 40 down‐gradient wells, screened primarily at the OMB and within a 7 × 12 m well field. We maintained an elevated head in the injection pit for 8.5 hr to simulate a storm. Tracer breakthrough velocities ranged from <0.12 to 0.93 m hr^?1, with the highest velocities in lichen‐covered soils. After 12 hr and cessation of the elevated head, the tracer coalesced and was only detected in thick mosses at a trough in the OMB. By 24 hr, approximately 17% of the tracer mass could be accounted for. The majority of the mass loss occurred between 4 and 12 hr, while the tracer was in contact with lichen‐covered soils, which is consistent with tracer transport into deeper flow paths via preferential flow through discrete gravelly areas. Slow breakthroughs suggest that storage and exchange also occurred in shallow soils, likely related to saturation and drainage in fine‐grained mineral soils caused by the elevated hydraulic head. These findings highlight the complex nature of storage and transmission of water and solutes from boreal hillslopes to streams and are particularly relevant given rapid changes to boreal environments related to climate change, thawing permafrost and increasing fire severity.     

Groundwater recharge and flow in a small mountain catchment in northern Ethiopia

Abstract

The hydrodynamic behaviour of a sloped phreatic aquifer in the Tigray Highlands in northern Ethiopia is described. The aquifer is situated in the soils of a plateau on top of a basalt sequence and lies on steep slopes; the latter lead to hydraulic gradients that can cause high discharge fluxes. Distinct wet and dry seasons characterize the climate of the Tigray Highlands and recharge is absent during the dry season. Because of the fertile vertisols that have developed, the plateau is heavily cultivated and thus has great local economic, and hence social, importance. Water for land irrigation is almost exclusively delivered by rainfall, which is largely restricted to the period June—September. During the dry season, the water table drops dramatically and the aquifer drains nearly completely, under the strong gravity-driven, sustained discharges. This study strives to give insights into recharge and discharge mechanisms of the aquifer, in order to improve the effectiveness of the implemented water conservation measures.