Subsurface Water Flow Paths Research Articles

Subsurface flow paths in a chronosequence of calcareous soils: impact of soil age and rainfall intensities on preferential flow occurrence

Abstract. Soil hydrologic processes play an important role in the hydro-pedo-geomorphological feedback cycle of landscape evolution. Soil properties and subsurface flow paths both change over time, but due to a lack of observations, subsurface water flow paths are often not properly represented in soil and landscape evolution models. We investigated the evolution of subsurface flow paths across a soil chronosequence in the calcareous glacier forefield at the Griessfirn glacier in the Swiss Alps. Young soils developed from calcareous parent material usually have a high pH value, which likely affects vegetation development and pedogenesis and thus the evolution of subsurface flow paths. We chose four glacial moraines of different ages (110, 160, 4 900, and 13 500 years) and conducted sprinkling experiments with the dye tracer Brilliant Blue on three plots at each moraine. Each plot was divided into three equal subplots, and dyed water was applied with three different irrigation intensities (20, 40, and 60 mm h−1) and an irrigation amount of 40 mm. Subsequent excavation of soil profiles enabled the tracing of subsurface flow paths. A change in flow types with increasing moraine age was observed from a rather homogeneous matrix flow at 110 and 160 years to heterogeneous matrix and finger-shaped flow at 4 900 and 13 500 years. However, the proportion of preferential flow paths is not necessarily directly related to the moraine age but rather to soil properties such as texture, soil layering, organic matter content, and vegetation characteristics such as root length density and biomass. Irrigation intensity had an effect on the number of finger-shaped flow paths at the two old moraines. We also found that flow paths in this calcareous material evolved differently compared to a previous study in siliceous material, which emphasizes the importance of parent material for flow path evolution. Our study provides a rare systematic dataset and observations on the evolution of vertical subsurface flow paths in calcareous soils, which is useful to improve their representation in the context of landscape evolution modeling.

Topology and spatial-pressure-distribution reconstruction of an englacial channel

Abstract. Information about glacier hydrology is important for understanding glacier and ice sheet dynamics. However, our knowledge about water pathways and pressure remains limited, as in situ observations are sparse and methods for direct area-wide observations are limited due to the extreme and hard-to-access nature of the environment. In this paper, we present a method that allows for in situ data collection in englacial channels using sensing drifters. Furthermore, we demonstrate a model that takes the collected data and reconstructs the planar subsurface water flow paths providing spatial reference to the continuous water pressure measurements. We showcase this method by reconstructing the 2D topology and the water pressure distribution of a free-flowing englacial channel in Austre Brøggerbreen (Svalbard). The approach uses inertial measurements from submersible sensing drifters and reconstructs the water flow path between given start and end coordinates. Validation of the method was done on a separate supraglacial channel, showing an average error of 3.9 m and the total channel length error of 29 m (6.5 %). At the englacial channel, the average error is 12.1 m; the length error is 107 m (11.6 %); and the water pressure standard deviation is 3.4 hPa (0.3 %). Our method allows for mapping of subsurface water flow paths and spatially referencing the pressure distribution within. Further, our method would be extendable to the reconstruction of other, previously underexplored subsurface fluid flow paths such as pipelines or karst caves.

A new approach for estimating and delineating within-field crop water stress zones with satellite imagery

ABSTRACT As water scarcity and demand for food both grow, viable solutions to efficiently irrigate crops are needed. The recent advent of variable rate irrigation systems provides technology for growers to water their fields accounting for the spatially heterogeneous demand of their crops for water. Some zones within a field may require less or more water based on varying solar radiation, soil type, or surface or subsurface water flow paths. Since these processes are heterogeneous and complex, they are challenging to model or measure with the spatial resolution necessary to implement successful variable rate irrigation. In this paper, we present a new approach that uses high-resolution (20 to 100 m) satellite imagery obtained from Sentinel-2 and Landsat 8 to estimate zones of crop water stress. Our approach integrates two independent measures of crop water availability – Normalized Difference Water Index, and Evapotranspiration modelled using remotely sensed land surface temperature – with k-means clustering to delineate zones in an automated fashion. Our study area in southern Idaho has been set up experimentally to grow crops that are water-stressed and not water-stressed, thus providing a rigorous design to test and validate our method. The results of our method agree well with our independent suite of crop water stress measures, such as yield, variable irrigation rates, and matric potential. This method can be used to develop data-driven zones for variable rate irrigation, which has the potential to improve the efficiency of water use without sacrificing yield.

Tracing silicate weathering processes in the permafrost-dominated Lena River watershed using lithium isotopes

Increasing global temperatures are causing widespread changes in the Arctic, including permafrost thawing and altered freshwater inputs and trace metal and carbon fluxes into the ocean and atmosphere. Changes in the permafrost active layer thickness can affect subsurface water flow-paths and water-rock interaction times, and hence weathering processes. Riverine lithium isotope ratios (reported as δ7Li) are tracers of silicate weathering that are unaffected by biological uptake, redox, carbonate weathering and primary lithology. Here we use Li isotopes to examine silicate weathering processes in one of the largest Russian Arctic rivers: the Lena River in eastern Siberia. The Lena River watershed is a large multi-lithological catchment, underlain by continuous permafrost. An extensive dataset of dissolved Li isotopic compositions of waters from the Lena River main channel, two main tributaries (the Aldan and Viliui Rivers) and a range of smaller sub-tributaries are presented from the post-spring flood/early-summer period at the onset of active layer development and enhanced water-rock interactions. The Lena River main channel (average δ7Lidiss ∼ 19‰) has a slightly lower isotopic composition than the mean global average of 23‰ (Huh et al., 1998a). The greatest range of [Li] and δ7Lidiss are observed in catchments draining the south-facing slopes of the Verkhoyansk Mountain Range. South-facing slopes in high-latitude, permafrost-dominated regions are typically characterised by increased summer insolation and higher daytime temperatures relative to other slope aspects. The increased solar radiation on south-facing catchments promotes repeated freeze-thaw cycles, and contributes to more rapid melting of snow cover, warmer soils, and increased active layer thaw depths. The greater variability in δ7Li and [Li] in the south-facing rivers likely reflect the greater infiltration of melt water and enhanced water-rock interactions within the active layer.A similar magnitude of isotopic fractionation is observed between the low-lying regions of the Central Siberian Plateau (and catchments draining into the Viliui River), and catchments draining the Verkhoyansk Mountain Range into the Aldan River. This is in contrast to global rivers in non-permafrost terrains that drain high elevations or areas of rapid uplift, where high degrees of physical erosion promote dissolution of freshly exposed primary rock typically yielding low δ7Lidiss, and low-lying regions exhibit high riverine δ7Li values resulting from greater water-rock interaction and formation of secondary mineral that fractionates Li isotopes. Overall, the range of Li concentrations and δ7Lidiss observed within the Lena River catchment are comparable to global rivers located in temperate and tropical regions. This suggests that cryogenic weathering features specific to permafrost regions (such as the continual exposure of fresh primary minerals due to seasonal freeze-thaw cycles, frost shattering and salt weathering), and climate (temperature and runoff), are not a dominant control on δ7Li variations. Despite vastly different climatic and weathering regimes, the same range of riverine δ7Li values globally suggests that the same processes govern Li geochemistry – that is, the balance between primary silicate mineral dissolution and the formation (or exchange with) secondary minerals. This has implications for the use of δ7Li as a palaeo-weathering tracer for interpreting changes in past weathering regimes.

Effect of hydrogeomorphic setting on calcareous fen hydrology

AbstractCalcareous fens are species‐rich peatlands that are dependent on minerotrophic water sources for wetland functioning, with current conceptual models suggesting the water source is ubiquitously groundwater upwelling. By quantifying the water balance and subsurface water flow paths and fluxes over 3 growing seasons for calcareous fens in 3 different hydrogeomorphic settings (Riparian, Trough, and Basin), we show evidence that challenges this conceptual model. The Riparian Fen received an order of magnitude more water inputs than the Trough or Basin Fens and was dominated by stream recharge inputs and groundwater outputs. Precipitation and evaporation dominated the water balance of the Trough Fen whereas only the Basin Fen received sizeable groundwater inputs. Indeed, subsurface water fluxes were low at all fens due to weak hydraulic gradients and low saturated hydraulic conductivity in some areas of each wetland, though variations in growing season precipitation led to subsurface flow reversals in all 3 fens. Our results demonstrate the importance of understanding landscape position, or hydrogeomorphic setting, on calcareous fen hydrology for improving conservation, management, and restoration efforts of these important ecosystems.

Enhancement of carbon and nitrogen removal by helophytes along subsurface water flowpaths receiving treated wastewater

Wastewater treatment plant (WWTP) effluents are sources of dissolved organic carbon (DOC) and inorganic nitrogen (DIN) to receiving streams, which can eventually become saturated by excess of DIN. Aquatic plants (i.e., helophytes) can modify subsurface water flowpaths as well as assimilate nutrients and enhance microbial activity in the rhizosphere, yet their ability to increase DIN transformation and removal in WWTP-influenced streams is poorly understood. We examined the influence of helophytes on DIN removal along subsurface water flowpaths and how this was associated with DOC removal and labile C availability. To do so, we used a set of 12 flow-through flumes fed with water from a WWTP effluent. The flumes contained solely sediments or sediments with helophytes. Presence of helophytes in the flumes enhanced both DIN and DOC removal. Experimental addition of a labile C source into the flumes resulted in a high removal of the added C within the first meter of the flumes. Yet, no concomitant increases in DIN removal were observed. Moreover, results from laboratory assays showed significant increases in the potential denitrifying enzyme activity of sediment biofilms from the flumes when labile C was added; suggesting denitrification was limited by C quality. Together these results suggest that lack of DIN removal response to the labile C addition in flumes was likely because potential increases in denitrification by biofilms from sediments were counterbalanced by high rates of mineralization of dissolved organic matter. Our results highlight that helophytes can enhance DIN removal in streams receiving inputs from WWTP effluents; and thus, they can become a relevant bioremediation tool in WWTP-influenced streams. However, results also suggest that the quality of DOC from the WWTP effluent can influence the N removal capacity of these systems.

Water Flow Path Characterization in Shallow Vadose Zone Using Tensiometers

In this project, we will present the findings of a study using Tensiometer systems designed to investigate the water flow path pattern in shallow vadose zone. The purpose of this paper is to evaluate water flow path in shallow vadose zone and to calculate the infiltration rate and hydraulic conductivity of a soil using Tensiometer. We have measured the subsurface water flow paths in sandy clay loam soil following infiltration experiment using Tensiometers. The matric potential and hydraulic conductivity measurements show that subsequent infiltration and water movement in unsaturated (vadose) zone are vertical, but it can have large lateral component under steady condition. This shows that water moves generally from high water content to lower water content region. Average pressure head for the percolation test conduction locations EB and HB was -30 and -80 cm respectively. Hysteresis produces another interesting situation when the soil is drained. We found that the wetter portion of the soil in vadose zone could be at a lower potential (head) than the dryer portions, resulting in lateral driving force for a preferential flow of water from the dryer to the wet soil. The infiltration rate for the 5 cm ponded water was calculated at 5.45 cm/hr. The infiltration rate curve shows that the rate of infiltration decreases with the time. When infiltration first starts, the wetting front is steep and very close to the surface. Similarly, due to the pressure head gradient, large value for infiltration is recorded. Under these conditions, we believe that the gradient in pressure head is responsible for the rapid movement of water into the dry soil. The effect of gravity is less on water during the initial stages of infiltration; however, it is more effective for preferential flow pattern. In the latter infiltration event, the wetting front has moved deeper into the soil. As a result, the pressure head gradient at the surface is much smaller and consequently has little effect. When the pressure head is equal zero, infiltration rate approaches almost to the lowest level. We also observed that even a minor change in soil-water pressure due to slope could change both direction and magnitude of water flux.

Hillslope runoff generation influenced by layered subsurface in a headwater catchment in Ore Mountains, Germany

The Ore Mountains were one of the important flood source areas for several heavy floods over the last years. Reducing damages caused by floods demands sufficient information on the runoff generation processes in the catchments. The aim of this study is to provide insights into prevailing flow pathways, retention times and threshold behavior of a representative hillslope catchment with layered subsurface in the Ore Mountains. The study site is a forested headwater with gneiss as bedrock. We used hydrometrical methods, soil temperature data and sprinkler experiments. Results indicate that the hydraulic anisotropic structure of the deepest layer in 0.9–1.7 m depth is the major controlling factor for subsurface water flow paths. On one hand, this layer acts as an aquitard for seeping water because of its high bulk density. On the other hand, water within the layer is able to flow laterally because of the sandy texture and coarse clasts oriented parallel to the slope. Moreover, three pre-moisture controlled types of runoff processes were addressed. With low antecedent soil moisture, saturation overland flow dominates in the spring bog. With intermediate or high pre-moisture, interflow is generated. The measured runoff coefficients increase in a nonlinear manner with rising pre-moisture. A soil water tension threshold value near field capacity is the tipping point for nonlinear runoff response. These findings emphasize the impact of the layered structure of the subsurface and of antecedent soil moisture for runoff generation in low mountain ranges and may be useful for establishing flood warning systems.

Hillslope permeability architecture controls on subsurface transit time distribution and flow paths

Summary Defining the catchment transit time distribution remains a challenge. Here, we used a new semi-analytical physically-based integrated subsurface flow and advective–dispersive particle movement model to assess the subsurface controls on subsurface water flow paths and transit time distributions. First, we tested the efficacy of the new model for simulation of the observed groundwater dynamics at the well-studied S-transect hillslope (Vastrabacken sub-catchment, Sweden). This system, like many others, is characterized by exponential decline in saturated hydraulic conductivity and porosity with soil depth. The model performed well relative to a tracer-based estimate of transit time distribution as well as observed groundwater depth–discharge relationship within 30 m of the stream. Second, we used the model to assess the effect of changes in the subsurface permeability architecture on flow pathlines and transit time distribution in a set of virtual experiments. Vertical patterns of saturated hydraulic conductivity and porosity with soil depth significantly influenced hillslope transit time distribution. Increasing infiltration rates significantly decreased mean groundwater age, but not the distribution of transit times relative to mean groundwater age. The location of hillslope hydrologic boundaries, including the groundwater divide and no-flow boundary underlying the hillslope, changed the transit time distribution less markedly. These results can guide future decisions on the degree of complexity that is warranted in a physically-based rainfall–runoff model to efficiently and explicitly estimate time invariant subsurface pathlines and transit time distribution.

Invited perspective: Why I am an optimist

I address a range of topics that provide the sociopolitical‐technological setting for my professional life. I discuss some influential features of post–World War II world geopolitics, landmark technological developments of that era, and the resulting follow‐up technologies that have made it possible to approach various problems in hydrology and water resources. I next address societal needs that have driven developments in hydrology and water resources engineering and follow with a discussion of the modern foundations of our science and what I think are the principal issues in hydrology. I pose three community challenges that when accomplished should advance hydrologic science: data network needs for improving the water budgets at all scales, characterizing subsurface water flow paths, and the information archiving and mining needs from instruments that will generate substantially richer data detail than have been used for most hydrologic work to the present. I then discuss several hydrologic and water resource risk‐based decision issues that matter to society to illustrate how such risks have been addressed successfully in the past. I conclude with a long‐term community “grand challenge,” the coupled modeling of the ocean‐atmosphere‐landform hydrologic cycle for the purpose of long–lead time hydrologic prediction.

Soil Moisture and Metolachlor Volatilization Observations over Three Years

A 3-yr study was conducted to focus on the impact of surface soil water content on metolachlor (2-chloro-N-(2-ethyl-6-methylphenyl)-N-(2-methoxy-1-methylethyl) acetamide) volatilization from a field with different surface soil water regimes created by subsurface water flow paths. Metolachlor vapor fluxes were measured at two locations within the field where local meteorological and soil conditions were relatively constant, except for surface soil water content, which differed significantly. Surface soil water content at the two sites differed in response to the presence of subsurface flow pathways. Detailed soil moisture observations over the duration of the study showed that for the first 2 yr (2004 and 2005), surface soil water contents at the dry location (V1) were nearly half those at the wetter location (V2). Cumulative metolachlor vapor fluxes during 2004 and 2005 at V1 were also about half that at V2. In the third year (2006), early-season drought conditions rendered the soil water content at the two locations to be nearly identical, resulting in similar metolachlor volatilization losses. Analysis of infrared soil surface temperatures suggests a correlation between surface soil temperatures and metolachlor volatilization when soils are wet (2004 and 2005) but not when the soils are dry (2006). Field-averaged metolachlor volatilization losses were highly correlated with increasing surface soil water contents (r(2) = 0.995).

The role of subsurface water flow paths on hillslope hydrological processes, landslides and landform development in steep mountains of Japan

AbstractHydrological monitoring was conducted in high‐relief watersheds in the Japan Alps to investigate the relationship between hillslope hydrological processes and landform evolution in steep granite and shale mountains. In the Koshibu watershed, underlain by Mesozoic shale, the drainage density and frequency was significantly lower than in the Yotagiri watershed underlain by granite. Drainage micro‐morphology analysis showed that hillslopes in the watersheds K1 and K6 (Koshibu basin) are mostly combinations of talus and bedrock exposures. In contrast, watershed Y1 (Yotagiri basin) is composed of several zero‐order streams with hollows. Infinite slope stability analysis indicates that the regolith shear strength in the K6 watershed (Koshibu basin) is lower than that of the Y1 hillslope, but groundwater levels were higher in the Y1 hillslope than in the K6 hillslope during storm events. These data suggest that, although the shear strength of the soil is stronger in the Yotagiri watershed, the slopes are unstable because of the groundwater conditions, whereas deep‐seated landslides may occur episodically in the Koshibu watershed associated with extreme storms and very high antecedent soil moisture. These differences would strongly contribute to the different observed hillslope processes and drainage characteristics. Copyright © 2004 John Wiley & Sons, Ltd.

Hydrology of an impounded lotic wetland—subsurface hydrology

Although the importance of hydrology to wetland ecosystem productivity is recognized, there are few detailed analyses of shallow subsurface water flow (<1 m depth) through wetland ecosystems. Shallow subsurface water flow paths and average linear velocity were quantified with nested lysimeters in a riverine wetland located in Hale County, AL. Water-table elevation was spatially and temporally variable and was directly correlated to precipitation at the wetland. In the beginning portion of the study (May to November 1996), average monthly precipitation was 73 mm and water-table elevation, decreased across much of the wetland. From December 1996 through July 1997, average monthly precipitation increased to 110 mm with a concomitant increase in water-table elevation throughout the wetland. During periods of low average monthly precipitation, wetland subsurface water was recharged from the hyporheic zone of the main inflow stream (northern side of the wetland). With increased precipitation, subsurface water recharge shifted from the stream hyporheos to recharge from an upgradient area (entire northern portion of the wetland) dominated by alder, Alnus serrulata. Horizontal and vertical average linear velocity was greatest in the upper 20 cm of the wetland compared to lower sediment depths. Vertical average linear velocity was similar among most sample intervals and was not correlated to antecedent precipitation. In laboratory mesocosm experiments, evapotranspiration by the emergent macrophyte Juncus effusus increased the rate of subsurface water flow in mesocosms with sand sediments, while there was no macrophyte effect in mesocosms with predominantly clay sediments.

Ba/Sr, Ca/Sr and 87Sr/ 86Sr ratios in soil water and groundwater: implications for relative contributions to stream water discharge

Barium/Sr and Ca/Sr ratios have been used to model the relative importance of different sources of stream water. Major and trace element concentrations together with 87Sr/ 86Sr ratios were measured in precipitation, soil water, groundwater and stream water in a small (9.4 km 2) catchment in northern Sweden. The study catchment is drained by a first order stream and mainly covered with podzolized Quaternary till of granitic composition. It is underlain by a 1.8 Ga granite. A model with mixing equations used in an iterative mode was developed in order to separate the stream water into 3 subsurface components: soil water, shallow groundwater, and deep groundwater. Contributions from precipitation are thus not included in the model. This source may be significant for the stream water generation, but it does not interfere with the calculations of the relative contributions from the subsurface components. The results show that the deep groundwater constitutes between 5 and 20% of the subsurface water discharge into the stream water. The highest values of the deep groundwater fraction occur during base flow. Soil water dominates during snowmelt seasons, whereas during base flow it is the least important fraction. Soil water accounts for 10–100% of the subsurface water discharge into the stream water. Shallow groundwater accounts for up to 80% of the subsurface water discharge with the lowest values at peak discharge during snowmelt seasons and the highest values during base flow. The validity of the model was tested by comparing the measured 87Sr/ 86Sr ratios in the stream water with the 87Sr/ 86Sr ratios predicted by the model. There was a systematic difference between the measured and modelled 87Sr/ 86Sr ratios which suggests that the fraction of soil water is overestimated by the model, especially during spring flood. As a consequence of this overestimation of soil water the amount of shallow groundwater is probably underestimated during this period. However, it is concluded that the differences between measured and predicted values are relatively small, and that element ratios are potentially effective tracers for different subsurface water flowpaths in catchments.

Tracing Subsurface Water Flowpaths by Means of Dissolved Ba/Sr, Ca/Sr and 87Sr/86Sr Ratios

Tracing Subsurface Water Flowpaths by Means of Dissolved Ba/Sr, Ca/Sr and 87Sr/86Sr Ratios