Search USGS⌕ Search

Geology topics

Allen M. Shapiro

Publications and source records attributed to Allen M. Shapiro.

At least 19 recordsLinked to original sources
USGS2023

Benefits and cautions in data assimilation strategies: An example of modeling groundwater recharge

Assimilating recent observations improves model outcomes for real-time assessments of groundwater processes. This is demonstrated in estimating time-varying recharge to a shallow fractured-rock aquifer in response to precipitation. Results from estimating the time-varying water-table altitude ( h ) and recharge, and their error covariances, are compared for forecasting, filtering, and fixed-lag smoothing (FLS), which are implemented using the Kalman Filter as applied to a data-driven, mechanistic model of recharge. Forecasting uses past observations to predict future states and is the current paradigm in most groundwater modeling investigations; filtering assimilates observations up to the current time to estimate current states; and FLS estimates states following a time lag over which additional observations are collected. Results for forecasting yield a large error covariance relative to the magnitude of the expected recharge. With assimilating recent observations of h , filtering and FLS produce estimates of recharge that better represent time-varying observations of h and reduce uncertainty in comparison to forecasting. Although model outcomes from applying data assimilation through filtering or FLS reduce model uncertainty, they are not necessarily mass conservative, whereas forecasting outcomes are mass conservative. Mass conservative outcomes from forecasting are not necessarily more accurate, because process errors are inherent in any model. Improvements in estimating real-time groundwater conditions that better represent observations need to be weighed for the model application against outcomes with inherent process deficiencies. Results from data assimilation strategies discussed in this investigation are anticipated to be relevant to other groundwater processes models where system states are sensitive to system inputs.

Groundwater↗
USGS2022

Mapping areas of groundwater susceptible to transient contamination events from rapid infiltration into shallow fractured-rock aquifers in agricultural regions of the conterminous United States

Current time-invariant groundwater vulnerability assessments may not capture intermittent contamination events in landscape areas that experience rapid infiltration following precipitation or snowmelt. Occurrences of rapid infiltration and intermittent degradation of groundwater quality are frequently reported in fractured-rock aquifers. This investigation identifies landscape areas underlain by fractured rock within the conterminous United States (CONUS) that may be susceptible to rapid infiltration and where groundwater is a principal source of water supply to the population. Our analysis shows that approximately 27 percent of the CONUS, corresponding to a population of approximately 150 million people, is both underlain by fractured rock and denoted as an area of significant groundwater use. The results of this survey identified shallow fractured-rock aquifers underlying glacial sediments in the upper Midwest and northeastern United States as areas that may be subject to rapid infiltration events. Additionally, aquifers associated with the early Mesozoic basins located in the northeastern and mid-Atlantic United States and bands of carbonate aquifers in the southeastern United States show high susceptibility to rapid infiltration. Index values used in this investigation indicate isolated areas in the western half of the United States also show high susceptibility to rapid infiltration. The isolated areas in Oklahoma, Texas, Arkansas, and southwestern Missouri correspond to karst regions of carbonate aquifers. The isolated areas showing high susceptibility to rapid infiltration and contamination from agricultural sources are locations where more detailed investigations of transient contamination events are warranted. This survey also addresses the potential for contaminant longevity in fractured-rock aquifers stemming from intermittent contamination events. Contaminants that can dissolve into the groundwater following infiltration may be introduced into fractures, and the dissolved constituents can diffuse from fractures into the porosity of the adjacent rock matrix. These constituents can then diffuse back into permeable fractures and adversely affect groundwater quality at downgradient locations over an extended time frame. Rock types with larger matrix porosities have the capacity to retain and then release larger quantities of dissolved constituents, resulting in longer residence times for dissolved groundwater contaminants. The magnitude of the dissolved contaminant concentration infiltrating to the water table will also dictate whether the contaminant concentration in the groundwater exceeds limits for human consumption over the duration of a contamination event. In general, sedimentary- and carbonate-rock aquifers have larger matrix porosities in comparison to igneous- and metamorphic-rock aquifers, and thus, they are more susceptible to longer contaminant residence times. Aquifers composed of sedimentary or carbonate rock constitute approximately 51 percent of the CONUS, and 19 percent of the CONUS is associated with sedimentary- or carbonate-rock aquifers that are of significance for groundwater use. Depending on the contaminants of concern and the concentration of the contaminants introduced into the groundwater from infiltrating water, it would be beneficial for investigations of susceptibility to rapid infiltration to also consider the potential for contaminant longevity. This investigation identifies areas of rapid infiltration into fractured rock using index values applied to the attributes (1) depth to the water table, (2) depth to bedrock, and (3) percentage of sand in soil, where larger index values indicate a greater susceptibility to rapid infiltration. These attributes are selected as the most likely factors that affect rapid infiltration to the water table. The combination of depth to water table and depth to bedrock highlight those aquifer settings that are characterized as shallow fractured-rock aquifers, where the water table may reside either in the bedrock or in overlying unconsolidated geologic materials. In addition, we consider the percentage of agricultural use as a land-use attribute when formulating an index of susceptibility to rapid infiltration and contamination. Agricultural areas are well recognized as nonpoint sources of contaminants that can affect groundwater quality because of seasonal amendments applied to the land surface. Rural agricultural areas are also characterized by septic tanks and leach fields for onsite treatment of wastewater, which may also be a source of contamination that may be introduced into the groundwater following precipitation or snowmelt events.

Conterminous United States↗
USGS2022

Application of recursive estimation to heat tracing for groundwater/surface-water exchange

We present and demonstrate a recursive-estimation framework to infer groundwater/surface-water exchange based on temperature time series collected at different vertical depths below the sediment/water interface. We formulate the heat-transport problem as a state-space model (SSM), in which the spatial derivatives in the convection/conduction equation are approximated using finite differences. The SSM is calibrated to estimate time-varying specific discharge using the Extended Kalman Filter (EKF) and Extended Rauch-Tung-Striebel Smoother (ERTSS). Whereas the EKF is suited to real-time (“online”) applications and uses only the past and current measurements for estimation (filtering), the ERTSS is intended for near-real time or batch-processing (“offline”) applications and uses a window of data for batch estimation (smoothing). The two algorithms are demonstrated with synthetic and field-experimental data and are shown to be efficient and rapid for the estimation of time-varying flux over seasonal periods; further, the recursive approaches are effective in the presence of rapidly changing flux and (or) nonperiodic thermal boundary conditions, both of which are problematic for existing approaches to heat tracing of time-varying groundwater/surface-water exchange.

Water Resources Research↗
USGS2022

Incorporating snowmelt into daily estimates of recharge using a state-space model of infiltration

A state-space model (SSM) of infiltration estimates daily groundwater recharge using time-series of groundwater-level altitude and meteorological inputs (liquid precipitation, snowmelt, and evapotranspiration). The model includes diffuse and preferential flow through the unsaturated zone, where preferential flow is a function of liquid precipitation and snowmelt rates and a threshold rate, above which there is direct recharge to the water table. Model parameters are estimated over seasonal periods and the SSM is coupled with the Kalman Filter (KF) to assimilate recent observations (hydraulic head) and meteorological inputs into recharge estimates. The approach can take advantage of real-time hydrologic and meteorological data to deliver real-time recharge estimates. The model is demonstrated on daily observations from two bedrock wells in carbonate aquifers of northwestern New York (USA) between 2013 and 2018. Meteorological inputs for liquid precipitation and snowmelt are compiled from SNODAS (2021). Results for recharge during winter and spring seasons show preferential flow events to the water table from liquid precipitation, snowmelt, or a combination of the two. Recharge estimates summed annually are consistent with previous estimates of recharge reported from groundwater flow and surface-process models. Results from the SSM and KF point to errors in meteorological inputs, such as the snowmelt rate, that are not compatible with hydraulic head observations. Whereas liquid and solid precipitation are measured at discrete stations and extrapolated to 1-km 2 grid cells, snowmelt is a meteorological modeled outcome that may not represent conditions in the vicinity of monitoring well locations.

Groundwater↗
USGS2022

Contributing areas to domestic wells in dipping sedimentary rocks under extreme recharge events

We use particle tracking to determine contributing areas (CAs) to wells for transient flow models that simulate cyclic domestic pumping and extreme recharge events in a small synthetic watershed underlain by dipping sedimentary rocks. The CAs consist of strike-oriented bands at locations where the water table intersects high-hydraulic conductivity beds, and from which groundwater flows to the pumping well. Factors that affect the size and location of the CAs include topographic flow directions, rock dip direction, cross-bed fracture density, and position of the well relative to streams. For an effective fracture porosity (n e ) of 10 -4 , the fastest advective travel times from CAs to wells are only a few hours. These results indicate that wells in this type of geologic setting can be highly vulnerable to contaminants or pathogens flushed into the subsurface during extreme recharge events. Increasing n e to 10 -3 results in modestly smaller CAs and delayed well vulnerability due to slower travel times. CAs determined for steady-state models of the same setting, but with long-term average recharge and pumping rates, are smaller than CAs in the models with extreme recharge. Also, the earliest-arriving particles arrive at the wells later in the steady-state models than in the extreme-recharge models. The results highlight the importance of characterizing geologic structure, simulating plausible effective porosities, and simulating pumping and recharge transience when determining CAs in fractured rock aquifers to assess well vulnerability under extreme precipitation events.

Groundwater↗
USGS2021

Estimating and forecasting time-varying groundwater recharge in fractured rock: A state-space formulation with preferential and diffuse flow to the water table

Rapid infiltration following precipitation may result in groundwater contamination from surface contaminants or pathogens. In fractured rock, contaminants can migrate rapidly to points of groundwater withdrawals. In contrast to the temporal availability of groundwater quality chemical indicators, meteorological and groundwater level observations are available in real-time to estimate time-varying recharge, which can act as a surrogate to identify periods of rapid infiltration that may indicate contamination susceptibility. Estimating recharge using methods, such as base-flow recession, unsaturated infiltration models, or Water-Table Fluctuations (WTF), cannot capitalize on currently available technologies and telecommunication infrastructure to conduct real-time recharge estimation at scales relevant to characterizing rapid infiltration. We present a linear, physics-based State-Space (SS) model of one-dimensional infiltration to estimate recharge, which includes preferential and diffuse-flow to the water table. The model can take advantage of real-time data for water-table altitude, precipitation, and evapotranspiration. Model parameters are calibrated over an observation period, and the Kalman Filter (KF) is subsequently applied to continuously update the observed (water-table altitude) and unobserved (groundwater recharge) system states and predict future states as new data become available. The SS/KF algorithm is demonstrated at the Masser Groundwater Recharge Site in Pennsylvania, USA and comparisons are made with recharge estimates from WTF methods. Model results indicate that the frequency of observations (daily versus sub-daily) dictates the allocation between preferential and diffuse flow. Additionally, because infiltration processes encompass many nonlinearities, model parameters estimated from observation periods need to be updated at least seasonally to account for changing recharge conditions.

Water Resources Research↗
USGS2019

The complex spatial distribution of trichloroethene and the probability of NAPL occurrence in the rock matrix of a mudstone aquifer

Methanol extractions for chloroethene analyses are conducted on rock samples from seven closely spaced coreholes in a mudstone aquifer that was subject to releases of the nonaqueous phase liquid (NAPL) form of trichloroethene (TCE) between the 1950's and 1990's. Although TCE concentration in the rock matrix over the length of coreholes is dictated by proximity to subhorizontal bedding planefractures, elevated TCE concentrations in the rock matrix are not continuous along the most permeable bedding plane fractures. A complex configuration of subvertical and subhorizontal fractures appears to be responsible for the TCE distribution from prior TCE releases at land surface. Phase partitioning calculations of TCE in the rock matrix show that most TCE is adsorbed to solid surfaces because of the large fraction of organic carbon ( f oc ) in the mudstone. Large TCE content in some cores indicate the likely presence of the NAPL form of TCE in the rock matrix. Using average values of porosity ( n ) and f oc in phase partitioning calculations identifies a number of locations of possible NAPL occurrence in the rock matrix. Samples of mudstone analyzed for n and f oc show variability in these properties over several orders of magnitude. Accounting for this variability in phase partitioning calculations identifies a probability of NAPL occurrence, P NAPL . The spatial variability of P NAPL along coreholes identifies a configuration that may be attributed to a TCE source zone that has evolved after emplacement due to NAPL dissolution, adsorption, and matrix diffusion.

New Jersey↗
USGS2018

Variability of organic carbon content and the retention and release of trichloroethene in the rock matrix of a mudstone aquifer

Contaminants diffusing from fractures into the immobile porosity of the rock matrix are subject to prolonged residence times. Organic contaminants can adsorb onto organic carbonaceous materials in the matrix extending contaminant retention. An investigation of spatial variability of the fraction of organic carbon ( f oc ) is conducted on samples of rock core from seven closely spaced boreholes in a mudstone aquifer contaminated with trichloroethene (TCE). A total of 378 samples were analyzed at depths between 14 and 36 m below land surface. Mudstone units associated with deep water deposition have the largest f oc , with a maximum value of 0.0396, and units associated with shallow water deposition have the smallest f oc . Even though f oc correlates with depositional conditions, f oc still varies over more than an order of magnitude in continuous mudstone layers between boreholes, and there is large variability in f oc over short distances perpendicular to bedding. Simulations of diffusion and linear equilibrium adsorption of TCE using spatially variable f oc in the rock matrix show order of magnitude variability in the adsorbed TCE over short distances in the matrix and residence times extending to hundreds of years following remediation in adjacent fractures. Simulations using average values of f oc do not capture the range of TCE mass that can be retained in a rock matrix characterized by spatially variable f oc . Bounds on TCE mass within the rock matrix can be obtained by simulations with spatially uniform values of f oc equal to the maximum and minimum values of f oc for a given mudstone unit.

New Jersey↗
USGS2017

Bioremediation in fractured rock: 1. Modeling to inform design, monitoring, and expectations

Field characterization of a trichloroethene (TCE) source area in fractured mudstones produced a detailed understanding of the geology, contaminant distribution in fractures and the rock matrix, and hydraulic and transport properties. Groundwater flow and chemical transport modeling that synthesized the field characterization information proved critical for designing bioremediation of the source area. The planned bioremediation involved injecting emulsified vegetable oil and bacteria to enhance the naturally occurring biodegradation of TCE. The flow and transport modeling showed that injection will spread amendments widely over a zone of lower‐permeability fractures, with long residence times expected because of small velocities after injection and sorption of emulsified vegetable oil onto solids. Amendments transported out of this zone will be diluted by groundwater flux from other areas, limiting bioremediation effectiveness downgradient. At nearby pumping wells, further dilution is expected to make bioremediation effects undetectable in the pumped water. The results emphasize that in fracture‐dominated flow regimes, the extent of injected amendments cannot be conceptualized using simple homogeneous models of groundwater flow commonly adopted to design injections in unconsolidated porous media (e.g., radial diverging or dipole flow regimes). Instead, it is important to synthesize site characterization information using a groundwater flow model that includes discrete features representing high‐ and low‐permeability fractures. This type of model accounts for the highly heterogeneous hydraulic conductivity and groundwater fluxes in fractured‐rock aquifers, and facilitates designing injection strategies that target specific volumes of the aquifer and maximize the distribution of amendments over these volumes.

Groundwater↗
USGS2017

Bioremediation in fractured rock: 2. Mobilization of chloroethene compounds from the rock matrix

A mass balance is formulated to evaluate the mobilization of chlorinated ethene compounds (CE) from the rock matrix of a fractured mudstone aquifer under pre- and postbioremediation conditions. The analysis relies on a sparse number of monitoring locations and is constrained by a detailed description of the groundwater flow regime. Groundwater flow modeling developed under the site characterization identified groundwater fluxes to formulate the CE mass balance in the rock volume exposed to the injected remediation amendments. Differences in the CE fluxes into and out of the rock volume identify the total CE mobilized from diffusion, desorption, and nonaqueous phase liquid dissolution under pre- and postinjection conditions. The initial CE mass in the rock matrix prior to remediation is estimated using analyses of CE in rock core. The CE mass mobilized per year under preinjection conditions is small relative to the total CE mass in the rock, indicating that current pump-and-treat and natural attenuation conditions are likely to require hundreds of years to achieve groundwater concentrations that meet regulatory guidelines. The postinjection CE mobilization rate increased by approximately an order of magnitude over the 5 years of monitoring after the amendment injection. This rate is likely to decrease and additional remediation applications over several decades would still be needed to reduce CE mass in the rock matrix to levels where groundwater concentrations in fractures achieve regulatory standards.

New Jersey↗
USGS2017

Porosity and pore size distribution in a sedimentary rock: Implications for the distribution of chlorinated solvents

Characterizing properties of the rock matrix that control retention and release of chlorinated solvents is essential in evaluating the extent of contamination and the application of remediation technologies in fractured rock. Core samples from seven closely spaced boreholes in a mudstone subject to trichloroethene (TCE) contamination were analyzed using Mercury Intrusion Porosimetry to investigate porosity and pore size distribution as a function of mudstone characteristics, and depth and lateral extent in the aquifer; organic carbon content was also evaluated to identify the potential for adsorption. Porosity and retardation factor varied over two orders of magnitude, with the largest porosities and largest retardation factors associated with carbon-rich mudstone layers. Larger porosities were also measured in the shallow rock that has been subject to enhanced groundwater flow. Porosity also varied over more than an order of magnitude in spatially continuous mudstone layers. The analyses of the rock cores indicated that the largest pore diameters may be accessible to entry of the nonaqueous form of TCE. Although the porosity associated with the largest pore diameters is small (~ 0.1%), that volume of TCE can significantly affect the total TCE that is retained in the rock matrix. The dimensions of the largest pore diameters may also be accessible to microbes responsible for reductive dechlorination; however, the small percentage of the pore space that can accommodate microbes may limit the extent of reductive dechlorination in the rock matrix.

Journal of Contaminant Hydrology↗
USGS2015

Interpretation of hydraulic conductivity in a fractured-rock aquifer over increasingly larger length dimensions

A comparison of the hydraulic conductivity over increasingly larger volumes of crystalline rock was conducted in the Piedmont physiographic region near Bethesda, Maryland, USA. Fluid-injection tests were conducted on intervals of boreholes isolating closely spaced fractures. Single-hole tests were conducted by pumping in open boreholes for approximately 30 min, and an interference test was conducted by pumping a single borehole over 3 days while monitoring nearby boreholes. An estimate of the hydraulic conductivity of the rock over hundreds of meters was inferred from simulating groundwater inflow into a kilometer-long section of a Washington Metropolitan Area Transit Authority tunnel in the study area, and a groundwater modeling investigation over the Rock Creek watershed provided an estimate of the hydraulic conductivity over kilometers. The majority of groundwater flow is confined to relatively few fractures at a given location. Boreholes installed to depths of approximately 50 m have one or two highly transmissive fractures; the transmissivity of the remaining fractures ranges over five orders of magnitude. Estimates of hydraulic conductivity over increasingly larger rock volumes varied by less than half an order of magnitude. While many investigations point to increasing hydraulic conductivity as a function of the measurement scale, a comparison with selected investigations shows that the effective hydraulic conductivity estimated over larger volumes of rock can either increase, decrease, or remain stable as a function of the measurement scale. Caution needs to be exhibited in characterizing effective hydraulic properties in fractured rock for the purposes of groundwater management.

Hydrogeology Journal↗
USGS2015

Imaging pathways in fractured rock using three-dimensional electrical resistivity tomography

Major challenges exist in delineating bedrock fracture zones because these cause abrupt changes in geological and hydrogeological properties over small distances. Borehole observations cannot sufficiently capture heterogeneity in these systems. Geophysical techniques offer the potential to image properties and processes in between boreholes. We used three-dimensional cross borehole electrical resistivity tomography (ERT) in a 9 m (diameter) × 15 m well field to capture high-resolution flow and transport processes in a fractured mudstone contaminated by chlorinated solvents, primarily trichloroethylene. Conductive (sodium bromide) and resistive (deionized water) injections were monitored in seven boreholes. Electrode arrays with isolation packers and fluid sampling ports were designed to enable acquisition of ERT measurements during pulsed tracer injections. Fracture zone locations and hydraulic pathways inferred from hydraulic head drawdown data were compared with electrical conductivity distributions from ERT measurements. Static ERT imaging has limited resolution to decipher individual fractures; however, these images showed alternating conductive and resistive zones, consistent with alternating laminated and massive mudstone units at the site. Tracer evolution and migration was clearly revealed in time-lapse ERT images and supported by in situ borehole vertical apparent conductivity profiles collected during the pulsed tracer test. While water samples provided important local information at the extraction borehole, ERT delineated tracer migration over spatial scales capturing the primary hydrogeological heterogeneity controlling flow and transport. The fate of these tracer injections at this scale could not have been quantified using borehole logging and/or borehole sampling methods alone.

New Jersey↗
USGS2014

Integration of stable carbon isotope, microbial community, dissolved hydrogen gas, and 2 H H 2 O tracer data to assess bioaugmentation for chlorinated ethene degradation in fractured rocks

An in situ bioaugmentation (BA) experiment was conducted to understand processes controlling microbial dechlorination of trichloroethene (TCE) in groundwater at the Naval Air Warfare Center (NAWC), West Trenton, NJ. In the BA experiment, an electron donor (emulsified vegetable oil and sodium lactate) and a chloro-respiring microbial consortium were injected into a well in fractured mudstone of Triassic age. Water enriched in 2 H was also injected as a tracer of the BA solution, to monitor advective transport processes. The changes in concentration and the δ 13 C of TCE, cis-dichloroethene (cis-DCE), and vinyl chloride (VC); the δ 2 H of water; changes in the abundance of the microbial communities; and the concentration of dissolved H2 gas compared to pre- test conditions, provided multiple lines of evidence that enhanced biodegradation occurred in the injection well and in two downgradient wells. For those wells where the biodegradation was stimulated intensively, the sum of the molar chlorinated ethene (CE) concentrations in post-BA water was higher than that of the sum of the pre-BA background molar CE concentrations. The concentration ratios of TCE/(cis-DCE + VC) indicated that the increase in molar CE concentration may result from additional TCE mobilized from the rock matrix in response to the oil injection or due to desorption/diffusion. The stable carbon isotope mass-balance calculations show that the weighted average 13 C isotope of the CEs was enriched for around a year compared to the background value in a two year monitoring period, an effective indication that dechlorination of VC was occurring. Insights gained from this study can be applied to efforts to use BA in other fractured rock systems. The study demonstrates that a BA approach can substantially enhance in situ bioremediation not only in fractures connected to the injection well, but also in the rock matrix around the well due to processes such as diffusion and desorption. Because the effect of the BA was intensive only in wells where an amendment was distributed during injection, it is necessary to adequately distribute the amendments throughout the fractured rock to achieve substantial bioremediation. The slowdown in BA effect after a year is due to some extend to the decrease abundant of appropriate microbes, but more likely the decreased concentration of electron donor.

New Jersey↗
USGS2013

Physical, chemical, and isotopic data from groundwater in the watershed of Mirror Lake, and in the vicinity of Hubbard Brook, near West Thornton, New Hampshire, 1983 to 1997

Research on the hydrogeologic setting of Mirror Lake near West Thornton, New Hampshire (43° 56.5’ N, 71° 41.5’ W), includes the study of the physical, chemical, and isotopic characteristics of groundwater in the vicinity of the lake and nearby Hubbard Brook. Presented here are those physical, chemical, and isotopic data for the period 1983 to 1997. Data were collected from observation wells installed in glacial drift and bedrock, as well as from one domestic well in the general area of the lake and Hubbard Brook. Also presented are data for Mirror Lake for August 1, 1993, to place groundwater data in context with chemical and isotopic characteristics of the lake.

New Hampshire↗