Search for: All records

Analytical solutions for the three‐dimensional groundwater flow equation have been widely used to gain insight about subsurface flow structure and as an alternative to computationally expensive numerical models. Of particular interest are solutions that decompose prescribed hydraulic head boundaries (e.g., Dirichlet boundary condition) into a collection of harmonic functions. Previous studies estimate the frequencies and amplitudes of these harmonics with a least‐square approach where the amplitudes are fitted given a pre‐assigned set of frequencies. In these studies, an ad hoc and structured discretization of the frequency domain is typically used, excluding dominant frequencies while assigning importance to spurious frequencies, with significant consequences for estimating the fluxes and residence times. This study demonstrates the advantages of using a pre‐assigned frequency spectrum that targets the dominant frequencies based on rigorous statistical analysis with predefined significance levels. The new approach is tested for three hydrologic conceptualizations: (a) a synthetic periodic basin, (b) synthetic bedforms, and (c) a natural mountainous watershed. The performance of the frequency spectrum selection is compared with exact analytical or approximate numerical solutions. We found that the new approach better describes the fluxes and residence times for Dirichlet boundaries with well‐defined characteristics spatial scales (e.g., periodic basins and bedforms). For more complex scenarios, such as natural mountainous watersheds, both pre‐assigned frequency spectrums present similar performance. The spectral solutions presented here can play a central role in developing reduced‐complexity models for assessing regional water and solute fluxes within mountain watersheds and hyporheic zones.

 

Storm direction modulates a hydrograph's magnitude and duration, thus having a potentially large effect on local flood risk. However, how changes in the preferential storm direction affect the probability distribution of peak flows remains unknown. We address this question with a novel Monte Carlo approach where stochastically transposed storms drive hydrologic simulations over medium and mesoscale watersheds in the Midwestern United States. Systematic rotations of these watersheds are used to emulate changes in the preferential storm direction. We found that the peak flow distribution impacts are scale‐dependent, with larger changes observed in the mesoscale watershed than in the medium‐scale watershed. We attribute this to the high diversity of storm patterns and the storms' scale relative to watershed size. This study highlights the potential of the proposed stochastic framework to address fundamental questions about hydrologic extremes when our ability to observe these events in nature is hindered by technical constraints and short time records.

 

Hyporheic exchange is a crucial control of the type and rates of streambed biogeochemical processes, including metabolism, respiration, nutrient turnover, and the transformation of pollutants. Previous work has shown that increasing discharge during an individual peak flow event strengthens biogeochemical turnover by enhancing the exchange of water and dissolved solutes. However, due to the nonsteady nature of the exchange process, successive peak flow events do not exhibit proportional variations in residence time and turnover, and in some cases, can reduce the hyporheic zones' biogeochemical potential. Here, we used a process‐based model to explore the role of successive peak flow events on the flow and transport characteristics of bedform‐induced hyporheic exchange. We conducted a systematic analysis of the impacts of the events' magnitude, duration, and time between peaks in the hyporheic zone's fluxes, penetration, and residence times. The relative contribution of each event to the transport of solutes across the sediment‐water interface was inferred from transport simulations of a conservative solute. In addition to temporal variations in the hyporheic flow field, our results demonstrate that the separation between two events determines the temporal evolution of residence time and that event time lags longer than the memory of the system result in successive events that can be treated independently. This study highlights the importance of discharge variability in the dynamics of hyporheic exchange and its potential implications for biogeochemical transformations and fate of contaminants along river corridors.

 

Coupled groundwater flow and heat transport within hyporheic zones extensively affect water, energy, and solute exchange with surrounding sediments. The local and cumulative implications of this tightly coupled process strongly depend on characteristics of drivers (i.e., discharge and temperature of the water column) and modulators (i.e., hydraulic and thermal properties of the sediment). With this in mind, we perform a systematic numerical analysis of hyporheic responses to understand how the temporal variability of river discharge and temperature affect flow and heat transport within hyporheic zones. We identify typical time series of river discharge and temperature from gauging stations along the headwater region of Mississippi River Basin, which are characterized by different degrees of flow alteration, to drive a physics‐based model of the hyporheic exchange process. Our modeling results indicate that coupled groundwater flow and heat transport significantly affects the dynamic response of hyporheic zones, resulting in substantial differences in exchange rates and characteristic time scales of hyporheic exchange processes. We also find that the hyporheic zone dampens river temperature fluctuations increasingly with higher frequency of temperature fluctuations. This dampening effect depends on the system transport time scale and characteristics of river discharge and temperature variability. Furthermore, our results reveal that the flow alteration reduces the potential of hyporheic zones to act as a temperature buffer and hinders denitrification within hyporheic zones. These results have significant implications for understanding the drivers of local variability in hyporheic exchange and the implications for the development of thermal refugia and ecosystem functioning in hyporheic zones.

 

Over the past 35 years the Buckman wellfield near Santa Fe, New Mexico, experienced production well drawdowns in excess of 180 m, resulting in ground subsidence and surface cracks. Increased reliance on surface water diversions since 2011 has reduced pumping and yielded water level recovery. To characterize the impact of wellfield management decisions on the aquifer system, we reconstruct the surface deformation history through the European Remote Sensing Satellite, Advanced Land Observing Satellite, and Sentinel‐1 Interferometric Synthetic Aperture Radar (InSAR) time series analysis during episodes of drawdown (1993–2000), recovery (2007–2010), and modern management (2015–2018) in discontinuous observations over a 25‐year period. The observed deformation generally reflects changes in hydraulic head. However, at times during the wellfield recovery, the deformation signal is complex, with patterns of uplift and subsidence suggesting a compartmentalized aquifer system. Recent records of locally high geothermal gradients and an overall warming of the system (~0.5°C during the water level recovery) obtained from repeat temperature measurements between 2013 and 2018 constrain a conceptual model of convective heat transfer that requires a vertical permeable zone near an observed fault. To reproduce observed temperature patterns at monitoring wells, high basal heat flow and convective cooling associated with downward flow of water from cool shallow aquifers during the drawdown period is necessary. The fault, however, appears to die out southward or may be locally permeable, as conceptual cross‐sectional hydrologic modeling reproduces the surface deformation without such a structure. Our work demonstrates the importance of incorporating well‐constrained stratigraphy and structure when modeling near‐surface deformation induced by, for instance, groundwater production.

 

Physics‐based distributed hydrological models that include groundwater are widely used to understand and predict physical and biogeochemical processes within watersheds. Typically, due to computational limitations, watershed modelers minimize the number of elements used in domain discretization, smoothing or even ignoring critical topographic features. We use an idealized model to investigate the implications of mesh refinement along streams and ridges for modeling three‐dimensional groundwater flow and transport in mountainous watersheds. For varying degrees of topographic complexity level (TCL), which increases with the level of mesh refinement, and geological heterogeneity, we estimate and compare steady state baseflow discharge, mean age, and concentration of subsurface weathering products. Results show that ignoring lower‐order streams or ridges diminishes flow through local flow paths and biases higher the contribution of intermediate and regional flow paths, and biases baseflow older. The magnitude of the bias increases for systems where permeability rapidly decreases with depth and is dominated by shallow flow paths. Based on a simple geochemical model, the concentration of weathering products is less sensitive to the TCL, partially due to the thermodynamic constraints on chemical reactions. Our idealized model also reproduces the observed emergent scaling relationship between the groundwater contribution to streamflow and drainage area, and finds that this scaling relationship is not sensitive to mesh TCL. The bias effects have important implications for the use of hydrological models in the interpretation of environmental tracer data and the prediction of biogeochemical evolution of stream water in mountainous watersheds.

 

The Critical Zone encompasses the biosphere and its heterogeneities, with an extremely high differentiation of properties and processes within each compartment from bedrock to canopy, and across terrestrial and aquatic interfaces. Given this complexity, a comprehensive areal characterization of the critical zone environment at multiple temporal resolutions is needed but not always possible, and failing which the ecosystem fluxes, exchange rates and biogeochemical functioning may be under- or over-predicted. The hot spots hot moments (HSHMs) concept provides an opportunity to identify the dominant controls on carbon, nutrients, water and energy exchanges. Hot spots are regions or sites that show disproportionately high reaction rates relative to surrounding area, while hot moments are defined as times that show disproportionately high reaction rates relative to longer intervening time periods (McClain et al. 2003). 

Biogeochemical processes are often spatially discrete (hot spots) and temporally isolated (hot moments) due to variability in controlling factors like hydrologic fluxes, lithological characteristics, bio-geomorphic features, and external forcing. Although these hot spots and hot moments (HSHMs) account for a high percentage of carbon, nitrogen and nutrient cycling within the Critical Zone, the ability to identify and incorporate them into reactive transport models remains a significant challenge. This chapter provides an overview of the hot spots hot moments (HSHMs) concepts, where past work has largely focused on carbon and nitrogen dynamics within riverine systems. This work is summarized in the context of process-based and data-driven modeling approaches, including a brief description of recent research that casts a wider net to incorporate Hg, Fe and other Critical Zone elements, and focuses on interdisciplinary approaches and concepts. The broader goal of this chapter is to provide an overview of the gaps in our current understanding of HSHMs, and the opportunities therein, while specifically focusing on the underlying parameters and processes leading to their prognostic and diagnostic representation in reactive transport models. 

Abstract. Groundwater table dynamics extensively modify the volume of the hyporheic zoneand the rate of hyporheic exchange processes. Understanding the effects ofdaily groundwater table fluctuations on the tightly coupled flow and heattransport within hyporheic zones is crucial for water resourcesmanagement. With this aim in mind, a physically based model is used to explorehyporheic responses to varying groundwater table fluctuationscenarios. The effects of different timing and amplitude of groundwater tabledaily drawdowns under gaining and losing conditions are explored in hyporheiczones influenced by natural flood events and diel river temperaturefluctuations. We find that both diel river temperature fluctuations and dailygroundwater table drawdowns play important roles in determining thespatiotemporal variability of hyporheic exchange rates, temperature ofexfiltrating hyporheic fluxes, mean residence times, and hyporheicdenitrification potentials. Groundwater table dynamics present substantiallydistinct impacts on hyporheic exchange under gaining or losing conditions. Thetiming of groundwater table drawdown has a direct influence on hyporheicexchange rates and hyporheic buffering capacity on thermaldisturbances. Consequently, the selection of aquifer pumping regimes hassignificant impacts on the dispersal of pollutants in the aquifer and thermalheterogeneity in the sediment.