Search

To address groundwater issues, it is often necessary to develop geological and hydrogeological models. Combining geological, geophysical and hydrogeological data available on a site to build such models is often a challenge. This paper presents a methodology to integrate such data within a geologically consistent model with robust error estimation. The methodology combines the Ensemble Smoother with Multiple Data Assimilation (ESMDA) algorithm with a hierarchical geological modeling approach (ArchPy). Geophysical and hydrogeological field data are jointly assimilated in a stochastic ESMDA framework. To speed up the inversion process, forward responses are computed in lower‐dimensional spaces relevant to each physical problem. By doing so, the final models take into account multiple data sources and regional conceptual geological knowledge. This study illustrates the applicability of this novel approach using actual data from the upper Aare Valley, Switzerland. The results of cross‐validation show that the combination of different data types, each sensitive to different spatial dimensions, enhances the quality of the model within a reasonable computing time. The proposed methodology allows the automatic generation of groundwater models with robust uncertainty estimation and could be applied to a wide variety of hydrogeological issues. Plain Language Summary: When dealing with groundwater, it is necessary to develop underground models. However, taking into account all the different data types on a site is time‐consuming, and there is a lack of uncertainty quantification. In this study, we develop an approach that automatically combines different types of data, including boreholes, geophysical data, and hydrogeological measurements. All data are combined using a stochastic algorithm and produce an ensemble of plausible and data‐compatible models. These models can be used, for example, to forecast groundwater availability, pollutant distribution, or the effect of climate change on groundwater. Key Points: A methodology is proposed to assimilate geological, hydrogeological, and geophysical data in consistent stochastic modelsThe methodology combines a hierarchical geological modeling technique with the ensemble smoother with multiple data assimilationThe applicability of the methodology is demonstrated using actual field data from the upper Aare valley aquifer in Switzerland [ABSTRACT FROM AUTHOR]

Hydraulic connections between a river and an adjacent aquifer are controlled by the river resistance and aquifer diffusivity. In this paper, we derive a spectral solution linking the power spectrum of river stage fluctuations to that of the hydraulic head of a confined aquifer by means of a physical scaling factor. The physical scaling factor represents an algebraic expression of the river resistance and aquifer diffusivity and is included in an exact spectral solution derived herein. Statistical measures of the aquifer diffusivity and river resistance are provided by fitting the solution versus observed groundwater hydraulic head obtained at several distances and/or frequencies. At a study site in the middle reach of the Yangtze River and downstream of the Three Gorges Dam in China, we find systematic damping of the hydraulic head variations with distance from the river, which follows a fractal pattern driven by the river stage. In general, the estimated parameters are consistent with results reported in the literature, which supports the validity of the proposed spectral approach, although the paper discusses advantages and limitations due to application conditions. Plain Language Summary: The characterization of hydraulic connections between a river and an adjacent aquifer requires knowledge of two critical parameters: the river resistance and the aquifer diffusivity. In this paper, we derive a spectral solution to link the power spectrum of river stage fluctuations to that of the groundwater hydraulic head in a confined aquifer and show how such solutions can be used to simultaneously assess the two critical parameters. Parameters estimated from the field case located in the middle reach of the Yangtze River and downstream of the Three Gorges Dam in China are consistent with results reported in the literature. Application conditions, advantages, and limitations of this spectral approach are also discussed. Key Points: The power spectrum of the hydraulic head is linked to that of the river stage by a scaling factorRiver resistance and aquifer diffusivity are analyzed in spectral formsTheoretical developments are applied to a field case study of the Yangtze River, China [ABSTRACT FROM AUTHOR]

Low‐elevation coastal areas are increasingly vulnerable to seawater flooding as sea levels rise and the frequency and intensity of large storms increase with climate change. Seawater flooding can lead to the salinization of fresh coastal aquifers by vertical saltwater intrusion (SWI). Vertical SWI is often overlooked in coastal zone threat assessments despite the risk it poses to critical freshwater resources and salt‐intolerant ecosystems that sustain coastal populations. This review synthesizes field and modeling approaches for investigating vertical SWI and the practical and theoretical understanding of salinization and flushing processes obtained from prior studies. The synthesis explores complex vertical SWI dynamics that are influenced by density‐dependent flow and oceanic, hydrologic, geologic, climatic, and anthropogenic forcings acting on coastal aquifers across spatial and temporal scales. Key knowledge gaps, management challenges, and research opportunities are identified to help advance our understanding of the vulnerability of fresh coastal groundwater. Past modeling studies often focus on idealized aquifer systems, and thus future work could consider more diverse geologic, climatic, and topographic environments. Concurrent field and modeling programs should be sustained over time to capture interactions between physical processes, repeated salinization and flushing events, and delayed aquifer responses. Finally, this review highlights the need for improved coordination and knowledge translation across disciplines (e.g., coastal engineering, hydrogeology, oceanography, social science) to gain a more holistic understanding of vertical SWI. There also needs to be more education of communities, policy makers, and managers to motivate societal action to address coastal groundwater vulnerability in a changing climate. Plain Language Summary: Along the world's coastlines, seawater flooding on the land surface by storms and high tide events can drive the downward infiltration of seawater into coastal aquifers that are an important source of fresh groundwater. Rising sea levels and the intensification of coastal storms are expected to exacerbate seawater flooding events globally. The salinization of fresh groundwater along the coast following floods is often "out of sight, out of mind," but this process threatens freshwater resources that sustain dense coastal human populations and valuable, sensitive ecosystems. This paper reviews prior studies that have applied monitoring and/or computer modeling methods to better understand processes and controls on vertical saltwater intrusion. Our review details how surface flood dynamics, groundwater processes, environmental characteristics, and human activities combine to create complex patterns of salinization and flushing over space and time. From our synthesis we identify future research needs, including long‐term, globally dispersed monitoring and collaboration across scientific disciplines. We further emphasize the need for communication with local stakeholders about this topic to motivate informed management and research initiatives to understand and reduce the vulnerability of coastal fresh groundwater resources in a changing climate. Key Points: Vertical saltwater intrusion following coastal floods threatens fresh groundwater in an age of rising seas and intensifying stormsWe review salinization and flushing processes identified with field and modeling methods and identify future research opportunitiesKnowledge gaps can be addressed through multidisciplinary studies to advance science and inform water resources management and policy [ABSTRACT FROM AUTHOR]

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance‐based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain‐scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. These results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems. Plain Language Summary: Models that simulate the movement of groundwater and contaminants in aquifers, known as reactive transport models (RTMs), are complex. Multiple competing processes, including the physics of groundwater flow and the chemistry of microbial interactions, make the results of such models difficult to understand without additional statistical analyses. Sensitivity analysis, a technique for calculating the effect of each input variable on model output, is one potential tool for interpreting complex models, but it is rarely performed on RTMs due to the computational resources required. In this paper, we show that distance‐based generalized sensitivity analysis (DGSA) can be used to interpret reactive transport models and has several advantages over other techniques. As a case study, we use DGSA to measure the effect of 17 input variables on a model that simulates iron cycling within a shallow aquifer. Results show that groundwater flow rates control dissolved iron concentrations and are generally more influential than geochemical reaction rates. Though we focus on iron, the findings are relevant for other microbially driven reactions, which can control the mobility of many nutrients, contaminants, and metals in groundwater. Key Points: Global sensitivity analysis is a valuable tool for improving process understanding in subsurface biogeochemical modelsParameter rankings from local and global sensitivity analyses can differ when applied to reactive transport modelsIn certain settings, sediment permeability exerts a more dominant control on redox cycling than do geochemical reaction rates [ABSTRACT FROM AUTHOR]

The slug test is a common field technique for obtaining local hydraulic parameters near wells, applied, for example, for the hydrogeological investigation at contaminated sites. Although many slug test models have been developed for interpretation of measurements, only a few of them have considered heterogeneous conditions, and water column inertial effects are usually neglected. In this paper, we propose a novel three‐dimensional slug test model (3DHIM) for application in heterogeneous aquifers, considering inertial effects associated with skin effects and linear friction forces. After comparison with existing analytical and numerical solutions of slug tests, the model is applied to an aquifer analog to simulate a series of slug tests. The results from single‐well slug tests show that the well geometry (i.e., the well radius, well depth, and screen length) has an impact on the water level response. For cross‐well slug tests, the results indicate that the water level fluctuations not only include information on the hydraulic signal propagation process in the aquifer but also on well characteristics, such as wellbore storage and inertial effects. These effects cause a phase shift and amplitude change of the water level fluctuation. As the observation and test wells have a good hydraulic connection and similar well geometry, the water level amplitude could be amplified relative to aquifer pressure at the measured position. Therefore, we suggest considering wellbore storage and in‐well inertial effects in slug test‐based subsurface investigations, otherwise the parameter estimates based on well water levels may include errors, particularly in highly permeable layers. Key Points: We establish a three‐dimensional slug test model that accounts for aquifer heterogeneity and inertial effects in wellsThe new 3DHIM model is demonstrated in an aquifer analog studyThe impact of aquifer heterogeneity, in‐well inertial, and wellbore storage effects on simulation results is evaluated [ABSTRACT FROM AUTHOR]

This paper introduces the special section on 'new modeling approaches and novel experimental technologies for improved understanding of process dynamics at aquifer-surface water interfaces.' It is contextualizing the framework for the 27 research papers of the special section by firth identifying research gaps and imminent challenges for ecohydrological research at aquifer-surface water interfaces and then discussing the specific paper contributions on (i) new developments in temperature/heat tracing at GW-SW interfaces, (ii) new methods to capture the temporal and spatial variability of groundwater-surface water exchange, (iii) new approaches in modeling aquifer-river exchange flow, and (iv) new concepts and advanced theory of groundwater-surface water exchange. [ABSTRACT FROM AUTHOR]

Under transient flow conditions, the groundwater age distribution (GAD) is affected by both the physical characteristics of the aquifer, such as the hydraulic conductivity field, and the past temporal flow variations. Therefore, measuring the concentration of single or multiple tracers in a water sample collected at one time does not provide sufficient information on the transient nature of groundwater flows. However, obtaining multiple tracers at different times may contain such information. To find the temporal relationship between the GADs at a certain location, and to gain understanding about the current characteristics of an aquifer using environmental tracers, it is useful to separate the effects of temporal flow variations from that of the aquifer structure. In this paper, a new approach is presented, aiming at establishing a link between a reference residence time distribution based on steady-state flow and the transient age distributions at different times. Such information allows us to look collectively at sets of tracers measured at different times, when inferring aquifer structure and temporal trends in flow. The potential of using multiple tracers collected at different times to infer the GAD and its temporal variations is shown through two demonstration cases where measurements of CFCs, SF6, and 85Kr in hypothetical water samples are used to infer the overall temporal flow pattern in the aquifer. Bayesian inverse modeling via Markov Chain Monte Carlo method was used to estimate lumped-parameter groundwater ages and temporal flow patterns. [ABSTRACT FROM AUTHOR]

Using a first-order cross-correlation analysis, this paper investigates the relationship between observed heads and hydraulic properties in the saturated and vadose zones at different times and locations of three-dimensional unconfined aquifers during pumping tests. Cross-correlation analysis is a weighted sensitivity analysis casted into a stochastic framework. It determines the relative impact of each parameter with respect to others in time and space on the observed heads according to uncertainty or spatial variability of each parameter. It reveals the information content in measured drawdowns about heterogeneity during a pumping test in an unconfined aquifer, which is critical for aquifer parameter estimation. Based on a synthetic, numerical example, our cross-correlation analysis reveals that heads in the saturated zone at late times carry the greatest nonsymmetrically weighted information content about the hydraulic conductivity (KS) distribution within the cone of depression. On the other hand, heads in the saturated zone at early times contain the most information about the specific storage (SS) heterogeneity in a narrow region between the observation and pumping locations. During intermediate and late times, heads in the saturated zone largely reflect the effects of saturated water content (θS) and pore-size parameter (α) in the thin unsaturated region near the water table above the pumping and observation locations. At last, heads in the vadose zone at late times carry the greatest information about θS and α around the observation point. [ABSTRACT FROM AUTHOR]

Hydrologists routinely analyze pumping test data using conventional interpretation methods that are based on the assumption of homogeneity and that, consequently, yield single estimates of representative flow parameters. However, natural subsurface formations are intrinsically heterogeneous, and hence, the flow parameters influencing the drawdown vary as the cone of depression expands in time. In this paper a novel procedure for the analysis of pumping tests in heterogeneous confined aquifers is developed. We assume that a given heterogeneous aquifer can be represented by a homogeneous system whose flow parameters evolve in time as the pumping test progresses. At any point in time, the interpreted flow parameters are estimated using the ratio of the drawdown and its derivative observed at that particular time. The procedure is repeated for all times, yielding time-dependent estimates of transmissivity Ti( t) and storativity, Si( t). Based on the analysis of the sensitivity of drawdown to inhomogeneities in the T field, the time-dependent interpreted transmissivity values are found to be a good estimate of Tg( r), the geometric mean of the transmissivity values encompassed within a progressively increasing radius r from the well. The procedure is illustrated for Gaussian heterogeneous fields with ln( T) variances up to a value of 2. The impact of the separation distance between the pumping well and observation point on data interpretation is discussed. The results show that information about the spatial variability of the transmissivity field can be inferred from time-drawdown data collected at a single observation point. [ABSTRACT FROM AUTHOR]

Do observation wells in fractured porous aquifers measure water head in the fracture network, water head in the matrix, or some combination of both? This question necessarily arises when calibrating dual‐continuum hydrogeological models against on‐field data. One can assume that observation wells measure fracture water head, because matrix permeability is negligible compared to fracture permeability. Nevertheless, this reasoning is invalid for wells poorly‐connected to the fractures. Yet, the possibility of such a poor connection at given depths has never been implemented in a physics‐based manner when comparing matrix and fracture water heads simulated by dual‐continuum models to on‐field data. To fill this knowledge gap, a physically based, easy to calibrate, open‐source postprocessing tool, POWeR‐FADS (Program for Observation Well Representation in Fractured Aquifer Dual‐continuum Simulations), available at https://github.com/BJeannot1/POWeR-FADS, has been developed. It introduces as parameters well geometry and the altitude of lowest interception of the fractures by the well. From these, POWeR‐FADS nonintrusively postprocesses time series of matrix and fracture water heads at the well, as simulated by any planar, bidimensional dual‐continuum hydrogeological model, to calculate water exchanges involving the observation well and thus the evolution of water level in the well. Synthetic test cases show that POWeR‐FADS makes it possible to simulate peculiar behaviors that are similar to patterns actually observed by the authors in on‐site observation wells of a fractured porous aquifer, like "floors" in observed water levels, delayed but sharp rises at the beginning of recharge events, or inflexion points accelerating the drawdown velocity during the recession phase. Key Points: A program is made to simulate the water level in observation wells from matrix and fracture water heads computed by dual‐continuum modelsFor wells poorly connected to the fractures, the code simulates patterns that could not be modeled by a dual‐continuum model aloneThis postprocessing tool is easy to use, as its most sensitive parameter can be deduced from drill cores or observation well data [ABSTRACT FROM AUTHOR]

Unsaturated flow influences both the seawater extent under steady‐state conditions and the propagation of tides in coastal aquifers. However, its effects on salt distributions in tidally influenced coastal aquifers are little investigated. The present study used numerical simulations and data from laboratory experiments to analyze the effects of unsaturated flow on density‐dependent solute transport in coastal unconfined aquifers. The effects of the inland boundary condition (i.e., constant‐head or constant‐flux) were tested. Compared to a stable sea level, the results show that unsaturated flow has a more pronounced influence on salt distributions in coastal unconfined aquifers when tides are considered, regardless of the type of inland boundary condition. Neglect of unsaturated flow effects leads to expansion of the upper saline plume (USP), shrinkage of the saltwater wedge (seaward movement of saltwater wedge), and overestimation of water and salt exchange across the aquifer‐ocean interface. This is caused by a lower head in the nearshore area during high‐tide periods with the unsaturated zone effects removed. Thus, without the unsaturated zone, stronger head gradients within the nearshore aquifer occur at high tide, leading to stronger tidally driven seawater infiltration and hence a larger USP. Counterintuitively, ignoring unsaturated flow effects leads to greater average inland head over a tidal period, which shifts the saltwater wedge seaward. It is concluded that unsaturated zone effects should not be neglected for modeling tide‐affected seawater intrusion, especially if quantification of near‐shore conditions is important. Key Points: Tidal fluctuations intensify unsaturated flow effects on salt distributions in coastal unconfined aquifersNeglecting unsaturated flow effects leads to a more extensive upper saline plume and shrinking of the saltwater wedgeNeglecting unsaturated flow effects intensifies the water and salt exchanges across the aquifer‐ocean interface [ABSTRACT FROM AUTHOR]

Among the impacts of earthquakes on aquifers, permeability change is one of the most challenging to quantify, since techniques to measure permeability evolution are scarce. The study of tidal response of boreholes is one of the most promising, yet complex to use in practice. We used 14 years of piezometric level measurements and two concurrent source signals, earth tidal strain and barometric pressure, for which we separated the respective contribution in a state‐of‐the‐art tidal analysis. We developed a new general analytical hydrogeological model, based on geological observations of a confined bi‐layer aquifer. It is able to match combined observations of earth and barometric tide phase lags which could not be explained by existing models. We demonstrate that its relative complexity can be overcome thanks to the results of tidal analysis, yielding a simpler model adapted to the Fond Lahaye site of the Martinique Island. The resulting evolution of diffusivity and loading efficiency, was validated independently with several pumping tests occurring all along the studied period. The transient diffusivity increases and decreases indicate which earthquakes impacted the aquifer, enabling to establish an empirical magnitude‐distance relationship criterion. This criterion confirms the suspected dependence on dynamic stresses, which decrease as the square of the hypocentral distance. Additionally, we investigate two other factors of diffusivity changes: heavy rainfall events and aquifer withdrawals, which demonstrates the sensitivity of volcanic aquifers properties to environmental and anthropogenic influence. Key Points: A model for the tidal response of bi‐layer aquifers was developed, focusing on pressure diffusion and exchange between two confined layersUsing 14 years of hourly piezometric data, the temporal evolution of aquifer diffusivity was inferred and confirmed independently by pumping testsThe change of permeability is shown to be due to earthquakes, extreme rainfall events and aquifer withdrawals in upstream borehole [ABSTRACT FROM AUTHOR]

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. Plain Language Summary: The findings of this study provide insight into the vulnerability of deltaic aquifers to three contamination processes: (a) widespread contaminant transport from the land surface, (b) river water infiltration, and (c) seawater intrusion. We consider how contamination is affected by the location of contaminants and the processes associated with the accumulation of sediments in deltas. Our work shows that vulnerability to contamination depends on how the aquifer is deposited. The results also demonstrate that the distribution of sandy channels preserved in the subsurface, as well as rivers on the surface, controls vertical contaminant transport. We find that these effects vary from upstream to downstream in the delta because of spatial differences in depositional processes. These findings will help to improve predictions of groundwater contamination and manage groundwater development in deltas around the world. Key Points: In deltaic aquifers, contaminant distribution is controlled by the sand fraction and channel stacking patternsRivers on the surface have a strong influence on vertical contaminant transport in aquifersDue to differences in channel distribution, vertical contamination varies nonlinearly from upstream to downstream regions [ABSTRACT FROM AUTHOR]

Large‐scale groundwater models are required to estimate groundwater availability and to inform water management strategies on the national scale. However, parameterization of large‐scale groundwater models covering areas of major river basins and more is challenging due to the lack of observational data and the mismatch between the scales of modeling and measurements. In this work, we propose to bridge the scale gap and derive regional hydraulic parameters by spectral analysis of groundwater level fluctuations. We hypothesize that specific locations in aquifers can reveal regional parameters of the hydraulic system. We first generate ensembles of synthetic but realistic aquifers which systematically differ in complexity. Applying Liang and Zhang's (2013), https://doi.org/10.1016/j.jhydrol.2012.11.044, semi‐analytical solution for the spectrum of hydraulic head time series, we identify for each ensemble member and at different locations representative aquifer parameters. Next, we extend our study to investigate the use of spectral analysis in more complex numerical models and in real settings. Our analyses indicate that the variance of inferred effective transmissivity and storativity values for stochastic aquifer ensembles is small for observation points which are far away from the Dirichlet boundary. Moreover, the head time series has to cover a period which is roughly 10 times as long as the characteristic time of the aquifer. In deterministic aquifer models we infer equivalent, regionally valid parameters. A sensitivity analysis further reveals that as long as the aquifer length and the position of the groundwater measurement location is roughly known, the parameters can be robustly estimated. Plain Language Summary: We build large‐scale (regional) computer models of the subsurface flow conditions in order to quantify the long‐term shift in groundwater storage and response on the national level under changing climatic conditions and increasing human water demands. These models must be fed with hydrogeological parameters obtained from subsurface observation wells, drilling logs, and hydraulic tests in conjunction with (hydro)geological and geostatistical methods. In some regions these wells are sparsely distributed and derived parameters are representative only for small areas. We hypothesize that groundwater level records can reveal regional aquifer information when analyzed in the spectral domain. In order to bridge that scale gap and because groundwater level time series are generally available, we propose to infer regional parameters by analyzing the frequency content (spectrum) of long groundwater level time series. The required parameters were determined using mathematical formulations of the theoretical spectrum for simplified settings. We tested the methodology in computer models with limited complexity and found that the groundwater level time series indeed contain regional information if the time of observation is sufficiently long. Lastly, we apply the spectral analysis to real groundwater data to test the capability of the method to infer regional aquifer parameters in real aquifers. Key Points: We successfully tested the spectral analysis of groundwater level fluctuations in numerical models and obtained regional aquifer parametersIn a sensitivity analysis of the spectral analysis using field data, the storativity and the response times could be robustly estimatedThe application of the suggested methodology to the field data from a catchment in central Germany produced plausible results [ABSTRACT FROM AUTHOR]

Microbial pollution of aquifers is a persistent water quality problem globally which poses significant risks to public health. Karst aquifer systems are exceptionally vulnerable to pollution from fecal contamination sources as a result of rapid recharge of water from the surface via discrete pathways linked to highly conductive, solutionally enlarged conduits alongside strong aquifer heterogeneity. Consequently, rapid changes in microbial water quality, which are difficult to monitor with expensive and time‐consuming conventional microbiological methods, are a major concern in karst environments. This study examined flow cytometric (FCM) fingerprinting of bacterial cells in groundwater together with fecal indicator bacteria (FIB) at nine separate karst springs of varying catchment size over a 14 month period in order to assess whether such a technique can provide faster and more descriptive information about microbial pollution through such karst aquifer systems. Moreover, the data have also been evaluated with respect to the potential of using turbidity as an easy‐to‐measure proxy indicator of microbial pollution in a novel way. We argue that FCM provides additional data from which enhanced insights into fecal pollution sources and its fate and transport in such karst catchments can be gained. We also present valuable new information on the potential and limitations of turbidity as an indicator of fecal groundwater contamination in karst. FCM has the potential to become a more widely used tool in the field of contaminant hydrogeology. Key Points: Flow cytometry can provide enhanced insights into the fate and transport of fecal pollution in karst aquifersStable bacterial communities in karst aquifers can limit the potential of flow cytometry to predict changes in fecal indicator bacteriaEnterococci should be analyzed with Escherichiacoli if turbidity is to be evaluated as an inexpensive groundwater quality indicator in karst springs [ABSTRACT FROM AUTHOR]

Uncertainty in model predictions is caused to a large extent by the uncertainty in model parameters, while the identification of model parameters is demanding because of the inherent heterogeneity of the aquifer. A variety of inverse methods has been proposed for parameter identification. In this paper we present a novel inverse method to constrain the model parameters (hydraulic conductivities) to the observed state data (hydraulic heads). In the method proposed we build a conditioning pattern consisting of simulated model parameters and observed flow data. The unknown parameter values are simulated by pattern searching through an ensemble of realizations rather than optimizing an objective function. The model parameters do not necessarily follow a multi-Gaussian distribution, and the nonlinear relationship between the parameter and the response is captured by the multipoint pattern matching. The algorithm is evaluated in two synthetic bimodal aquifers. The proposed method is able to reproduce the main structure of the reference fields, and the performance of the updated model in predicting flow and transport is improved compared with that of the prior model. [ABSTRACT FROM AUTHOR]

In this study, we developed two‐dimensional hydrodynamic age expressions for the free‐surface aquifer problem required for radiocarbon data interpretation in terms of recharge. Time‐varying recharge is accounted for in fully transient or pseudo‐steady state alternative analytical expressions applying to weakly or highly hydraulically reactive unconfined aquifers, respectively. These expressions extend the well‐known model proposed by Vogel in 1967 and represent a convenient alternative approach to the more complex numerical solution of a particle tracking problem considering climate variability. An application of the developed expressions was carried out to assess past (40 Kyr) and present recharge rates for the North Western Sahara Aquifer System (NWSAS) and the Djeffara plain (Tunisia‐Libya). Using the available 14C data at different outcrops, we obtained local recharge values ranging from 0.2 ± 1.7 to 4.5 ± 3.2 mm yr−1. Here, it is shown that using the standard purely steady‐state model for radiocarbon interpretation leads to a present recharge overestimation of up to 360%. A linear relation between the present recharge and the average annual rainfall was obtained and used to identify the regional distribution of the recharge and domain‐averaged values. Global values of 1.6 ± 2.3 mm yr−1 for the NWSAS and 3.5 ± 2.3 mm yr−1 for the Djeffara are in excellent agreement with the values obtained by hydrogeological models and derived using Satellite gravity data. The relation between recharge and rainfall obtained here can be used as a first‐order estimate for hydrogeological modeling or water management analysis for aquifers in the Saharan area. Key Points: Analytical expressions for Radiocarbon distribution in transient state and unconfined aquifersDetermination of recharge rates of North Africa aquifers by using radiocarbon tracers accounting for past climate variabilityThe use of Vogel's standard solution for North Africa aquifers leads to an average overestimation of present‐day recharge of about 70% [ABSTRACT FROM AUTHOR]

As shallow aquifers become depleted or contaminated worldwide, use of deep aquifers will likely increase to meet growing water demands despite expensive drilling costs and well maintenance, and limited recharge to many of these deep systems. We discuss depletion of the Cambrian‐Ordovician sandstone aquifer, a deep bedrock aquifer that has been a major source of water in the Midwestern US for over 150 years. Using a participatory groundwater modeling approach with active stakeholder engagement, we developed a conceptual framework to assess risk to different sandstone units of the aquifer based on simulated heads and production well behavior (pumping levels and specific capacities). Over 320 m of static head decline has occurred since predevelopment in the deepest sandstone unit. Shallower units are actively being dewatered. Drawdown and head separation between sandstone units is even greater under pumping conditions and is most acute near a regional fault zone. Current depletion is prompting massive infrastructure changes to mitigate anticipated water supply shortages. Our results show that even with a substantial reduction in water use by 35% in 2030, many areas of the aquifer remain at High or Severe risk out to 2070. Head decline, as opposed to storage loss, is a critical metric for understanding depletion of deep bedrock aquifers. The participatory modeling process made clear that pumping levels at production wells should not be ignored in evaluations of aquifer sustainability and risk. Our conceptual risk framework could be applied to other deep bedrock aquifers undergoing depletion. Key Points: Over 320 m of head decline has occurred in the deep sandstone aquifer system underlying the United States' third largest metropolitan areaParticipatory modeling allows assessment of risk to groundwater supplies under pumping conditionsAssessment of deep aquifer sustainability requires understanding of multi‐aquifer wells, geologic barriers, and production well behavior [ABSTRACT FROM AUTHOR]

The depth of active groundwater circulation is a fundamental control on stream flows and chemistry in mountain watersheds, yet it remains challenging to characterize and is rarely well constrained. We collected hydraulic conductivity, hydraulic head, temperature, chemical, noble gas, and 3H/3He groundwater age data from discrete levels in two boreholes 46 and 81 m deep in an alpine watershed, in combination with chemical and age data from shallow groundwater discharge, to discern groundwater flow rates at different depths and directly observe active and inactive groundwater. Vertical head gradients are steep (average of 0.4) and thermal profiles are consistent with typical linear conductive continental geotherms. Groundwater deeper than ∼20 m is distinct from shallow groundwater and creek water in its chemistry, noble gas signature, and age (dominantly >65 years compared to <9 years). Together these results suggest low vertical groundwater flow velocities and a relatively shallow active circulation depth of ∼20 m. This hypothesis is tested with a simple 2‐D numerical fluid flow and heat transport model representing a hillslope transect through the two boreholes. The modeling indicates that the subhorizontally bedded sedimentary rocks underlying the basin are highly anisotropic with low vertical hydraulic conductivity, and at most ∼10% of bedrock recharge (equivalent to <2% of stream baseflow) flows below a depth of 20 m. The study demonstrates the considerable value of discrete‐depth hydrogeologic, chemical, and age data for determining active circulation depth, and illustrates an approach for maximizing the utility of individual boreholes drilled for mountain bedrock aquifer characterization. Key Points: The depth of active groundwater circulation is a fundamental characteristic of mountain watersheds, yet is rarely well constrainedHydrogeologic, groundwater chemistry, and 3H/3He age data were collected at discrete depths from two boreholes 46 and 81 m deepBorehole data combined with shallow groundwater data and numerical modeling indicate an active circulation depth of ∼20 m [ABSTRACT FROM AUTHOR]

GRACE satellite data are widely used to estimate groundwater storage (GWS) changes in aquifers globally; however, comparisons with GW monitoring and modeling data are limited. Here we compared GWS changes from GRACE over 15 yr (2002–2017) in 14 major U.S. aquifers with groundwater‐level (GWL) monitoring data in ~23,000 wells and with regional and global hydrologic and land surface models. Results show declining GWS trends from GRACE data in the six southwestern and south‐central U.S. aquifers, totaling −90 km3 over 15 yr, related to long‐term (5–15 yr) droughts, and exceeding Lake Mead volume by ~2.5×. GWS trends in most remaining aquifers were stable or slightly rising. GRACE‐derived GWS changes agree with GWL monitoring data in most aquifers (correlation coefficients, R = 0.52–0.95), showing that GRACE satellites capture groundwater (GW) dynamics. Regional GW models (eight models) generally show similar or greater GWS trends than those from GRACE. Large discrepancies in the Mississippi Embayment aquifer, with modeled GWS decline approximately four times that of GRACE, may reflect uncertainties in model storage parameters, stream capture, pumpage, and/or recharge rates. Global hydrologic models (2003–2014), which include GW pumping, generally overestimate GRACE GWS depletion (total: approximately −172 to −186 km3) in heavily exploited aquifers in southwestern and south‐central U.S. by ~2.4× (GRACE: −74 km3), underscoring needed modeling improvements relative to anthropogenic impacts. Global land surface models tend to track GRACE GWS dynamics better than global hydrologic models. Intercomparing remote sensing, monitoring, and modeling data underscores the importance of considering all data sources to constrain GWS uncertainties. Plain Language Summary: The major U.S. aquifers provide an ideal system to assess GRACE (Gravity Recovery and Climate Experiment) satellite data. We compared GRACE groundwater storage anomalies (GWSAs) with groundwater level anomalies (GWLAs) from ~23,000 wells and with groundwater storage (GWS) from regional and global models in 14 major U.S. aquifers. Results show large GWS declines from GRACE in southwestern (Central Valley and Arizona Alluvial Basins) and south‐central (Central and Southern High Plains and Texas) aquifers from multiyear droughts (5–15 yr). In contrast, GWS trends in aquifers throughout the rest of the U.S. showed mostly stable or rising values. Time series of GRACE GWSAs compared favorably with GWLAs from most aquifers, suggesting that GRACE data track groundwater (GW) dynamics. Regional GW models show similar or greater declines in GWS compared with GRACE data, with the largest discrepancy of a factor of four times greater modeled depletion in the Mississippi Embayment. Global hydrologic models show minimal storage dynamics but greatly overestimated GWS declines by ~2.4× in southwestern and south‐central aquifers with intensive irrigation compared with GRACE data. In contrast, global land surface models show similar GW dynamics to GRACE data but underestimated GWS declines in heavily exploited aquifers because these land surface models do not include human intervention. Key Points: GRACE satellites track GW storage in U.S. aquifers based on good agreement between GRACE and in situ GW level monitoring (23,000 wells)Regional models capture GRACE GW storage in most aquifers, except the Mississippi Embayment, which overestimates GRACE GWS trends by approximately four timesGlobal hydrologic models overestimate GW depletion compared with GRACE by ~2.4× in heavily depleted aquifers in the SW and SC United States [ABSTRACT FROM AUTHOR]

Riparian lowlands act as interfaces between the streams and upland areas. This study identified and quantified local flow paths in four transects of a 26 ha Danish riparian lowland in a glacial till landscape. The riparian lowland was fed by drain water from the upland agricultural drainage catchments. Precipitation, stream stage, and drainage discharge into the riparian lowland were measured continuously, while groundwater hydraulic heads were measured in piezometer pipes twice per month. A water balance model was developed to quantify water fluxes leaving the riparian lowland area via evapotranspiration, leakage to a deeper aquifer, and via groundwater flow, drain flow, and overland flow to the adjacent stream. Overland flow originating from the tile drains was the main flow path in all four transects, and also fluxes to the stream via groundwater or lowland tile drains were significant in some subareas. The presence of a secondary tile drainage network within parts of the riparian lowland reduced overland flow and increased interaction with the riparian lowland soils. Area‐normalized fluxes varied greatly between transects, largely reflecting variations in hydraulic loading rate (ratio of water input rate to the area of the receiving riparian lowland). This, combined with the significance of groundwater flow and riparian lowland tile drain flow in some of the investigated transects, revealed a heterogeneous distribution of flow paths within the small headwater lowland. The rate of overland flow was highly correlated to the hydraulic loading rate, which in turn was dominated by the drainage discharge rate at the hillslope boundary. Key Points: Overland flow, bypassing the riparian lowland, is produced when the hydraulic loading exceeds riparian lowland infiltrationMagnitudes of overland flow are correlated with the ratio of catchment area to riparian lowland areaRiparian lowland drains reduce overland flow and increase water residence times [ABSTRACT FROM AUTHOR]

Mountain‐block recharge (MBR) is the subsurface inflow of groundwater to lowland aquifers from adjacent mountains. MBR can be a major component of recharge but remains difficult to characterize and quantify due to limited hydrogeologic, climatic, and other data in the mountain block and at the mountain front. The number of MBR‐related studies has increased dramatically in the 15 years since the last review of the topic was conducted by Wilson and Guan (2004), generating important advancements. We review this recent body of literature, summarize current understanding of factors controlling MBR, and provide recommendations for future research priorities. Prior to 2004, most MBR studies were performed in the southwestern United States. Since then, numerous studies have detected and quantified MBR in basins around the world, typically estimating MBR to be 5–50% of basin‐fill aquifer recharge. Theoretical studies using generic numerical modeling domains have revealed fundamental hydrogeologic and topographic controls on the amount of MBR and where it originates within the mountain block. Several mountain‐focused hydrogeologic studies have confirmed the widespread existence of mountain bedrock aquifers hosting considerable groundwater flow and, in some cases, identified the occurrence of interbasin flow leaving headwater catchments in the subsurface—both of which are required for MBR to occur. Future MBR research should focus on the collection of high‐priority data (e.g., subsurface data near the mountain front and within the mountain block) and the development of sophisticated coupled models calibrated to multiple data types to best constrain MBR and predict how it may change in response to climate warming. Key Points: Mountain‐block recharge confirmed as important source of recharge to basin aquifers in a variety of climatic and geologic settings globallyRecent work advanced understanding of fundamental controls on mountain‐block recharge and somewhat improved methods for characterizationFuture research should aim to acquire subsurface data in mountain block and at the mountain front [ABSTRACT FROM AUTHOR]

Hydrogeological field studies rely often on a single conceptual representation of the subsurface. This is problematic since the impact of a poorly chosen conceptual model on predictions might be significantly larger than the one caused by parameter uncertainty. Furthermore, conceptual models often need to incorporate geological concepts and patterns in order to provide meaningful uncertainty quantification and predictions. Consequently, several geologically realistic conceptual models should ideally be considered and evaluated in terms of their relative merits. Here, we propose a full Bayesian methodology based on Markov chain Monte Carlo to enable model selection among 2‐D conceptual models that are sampled using training images and concepts from multiple‐point statistics. More precisely, power posteriors for the different conceptual subsurface models are sampled using sequential geostatistical resampling and Graph Cuts. To demonstrate the methodology, we compare and rank five alternative conceptual geological models that have been proposed in the literature to describe aquifer heterogeneity at the MAcroDispersion Experiment site in Mississippi, USA. We consider a small‐scale tracer test for which the spatial distribution of hydraulic conductivity impacts multilevel solute concentration data observed along a 2‐D transect. The thermodynamic integration and the stepping‐stone sampling methods were used to compute the evidence and associated Bayes factors using the computed power posteriors. We find that both methods are compatible with multiple‐point statistics‐based inversions and provide a consistent ranking of the competing conceptual models considered. Key Points: A full Bayesian method is proposed for model selection among geologically realistic conceptual modelsThermodynamic integration and stepping‐stone sampling provide consistent ranking of conceptual modelsThe method is applied to solute concentration data collected during a tracer test at the MADE site [ABSTRACT FROM AUTHOR]

Hydrologic responses to earthquakes such as streamflow increases, water‐level changes, and changes in geyser eruption frequency often reflect changes in permeability caused by seismic waves. The dynamic nature of permeability, as revealed by coseismic hydrologic phenomena, holds implications for groundwater systems, geothermal resources, mineral resources, and geologic hazards. Analysis of water‐level responses to solid Earth tides and changes in atmospheric pressure provides a passive way to continuously monitor changes in permeability and storage properties in tectonically active regions. Plain Language Summary: Earthquakes affect groundwater by changing properties of aquifers and their ability to transmit fluids. Monitoring how water levels in wells respond to tides and atmospheric pressure provides an inexpensive way to monitor changes in aquifers. Properties of aquifers are not static but change over time. Key Points: Earthquakes affect groundwater by changing properties of aquifers and their ability to transmit fluidsMonitoring how water levels in wells respond to tides and atmospheric pressure provides an inexpensive way to monitor changes in aquifersProperties of aquifers are not static but change over time [ABSTRACT FROM AUTHOR]

Earthquakes have been inferred to induce hydrological changes in aquifers on the basis of either changes to well water‐levels or tidal behavior, but the relationship between these changes remains unclear. Here, changes in tidal behavior and water‐levels are quantified using a hydrological network monitoring gravel aquifers in Canterbury, New Zealand, in response to nine earthquakes (of magnitudes Mw 5.4 to 7.8) that occurred between 2008 and 2015. Of the 161 wells analyzed, only 35 contain water‐level fluctuations associated with "Earth + Ocean" (7) or "Ocean" (28) tides. Permeability reduction manifest as changes in tidal behavior and increased water‐levels in the near field of the Canterbury earthquake sequence of 2010–2011 support the hypothesis of shear‐induced consolidation. However, tidal behavior and water‐level changes rarely occurred simultaneously (~2%). Water‐level changes that occurred with no change in tidal behavior reequilibrated at a new postseismic level more quickly (on timescales of ~50 min) than when a change in tidal behavior occurred (~240 min to 10 days). Water‐level changes were more than likely to occur above a peak dynamic stress of ~50 kPa and were more than likely to not occur below ~10 kPa. The minimum peak dynamic stress required for a tidal behavior change to occur was ~0.2 to 100 kPa. Key Points: Tidal behavior and water‐level changes were quantified in 35 wells in response to nine earthquakes larger than Mw 5.4Earthquake‐induced water‐level and tidal behavior changes rarely occurred simultaneously in the same monitoring well (~2%)Permeability change thresholds correspond to peak dynamic stresses of ~0.2 to 100 kPa [ABSTRACT FROM AUTHOR]

We address the problem of stochastic simulation of soil particle-size curves (PSCs) in heterogeneous aquifer systems. Unlike traditional approaches that focus solely on a few selected features of PSCs (e.g., selected quantiles), our approach considers the entire particle-size curves and can optionally include conditioning on available data. We rely on our prior work to model PSCs as cumulative distribution functions and interpret their density functions as functional compositions. We thus approximate the latter through an expansion over an appropriate basis of functions. This enables us to (a) effectively deal with the data dimensionality and constraints and (b) to develop a simulation method for PSCs based upon a suitable and well defined projection procedure. The new theoretical framework allows representing and reproducing the complete information content embedded in PSC data. As a first field application, we demonstrate the quality of unconditional and conditional simulations obtained with our methodology by considering a set of particle-size curves collected within a shallow alluvial aquifer in the Neckar river valley, Germany. [ABSTRACT FROM AUTHOR]

We analyze the probability distribution of the hazard attenuation factor for a noncarcinogenic reactive compound captured by a well in heterogeneous porous formations. The hazard attenuation factor is defined as the ratio between the hazard index HI at a detection well and at the source. Heterogeneity of the aquifer is represented through the multi-indicator model (a collection of blocks of independent permeability) while flow and transport are solved by the means of the self-consistent approach that is able to deal with any degree of heterogeneity. Due to formation heterogeneity, HI is a random variable and similar for hazard attenuation index. The latter can be fully characterized by its cumulative distribution function (CDF), which in turn can be related to the statistics of the travel time of solute particles, from the source to the detection well. The approach is applied to the case of a solute which undergoes decay and a well with a screen much smaller than the correlation scale of hydraulic conductivity. The results show that the probability of exceeding a given acceptable threshold of the hazard index is significantly affected by the level of heterogeneity comparable to the one observed for the MADE site, and the distance between the source and the well. [ABSTRACT FROM AUTHOR]

A sequential procedure of hydraulic tomographical inversion is applied to characterize at high resolution the spatial heterogeneity of hydraulic conductivity and specific storage at the field test site Stegemühle, Germany. The shallow aquifer at this site is examined by five short-term multilevel pumping tests with 30 pumping-observation pairs between two wells. Utilizing travel time diagnostics of the recorded pressure response curves, fast eikonal-based inversion is shown to deliver insight into the sedimentary structures. Thus, the structural information from the generated travel time tomogram is exploited to constrain full calibration of the pressure response curves. Based on lateral extrapolation from the measured inter-well profile, a three-dimensional reconstruction of the aquifer is obtained. It is demonstrated that calibration of spatially variable specific storage in addition to hydraulic conductivity can improve the fitting of the model while the structural features are only slightly changed. At the field site, two tracer tests with uranine and sodium-naphthionate were also performed and their concentrations were monitored for 2 months. The measured tracer breakthrough curves are employed for independent validation of the hydraulic tomographical reconstruction. It is demonstrated that major features of the observed solute transport can be reproduced, and structures relevant for macrodispersive tracer spreading could be resolved. However, for the mildly heterogeneous aquifer, the tracer breakthrough curves can also be approximated by a simplified homogeneous model with higher dispersivity. Therefore, improved validation results that capture specific characteristics of the breakthrough curves would require additional hydraulic measurements. [ABSTRACT FROM AUTHOR]

Aquifer hydraulic properties such as hydraulic conductivity ( K) are ubiquitously heterogeneous and typically only a statistical characterization can be sought. Additionally, statistical anisotropy at typical characterization scales is the rule. Thus, regardless of the processes governing solute transport at the local (pore) scale, transport becomes non-Fickian. Mass-transfer models provide an efficient tool that reproduces observed anomalous transport; in some cases though, these models lack predictability as model parameters cannot readily be connected to the physical properties of aquifers. In this study, we focus on a multirate mass-transfer model (MRMT), and in particular the apparent capacity coefficient ( β), which is a strong indicator of the potential of immobile zones to capture moving solute. We aim to find if the choice of an apparent β can be phenomenologically related to measures of statistical anisotropy. We analyzed an ensemble of random simulations of three-dimensional log-transformed multi-Gaussian permeability fields with stationary anisotropic correlation under convergent flow conditions. It was found that apparent β also displays an anisotropic behavior, physically controlled by the aquifer directional connectivity, which in turn is controlled by the anisotropic correlation model. A high hydraulic connectivity results in large β values. These results provide new insights into the practical use of mass-transfer models for predictive purposes. [ABSTRACT FROM AUTHOR]

Calibration of subsurface flow models is an essential step for managing ground water aquifers, designing of contaminant remediation plans, and maximizing recovery from hydrocarbon reservoirs. We investigate an efficient sampling algorithm known as nested sampling (NS), which can simultaneously sample the posterior distribution for uncertainty quantification, and estimate the Bayesian evidence for model selection. Model selection statistics, such as the Bayesian evidence, are needed to choose or assign different weights to different models of different levels of complexities. In this work, we report the first successful application of nested sampling for calibration of several nonlinear subsurface flow problems. The estimated Bayesian evidence by the NS algorithm is used to weight different parameterizations of the subsurface flow models (prior model selection). The results of the numerical evaluation implicitly enforced Occam's razor where simpler models with fewer number of parameters are favored over complex models. The proper level of model complexity was automatically determined based on the information content of the calibration data and the data mismatch of the calibrated model. [ABSTRACT FROM AUTHOR]

Numerical solutions to the nonlinear Boussinesq equation, applied to a steeply sloping aquifer and assuming uniform hydraulic conductivity, indicate that late-time recession discharge decreases nearly linearly in time. When recession discharge is characterized by −d Q/d t = aQb, this is equivalent to constant d Q/d t or b = 0. This result suggests that a previously reported exponential decrease with time ( b = 1) of modeled recession discharge from a similar sloping aquifer represented by the same equation appears to be an artifact of the numerical solution scheme and its interpretation. Because the linearly decreasing recession discharge ( b = 0) is not known from field studies, these findings challenge the application of a nonlinear Boussinesq framework assuming uniform conductivity and geometric similarity to infer hydraulic properties of sloping aquifers from observations of streamflow. This finding also questions the validity of the physical interpretation of the exponential decline in late time resulting from the commonly used linearized form of the Boussinesq equation, opposed to the full nonlinear equation, when applied under these conditions. For this reason, application of the linearized equation to infer hydraulic properties of sloping aquifers is also challenged, even if the observed recession is consistent with that of the linearized Boussinesq equation. [ABSTRACT FROM AUTHOR]

Subsurface aquifer characterization often involves high parameter dimensionality and requires tremendous computational resources if employing a full Bayesian approach. Ensemble-based data assimilation techniques, including filtering and smoothing, are computationally efficient alternatives. Despite the increasing use of ensemble-based methods in assimilating flow and transport related data for subsurface aquifer characterization, most applications have been limited to synthetic studies or two-dimensional problems. In this study, we applied ensemble-based techniques adapted for parameter estimation, including the p-space ensemble Kalman filter and ensemble smoother, for assimilating field tracer experimental data obtained from the Integrated Field Research Challenge (IFRC) site at the Hanford 300 Area. The forward problem was simulated using the massively parallel three-dimensional flow and transport code PFLOTRAN to effectively deal with the highly transient flow boundary conditions at the site and to meet the computational demands of ensemble-based methods. This study demonstrates the effectiveness of ensemble-based methods for characterizing a heterogeneous aquifer by assimilating experimental tracer data, with refined prior information obtained from assimilating other types of data available at the site. It is demonstrated that high-performance computing enables the use of increasingly mechanistic nonlinear forward simulations for a complex system within the data assimilation framework with reasonable turnaround time. [ABSTRACT FROM AUTHOR]

The response of borehole water levels to barometric pressure changes in semiconfined aquifers can be used to determine barometric response functions from which aquifer and confining layer properties can be obtained. Following earlier work on barometric response functions and aquifer confinement, we explore the barometric response function as a tool to improve the assessment of groundwater vulnerability in semiconfined aquifers, illustrated through records from two contrasting boreholes in the semiconfined Chalk Aquifer, East Yorkshire, UK. After removal of recharge and Earth tide influences on the water level signal, barometric response functions were estimated and aquifer and confining layer properties determined through an analytical model of borehole water level response to barometric pressure. A link between the thickness and vertical diffusivity of the confining layer determined from the barometric response function, and groundwater vulnerability is proposed. The amplitude spectrum for barometric pressure and instrument resolution favor determination of the barometric response function at frequencies to which confining layer diffusivities are most sensitive. Numerical modeling indicates that while the high frequency response reflects confining layer properties in the immediate vicinity of the borehole, the low frequency response reflects vertical, high diffusivity pathways though the confining layer some hundreds of meters distant. A characteristic time scale parameter, based on vertical diffusivities and thicknesses of the saturated and unsaturated confining layer, is introduced as a measure of semiconfined aquifer vulnerability. The study demonstrates that the barometric response function has potential as a tool for quantitative aquifer vulnerability assessment in semiconfined aquifers. [ABSTRACT FROM AUTHOR]

Although mechanistic reaction networks have been developed to quantify the biogeochemical evolution of subsurface systems associated with bioremediation, it is difficult in practice to quantify the onset and distribution of these transitions at the field scale using commonly collected wellbore datasets. As an alternative approach to the mechanistic methods, we develop a data-driven, statistical model to identify biogeochemical transitions using various time-lapse aqueous geochemical data (e.g., Fe(II), sulfate, sulfide, acetate, and uranium concentrations) and induced polarization (IP) data. We assume that the biogeochemical transitions can be classified as several dominant states that correspond to redox transitions and test the method at a uranium-contaminated site. The relationships between the geophysical observations and geochemical time series vary depending upon the unknown underlying redox status, which is modeled as a hidden Markov random field. We estimate unknown parameters by maximizing the joint likelihood function using the maximization-expectation algorithm. The case study results show that when considered together aqueous geochemical data and IP imaginary conductivity provide a key diagnostic signature of biogeochemical stages. The developed method provides useful information for evaluating the effectiveness of bioremediation, such as the probability of being in specific redox stages following biostimulation where desirable pathways (e.g., uranium removal) are more highly favored. The use of geophysical data in the approach advances the possibility of using noninvasive methods to monitor critical biogeochemical system stages and transitions remotely and over field relevant scales (e.g., from square meters to several hectares). [ABSTRACT FROM AUTHOR]