We present a new field measurement and numerical interpretation method (combined termed "test") to parameterize the diffusion of trichloroethene (TCE) and its biodegradation products (DPs) from the matrix of sedimentary rock. The method uses a dual-packer system to interrogate a low-permeability section of the rock matrix adjacent to a previously contaminated borehole and uses the borehole monitoring history to establish the pretest condition. TCE and its DPs are removed from the groundwater between the packers at the onset of the testing. The parameters estimated by fitting a radial diffusion model to the concentration history and borehole concentration data, also termed back diffusion, are the tortuosity factor and sorption coefficients of TCE and DPs in the rock matrix and the TCE and DP biodegradation rate coefficients in the borehole. We demonstrate the equipment design and the interpretive method using a borehole accessing the gray mudstone at a TCE contaminated site in the Newark Basin. In this test, both nonreactive (bromide) and reactive (trichlorofluoroethene) tracers are used to constrain the estimated parameters; however, the bromide tracer was not needed to estimate the parameters in this test. The parameters estimated from the field test are consistent with values measured independently in laboratory experiments using field samples of similar lithology. From the interpretation, we compute the TCE and DP concentration distributions in the rock matrix prior to the test to illustrate how the results can be used to enhance understanding of contaminant distribution in the rock matrix.

First posted February 6, 2020 For additional information, contact: Director, Pennsylvania Water Science CenterU.S. Geological Survey215 Limekiln RoadNew Cumberland, PA 17070 In 2014, groundwater samples from residential and public supply wells in the vicinity of two former U.S. Navy bases at Willow Grove and Warminster, and an active Air National Guard Station at Horsham, Bucks and Montgomery Counties, Pennsylvania, were found to have concentrations of perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS), which are per- and polyfluoroalkyl substances (PFAS), above U.S. Environmental Protection Agency (EPA) provisional health advisory (HA) levels for drinking water. Five supply wells near the bases were shut down because of PFAS contamination. In 2016, after EPA established a Lifetime HA for PFAS in drinking water that is lower than the provisional HA in place in 2014, at least 13 additional supply wells near the bases were shut down because of PFAS contamination. At the request of the U.S. Navy, and in consultation with other Federal and State agencies and local stakeholders, the U.S. Geological Survey used historical and recent data on well withdrawals, recharge rates, aquifer properties, groundwater levels, and stream base flow to evaluate regional groundwater-flow paths from identified areas of PFAS groundwater contamination or potential PFAS sources at the bases. Groundwater withdrawals near the bases from public supply and other large wells decreased substantially from the 1990s to 2017, increasing the proportion of groundwater recharge that discharged to local streams. A preliminary groundwater-flow model, calibrated using 1,009 groundwater levels and 17 stream base flow estimates, simulated regional flow paths from the bases and showed that recharge at the bases discharged to withdrawal wells and local streams, generally within a mile or two of the bases. Supply and remediation wells at the bases captured some of the recharge on base areas of possible PFAS contamination, whereas other base recharge was simulated to flow to nearby public supply wells and streams, depending on water use and aquifer recharge conditions between 1999 and 2017. The locations of many residential wells near the bases that were identified by the Navy and Air National Guard as having elevated PFAS concentrations were generally consistent with the simulated flow paths from possible sources at the bases. However, there are some areas of observed PFAS contamination where no flow paths from base sources were simulated. Additionally, no data were available on PFAS concentrations in groundwater in some areas of simulated flow paths from base sources. Data and models used for this study are provided in this report and in digital data releases to support further investigations and model revisions.

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 plane fractures, 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 (foc) 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 foc in phase partitioning calculations identifies a number of locations of possible NAPL occurrence in the rock matrix. Samples of mudstone analyzed for n and foc show variability in these properties over several orders of magnitude. Accounting for this variability in phase partitioning calculations identifies a probability of NAPL occurrence, PNAPL. The spatial variability of PNAPL 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.

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.

First posted December 21, 2017 For additional information, contact: Director, Integrated Modeling and Prediction Division U.S. Geological Survey MS 415 National Center 12201 Sunrise Valley Drive Reston, VA 20192Contact Pubs Warehouse Generating an informative display of groundwater withdrawals can sometimes be difficult because the symbols for closely spaced wells can overlap. An alternative method for displaying groundwater withdrawals is to generate a “footprint” of the withdrawals. WellFootprint version 1.0 implements the Footprint algorithm with two optional variations that can speed up the footprint calculation. ModelMuse has been modified in order to generate the input for WellFootprint and to read and graphically display the output from WellFootprint.

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.

First posted June 6, 2017 For additional information, contact: Director, Pennsylvania Water Science Center U.S. Geological Survey 215 Limekiln Road New Cumberland, PA 17070 A previously developed regional groundwater flow model was used to simulate the effects of changes in pumping rates on groundwater-flow paths and extent of recharge discharging to wells for a contaminated fractured bedrock aquifer in southeastern Pennsylvania. Groundwater in the vicinity of the North Penn Area 7 Superfund site, Montgomery County, Pennsylvania, was found to be contaminated with organic compounds, such as trichloroethylene (TCE), in 1979. At the time contamination was discovered, groundwater from the underlying fractured bedrock (shale) aquifer was the main source of supply for public drinking water and industrial use. As part of technical support to the U.S. Environmental Protection Agency (EPA) during the Remedial Investigation of the North Penn Area 7 Superfund site from 2000 to 2005, the U.S. Geological Survey (USGS) developed a model of regional groundwater flow to describe changes in groundwater flow and contaminant directions as a result of changes in pumping. Subsequently, large decreases in TCE concentrations (as much as 400 micrograms per liter) were measured in groundwater samples collected by the EPA from selected wells in 2010 compared to 2005‒06 concentrations.To provide insight on the fate of potentially contaminated groundwater during the period of generally decreasing pumping rates from 1990 to 2010, steady-state simulations were run using the previously developed groundwater-flow model for two conditions prior to extensive remediation, 1990 and 2000, two conditions subsequent to some remediation 2005 and 2010, and a No Pumping case, representing pre-development or cessation of pumping conditions. The model was used to (1) quantify the amount of recharge, including potentially contaminated recharge from sources near the land surface, that discharged to wells or streams and (2) delineate the areas contributing recharge that discharged to wells or streams for the five conditions.In all simulations, groundwater divides differed from surface-water divides, partly because of differences in stream elevations and because of geologic structure and pumping. In the 1990 and 2000 simulations, all recharge in and near the vicinity of North Penn Area 7 discharged to wells, but in the 2005 and 2010 simulations some recharge in this area discharged to streams, indicating possible discharge of contaminated groundwater from North Penn Area 7 sources to streams. As the amount of groundwater withdrawals by wells has declined since 1990, the area contributing recharge to wells in the vicinity of North Penn Area 7 has decreased.To determine the effect of changes in pumping on flow paths and possible flow-path-related contributions to the observed changes in spatial distribution of contaminants in groundwater from 2005 to 2010, the USGS conducted simulations using the previously developed regional groundwater-flow model using reported pumping and estimated recharge rates for 2005 and 2010. Flow paths from recharge at known contaminant source areas to discharge locations at wells or streams were simulated under steady-state conditions for the two periods. Simulated groundwater-flow paths shifted only slightly from 2005 to 2010 as a result of changes in pumping rates. These slight changes in groundwater-flow paths from known sources of contamination are not coincident with the spatial distribution of observed changes in TCE concentrations from 2005 to 2010, indicating that the decreases of TCE concentrations may be a result of other processes, such as source removal or degradation. Results of the simulations and the absence of increases in TCE-degradation-product concentrations indicate that the decreases of TCE concentrations observed in 2010 may be at least partly related to contaminant-source removal by soil excavation completed in 2005, although additional data would be needed to confirm this preliminary explanation.

This data set presents results from two groundwater tracing experiments conducted in the mudstone aquifer underlying the former Naval Air Warfare Center (NAWC), West Trenton, NJ. In each test, a bromide solution was introduced into a hydraulically isolated section of borehole 36BR (denoted as 36BR-A); the hydraulically isolated section of the borehole isolated specific bedding plane parting fractures intersecting the borehole 36BR. In the first test, initiated in July 2007, water samples were collected from borehole 15BR; borehole 15BR is outfitted with a submersible pump operating continuously as part of the pump-and-treat operation at the NAWC to prevent off-site migration of groundwater contaminants. In the second tracer test, initiated in September 2008, groundwater samples were collected from a hydraulically isolated section of borehole 73BR (denoted as 73BR-D2); a peristaltic pump was used extract groundwater from the selected interval. In each test, the groundwater samples were analyzed for the concentration of bromide using ion chromatography.

A simple method is proposed to visualize the magnitude of groundwater withdrawals from wells relative to user-defined water-resource metrics. The map is solely an illustration of the withdrawal magnitudes, spatially centered on wells—it is not capture zones or source areas contributing recharge to wells. Common practice is to scale the size (area) of withdrawal well symbols proportional to pumping rate. Symbols are drawn large enough to be visible, but not so large that they overlap excessively. In contrast to such graphics-based symbol sizes, the proposed method uses a depth-rate index (length per time) to visualize the well withdrawal rates by volumetrically consistent areas, called “footprints”. The area of each individual well’s footprint is the withdrawal rate divided by the depth-rate index. For example, the groundwater recharge rate could be used as a depth-rate index to show how large withdrawals are relative to that recharge. To account for the interference of nearby wells, composite footprints are computed by iterative nearest-neighbor distribution of excess withdrawals on a computational and display grid having uniform square cells. The map shows circular footprints at individual isolated wells and merged footprint areas where wells’ individual footprints overlap. Examples are presented for depth-rate indexes corresponding to recharge, to spatially variable stream baseflow (normalized by basin area), and to the average rate of water-table decline (scaled by specific yield). These depth-rate indexes are water-resource metrics, and the footprints visualize the magnitude of withdrawals relative to these metrics.

Synthesis of rock-core sampling and chlorinated volatile organic compound (CVOC) analysis at five coreholes, with hydraulic and water-quality monitoring and a detailed hydrogeologic framework, was used to characterize the fine-scale distribution of CVOCs in dipping, fractured mudstones of the Lockatong Formation of Triassic age, of the Newark Basin in West Trenton, New Jersey. From these results, a refined conceptual model for more than 55years of migration of CVOCs and depth- and strata-dependent rock-matrix contamination was developed. Industrial use of trichloroethene (TCE) at the former Naval Air Warfare Center (NAWC) from 1953 to 1995 resulted in dense non-aqueous phase liquid (DNAPL) TCE and dissolved TCE and related breakdown products, including other CVOCs, in underlying mudstones. Shallow highly weathered and fractured strata overlie unweathered, gently dipping, fractured strata that become progressively less fractured with depth. The unweathered lithology includes black highly fractured (fissile) carbon-rich strata, gray mildly fractured thinly layered (laminated) strata, and light-gray weakly fractured massive strata. CVOC concentrations in water samples pumped from the shallow weathered and highly fractured strata remain elevated near residual DNAPL TCE, but dilution by uncontaminated recharge, and other natural and engineered attenuation processes, have substantially reduced concentrations along flow paths removed from sources and residual DNAPL. CVOCs also were detected in most rock-core samples in source areas in shallow wells. In many locations, lower aqueous concentrations, compared to rock core concentrations, suggest that CVOCs are presently back-diffusing from the rock matrix. Below the weathered and highly fractured strata, and to depths of at least 50 meters (m), groundwater flow and contaminant transport is primarily in bedding-plane-oriented fractures in thin fissile high-carbon strata, and in fractured, laminated strata of the gently dipping mudstones. Despite more than 18 years of pump and treat (P&T) remediation, and natural attenuation processes, CVOC concentrations in aqueous samples pumped from these deeper strata remain elevated in isolated intervals. DNAPL was detected in one borehole during coring at a depth of 27 m. In contrast to core samples from the weathered zone, concentrations in core samples from deeper unweathered and unfractured strata are typically below detection. However, high CVOC concentrations were found in isolated samples from fissile black carbon-rich strata and fractured gray laminated strata. Aqueous-phase concentrations were correspondingly high in samples pumped from these strata via short-interval wells or packer-isolated zones in long boreholes. A refined conceptual site model considers that prior to P&T remediation groundwater flow was primarily subhorizontal in the higher-permeability near surface strata, and the bulk of contaminant mass was shallow. CVOCs diffused into these fractured and weathered mudstones. DNAPL and high concentrations of CVOCs migrated slowly down in deeper unweathered strata, primarily along isolated dipping bedding-plane fractures. After P&T began in 1995, using wells open to both shallow and deep strata, downward transport of dissolved CVOCs accelerated. Diffusion of TCE and other CVOCs from deeper fractures penetrated only a few centimeters into the unweathered rock matrix, likely due to sorption of CVOCs on rock organic carbon. Remediation in the deep, unweathered strata may benefit from the relatively limited migration of CVOCs into the rock matrix. Synthesis of rock core sampling from closely spaced boreholes with geophysical logging and hydraulic testing improves understanding of the controls on CVOC delineation and informs remediation design and monitoring.

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 2H was also injected as a tracer of the BA solution, to monitor advective transport processes. The changes in concentration and the δ13C of TCE, cis-dichloroethene (cis-DCE), and vinyl chloride (VC); the δ2H 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 13C 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.

Streams crossing underground coal mines may lose flow, whereas abandoned mine drainage (AMD) restores flow downstream. During 2005–2012, discharge from the Pine Knot Mine Tunnel, the largest AMD source in the upper Schuylkill River Basin, had near‐neutral pH and elevated concentrations of iron, manganese and sulphate. Discharge from the tunnel responded rapidly to recharge but exhibited a prolonged recession compared with nearby streams, consistent with rapid infiltration of surface water and slow release of groundwater from the mine complex. Dissolved iron was attenuated downstream by oxidation and precipitation, whereas dissolved CO2 degassed and pH increased. During high flow conditions, the AMD and downstream waters exhibited decreased pH, iron and sulphate with increased acidity that were modelled by mixing net‐alkaline AMD with recharge or run‐off having low ionic strength and low pH. Attenuation of dissolved iron within the river was least effective during high flow conditions because of decreased transport time coupled with inhibitory effects of low pH on oxidation kinetics.

Groundwater in the vicinity of several industrial facilities in Upper Gwynedd Township and vicinity, Montgomery County, in southeast Pennsylvania has been shown to be contaminated with volatile organic compounds (VOCs), the most common of which is the solvent trichloroethylene (TCE). The 2-square-mile area was placed on the National Priorities List as the North Penn Area 7 Superfund site by the U.S. Environmental Protection Agency (USEPA) in 1989. The U.S. Geological Survey (USGS) conducted geophysical logging, aquifer testing, and water-level monitoring, and measured streamflows in and near North Penn Area 7 from fall 2000 through fall 2006 in a technical assistance study for the USEPA to develop an understanding of the hydrogeologic framework in the area as part of the USEPA Remedial Investigation. In addition, the USGS developed a groundwater-flow computer model based on the hydrogeologic framework to simulate regional groundwater flow and to estimate directions of groundwater flow and pathways of groundwater contaminants. The study area is underlain by Triassic- and Jurassic-age sandstones and shales of the Lockatong Formation and Brunswick Group in the Mesozoic Newark Basin. Regionally, these rocks strike northeast and dip to the northwest. The sequence of rocks form a fractured-sedimentary-rock aquifer that acts as a set of confined to partially confined layers of differing permeabilities. Depth to competent bedrock typically is less than 20 ft below land surface. The aquifer layers are recharged locally by precipitation and discharge locally to streams. The general configuration of the potentiometric surface in the aquifer is similar to topography, except in areas affected by pumping. The headwaters of Wissahickon Creek are nearby, and the stream flows southwest, parallel to strike, to bisect North Penn Area 7. Groundwater is pumped in the vicinity of North Penn Area 7 for industrial use, public supply, and residential supply. Results of field investigations by USGS at the site and results from other studies support, and are consistent with, a conceptual model of a layered leaky aquifer where the dip of the beds has a strong control on hydraulic connections in the groundwater system. Connections within and (or) parallel to bedding tend to be greater than across bedding. Transmissivities of aquifer intervals isolated by packers ranged over three orders of magnitude [from about 2.8 to 2,290 square feet per day (ft2/d) or 0.26 to 213 square meters per day (m2/d)], did not appear to differ much by mapped geologic unit, but showed some relation to depth being relatively smaller in the shallowest and deepest intervals (0 to 50 ft and more than 250 ft below land surface, respectively) compared to the intermediate depth intervals (50 to 250 ft below land surface) tested. Transmissivities estimated from multiple-observation well aquifer tests ranged from about 700 to 2,300 ft2/d (65 to 214 m2/d). Results of chemical analyses of water from isolated intervals or monitoring wells open to short sections of the aquifer show vertical differences in concentrations; chloride and silica concentrations generally were greater in shallow intervals than in deeper intervals. Chloride concentrations greater than 100 milligrams per liter (mg/L), combined with distinctive chloride/bromide ratios, indicate a different source of chloride in the western part of North Penn Area 7 than elsewhere in the site. Groundwater flow at a regional scale under steady-state conditions was simulated by use of a numerical model (MODFLOW-2000) for North Penn Area 7 with different layers representing saprolite/highly weathered rock near the surface and unweathered competent bedrock. The sedimentary formations that underlie the study area were modeled using dipping model layers for intermediate and deep zones of unweathered, fractured rock. Horizontal cell model size was 100 meters (m) by 100 meters (328 ft by 328 ft), and model layer thickness ranged from 6 m (19.7 ft) representing shallow weathered rock and saprolite up to 200 m (656 ft) representing deeper dipping bedrock. The model did not include detailed structure to account for local-scale differences in hydraulic properties, with the result that local-scale groundwater flow may not be well simulated. Additional detailed multi-well aquifer tests would be needed to establish the extent of interconnection between intervals at the local scale to address remediation of contamination at each source area. This regional groundwater-flow model was calibrated against measured groundwater levels (1996, 2000, and 2005) and base flow estimated from selected streamflow measurements by use of nonlinear-regression parameter-estimation algorithms to determine hydraulic conductivity and anisotropy of hydraulic conductivity, streambed hydraulic conductivity, and recharge during calibration periods. Results of the simulation using the calibrated regional model indicate that the aquifer appears to be anisotropic where hydraulic conductivity is greatest parallel to the orientation of bedding of the formations underlying the area and least in the cross-bed direction. The maximum hydraulic conductivity is aligned with the average regional strike of the formations, which is “subhorizontal” in the model because the altitudes of the beds and model cells vary in the strike, as well as dip, direction. Estimated subhorizontal hydraulic conductivities (in strike direction parallel to dipping beds) range from 0.001 to 1.67 meters per day (0.0032 to 5.5 feet per day). The ratio of minimum (dip direction) to maximum (strike direction) subhorizontal hydraulic conductivity ranges from 1/3.1 to 1/8.6, and the ratio of vertical to horizontal hydraulic conductivity ranges from 1/1 to 1/478. However, limited available field data precluded rigorous calibration of vertical anisotropy in the model. Estimated recharge rates corresponding to calibration periods in 1996, 2000, and 2005 are 150, 109, and 124 millimeters per year (5.9, 4.3, and 4.9 inches per year), respectively. The calibrated groundwater-flow model was used to simulate groundwater flow under steady-state conditions during periods of relatively high withdrawals (pumpage) (1990) and relatively low withdrawals (2000 and 2005). Groundwater-flow paths originating from recharge areas near known areas of soil contamination (sources) were simulated. Pumped industrial and production wells captured more groundwater from several of these sources during 1990 than after 1990 when pumping declined or ceased and greater amounts of contaminated groundwater moved away from North Penn Area 7 Superfund site to surrounding areas. Uncertainty in simulated groundwater-flow paths from contaminant sources and contributing areas, resulting from uncertainty in estimated hydraulic properties of the model, was illustrated through Monte Carlo simulations. The effect of uncertainty in the vertical anisotropy was not included in the Monte Carlo simulations. Contributing areas indicating the general configuration of groundwater flow towards production well MG-202 (L-22) in the study area also were simulated for the different time periods; as simulated, the flow paths do not pass through any identified contaminant source in North Penn Area 7. However, contributing areas to wells, such as MG-202, located near many pumped wells are particularly complex and, in some cases, include areas that contribute flow to streams that subsequently recharge the aquifer through stream loss. In these cases, water-quality constituents, including contaminants that are present in surface water may be drawn into the aquifer to nearby pumped wells. Results of a simulated shutdown of well MG-202 under steady-state 2005 conditions showed that the area contributing recharge for nearby production well MG-76 (L-17), when MG-202 is not pumping, shifts downstream and is similar to the area contributing recharge for MG-202 when both wells are pumping. Concentrations of constituents in groundwater samples collected in fall 2005 or spring 2006 were compared to simulated groundwater-flow paths for the year 2005 to provide a qualitative assessment of model results. The observed spatial distribution of selected constituents, including TCE, CFC-11, and CFC-113 in groundwater in 2005 and the chloride/bromide mass ratios in 2006, generally were consistent with the model results of the simulated 2005 groundwater-flow paths at North Penn Area 7, indicating the presence of several separate sources of contaminants within North Penn Area 7.

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.

Changes in groundwater levels and salinity in six groundwater basins in Jordan were characterized by using linear trends fit to well-monitoring data collected from 1960 to early 2011. On the basis of data for 117 wells, groundwater levels in the six basins were declining, on average about -1 meter per year (m/yr), in 2010. The highest average rate of decline, -1.9 m/yr, occurred in the Jordan Side Valleys basin, and on average no decline occurred in the Hammad basin. The highest rate of decline for an individual well was -9 m/yr. Aquifer saturated thickness, a measure of water storage, was forecast for year 2030 by using linear extrapolation of the groundwater-level trend in 2010. From 30 to 40 percent of the saturated thickness, on average, was forecast to be depleted by 2030. Five percent of the wells evaluated were forecast to have zero saturated thickness by 2030. Electrical conductivity was used as a surrogate for salinity (total dissolved solids). Salinity trends in groundwater were much more variable and less linear than groundwater-level trends. The long-term linear salinity trend at most of the 205 wells evaluated was not increasing, although salinity trends are increasing in some areas. The salinity in about 58 percent of the wells in the Amman-Zarqa basin was substantially increasing, and the salinity in Hammad basin showed a long-term increasing trend. Salinity increases were not always observed in areas with groundwater-level declines. The highest rates of salinity increase were observed in regional discharge areas near groundwater pumping centers.

Significant microbial reductive dechlorination of [1,2 C-14] cis-dichloroethene (DCE) was observed in anoxic microcosms prepared with unamended, fractured rock aquifer materials, which were colonized in situ at multiple depths in two boreholes at the Naval Air Warfare Center (NAWC) in West Trenton, New Jersey. The lack of significant reductive dechlorination in corresponding water-only treatments indicated that chlororespiration activity in unamended, fractured rock treat-ments was primarily associated with colonized core material. In these unamended fractured rock microcosms, activity was highest in the shallow zones and generally decreased with increasing depth. Electron-donor amendment (biostimulation) enhanced chlororespiration in some but not all treatments. In contrast, combining electron-donor amendment with KB1 amendment (bioaug-mentation) enhanced chlororespiration in all treatments and substantially reduced the variability in chlororespiration activity both within and between treatments. These results indicate (1) that a potential for chlororespiration-based bioremediation exists at NAWC Trenton but is limited under nonengineered conditions, (2) that the limitation on chlororespiration activity is not entirely due to electron-donor availability, and (3) that a bioaugmentation approach can substantially enhance in situ bioremediation if the requisite amendments can be adequately distributed throughout the fractured rock matrix. (C) 2012 Wiley Periodicals, Inc.

There is a need to better characterize contaminated aquifers in fractured rock and to determine the fate of remediation injections away from boreholes. A synthetic cross borehole electrical resistivity study was conducted within a fractured rock setting representing the Naval Air Warfare Center (NAWC), New Jersey, which is contaminated with trichloroethylene (TCE), cis-1,2-dichloroethene (DCE), and vinyl chloride (VC). Four boreholes and two fracture seams, assumed to be the dominant electrical and hydraulically conductive pathways, were explicitly modeled within an unstructured tetrahedral mesh. Parallel processing was used to handle the large number of mesh elements associated with modeling such small-scale features. To simulate field conditions, 5% random noise was added to all synthetic datasets. Five 3D inversions were performed on the synthetic data, each with different regularization constraints. We began with an uninformed smoothness-constrained inversion, to which a priori information was incrementally added. We found major improvements when (1) smoothness regularization constraints were relaxed (or disconnected) along boreholes and fractures, (2) a homogeneous conductivity was assumed along boreholes, and (3) borehole and packer conductivity constraints were applied. Incorporating this information has the potential to improve imaging of fracture location, orientation, and changes in electrical conductivity distribution due to remediation injections. INTRODUCTION Delineation of fracture networks in a contaminated fractured rock aquifer is essential to successful remediation. Borehole logging and point sampling can provide information near the well; however, information between these discrete locations is limited and spatial interpolation of such data is highly uncertain. Cross borehole fracture orientation is particularly important in remediation efforts that include amendment injections where deliveries need to be monitored to ensure treatment of target areas. Thus, there is a need for techniques to image fractures and the delivery of amendments between boreholes. Using geophysical methods, in particular electrical resistivity, with standard practices can potentially provide spatially continuous information at a reduced