Quantifying effects of soil heterogeneity on groundwater pollution in four sites in the USA

Four sites located in the north-eastern region of the United States of America have been chosen to investigate the impacts of soil heterogeneity in the transport of solutes (bromide and chloride) through the vadose zone (the zone in the soil that lies below the root zone and above the permanent saturated groundwater). A recently proposed mathematical model based on the cumulative beta distribution has been deployed to compare and contrast the regions' heterogeneity from multiple sample percolation experiments. Significant differences in patterns of solute leaching were observed even over a small spatial scale, indicating that traditional sampling methods for solute transport, for example the gravity pan or suction lysimeters, or more recent inventions such as the multiple sample percolation systems may not be effective in estimating solute fluxes in soils when a significant degree of soil heterogeneity is present. Consequently, ignoring soil heterogeneity in solute transport studies will likely result in under- or overprediction of leached fluxes and potentially lead to serious pollution of soils and/or groundwater. The cumulative beta distribution technique is found to be a versatile and simple technique of gaining valuable information regarding soil heterogeneity effects on solute transport. It is also an excellent tool for guiding future decisions of experimental designs particularly in regard to the number of samples within one site and the number of sampling locations between sites required to obtain a representative estimate of field solute or drainage flux.

Theory and application of inverse problems in groundwater: numerical, laboratory and field studies.

Aquifers are a vital water resource whose quality characteristics must be safeguarded or, if damaged, restored. The extent and complexity of aquifer contamination is related to characteristics of the porous medium, the influence of boundary conditions, and the biological, chemical and physical processes. After the nineties, the efforts of the scientists have been increased exponentially in order to find an efficient way for estimating the hydraulic parameters of the aquifers, and thus, recover the contaminant source position and its release history. To simplify and understand the influence of these various factors on aquifer phenomena, it is common for researchers to use numerical and controlled experiments. This work presents some of these methods, applying and comparing them on data collected during laboratory, field and numerical tests. The work is structured in four parts which present the results and the conclusions of the specific objectives.

Concepts and methods of quantitative assessment of risk to groundwater resources

This work presents a two-dimensional approach of risk assessment method based on the quantification of the probability of the occurrence of contaminant source terms, as well as the assessment of the resultant impacts. The risk is calculated using Monte Carlo simulation methods whereby synthetic contaminant source terms were generated to the same distribution as historically occurring pollution events or a priori potential probability distribution. The spatial and temporal distributions of the generated contaminant concentrations at pre-defined monitoring points within the aquifer were then simulated from repeated realisations using integrated mathematical models. The number of times when user defined ranges of concentration magnitudes were exceeded is quantified as risk. The utilities of the method were demonstrated using hypothetical scenarios, and the risk of pollution from a number of sources all occurring by chance together was evaluated. The results are presented in the form of charts and spatial maps. The generated risk maps show the risk of pollution at each observation borehole, as well as the trends within the study area. This capability to generate synthetic pollution events from numerous potential sources of pollution based on historical frequency of their occurrence proved to be a great asset to the method, and a large benefit over the contemporary methods.

Trichloroethylene fate and transport studies and biodegradation kinetics in the saturated zone

Permeable reactive barriers (PRB) are constructed from soil solid amendments to support the growth of bacteria that are capable of degrading organic contaminants. The objective of this study was to identify low-cost soil solid amendments that could retard the movement of trichloroethylene (TCE) while serving as long-lived carbon sources to foster its biodegradation in shallow groundwater through the use of a PRB. The natural amendments high in organic carbon content such as eucalyptus mulch, compost, wetland peat, organic humus were compared based on their geophysical characteristics, such as pHw, porosity and total organic carbon (TOC), and as well as TCE sorption potentials. The pHw values were within neutral range except for pine bark mulch and wetland peat. All other geophysical characteristics of the amendments showed suitability for use in a PRB. While the Freundlich model showed better fit for compost and pine bark mulch, the linear sorption model was adequate for eucalyptus mulch, wetland peat and Everglades muck within the concentration range studied (0.2-0.8 mg/L TCE). According to these results, two composts and eucalyptus mulch were selected for laboratory column experiments to evaluate their effectiveness at creating and maintaining conditions suitable for TCE anaerobic dechlorination. The columns were monitored for pH, ORP, TCE degradation, longevity of nutrients and soluble TOC to support TCE dechlorination. Native bacteria in the columns had the ability to convert TCE to DCEs; however, the inoculation with the TCE-degrading culture greatly increased the rate of biodegradation. This caused a significant increase in by-product concentration, mostly in the form of DCEs and VC followed by a slow degradation to ethylene. Of the tested amendments eucalyptus mulch was the most effective at supporting the TCE dechlorination. The experimental results of TCE sequential dechlorination took place in eucalyptus mulch and commercial compost from Savannah River Site columns were then simulated using the Hydrus-1D model. The simulations showed good fit with the experimental data. The results suggested that sorption and degradation were the dominant fate and transport mechanisms for TCE and DCEs in the column, supporting the use of these amendments in a permeable reactive barrier to remediate the TCE.

Quantifying the effects of temperature and concentration on variable-density flow in numerical modeling of groundwater systems: Implications for predictive uncertainty and data collection

Groundwater systems of different densities are often mathematically modeled to understand and predict environmental behavior such as seawater intrusion or submarine groundwater discharge. Additional data collection may be justified if it will cost-effectively aid in reducing the uncertainty of a model's prediction. The collection of salinity, as well as, temperature data could aid in reducing predictive uncertainty in a variable-density model. However, before numerical models can be created, rigorous testing of the modeling code needs to be completed. This research documents the benchmark testing of a new modeling code, SEAWAT Version 4. The benchmark problems include various combinations of density-dependent flow resulting from variations in concentration and temperature. The verified code, SEAWAT, was then applied to two different hydrological analyses to explore the capacity of a variable-density model to guide data collection. ^ The first analysis tested a linear method to guide data collection by quantifying the contribution of different data types and locations toward reducing predictive uncertainty in a nonlinear variable-density flow and transport model. The relative contributions of temperature and concentration measurements, at different locations within a simulated carbonate platform, for predicting movement of the saltwater interface were assessed. Results from the method showed that concentration data had greater worth than temperature data in reducing predictive uncertainty in this case. Results also indicated that a linear method could be used to quantify data worth in a nonlinear model. ^ The second hydrological analysis utilized a model to identify the transient response of the salinity, temperature, age, and amount of submarine groundwater discharge to changes in tidal ocean stage, seasonal temperature variations, and different types of geology. The model was compared to multiple kinds of data to (1) calibrate and verify the model, and (2) explore the potential for the model to be used to guide the collection of data using techniques such as electromagnetic resistivity, thermal imagery, and seepage meters. Results indicated that the model can be used to give insight to submarine groundwater discharge and be used to guide data collection. ^

Use of tritium and helium to define groundwater flow conditions in Everglades National Park

The concentrations of tritium (3H) and helium isotopes (3He and4He) were used as tracers of groundwater flow in the surficial aquifer system (SAS) beneath Everglades National Park (ENP), south Florida. From ages determined by 3H/3He dating techniques, groundwater within the upper 28 m originated within the last 30 years. Below 28 m, waters originated prior to 30 years before present with evidence of mixing at the interface. Interannual variation of the 3H/3He ages within the upper 28 m was significant throughout the 3 year investigation, corresponding with varying hydrologic conditions. In the region of Taylor Slough Bridge, younger groundwater was consistently detected below older groundwater in the Biscayne Aquifer, suggesting preferential flow to the lower part of the aquifer. An increase in 4He with depth in the SAS indicated that radiogenic 4He produced in the underlying Hawthorn Group migrates into the SAS by diffusion. Higher Δ4He values in brackish groundwaters compared to fresh waters from similar depths suggested a possible enhanced vertical transport of4He in the seawater mixing zone. Groundwater salinity measurements indicated the presence of a wide (6–28 km) seawater mixing zone. Comparison of groundwater levels with surface water levels in this zone indicated the potential for brackish groundwater discharge to the overlying Everglades surface water.

Geochemical determination of the fate and transport of injected fresh wastewater to a deep saline aquifer

Deep well injection into non-potable saline aquifers of treated domestic wastewater has been used in Florida for decades as a safe and effective alternative to ocean outfall disposal. The objectives of this study were to determine the fate and transport of injected wastewater at two deep well injection sites in Miami Dade County, Florida, USA. Detection of ammonium in the Middle Confining units of the Floridan aquifer above the injection zone at both sites has been interpreted as evidence of upward migration of injected wastewater, posing a risk to underground sources of drinking water. Historical water quality data, including ammonia, chloride, temperature, and pH from existing monitoring wells at both sites from 1983 to 2008, major ions collected monthly from 2006 and 2008, and a synoptic sampling event for stable isotopes, tritium, and dissolved gases in 2008, were used to determine the source of ammonium in groundwater and possible migration pathways. Geochemical modeling was used to determine possible effects of injected wastewater on native water and aquifer matrix geochemistry. Injected wastewater was determined to be the source of elevated ammonium concentrations above ambient water levels, based on the results of major ion concentrations, tritium, dissolved noble gases and 15N isotopes analyses. Various possible fluid migration pathways were identified at the sites. Data for the south site suggest buoyancy-driven vertical pathways to overlying aquifers bypassing the confining units, with little mixing of injected wastewater with native water as it migrated upward. Once it is introduced into an aquifer, the injectate appeared to migrate advectively with the regional groundwater flow. Geochemical modeling indicated that CO 2-enriched injected wastewater allowed for carbonate dissolution along the vertical pathways, enhancing permeability along these flowpaths. At the north site, diffusive upward flow through the confining units or offsite vertical pathways were determined to be possible, however no evidence was detected for any on-site confining unit bypass pathway. No evidence was observed at either site of injected wastewater migration to the Upper Floridan aquifer, which is used as a municipal water supply and for aquifer storage and recovery.

Arsenic transport by zebrafish aquaglyceroporins

Background Arsenic is one of the most ubiquitous toxins and endangers the health of tens of millions of humans worldwide. It is a mainly a water-borne contaminant. Inorganic trivalent arsenic (AsIII) is one of the major species that exists environmentally. The transport of AsIII has been studied in microbes, plants and mammals. Members of the aquaglyceroporin family have been shown to actively conduct AsIII and its organic metabolite, monomethylarsenite (MAsIII). However, the transport of AsIII and MAsIII in in any fish species has not been characterized. Results In this study, five members of the aquaglyceroporin family from zebrafish (Danio rerio) were cloned, and their ability to transport water, glycerol, and trivalent arsenicals (AsIII and MAsIII) and antimonite (SbIII) was investigated. Genes for at least seven aquaglyceroporins have been annotated in the zebrafish genome project. Here, five genes which are close homologues to human AQP3, AQP9 and AQP10 were cloned from a zebrafish cDNA preparation. These genes were namedaqp3, aqp3l, aqp9a, aqp9b and aqp10 according to their similarities to the corresponding human AQPs. Expression of aqp9a, aqp9b, aqp3, aqp3l and aqp10 in multiple zebrafish organs were examined by RT-PCR. Our results demonstrated that these aquaglyceroporins exhibited different tissue expression. They are all detected in more than one tissue. The ability of these five aquaglyceroporins to transport water, glycerol and the metalloids arsenic and antimony was examined following expression in oocytes from Xenopus leavis. Each of these channels showed substantial glycerol transport at equivalent rates. These aquaglyceroporins also facilitate uptake of inorganic AsIII, MAsIII and SbIII. Arsenic accumulation in fish larvae and in different tissues from adult zebrafish was studied following short-term arsenic exposure. The results showed that liver is the major organ of arsenic accumulation; other tissues such as gill, eye, heart, intestine muscle and skin also exhibited significant ability to accumulate arsenic. The zebrafish larvae also accumulate considerable amounts of arsenic. Conclusion This is the first molecular identification of fish arsenite transport systems and we propose that the extensive expression of the fish aquaglyceroporins and their ability to transport metalloids suggests that aquaglyceroporins are the major pathways for arsenic accumulation in a variety of zebrafish tissues. Uptake is one important step of arsenic metabolism. Our results will contribute to a new understanding of aquatic arsenic metabolism and will support the use of zebrafish as a new model system to study arsenic associated human diseases.

Influence of hydrobiogeochemistry on transport of chromium through manganese -containing sediments: Experimental and modeling approaches

Chromium (Cr) is a metal of particular environmental concern, owing to its toxicity and widespread occurrence in groundwater, soil, and soil solution. A combination of hydrological, geochemical, and microbiological processes governs the subsurface migration of Cr. Little effort has been devoted to examining how these biogeochemical reactions combine with hydrologic processes influence Cr migration. This study has focused on the complex problem of predicting the Cr transport in laboratory column experiments. A 1-D reactive transport model was developed and evaluated against data obtained from laboratory column experiments. ^ A series of dynamic laboratory column experiments were conducted under abiotic and biotic conditions. Cr(III) was injected into columns packed with β-MnO 2-coated sand at different initial concentrations, variable flow rates, and at two different pore water pH (3.0 and 4.0). In biotic anaerobic column experiments Cr(VI) along with lactate was injected into columns packed with quartz sand or β-MnO2-coated sand and bacteria, Shewanella alga Simidu (BrY-MT). A mathematical model was developed which included advection-dispersion equations for the movement of Cr(III), Cr(VI), dissolved oxygen, lactate, and biomass. The model included first-order rate laws governing the adsorption of each Cr species and lactate. The equations for transport and adsorption were coupled with nonlinear equations for rate-limited oxidation-reduction reactions along with dual-monod kinetic equations. Kinetic batch experiments were conducted to determine the reduction of Cr(VI) by BrY-MT in three different substrates. Results of the column experiments with Cr(III)-containing influent solutions demonstrate that β-MnO2 effectively catalyzes the oxidation of Cr(III) to Cr(VI). For a given influent concentration and pore water velocity, oxidation rates are higher, and hence effluent concentrations of Cr(VI) are greater, at pH 4 relative to pH 3. Reduction of Cr(VI) by BrY-MT was rapid (within one hour) in columns packed with quartz sand, whereas Cr(VI) reduction by BrY-MT was delayed (57 hours) in presence of β-MnO 2-coated sand. BrY-MT grown in BHIB (brain heart infusion broth) reduced maximum amount of Cr(VI) to Cr(III) followed by TSB (tryptic soy broth) and M9 (minimum media). The comparisons of data and model results from the column experiments show that the depths associated with Cr(III) oxidation and transport within sediments of shallow aquatic systems can strongly influence trends in surface water quality. The results of this study suggests that carefully performed, laboratory column experiments is a useful tool in determining the biotransformation of redox-sensitive metals even in the presence of strong oxidant, like β-MnO2. ^

The influence of bedrock hydrogeology on catchment-scale nitrate fate and transport in fractured aquifers

Characterising catchment scale biogeochemical processes controlling nitrate fate in groundwater constitutes a fundamental consideration when applying programmes of measures to reduce risks posed by diffuse agricultural pollutants to water quality. Combining hydrochemical analyses with nitrate isotopic data and physical hydrogeological measurements permitted characterisation of biogeochemical processes influencing nitrogen fate and transport in the groundwater in two fractured bedrock aquifers with contrasting hydrogeology but comparable nutrient loads. Hydrochemical and isotopic analyses of groundwater samples collected from moderately fractured, diffusely karstified limestone indicated nitrification controlled dissolved nitrogen fate and delivery to aquatic receptors. By contrast nitrate concentrations in groundwater were considerably lower in a low transmissivity highly lithified sandstone and pyrite-bearing shale unit with patchy subsoil cover. Geophysical and hydrochemical investigations showed shallower intervals contained hydraulically active fractures where denitrification was reflected through lower nitrogen levels and an isotopic enrichment ratio of 1.7 between δ15N and δ18O. Study findings highlight the influence of bedrock hydrogeological conditions on aqueous nitrogen mobility. Investigation results demonstrate that bedrock conditions need to be considered when implementing catchment management plans to reduce the impact of agricultural practices on the quality of groundwater and baseflow in receiving rivers.

EFFECTS OF BUBBLE-FACILITATED VOC TRANSPORT FROM LNAPL SMEAR ZONES ON VAPOR INTRUSION

Light non-aqueous phase liquid (LNAPL) sources can pose a significant threat to indoor air through vapour intrusion (VI). Most conceptual and numerical models of VI assume that the transport of volatile organic compounds (VOCs) is a diffusion-limited process. Recently, alternate conditions have been identified that could lead to faster transport, including the presence of preferential pathways and methanogenic gas production. In this study, an additional mechanism that could lead to faster transport was investigated: bubble-facilitated VOC transport from LNAPL smear zones. A laboratory investigation was preformed using pentane in one-dimensional laboratory columns and two-dimensional visualization experiments. Results of the column experiments showed that average VOC mass fluxes in the bubble-facilitated columns were over two orders of magnitude greater than in the diffusion-limited columns. In addition, the flux signal was intermittent, consistent with expectations of bubble-facilitated transport as bubbles expand, mobilize and are released to the vadose zone at various times during the test. The results from the visualization experiments showed gas fingers growing and mobilizing over time, which supports the findings of the column experiments. In conclusion, these results demonstrate the potential for bubble-facilitated VOC transport to affect mass transfer in LNAPL smear zones, and lead to increased indoor air concentrations by VI.

2013 survey of Iowa groundwater and evaluation of public well vulnerability classifications for contaminants of emerging concern (2015)

Studies in Iowa have long documented the vulnerability of wells with less than 50 feet (15 meters) of confining materials above the source aquifer to contamination from nitrate and various pesticides. Recent studies in Wisconsin have documented the occurrence of viruses in untreated groundwater, even in wells considered to have little vulnerability to contamination from near-surface activities. In addition, sensitive methods have become available for analyses of pharmaceuticals and pesticides. This study represents the first comprehensive examination of contaminants of emerging concern in Iowa’s groundwater conducted to date, and one of the first conducted in the United States. Raw groundwater samples were collected from 66 public supply wells during the spring of 2013, when the state was recovering from drought conditions. Samples were analyzed for 206 chemical and biological parameters; including 20 general water-quality parameters and major ions, 19 metals, 5 nutrients, 10 virus groups, 3 species of pathogenic bacteria, 5 microbial indicators, 108 pharmaceuticals, 35 pesticides and pesticide degradates, and tritium. The wells chosen for this study represent a diverse range of ages, depths, confining material thicknesses, pumping rates, and land use settings. The most commonly detected contaminant group was pesticide compounds, which were present in 41% of the samples. As many as 6 pesticide compounds were found together in a sample, most of which were chloroacetanilide degradates. While none of the measured concentrations of pesticide compounds exceeded current benchmark levels, several of these compounds are listed on the U.S. Environmental Protection Agency’s Contaminant Candidate List and could be subject to drinking water standards in the future. Despite heavy use in the past decade, glyphosate was not detected, and its metabolite, aminomethylphosphonic acid, was only detected in two of 60 wells tested (3%) at the detection limit of 0.02 μg/L. Pharmaceutical compounds were detected in 35% of 63 samples. Of the 14 pharmaceuticals detected, six had reported concentrations above the method reporting limit, with the maximum reported concentration of 826 ng/L for acetaminophen. Diphenhydramine was the only pharmaceutical to have two detections above the reporting limit, at 24.5 and 145 ng/L. Eight pharmaceuticals had confirmed detections at concentrations below the method reporting limit. Caffeine was the most frequently detected pharmaceutical compound (25%), followed by the caffeine metabolite, 1,7-dimethylxanthine (16%). Microorganisms were detected in 21% of the wells using quantitative polymerase chain reaction methodologies. The most frequently detected microorganism was the pepper mild mottle virus (PMMV), a plant pathogen found in human waste. PMMV was detected in 17% of samples at concentrations ranging from 0.4 to 6.38 gene copies per liter. GII norovirus, human polyomavirus, bovine polyomavirus, and Campylobacter were also detected, while adenovirus, enterovirus, GI norovirus, swine hepatitis E, Salmonella, and enterohemmorhagic E. coli were not detected. No correlations were found between viruses or pathogenic bacteria and microbial indicators. Wells with less than 50 feet (15 meters) of confining material were shown to have greater incidence of surface-related contaminants; however, significant relationships (p<0.05) between confining layer thickness and contaminants were only found for nitrate and herbicides.

Nitrate in groundwater in a recharge area of Guarany aquifer in Brazil.

The region of Ribeirão Preto City located in São Paulo State, southeastern Brazil, is an important sugarcane, soybean and corn producing area. This region is also an important recharge area (Espraiado) for groundwater of the Guarany aquifer, a water supply source for the city and region. It has an intercontinental extension that comprises areas of eight Brazilian states, as well as significant portions of other South American countries like Argentina, Uruguay, and Paraguay, with a total area of approximately 1,200,000 Km2. Due to the high permeability of some soils present in this region, the high mobility of the herbicides and fertilizers applied, and being a recharge area, it is important to investigate the potential transport of applied fertilizers to underlying aquifer. The cultivation sugar cane in this area demands the frequent use of nitrogen as fertilizer. This research was conducted to characterize the potential contamination of groundwater with nitrogen in the recharge area of groundwater. Seven groundwater sample points were selected in the Espraiado stream watershed, during the years of 2005 and 2006. Samples were collected during the months of March, July, and December of each year. Three replications were collected at each site. Groundwater was also collected during the same months from county groundwater wells located throughout the city. The following six wells were studied: Central, Palmares, Portinari, Recreio Internacional, São Sebastião, and São José. Nitrate water samples were analyzed by Cadmium Reduction Method. No significant amount of nitrate was found in the recharge, agricultural, area. However, nitrate levels were detected at concentrations higher than the Maximum Concentration Level (MCL) of 10mg/L in downtown, urban, well located away from agricultural sites with no history of fertilizer or nitrogen application.

Testing the effectiveness of phytoscreening to monitor shallow groundwater contamination by chlorinated ethenes

The research presented herein aims to investigate the strengths and weaknesses of a relatively new technique called phytoscreening. Parallel to the well-known phytoremediation, it consists of exploiting the absorbing potential of trees to delineate groundwater contamination plumes, especially for chlorinated ethenes (i.e., PCE, TCE, 1,2-cis DCE, and VC). The latter are prevalent contaminants in groundwater but their fate and transport in surface ecosystems, such as trees, are still poorly understood and subjected to high variability. Moreover, the analytical validity of tree-coring is still limited in many countries due to a lack of knowledge of its application opportunities. Tree-cores are extracted from trunks and generally analyzed by gas chromatography/mass spectrometry. A systematic review of former literature on phytoscreening for chlorinated ethenes is presented in this PhD thesis to evaluate the factors influencing the effectiveness of the technique. Besides, we tested the technique by probing eight sites contaminated by chlorinated ethenes in Italy (Emilia-Romagna) in different hydrogeological and seasonal settings. We coupled the technique with the assessment of gaseous-phase concentrations directly on-site, inserting detector tubes or a photoionization detector in the tree-holes left by the coring tool. Finally, we applied rank order statistic analysis on field data along with literature data to assess under which conditions phytoscreening should be applied to either screen or monitor environmental contamination issues. A relatively high correlation exists between tree-core and groundwater concentrations (Spearman’s ρ > 0.6), being higher for compounds with higher sorption, for sites with shallower and thinner aquifers, and when sampling specific tree types with standardized sampling and extraction protocols. These results indicate the opportunities for assessing the occurrence, type, and concentration of solvents directly from the stem of trees. This can reduce the costs of characterization surveys, allowing rapid identification of hotspots and plume direction and thus optimizing the drilling of boreholes.