David W. Hyndman

Impacts of the 2004 tsunami on groundwater resources in Sri Lanka

The 26 December 2004 tsunami caused widespread destruction and contamination of coastal aquifers across southern Asia. Seawater filled domestic open dug wells and also entered the aquifers via direct infiltration during the first flooding waves and later as ponded seawater infiltrated through the permeable sands that are typical of coastal aquifers. In Sri Lanka alone, it is estimated that over 40,000 drinking water wells were either destroyed or contaminated. From February through September 2005, a team of United States, Sri Lankan, and Danish water resource scientists and engineers surveyed the coastal groundwater resources of Sri Lanka to develop an understanding of the impacts of the tsunami and to provide recommendations for the future of coastal water resources in south Asia. In the tsunami-affected areas, seawater was found to have infiltrated and mixed with fresh groundwater lenses as indicated by the elevated groundwater salinity levels. Seawater infiltrated through the shallow vadose zone as well as entered aquifers directly through flooded open wells. Our preliminary transport analysis demonstrates that the intruded seawater has vertically mixed in the aquifers because of both forced and free convection. Widespread pumping of wells to remove seawater was effective in some areas, but overpumping has led to upconing of the saltwater interface and rising salinity. We estimate that groundwater recharge from several monsoon seasons will reduce salinity of many sandy Sri Lankan coastal aquifers. However, the continued sustainability of these small and fragile aquifers for potable water will be difficult because of the rapid growth of human activities that results in more intensive groundwater pumping and increased pollution. Long-term sustainability of coastal aquifers is also impacted by the decrease in sand replenishment of the beaches due to sand mining and erosion.

Inferring the relation between seismic slowness and hydraulic conductivity in heterogeneous aquifers

Cross-well seismic tomography can be used to develop high-resolution seismic slowness (1/velocity) estimates along planes through aquifers. Unfortunately, the relation between seismic slowness and hydraulic conductivity is poorly understood, resulting in poor characterization of hydraulic properties from seismic data. This relation is generally developed from laboratory measurements, but slowness values measured with very high frequencies in the lab are often poorly correlated with lower frequency cross-well and surface seismic slowness values. To address this problem, we developed an approach to infer the relation between slowness and hydraulic conductivity using field scale geophysical and hydrogeologic measurements. We first develop an a priori relation between the conductivity measurements and the cross-well slowness estimates. Multiple three- dimensional slowness realizations, conditioned on the cross-well estimates, are then generated and remapped into log conductivity fields using the a priori slowness to log conductivity relation. We simulate groundwater flow and tracer transport through these conductivity fields and calculate the residuals between measured and simulated concentration arrival time quantiles and drawdown. The slope and intercept of the relation between slowness and log hydraulic conductivity and the dispersivity are then estimated for each slowness realization to minimize the sum of these squared residuals. We demonstrate this approach for the Kesterson aquifer, California, where seismic tomography provided valuable information about aquifer properties. The groundwater flow and tracer transport simulations, through the estimated conductivity fields, yield reasonable fits to the observed tracer concentration histories for two multiple-well tracer tests (one of which was not used in the inversion) and to the measured drawdown. This approach provides estimates of seismic slowness and hydraulic conductivity, and information about the relation between slowness and log conductivity for a field site.

Temporal variations in parameters reflecting terminal-electron-accepting processes in an aquifer contaminated with waste fuel and chlorinated solvents

Abstract A fundamental issue in aquifer biogeochemistry is the means by which solute transport, geochemical processes, and microbiological activity combine to produce spatial and temporal variations in redox zonation. In this paper, we describe the temporal variability of TEAP conditions in shallow groundwater contaminated with both waste fuel and chlorinated solvents. TEAP parameters (including methane, dissolved iron, and dissolved hydrogen) were measured to characterize the contaminant plume over a 3-year period. We observed that concentrations of TEAP parameters changed on different time scales and appear to be related, in part, to recharge events. Changes in all TEAP parameters were observed on short time scales (months), and over a longer 3-year period. The results indicate that (1) interpretations of TEAP conditions in aquifers contaminated with a variety of organic chemicals, such as those with petroleum hydrocarbons and chlorinated solvents, must consider additional hydrogen-consuming reactions (e.g., dehalogenation); (2) interpretations must consider the roles of both in situ (at the sampling point) biogeochemical and solute transport processes; and (3) determinations of microbial communities are often necessary to confirm the interpretations made from geochemical and hydrogeological measurements on these processes.

Subsurface imaging of vegetation, climate, and root-zone moisture interactions

Changes in global climate and land use affect important prolesses from evapotranspiration and groundwater recharge to carbon storage and biochemical cycling. Near surface soil moisture is pivotal to understand the consequences of these changes. However, the dynamic interactions between vegetation and soil moisture remain largely unresolved because it is difficult to monitor and quantify subsurface hydrologic fluxes at relevant scales. Here we use electrical resistivity to monitor the influence of climate and vegetation on root-zone moisture, bridging the gap between remotely-sensed and in-situ point measurements. Our research quantifies large seasonal differences in root-zone moisture dynamics for a forest-grassland ecotone. We found large differences in effective rooting depth and moisture distributions for the two vegetation types. Our results highlight the likely impacts of land transformations on groun ter recharge, streamflow, and land-atmosphere exchanges.

Hydrological consequences of land-cover change: Quantifying the influence of plants on soil moisture with time-lapse electrical resistivity

Electrical resistivity of soils and sediments is strongly influenced by the presence of interstitial water. Taking advantage of this dependency, electrical-resistivity imaging (ERI) can be effectively utilized to estimate subsurface soil-moisture distributions. The ability to obtain spatially extensive data combined with time-lapse measurements provides further opportunities to understand links between land use and climate processes. In natural settings, spatial and temporal changes in temperature and porewater salinity influence the relationship between soil moisture and electrical resistivity. Apart from environmental factors, technical, theoretical, and methodological ambiguities may also interfere with accurate estimation of soil moisture from ERI data. We have examined several of these complicating factors using data from a two-year study at a forest-grassland ecotone, a boundary between neighboring but different plant communities.At this site, temperature variability accounts for approximately 20%–4...

Traveltime inversion for the geometry of aquifer lithologies

Crosswell traveltime tomography can provide detailed descriptions of the geometry and seismic slowness of lithologic zones in aquifers and reservoirs. Traditional tomographic inversions that estimate a smooth slowness field to match traveltime data, provide limited information about the dominant scale of subsurface heterogeneity. We demonstrate an alternative method, called the multiple population inversion (MPI), that co-inverts traveltimes between multiple well pairs to identify the spatial distribution of a small number of slowness populations. We also compare the MPI with the split inversion method (SIM) that was recently introduced to address the same problem. The lithologies and hydraulic parameters for these populations can then be determined from core data and hydraulic testing. The MPI iteratively assigns pixels to a small number of slowness populations based on the histogram of slowness residuals. By constraining the number of slowness values, this method is less susceptible to inversion artifacts, such as those related to slight variations in ray coverage, and can resolve finer scale sedimentary structures better than methods that smooth the slowness field. We demonstrate the MPI in two dimensions with a synthetic aquifer and in three dimensions with the Kesterson aquifer in the central valley of California. In both cases, the constrained inversion algorithm converges to an equal or smaller average traveltime residual than obtained with unconstrained-value tomography. The MPI accurately images the dominant lithologies of the synthetic aquifer and provides a geologically reasonable image of the Kesterson aquifer.

Examining the influence of heterogeneous porosity fields on conservative solute transport.

It is widely recognized that groundwater flow and solute transport in natural media are largely controlled by heterogeneities. In the last three decades, many studies have examined the effects of heterogeneous hydraulic conductivity fields on flow and transport processes, but there has been much less attention to the influence of heterogeneous porosity fields. In this study, we use porosity and particle size measurements from boreholes at the Boise Hydrogeophysical Research Site (BHRS) to evaluate the importance of characterizing the spatial structure of porosity and grain size data for solute transport modeling. Then we develop synthetic hydraulic conductivity fields based on relatively simple measurements of porosity from borehole logs and grain size distributions from core samples to examine and compare the characteristics of tracer transport through these fields with and without inclusion of porosity heterogeneity. In particular, we develop horizontal 2D realizations based on data from one of the less heterogeneous units at the BHRS to examine effects where spatial variations in hydraulic parameters are not large. The results indicate that the distributions of porosity and the derived hydraulic conductivity in the study unit resemble fractal normal and lognormal fields respectively. We numerically simulate solute transport in stochastic fields and find that spatial variations in porosity have significant effects on the spread of an injected tracer plume including a significant delay in simulated tracer concentration histories.

Improved methods for satellite‐based groundwater storage estimates: A decade of monitoring the high plains aquifer from space and ground observations

The impacts of climate extremes and water use on groundwater storage across large aquifers can be quantified using Gravity Recovery and Climate Experiment (GRACE) satellite monitoring. We present new methods to improve estimates of changes in groundwater storage by incorporating irrigation soil moisture corrections to common data assimilation products. These methods are demonstrated using data from the High Plains Aquifer (HPA) for 2003 to 2013. Accounting for the impacts of observed and inferred irrigation on soil moisture significantly improves estimates of groundwater storage changes as verified by interpolated measurements from ~10,000 HPA wells. The resulting estimates show persistent declines in groundwater storage across the HPA, more severe in the southern and central HPA than in the north. Groundwater levels declined by an average of approximately 276 ± 23 mm from 2003 to 2013, resulting in a storage loss of 125 ± 4.3 km3, based on the most accurate of the three methods developed here.

Hydrogeophysical Case Studies at the Local Scale: The Saturated Zone

Modern geophysical methods provide significant promise for estimating subsurface aquifer properties of the saturated zone in a minimally invasive manner. Mapping aquifer boundaries and internal stratification, estimating spatial distribution of hydrogeologic parameters, or monitoring tracer and contaminant plumes are examples of geophysical tools successfully applied. A general benefit of geophysical methods is the ability to collect high-resolution data in the horizontal dimension, where core/borehole data is nearly always limited. The complementary nature of core/borehole data and 2-D or 3-D geophysical data promises to help improve the accuracy and resolution of aquifer characterization at a variety of scales. However, one significant remaining difficulty is transforming geophysical parameters into flow and transport properties. A series of approaches and petrophysical models have been developed to help in this transformation (e.g., see Chapter 4 and Chapter 9 of this volume), but often complex and non-unique parameter relationships complicate data analysis and interpretation.

Water Level Declines in the High Plains Aquifer: Predevelopment to Resource Senescence.

A large imbalance between recharge and water withdrawal has caused vital regions of the High Plains Aquifer (HPA) to experience significant declines in storage. A new predevelopment map coupled with a synthesis of annual water levels demonstrates that aquifer storage has declined by approximately 410 km(3) since the 1930s, a 15% larger decline than previous estimates. If current rates of decline continue, much of the Southern High Plains and parts of the Central High Plains will have insufficient water for irrigation within the next 20 to 30 years, whereas most of the Northern High Plains will experience little change in storage. In the western parts of the Central and northern part of the Southern High Plains, saturated thickness has locally declined by more than 50%, and is currently declining at rates of 10% to 20% of initial thickness per decade. The most agriculturally productive portions of the High Plains will not support irrigated production within a matter of decades without significant changes in management.