Abstract Mountain System Recharge processes are significant natural recharge pathways in many arid and semi‐arid mountainous regions. However, Mountain System Recharge processes are often poorly understood and characterized in hydrologic models. Mountains are the primary water supply source to valley aquifers via lateral groundwater flow from the mountain block (Mountain Block Recharge) and focused recharge from mountain streams contributing to focused Mountain Front Recharge at the piedmont zone. Here, we present a multi‐tool isogeochemical approach to characterize mountain flow paths and Mountain System Recharge in the northern Tulare Basin, California. We used groundwater chemistry data to delineate hydrochemical facies and explain the chemical evolution of groundwater from the Sierra Nevada to the Central Valley aquifer. Stable isotopes and radiogenic groundwater tracers validated Mountain System Recharge processes by differentiating focused from diffuse recharge, and estimating apparent groundwater age, respectively. Novel application of End‐Member Mixing Analysis using conservative chemical components revealed three Mountain System Recharge end‐members: (a) evaporated Ca‐HCO3water type associated with focused Mountain Front Recharge, (b) non‐evaporated Ca‐HCO3and Na‐HCO3water types with short residence times associated with shallow Mountain Block Recharge, and (c) Na‐HCO3groundwater type with long residence time associated with deep Mountain Block Recharge. We quantified the contribution of each Mountain System Recharge process to the valley aquifer by calculating mixing ratios. Our results show that deep Mountain Block Recharge is a significant recharge component, representing 31%–53% of the valley groundwater. Greater hydraulic connectivity between the Sierra Nevada and Central Valley has significant implications for parameterizing groundwater flow models. Our framework is useful for understanding Mountain System Recharge processes in other snow‐dominated mountain watersheds.
more »
« less
Groundwater Circulation Within the Mountain Block: Combining Flow and Transport Models With Magnetotelluric Observations to Untangle Its Nested Nature
Abstract Mountains are vital water sources for humans and ecosystems, continuously replenishing lowland aquifers through surface runoff and mountain recharge. Quantifying these fluxes and their relative importance is essential for sustainable water resource management. However, our mechanistic understanding of the flow and transport processes determining the connection between the mountain block and the basin aquifer remains limited. Traditional conceptualizations assume groundwater circulation within the mountain block is predominantly shallow. This view neglects the role of deep groundwater flowpaths significantly contributing to the water, solute, and energy budgets. Overcoming these limitations requires a holistic characterization of the multiscale nature of groundwater flow along the mountain‐to‐valley continuum. As a proof‐of‐concept, we use a coupled groundwater flow and transport model to design a series of numerical experiments that explore the role of geology, topography, and weathering rates in groundwater circulation and their resulting resistivity patterns. Our results show that accumulating solutes near stagnation zones create contrasting electrical resistivity patterns that separate local, intermediate, and regional flow cells, presenting a target for magnetotelluric observations. To demonstrate the sensitivity of magnetotelluric data to features in our resistivity models, we use the MARE2DEM electromagnetic modeling code to perform forward and inverse simulations. This study highlights the potential of magnetotelluric surveys to image the resistivity structure resulting from multiscale groundwater circulation through relatively impervious crystalline basement rocks in mountainous terrains. This capability could change our understanding of the critical zone, offering a holistic perspective that includes deep groundwater circulation and its role in conveying solutes and energy.
more »
« less
- Award ID(s):
- 1830172
- PAR ID:
- 10498295
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 60
- Issue:
- 4
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The Denver Basin Aquifer System (DBAS) is an important groundwater resource for Front Range communities and is currently experiencing increasing demand as populations grow and surface water supplies remain limited. It is necessary to better constrain aquifer recharge mechanisms to enable sustainable use of this resource. In other sedimentary basin aquifer systems, mountain front recharge has been shown to be a significant contributor to local basin groundwater recharge. In the DBAS, inputs from the mountain block are poorly understood, and previous numerical models have treated large segments of the mountain-front boundary as impermeable. However, there exist potential connections between the mountain block and the DBAS, either by direct contact of permeable units, which would facilitate underflow recharge into the basin, or by surface water infiltration to the aquifer units where they outcrop near the mountain front. To observe spatial and temporal relationships between mountain block water and DBAS water, we use water stable isotopes and characterize the d2H and d18O of monthly precipitation, seasonal surface waters, and groundwaters in and around the Front Range and Denver Basin. The goal of this study is to determine if differences in the isotopic composition of waters across the Front Range permit the use of d18O and d2H as tracers of water flow between Front Range streams and groundwater and the DBAS. We analyzed the unique signature of mountain-block water to compare with DBAS water stable isotope data collected from Castle Rock Water municipal wells. Stable isotope ratios varied spatially and temporally, with the greatest temporal variance observed in precipitation. Streams showed great spatial variance, and less significant seasonal variance between the three seasonal sampling events conducted. Groundwaters showed very little temporal variance but had great spatial variance both between the aquifer units of the DBAS and between different locations within the mountain block crystalline aquifer. The lowest d2H and d18O ratios were measured in winter precipitation, winter streams, and groundwater samples collected from the high-elevation Front Range. Samples of DBAS groundwaters with the lowest d2H and d18O ratios indicate potential hydrogeologic connection to the mountain block. Interpreted mixing lines on a d-excess versus d18O plot support the potential DBAS-mountain block connection. The deepest aquifer units of the DBAS (Arapahoe and Laramie-Fox Hills) show the least relationship with meteoric or surface waters on both a d2H and d18O plot and the d-excess versus d18O plot and have higher d18O values than would be predicted based on their previously measured recharge ages and paleoclimate data from the region. Characterizing the spatial and temporal variations in water stable isotope signatures of the Front Range and DBAS region enhances understanding of the region’s hydrology and hydrogeology. Additionally, these results help to better inform models of aquifer recharge and promote sustainable use of the DBAS resource.more » « less
-
Abstract Coastal deltaic aquifers are vulnerable to degradation from seawater intrusion, geogenic and anthropogenic contamination, and groundwater abstraction. The distribution and transport of contaminants are highly dependent on the subsurface sedimentary architecture, such as the presence of channelized features that preferentially conduct flow. Surface deposition changes in response to sea‐level rise (SLR) and sediment supply, but it remains unclear how these surface changes affect the distribution and transport of groundwater solutes in aquifers. Here, we explore the influence of SLR and sediment supply on aquifer heterogeneity and resulting effects on contaminant transport. We use realizations of subsurface heterogeneity generated by a process‐based numerical model, DeltaRCM, which simulates the evolution of a deltaic aquifer with different input sand fractions and rates of SLR. We simulate groundwater flow and solute transport through these deposits in three contamination scenarios: (a) vertical transport from widespread contamination at the land surface, (b) vertical transport from river water infiltration, and (c) lateral seawater intrusion. The simulations show that the vulnerability of deltaic aquifers to seawater intrusion correlates to sand fraction, while vertical transport of contaminants, such as widespread shallow contamination and river water infiltration, is influenced by channel stacking patterns. This analysis provides new insights into the connection between the depositional system properties and vulnerability to different modes of groundwater contamination. It also illustrates how vulnerability may vary locally within a delta due to depositional differences. Results suggest that groundwater management strategies may be improved by considering surface features, location within the delta, and the external forcings during aquifer deposition.more » « less
-
Abstract Evaluating stream water chemistry patterns provides insight into catchment ecosystem and hydrologic processes. Spatially distributed patterns and controls of stream solutes are well‐established for high‐relief catchments where solute flow paths align with surface topography. However, the controls on solute patterns are poorly constrained for low‐relief catchments where hydrogeologic heterogeneities and river corridor features, like wetlands, may influence water and solute transport. Here, we provide a data set of solute patterns from 58 synoptic surveys across 28 sites and over 32 months in a low‐relief wetland‐rich catchment to determine the major surface and subsurface controls along with wetland influence across the catchment. In this low‐relief catchment, the expected wetland storage, processing, and transport of solutes is only apparent in solute patterns of the smallest subcatchments. Meanwhile, downstream seasonal and wetland influence on observed chemistry can be masked by large groundwater contributions to the main stream channel. These findings highlight the importance of incorporating variable groundwater contributions into catchment‐scale studies for low‐relief catchments, and that understanding the overall influence of wetlands on stream chemistry requires sampling across various spatial and temporal scales. Therefore, in low‐relief wetland‐rich catchments, given the mosaic of above and below ground controls on stream solutes, modeling efforts may need to include both surface and subsurface hydrological data and processes.more » « less
-
Abstract Groundwater mixing dynamics play a crucial role in the biogeochemical cycling of shallow wetlands. In this paper, we conducted groundwater simulations to investigate the combined effects of evaporation and local heterogeneity on mixing dynamics in shallow wetland sediments. The results show that evaporation causes groundwater and solutes to upwell from deep sediments to the surface. As the solute reaches the surface, evaporation enhances the accumulation of the solute near the surface, resulting in a higher solute concentration than in deep sediments. Mapping of flow topology reveals that local heterogeneity generates spatially varied mixing patterns mainly along preferential flow pathways. The upwelling of groundwater induced by surface evaporation through heterogeneous sediments is likely to create distinct mixing hotspots that differ spatially from those generated by lateral preferential flows driven by large‐scale hydraulic gradients, which enhances the overall mixing in the subsurface. These findings have strong implications for biogeochemical processing in wetlands.more » « less
An official website of the United States government
