Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Contaminant Hydrology”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3Linked to original sources
USGS2006

Recharge processes drive sulfate reduction in an alluvial aquifer contaminated with landfill leachate

Natural attenuation of contaminants in groundwater depends on an adequate supply of electron acceptors to stimulate biodegradation. In an alluvial aquifer contaminated with leachate from an unlined municipal landfill, the mechanism of recharge infiltration was investigated as a source of electron acceptors. Water samples were collected monthly at closely spaced intervals in the top 2 m of the saturated zone from a leachate-contaminated well and an uncontaminated well, and analyzed for δ 18 O, δ 2 H, non-volatile dissolved organic carbon (NVDOC), SO 4 2− , NO 3 − and Cl − . Monthly recharge amounts were quantified using the offset of the δ 18 O or δ 2 H from the local meteoric water line as a parameter to distinguish water types, as evaporation and methanogenesis caused isotopic enrichment in waters from different sources. Presence of dissolved SO 4 2− in the top 1 to 2 m of the saturated zone was associated with recharge; SO 4 2− averaged 2.2 mM, with maximum concentrations of 15 mM. Nitrate was observed near the water table at the contaminated site at concentrations up to 4.6 mM. Temporal monitoring of δ 2 H and SO 4 2− showed that vertical transport of recharge carried SO 4 2− to depths up to 1.75 m below the water table, supplying an additional electron acceptor to the predominantly methanogenic leachate plume. Measurements of δ 34 S in SO 4 2− indicated both SO 4 2− reduction and sulfide oxidation were occurring in the aquifer. Depth-integrated net SO 4 2− reduction rates, calculated using the natural Cl − gradient as a conservative tracer, ranged from 7.5 × 10 − 3 to 0.61 mM·d − 1 (over various depth intervals from 0.45 to 1.75 m). Sulfate reduction occurred at both the contaminated and uncontaminated sites; however, median SO 4 2− reduction rates were higher at the contaminated site. Although estimated SO 4 2− reduction rates are relatively high, significant decreases in NVDOC were not observed at the contaminated site. Organic compounds more labile than the leachate NVDOC may be present in the root zone, and SO 4 2− reduction may be coupled to methane oxidation. The results show that sulfur (and possibly nitrogen) redox processes within the top 2 m of the aquifer are directly related to recharge timing and seasonal water level changes in the aquifer. The results suggest that SO 4 2− reduction associated with the infiltration of recharge may be a significant factor affecting natural attenuation of contaminants in alluvial aquifers.

Journal of Contaminant Hydrology↗
USGS2006

Investigating the role of gas bubble formation and entrapment in contaminated aquifers: Reactive transport modelling

In many natural and contaminated aquifers, geochemical processes result in the production or consumption of dissolved gases. In cases where methanogenesis or denitrification occurs, the production of gases may result in the formation and growth of gas bubbles below the water table. Near the water table, entrapment of atmospheric gases during water table rise may provide a significant source of O 2 to waters otherwise depleted in O 2 . Furthermore, the presence of bubbles will affect the hydraulic conductivity of an aquifer, resulting in changes to the groundwater flow regime. The interactions between physical transport, biogeochemical processes, and gas bubble formation, entrapment and release is complex and requires suitable analysis tools. The objective of the present work is the development of a numerical model capable of quantitatively assessing these processes. The multicomponent reactive transport code MIN3P has been enhanced to simulate bubble growth and contraction due to in-situ gas production or consumption, bubble entrapment due to water table rise and subsequent re-equilibration of the bubble with ambient groundwater, and permeability changes due to trapped gas phase saturation. The resulting formulation allows for the investigation of complex geochemical systems where microbially mediated redox reactions both produce and consume gases as well as affect solution chemistry, alkalinity, and pH. The enhanced model has been used to simulate processes in a petroleum hydrocarbon contaminated aquifer where methanogenesis is an important redox process. The simulations are constrained by data from a crude oil spill site near Bemidji, MN. Our results suggest that permeability reduction in the methanogenic zone due to in-situ formation of gas bubbles, and dissolution of entrapped atmospheric bubbles near the water table, both work to attenuate the dissolved gas plume emanating from the source zone. Furthermore, the simulations demonstrate that under the given conditions more than 50% of all produced CH 4 partitions to the gas phase or is aerobically oxidised near the water table, suggesting that these processes should be accounted for when assessing the rate and extent of methanogenic degradation of hydrocarbons.

Journal of Contaminant Hydrology↗
USGS2006

Analysis of a mesoscale infiltration and water seepage test in unsaturated fractured rock: Spatial variabilities and discrete fracture patterns

A mesoscale (21??m in flow distance) infiltration and seepage test was recently conducted in a deep, unsaturated fractured rock system at the crossover point of two underground tunnels. Water was released from a 3??m ?? 4??m infiltration plot on the floor of an alcove in the upper tunnel, and seepage was collected from the ceiling of a niche in the lower tunnel. Significant temporal and (particularly) spatial variabilities were observed in both measured infiltration and seepage rates. To analyze the test results, a three-dimensional unsaturated flow model was used. A column-based scheme was developed to capture heterogeneous hydraulic properties reflected by these spatial variabilities observed. Fracture permeability and van Genuchten ?? parameter [van Genuchten, M.T., 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44, 892-898] were calibrated for each rock column in the upper and lower hydrogeologic units in the test bed. The calibrated fracture properties for the infiltration and seepage zone enabled a good match between simulated and measured (spatially varying) seepage rates. The numerical model was also able to capture the general trend of the highly transient seepage processes through a discrete fracture network. The calibrated properties and measured infiltration/seepage rates were further compared with mapped discrete fracture patterns at the top and bottom boundaries. The measured infiltration rates and calibrated fracture permeability of the upper unit were found to be partially controlled by the fracture patterns on the infiltration plot (as indicated by their positive correlations with fracture density). However, no correlation could be established between measured seepage rates and density of fractures mapped on the niche ceiling. This lack of correlation indicates the complexity of (preferential) unsaturated flow within the discrete fracture network. This also indicates that continuum-based modeling of unsaturated flow in fractured rock at mesoscale or a larger scale is not necessarily conditional explicitly on discrete fracture patterns. ?? 2006 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗
USGS2006

Tracer test with As(V) under variable redox conditions controlling arsenic transport in the presence of elevated ferrous iron concentrations

To study transport and reactions of arsenic under field conditions, a small-scale tracer test was performed in an anoxic, iron-reducing zone of a sandy aquifer at the USGS research site on Cape Cod, Massachusetts, USA. For four weeks, a stream of groundwater with added As(V) (6.7 μM) and bromide (1.6 mM), was injected in order to observe the reduction of As(V) to As(III). Breakthrough of bromide (Br − ), As(V), and As(III) as well as additional parameters characterizing the geochemical conditions was observed at various locations downstream of the injection well over a period of 104 days. After a short lag period, nitrate and dissolved oxygen from the injectate oxidized ferrous iron and As(V) became bound to the freshly formed hydrous iron oxides. Approximately one week after terminating the injection, anoxic conditions had been reestablished and increases in As(III) concentrations were observed within 1 m of the injection. During the observation period, As(III) and As(V) were transported to a distance of 4.5 m downgradient indicating significant retardation by sorption processes for both species. Sediment assays as well as elevated concentrations of hydrogen reflected the presence of As(V) reducing microorganisms. Thus, microbial As(V) reduction was thought to be one major process driving the release of As(III) during the tracer test in the Cape Cod aquifer.

Journal of Contaminant Hydrology↗
USGS2007

Insights into the use of time-lapse GPR data as observations for inverse multiphase flow simulations of DNAPL migration

Perchloroethylene (PCE) saturations determined from GPR surveys were used as observations for inversion of multiphase flow simulations of a PCE injection experiment (Borden 9??m cell), allowing for the estimation of optimal bulk intrinsic permeability values. The resulting fit statistics and analysis of residuals (observed minus simulated PCE saturations) were used to improve the conceptual model. These improvements included adjustment of the elevation of a permeability contrast, use of the van Genuchten versus Brooks-Corey capillary pressure-saturation curve, and a weighting scheme to account for greater measurement error with larger saturation values. A limitation in determining PCE saturations through one-dimensional GPR modeling is non-uniqueness when multiple GPR parameters are unknown (i.e., permittivity, depth, and gain function). Site knowledge, fixing the gain function, and multiphase flow simulations assisted in evaluating non-unique conceptual models of PCE saturation, where depth and layering were reinterpreted to provide alternate conceptual models. Remaining bias in the residuals is attributed to the violation of assumptions in the one-dimensional GPR interpretation (which assumes flat, infinite, horizontal layering) resulting from multidimensional influences that were not included in the conceptual model. While the limitations and errors in using GPR data as observations for inverse multiphase flow simulations are frustrating and difficult to quantify, simulation results indicate that the error and bias in the PCE saturation values are small enough to still provide reasonable optimal permeability values. The effort to improve model fit and reduce residual bias decreases simulation error even for an inversion based on biased observations and provides insight into alternate GPR data interpretations. Thus, this effort is warranted and provides information on bias in the observation data when this bias is otherwise difficult to assess. ?? 2006 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗
USGS2007

The chemical response of particle-associated contaminants in aquatic sediments to urbanization in New England, U.S.A.

Relations between urbanization and particle-associated contaminants in New England were evaluated using a combination of samples from sediment cores, streambed sediments, and suspended stream sediments. Concentrations of PAHs, PCBs, DDT, and seven trace metals (Cd, Cr, Cu, Hg, Ni, Pb, Zn) were correlated strongly with urbanization, with the strongest relations to percentage commercial, industrial, and transportation (CIT) land use. Average PAH and metal concentrations in the most urbanized watersheds were approximately 30 and 6 times the reference concentrations, respectively, in remote, undeveloped watersheds. One-quarter to one-half of sampling sites had concentrations of PAHs, Cu, Pb, or Zn above the probable effects concentration, a set of sediment quality guidelines for adverse effects to aquatic biota, and sediments were predicted to be toxic, on average, when CIT land use exceeded about 10%. Trends in metals in cores from urban watersheds were dominantly downward, whereas trends in PAHs in a suburban watershed were upward. A regional atmospheric-fallout gradient was indicated by as much as order-of-magnitude-greater concentrations and accumulation rates of contaminants in cores from an undeveloped reference lake in Boston compared to those from remote reference watersheds. Contaminant accumulation rates in the lakes with urbanization in their watersheds, however, were 1-3 orders of magnitude greater than those of reference lakes, which indicate the dominance of local sources and fluvial transport of contaminants to urban lakes. These analyses demonstrate the magnitude of urban contamination of aquatic systems and air sheds, and suggest that, despite reductions in contaminant emissions in urban settings, streams and lakes will decline in quality as urbanization of their watersheds takes place. ?? 2006 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗
USGS2007

Anthropogenic contaminants as tracers in an urbanizing karst aquifer

Karst aquifers are uniquely vulnerable to contamination. In the Barton Springs segment of the karstic Edwards aquifer (Texas, U.S.A.), urban contaminants such as pesticides and volatile organic compounds frequently are detected in spring base flow. To determine whether contaminant concentrations change in response to storms, and if they therefore might act as tracers of focused recharge, samples were collected from Barton Springs at closely spaced intervals following three storms. Two herbicides (atrazine and simazine), two insecticides (carbaryl and diazinon), and a solvent (tetrachloroethene) described breakthrough curves over a 1-week period following one or more storms. The breakthrough curves were decomposed into two to five log-normal subcurves, which were interpreted as representing pulses of contaminants moving through the aquifer. Each subcurve could be used in the same way as an artificial tracer to determine travel time to and recovery at the spring. The contaminants have several advantages over artificial tracers: they represent the actual compounds of interest, they are injected essentially simultaneously at several points, and they are injected under those conditions when transport is of the most interest, i.e., following storms. The response of storm discharge, specific conductance, and contaminant loading at the spring depended on initial aquifer flow conditions, which varied from very low (spring discharge of 0.48??m3/s) to high (spring discharge of 2.7??m3/s): concentrations and recovery were the highest when initial aquifer flow conditions were low. This behavior provides information about aquifer structure and the influence of aquifer flow condition on transport properties. ?? 2006 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗
USGS2007

Linking ground-water age and chemistry data along flow paths: Implications for trends and transformations of nitrate and pesticides

Tracer-based ground-water ages, along with the concentrations of pesticides, nitrogen species, and other redox-active constituents, were used to evaluate the trends and transformations of agricultural chemicals along flow paths in diverse hydrogeologic settings. A range of conditions affecting the transformation of nitrate and pesticides (e.g., thickness of unsaturated zone, redox conditions) was examined at study sites in Georgia, North Carolina, Wisconsin, and California. Deethylatrazine (DEA), a transformation product of atrazine, was typically present at concentrations higher than those of atrazine at study sites with thick unsaturated zones but not at sites with thin unsaturated zones. Furthermore, the fraction of atrazine plus DEA that was present as DEA did not increase as a function of ground-water age. These findings suggest that atrazine degradation occurs primarily in the unsaturated zone with little or no degradation in the saturated zone. Similar observations were also made for metolachlor and alachlor. The fraction of the initial nitrate concentration found as excess N2 (N2 derived from denitrification) increased with ground-water age only at the North Carolina site, where oxic conditions were generally limited to the top 5??m of saturated thickness. Historical trends in fluxes to ground water were evaluated by relating the times of recharge of ground-water samples, estimated using chlorofluorocarbon concentrations, with concentrations of the parent compound at the time of recharge, estimated by summing the molar concentrations of the parent compound and its transformation products in the age-dated sample. Using this approach, nitrate concentrations were estimated to have increased markedly from 1960 to the present at all study sites. Trends in concentrations of atrazine, metolachlor, alachlor, and their degradates were related to the timing of introduction and use of these compounds. Degradates, and to a lesser extent parent compounds, were detected in ground water dating back to the time these compounds were introduced.

Journal of Contaminant Hydrology↗
USGS2008

Nutrient dynamics as indicators of karst processes: Comparison of the Chalk aquifer (Normandy, France) and the Edwards aquifer (Texas, U.S.A.)

Karst aquifers display a range of geologic and geomorphic characteristics in a wide range of climatic and land-use settings; identification of transport dynamics representative of karst aquifers in general could help advance our understanding of these complex systems. To this end, nutrient, turbidity, and major ion dynamics in response to storms were compared at multiple sites in two karst aquifers with contrasting characteristics and settings: the Chalk aquifer (Eure Department, Normandy, France) and the Barton Springs segment of the Edwards Aquifer (Texas, U.S.A.). The Chalk aquifer is typified by high matrix porosity, thick surficial deposits (up to 30??m thick), and agricultural land use; the Barton Springs segment is typified by low matrix porosity, outcropping limestone, and urban land use. Following one to three storms, from 5 to 16 samples from springs and wells were analyzed for major ions, and specific conductance and turbidity were monitored continuously. Comparison of the chemographs indicated some generalized responses, including an increase in turbidity and potassium concentrations and a decrease in major ion and nitrate concentrations with infiltrating storm runoff. Factor analysis of major ions and turbidity revealed strikingly similar behavior of the chemical variables for the two aquifers: The first two factors, explaining more than 75% of the variability, illustrate that dynamics of most major ions (including nitrate) are opposed to those of turbidity and of potassium. The results demonstrate that potassium and nitrate are effective tracers of infiltrating storm runoff and resident ground water, respectively, and the similar results for these two highly contrasting aquifers suggest that the dynamics identified might be applicable to karst systems in general. ?? 2008 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗
USGS2008

Distribution and variability of redox zones controlling spatial variability of arsenic in the Mississippi River Valley alluvial aquifer, southeastern Arkansas

Twenty one of 118 irrigation water wells in the shallow (25-30??m thick) Mississippi River Valley alluvial aquifer in the Bayou Bartholomew watershed, southeastern Arkansas had arsenic (As) concentrations (< 0.5 to 77????g/L) exceeding 10????g/L. Sediment and groundwater samples were collected and analyzed from the sites of the highest, median, and lowest concentrations of As in groundwater in the alluvial aquifers located at Jefferson County, Arkansas. A traditional five-step sequential extraction was performed to differentiate the exchangeable, carbonate, amorphous Fe and Mn oxide, organic, and hot HNO3-leachable fraction of As and other compounds in sediments. The Chao reagent (0.25??M hydroxylamine hydrochloride in 0.25??M HCl) removes amorphous Fe and Mn oxides and oxyhydroxides (present as coatings on grains and amorphous minerals) by reductive dissolution and is a measure of reducible Fe and Mn in sediments. The hot HNO3 extraction removes mostly crystalline metal oxides and all other labile forms of As. Significant total As (20%) is complexed with amorphous Fe and Mn oxides in sediments. Arsenic abundance is not significant in carbonates or organic matter. Significant (40-70????g/kg) exchangeable As is only present at shallow depth (0-1??m below ground surface). Arsenic is positively correlated to Fe extracted by Chao reagent (r = 0.83) and hot HNO3 (r = 0.85). Arsenic extracted by Chao reagent decreases significantly with depth as compared to As extracted by hot HNO3. Fe (II)/Fe (the ratio of Fe concentration in the extracts of Chao reagent and hot HNO3) is positively correlated (r = 0.76) to As extracted from Chao reagent. Although Fe (II)/Fe increases with depth, the relative abundance of reducible Fe decreases noticeably with depth. The amount of reducible Fe, as well as As complexed to amorphous Fe and Mn oxides and oxyhydroxides decreases with depth. Possible explanations for the decrease in reducible Fe and its complexed As with depth include historic flushing of As and Fe from hydrous ferric oxides (HFO) by microbially-mediated reductive dissolution and aging of HFO to crystalline phases. Hydrogeochemical data suggests that the groundwater in the area falls in the mildly reducing (suboxic) to relatively highly reducing (anoxic) zone, and points to reductive dissolution of HFO as the dominant As release mechanism. Spatial variability of gypsum solubility and simultaneous SO42- reduction with co-precipitation of As and sulfide is an important limiting process controlling the concentration of As in groundwater in the area. ?? 2008 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗
USGS2009

Hydraulic anisotropy characterization of pneumatic-fractured sediments using azimuthal self potential gradient

The pneumatic fracturing technique is used to enhance the permeability and porosity of tight unconsolidated soils (e.g. clays), thereby improving the effectiveness of remediation treatments. Azimuthal self potential gradient (ASPG) surveys were performed on a compacted, unconsolidated clay block in order to evaluate their potential to delineate contaminant migration pathways in a mechanically-induced fracture network. Azimuthal resistivity (ARS) measurements were also made for comparative purposes. Following similar procedures to those used in the field, compressed kaolinite sediments were pneumatically fractured and the resulting fracture geometry characterized from strike analysis of visible fractures combined with strike data from optical borehole televiewer (BHTV) imaging. We subsequently injected a simulated treatment (electrolyte/dye) into the fractures. Both ASPG and ARS data exhibit anisotropic geoelectric signatures resulting from the fracturing. Self potentials observed during injection of electrolyte are consistent with electrokinetic theory and previous laboratory results on a fracture block model. Visual (polar plot) analysis and linear regression of cross plots show ASPG lobes are correlated with azimuths of high fracture strike density, evidence that the ASPG anisotropy is a proxy measure of hydraulic anisotropy created by the pneumatic fracturing. However, ARS data are uncorrelated with fracture strike maxima and resistivity anisotropy is probably dominated by enhanced surface conduction along azimuths of weak 'starter paths' formed from pulverization of the clay and increases in interfacial surface area. We find the magnitude of electrokinetic SP scales with the applied N2 gas pressure gradient (??PN2) for any particular hydraulically-active fracture set and that the positive lobe of the ASPG anomaly indicates the flow direction within the fracture network. These findings demonstrate the use of ASPG in characterizing the effectiveness of (1) pneumatic fracturing and (2) defining likely flow directions of remedial treatments in unconsolidated sediments and rock. ?? 2008 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗
USGS2009

Microbial characterization of nitrification in a shallow, nitrogen-contaminated aquifer, Cape Cod, Massachusetts and detection of a novel cluster associated with nitrifying Betaproteobacteria

Groundwater nitrification is a poorly characterized process affecting the speciation and transport of nitrogen. Cores from two sites in a plume of contamination were examined using culture-based and molecular techniques targeting nitrification processes. The first site, located beneath a sewage effluent infiltration bed, received treated effluent containing O 2 (> 300 µM) and NH 4 + (51–800 µM). The second site was 2.5 km down-gradient near the leading edge of the ammonium zone within the contaminant plume and featured vertical gradients of O 2 , NH 4 + , and NO 3 − (0–300, 0–500, and 100–200 µM with depth, respectively). Ammonia- and nitrite-oxidizers enumerated by the culture-based MPN method were low in abundance at both sites (1.8 to 350 g − 1 and 33 to 35,000 g − 1 , respectively). Potential nitrifying activity measured in core material in the laboratory was also very low, requiring several weeks for products to accumulate. Molecular analysis of aquifer DNA (nested PCR followed by cloning and 16S rDNA sequencing) detected primarily sequences associated with the Nitrosospira genus throughout the cores at the down-gradient site and a smaller proportion from the Nitrosomonas genus in the deeper anoxic, NH 4 + zone at the down-gradient site. Only a single Nitrosospira sequence was detected beneath the infiltration bed. Furthermore, the majority of Nitrosospira -associated sequences represent an unrecognized cluster. We conclude that an uncharacterized group associated with Nitrosospira dominate at the geochemically stable, down-gradient site, but found little evidence for Betaproteobacteria nitrifiers beneath the infiltration beds where geochemical conditions were more variable.

Journal of Contaminant Hydrology↗
USGS2009

Biogeochemistry at a wetland sediment-alluvial aquifer interface in a landfill leachate plume

The biogeochemistry at the interface between sediments in a seasonally ponded wetland (slough) and an alluvial aquifer contaminated with landfill leachate was investigated to evaluate factors that can effect natural attenuation of landfill leachate contaminants in areas of groundwater/surface-water interaction. The biogeochemistry at the wetland-alluvial aquifer interface differed greatly between dry and wet conditions. During dry conditions (low water table), vertically upward discharge was focused at the center of the slough from the fringe of a landfill-derived ammonium plume in the underlying aquifer, resulting in transport of relatively low concentrations of ammonium to the slough sediments with dilution and dispersion as the primary attenuation mechanism. In contrast, during wet conditions (high water table), leachate-contaminated groundwater discharged upward near the upgradient slough bank, where ammonium concentrations in the aquifer where high. Relatively high concentrations of ammonium and other leachate constituents also were transported laterally through the slough porewater to the downgradient bank in wet conditions. Concentrations of the leachate-associated constituents chloride, ammonium, non-volatile dissolved organic carbon, alkalinity, and ferrous iron more than doubled in the slough porewater on the upgradient bank during wet conditions. Chloride, non-volatile dissolved organic carbon (DOC), and bicarbonate acted conservatively during lateral transport in the aquifer and slough porewater, whereas ammonium and potassium were strongly attenuated. Nitrogen isotope variations in ammonium and the distribution of ammonium compared to other cations indicated that sorption was the primary attenuation mechanism for ammonium during lateral transport in the aquifer and the slough porewater. Ammonium attenuation was less efficient, however, in the slough porewater than in the aquifer and possibly occurred by a different sorption mechanism. A stoichiometrically balanced increase in magnesium concentration with decreasing ammonium and potassium concentrations indicated that cation exchange was the sorption mechanism in the slough porewater. Only a partial mass balance could be determined for cations exchanged for ammonium and potassium in the aquifer, indicating that some irreversible sorption may be occurring. Although wetlands commonly are expected to decrease fluxes of contaminants in riparian environments, enhanced attenuation of the leachate contaminants in the slough sediment porewater compared to the aquifer was not observed in this study. The lack of enhanced attenuation can be attributed to the fact that the anoxic plume, comprised largely of recalcitrant DOC and reduced inorganic constituents, interacted with anoxic slough sediments and porewaters, rather than encountering a change in redox conditions that could cause transformation reactions. Nevertheless, the attenuation processes in the narrow zone of groundwater/surface-water interaction were effective in reducing ammonium concentrations by a factor of about 3 during lateral transport across the slough and by a factor of 2 to 10 before release to the surface water. Slough porewater geochemistry also indicated that the slough could be a source of sulfate in dry conditions, potentially providing a terminal electron acceptor for natural attenuation of organic compounds in the leachate plume.

Journal of Contaminant Hydrology↗
USGS2009

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.

Journal of Contaminant Hydrology↗
USGS2010

Vadose zone attenuation of organic compounds at a crude oil spill site - Interactions between biogeochemical reactions and multicomponent gas transport

Contaminant attenuation processes in the vadose zone of a crude oil spill site near Bemidji, MN have been simulated with a reactive transport model that includes multicomponent gas transport, solute transport, and the most relevant biogeochemical reactions. Dissolution and volatilization of oil components, their aerobic and anaerobic degradation coupled with sequential electron acceptor consumption, ingress of atmospheric O 2 , and the release of CH 4 and CO 2 from the smear zone generated by the floating oil were considered. The focus of the simulations was to assess the dynamics between biodegradation and gas transport processes in the vadose zone, to evaluate the rates and contributions of different electron accepting processes towards vadose zone natural attenuation, and to provide an estimate of the historical mass loss. Concentration distributions of reactive (O 2 , CH 4 , and CO 2 ) and non-reactive (Ar and N 2 ) gases served as key constraints for the model calibration. Simulation results confirm that as of 2007, the main degradation pathway can be attributed to methanogenic degradation of organic compounds in the smear zone and the vadose zone resulting in a contaminant plume dominated by high CH 4 concentrations. In accordance with field observations, zones of volatilization and CH 4 generation are correlated to slightly elevated total gas pressures and low partial pressures of N 2 and Ar, while zones of aerobic CH 4 oxidation are characterized by slightly reduced gas pressures and elevated concentrations of N 2 and Ar. Diffusion is the most significant transport mechanism for gases in the vadose zone; however, the simulations also indicate that, despite very small pressure gradients, advection contributes up to 15% towards the net flux of CH 4 , and to a more limited extent to O 2 ingress. Model calibration strongly suggests that transfer of biogenically generated gases from the smear zone provides a major control on vadose zone gas distributions and vadose zone carbon balance. Overall, the model was successful in capturing the complex interactions between biogeochemical reactions and multicomponent gas transport processes. However, despite employing a process-based modeling approach, honoring observed parameter ranges, and generally obtaining good agreement between field observations and model simulations, accurate quantification of natural attenuation rates remains difficult. The modeling results are affected by uncertainties regarding gas phase saturations, tortuosities, and the magnitude of CH 4 and CO 2 flux from the smear zone. These findings highlight the need to better delineate gas fluxes at the model boundaries, which will help constrain contaminant degradation rates, and ultimately source zone longevity.

Minnesota↗
USGS2010

In situ measurements of volatile aromatic hydrocarbon biodegradation rates in groundwater

Benzene and alkylbenzene biodegradation rates and patterns were measured using an in situ microcosm in a crude-oil contaminated aquifer near Bemidji, Minnesota. Benzene-D6, toluene, ethylbenzene, o-, m- and p-xylenes and four pairs of C3- and C4-benzenes were added to an in situ microcosm and studied over a 3-year period. The microcosm allowed for a mass-balance approach and quantification of hydrocarbon biodegradation rates within a well-defined iron-reducing zone of the anoxic plume. Among the BTEX compounds, the apparent order of persistence is ethylbenzene > benzene > m,p-xylenes > o-xylene ≥ toluene. Threshold concentrations were observed for several compounds in the in situ microcosm, below which degradation was not observed, even after hundreds of days. In addition, long lag times were observed before the onset of degradation of benzene or ethylbenzene. The isomer-specific degradation patterns were compared to observations from a multi-year study conducted using data collected from monitoring wells along a flowpath in the contaminant plume. The data were fit with both first-order and Michaelis-Menten models. First-order kinetics provided a good fit for hydrocarbons with starting concentrations below 1 mg/L and Michaelis-Menten kinetics were a better fit when starting concentrations were above 1 mg/L, as was the case for benzene. The biodegradation rate data from this study were also compared to rates from other investigations reported in the literature.

Journal of Contaminant Hydrology↗
USGS2011

Mineralogy, morphology, and textural relationships in coatings on quartz grains in sediments in a quartz-sand aquifer

Mineralogical studies of coatings on quartz grains and bulk sediments from an aquifer on Western Cape Cod, Massachusetts, USA were carried out using a variety of transmission electron microscopy (TEM) techniques. Previous studies demonstrated that coatings on quartz grains control the adsorption properties of these sediments. Samples for TEM characterization were made by a gentle mechanical grinding method and focused ion beam (FIB) milling. The former method can make abundant electron-transparent coating assemblages for comprehensive and quantitative X-ray analysis and the latter technique protects the coating texture from being destroyed. Characterization of the samples from both a pristine area and an area heavily impacted by wastewater discharge shows similar coating textures and chemical compositions. Major constituents of the coating include Al-substituted goethite and illite/chlorite clays. Goethite is aggregated into well-crystallized domains through oriented attachment resulting in increased porosity. Illite/chlorite clays with various chemical compositions were observed to be mixed with goethite aggregates and aligned sub-parallel to the associated quartz surface. The uniform spatial distribution of wastewater-derived phosphorus throughout the coating from the wastewater-contaminated site suggests that all of the coating constituents, including those adjacent to the quartz surface, are accessible to groundwater solutes. Both TEM characterization and chemical extraction results indicate there is a significantly greater amount of amorphous iron oxide in samples from wastewater discharge area compared to those from the pristine region, which might reflect the impact of redox cycling of iron under the wastewater-discharge area. Coating compositions are consistent with the moderate metal and oxy-metalloid adsorption capacities, low but significant cation exchange capacities, and control of iron(III) solubility by goethite observed in reactive transport experimental and modeling studies conducted at the site.

Massachusetts↗
USGS2011

Methane oxidation in a crude oil contaminated aquifer: Delineation of aerobic reactions at the plume fringes

High resolution direct-push profiling over short vertical distances was used to investigate CH 4 attenuation in a petroleum contaminated aquifer near Bemidji, Minnesota. The contaminant plume was delineated using dissolved gases, redox sensitive components, major ions, carbon isotope ratios in CH 4 and CO 2 , and the presence of methanotrophic bacteria. Sharp redox gradients were observed near the water table. Shifts in δ 13 C CH4 from an average of − 57.6‰ (± 1.7‰) in the methanogenic zone to − 39.6‰ (± 8.7‰) at 105 m downgradient, strongly suggest CH 4 attenuation through microbially mediated degradation. In the downgradient zone the aerobic/anaerobic transition is up to 0.5 m below the water table suggesting that transport of O 2 across the water table is leading to aerobic degradation of CH 4 at this interface. Dissolved N 2 concentrations that exceeded those expected for water in equilibrium with the atmosphere indicated bubble entrapment followed by preferential stripping of O 2 through aerobic degradation of CH 4 or other hydrocarbons. Multivariate and cluster analysis were used to distinguish between areas of significant bubble entrapment and areas where other processes such as the infiltration of O 2 rich recharge water were important O 2 transport mechanisms.

Journal of Contaminant Hydrology↗