A field and reactive transport model study of arsenic in a basaltic rock aquifer

The use of geothermal energy as a source for electricity and district heating has increased over recent decades. Dissolved As can be an important constituent of the geothermal fluids brought to the Earth's surface. Here the field application of laboratory measured adsorption coefficients of aqueous As species on basaltic glass surfaces is discussed. The mobility of As species in the basaltic aquifer in the Nesjavellir geothermal system, Iceland was modelled by the one-dimensional (1D) reactive transport model PHREEQC ver. 2, constrained by a long time series of field measurements with the chemical composition of geothermal effluent fluids, pH, Eh and, occasionally, Fe- and As-dissolved species measurements. Di-, tri- and tetrathioarsenic species (As(OH)S22-, AsS₃H^2-, AsS33- and As(SH)4-) were the dominant form of dissolved As in geothermal waters exiting the power plant (2.556μM total As) but converted to some extent to arsenite (H₃AsO₃) and arsenate HAsO42- oxyanions coinciding with rapid oxidation of S2- to S2O32- and finally to SO42- during surface runoff before feeding into a basaltic lava field with a total As concentration of 0.882μM following dilution with other surface waters. A continuous 25-a data set monitoring groundwater chemistry along a cross section of warm springs on the Lake Thingvallavatn shoreline allowed calibration of the 1D model. Furthermore, a series of ground water wells located in the basaltic lava field, provided access along the line of flow of the geothermal effluent waters towards the lake. The conservative ion Cl^- moved through the basaltic lava field (4100m) in less than10a but As was retarded considerably due to surface reactions and has entered a groundwater well 850m down the flow path as arsenate in accordance to the prediction of the 1D model. The 1D model predicted a complete breakthrough of arsenate in the year 2100. In a reduced system arsenite should be retained for about 1ka. © 2011 Elsevier Ltd.

Characterizing Aquifer Heterogeneity Using Bacterial and Bacteriophage Tracers

Gravel aquifers act as important potable water sources in central western Europe yet they are subject to numerous contamination pressures. Compositional and textural heterogeneity makes protection zone delineation around groundwater supplies in these units challenging; artificial tracer testing aids characterization. This paper reappraises previous tracer test results in light of new geological and microbiological data. Comparative passive gradient testing, using a fluorescent solute (Uranine), virus (H40/1 bacteriophage), and comparably sized bacterial tracers Escherichia coli and Pseudomonas putida, was used to investigate a calcareous gravel aquifer’s ability to remove microbiological contaminants at a test site near Munich, Germany. Test results revealed E. coli relative recoveries could exceed those of H40/1 at monitoring wells 10 m and 20 m from an injection well by almost four times; P. putida recoveries varied by a factor of up to three between wells. Application of filtration theory suggested greater attenuation of H40/1 relative to similarly charged E. coli occurred due to differences in microorganism size, while estimated collision efficiencies appeared comparable. By contrast, more positively charged P. putida experienced greater attenuation at one monitoring point, while lower attenuation rates at the second location indicated the influence of geochemical heterogeneity. Test findings proved consistent with observations from nearby fresh outcrops that suggested thin open framework gravel beds dominated mass transport in the aquifer, while discrete intervals containing stained clasts reflect localized geochemical heterogeneity. Study results highlight the utility of reconciling outcrop observations with artificial tracer test responses, using microbiological tracers with well-defined properties, to characterize aquifer heterogeneity.

Fate and transport of volatile organic compounds in glacial till and groundwater at an industrial site in Northern Ireland

Volatile organic compound (VOC) contamination of subsurface geological material and groundwater was discovered on the Nortel Monkstown industrial site, Belfast, Northern Ireland. The objectives of this study were to (1) investigate the characteristics of the geological material and its influences on contaminated groundwater flow across the site using borehole logs and hydrological evaluations, and (2) identify the contaminants and examine their distribution in the subsurface geological material and groundwater using chemical analysis. This report focuses on the eastern car park (ECP) which was a former storage area associated with trichloroethene (TCE) degreasing operations. This is where the greatest amount of volatile organic compounds (VOCs), particularly TCE, were detected. The study site is on a complex deposit of clayey glacial till with discontinuous coarser grained lenses, mainly silts, sands and gravel, which occur at 0.45–7.82 m below ground level (bgl). The lenses overall form an elongated formation that acts as a small unconfined shallow aquifer. There is a continuous low permeable stiff clayey till layer beneath the lenses that performs as an aquitard to the groundwater. Highest concentrations of VOCs, mainly TCE, in the geological material and groundwater are in these coarser lenses at ~4.5–7 m bgl. Highest TCE measurements at 390,000 µg L-1 for groundwater and at 39,000 µg kg-1 at 5.7 m for geological material were in borehole GA19 in the coarse lens zone. It is assumed that TCE gained entrance to the subsurface near this borehole where the clayey till was thin to absent above coarse lenses which provided little retardation to the vertical migration of this dense non-aqueous phase liquid (DNAPL) into the groundwater. However, TCE is present in low concentrations in the geological material overlying the coarse lens zone. Additionally, VOCs appear to be associated with poorly drained layers and in peat

Anatomy and facies association of a drumlin in Co. Down, Northern Ireland, from seismic and electrical resistivity surveys

Seismic refraction and electrical resistivity geophysical techniques were used to reconstruct the internal architecture of a drumlin in Co. Down, Northern Ireland. Geophysical results were both validated and complemented by borehole drilling, ground water flow modelling, and geologic mapping. The geophysical anatomy of the drumlin consists of five successive layers with depth including; topsoil, partially saturated and saturated glacial tills, and weathered and more competent greywacke bedrock. There are numerous, often extensive inclusions of clay, sand, gravel, cobbles, and boulders within the topsoil and the till units. Together geophysical and geotechnical findings imply that the drumlin is part of the subglacial lodgement, melt-out, debris flow, sheet flow facies described by previous authors, and formed by re-sedimentation and streamlining of pre-existing sediments during deglaciation of the Late Devensian ice sheet. Seismic refraction imaging is particularly well suited to delineating layering within the drumlin, and is able to reconstruct depths to interfaces to within ± 0.5 m accuracy. Refraction imaging ascertained that the weathered bedrock layer is continuous and of substantial thickness, so that it acts as a basal aquifer which underdrains the bulk of the drumlin. Electrical resistivity imaging was found to be capable of delineating relative spatial changes in the moisture content of the till units, as well as mapping sedimentary inclusions within the till. The moisture content appeared to be elevated near the margins of the drumlin, which may infer a weakening of the drumlin slopes. Our findings advocate the use of seismic refraction and electrical resistivity methods in future sedimentological and geotechnical studies of internal drumlin architecture and drumlin formation, owing particularly to the superior, 3- D spatial coverage of these methods.

Comparative tracing experiments to investigate epikarst structural and compositional heterogeneity

Comparative tracer testing may be used to evaluate the vulnerability of groundwater to specific contaminants by comparing reactive tracer response to that of a simultaneously injected non-reactive “conservative” substance. Conversely, knowledge of tracer reaction with specific materials permits information about subsurface heterogeneity to be inferred. A series of tests completed in the vadose zone overlying a limestone aquifer employed a cocktail of particles along with reactive and non-reactive solute tracers to investigate transport rates between the ground surface and monitoring points approximately 10 m below ground. Short pulse tests revealed both solutes and particulate contaminants could travel at rates of over 10 m/h. Comparison of particle (microorganisms) and non-reactive solute tracer breakthrough revealed that particle tracers experience pore exclusion resulting in higher peak relative concentrations which arrive earlier than those of the solute. Prolonged tracer injection during subsequent experiments confirmed the response observed and illustrated that over 40 % of flow paths between injection and monitoring points were inaccessible to particles, but could allow solutes to pass through them. Similarly, the difference in response between various reactive tracers demonstrated tracers reached monitoring points via multiple flow paths and suggests geochemical heterogeneity plays an important role in influencing tracer behaviour. The results of this investigation highlight the complexity of water flow through the epikarst and the vulnerability of groundwater in karst aquifers to contamination when soil cover is thin to absent.

The typology of Irish hard rock aquifers based on an integrated hydrogeological and geophysical approach

Groundwater flow in hard-rock aquifers is strongly controlled by the characteristics and distribution of structural heterogeneity. A methodology for catchment-scale characterisation is presented, based on the integration of complementary, multi-scale hydrogeological, geophysical and geological approaches. This was applied to three contrasting catchments underlain by metamorphic rocks in the northern parts of Ireland (Republic of Ireland and Northern Ireland, UK). Cross-validated surface and borehole geophysical investigations confirm the discontinuous overburden, lithological compartmentalisation of the bedrock and important spatial variations of the weathered bedrock profiles at macro-scale. Fracture analysis suggests that the recent (Alpine) tectonic fabric exerts strong control on the internal aquifer structure at meso-scale, which is likely to impact on the anisotropy of aquifer properties. The combination of the interpretation of depth-specific hydraulic-test data with the structural information provided by geophysical tests allows characterisation of the hydrodynamic properties of the identified aquifer units. Regionally, the distribution of hydraulic conductivities can be described by inverse power laws specific to the aquifer litho-type. Observed groundwater flow directions reflect this multi-scale structure. The proposed integrated approach applies widely available investigative tools to identify key dominant structures controlling groundwater flow, characterising the aquifer type for each catchment and resolving the spatial distribution of relevant aquifer units and associated hydrodynamic parameters.

Detection of Mixing Dynamics During Pumping of a Flooded Coal Mine

In complex hydrogeological environments the effective management of groundwater quality problems by pump-and-treat operations can be most confidently achieved if the mixing dynamics induced within the aquifer by pumping are well understood. The utility of isotopic environmental tracers (C-, H-, O-, S-stable isotopic analyses and age indicators—14C, 3H) for this purpose is illustrated by the analysis of a pumping test in an abstraction borehole drilled into flooded, abandoned coal mineworkings at Deerplay (Lancashire, UK). Interpretation of the isotope data was undertaken conjunctively with that of major ion hydrochemistry, and interpreted in the context of the particular hydraulic setting of flooded mineworkings to identify the sources and mixing of water qualities in the groundwater system. Initial pumping showed breakdown of initial water quality stratification in the borehole, and gave evidence for distinctive isotopic signatures (d34S(SO4) ~= -1.6‰, d18O(SO4) ~= +15‰) associated with primary oxidation of pyrite in the zone of water table fluctuation—the first time this phenomenon has been successfully characterized by these isotopes in a flooded mine system. The overall aim of the test pumping—to replace an uncontrolled outflow from a mine entrance in an inconvenient location with a pumped discharge on a site where treatment could be provided—was swiftly achieved. Environmental tracing data illustrated the benefits of pumping as little as possible to attain this aim, as higher rates of pumping induced in-mixing of poorer quality waters from more distant old workings, and/or renewed pyrite oxidation in the shallow subsurface.

Arsenic mobilization from iron oxyhydroxides is regulated by organic matter carbon to nitrogen (C:N) ratio

Arsenic (As) is mobilized from delta and floodplain aquifer sediments throughout S.E. Asia via reductive dissolution of As bound to iron (Fe) oxyhydroxides. The reductive driving force is organic carbon, but its source and constitution is uncertain. Here batch incubation experiments were conducted to investigate the role of organic matter (OM) carbon:nitrogen (C:N) ratio on the mobilization of arsenic, Fe and N from As dosed, Fe oxyhydroxide coated sands. As mobilization into pore waters from the sand was strongly regulated by the C:N ratio of the OM, and also the concentration of OM present. The lower the C:N, the more As released. Fe and ammonium release were similarly dependent on the quality and quantity of OM, but Fe mobilization was more rapid and ammonium release slower than As suggesting that the mobilization of these 3 moieties although interdependent, were not directly linked. It was concluded that low C:N ratios for OM responsible for reducing aquifers were As in groundwater is observed were likely.

Codeposition of organic carbon and arsenic in Bengal Delta aquifers

We present data showing that arsenic (As) was codeposited with organic carbon (OC) in Bengal Delta sediments as As and OC concentrations are highly (p <0.001) positively correlated in core profiles collected from widely dispersed geographical sites with different sedimentary depositional histories. Analysis of modern day depositional environments revealed that the As-OC correlations observed in cores are due to As retention and high OC inputs in vegetated zones of the deltaic environment. We hypothesize that elevated concentrations of As occur in vegetated wetland sediments due to concentration and retention of arsenate in aerated root zones and animal burrows where copious iron(III) oxides are deposited. On burial of the sediment, degradation of organic carbon from plant and animal biomass detritus provides the reducing conditions to dissolve iron(III) oxides and release arsenite into the porewater. As tubewell abstracted aquifer water is an invaluable resource on which much of Southeast Asia is now dependent, this increased understanding of the processes responsible for As buildup and release will identify, through knowledge of the palaeosedimentary environment, which sediments are at most risk of having high arsenic concentrations in porewater. Our data allow the development of a new unifying hypothesis of how As is mobilized into groundwaters in river flood plains and deltas of Southeast Asia, namely that in these highly biologically productive environments, As and OC are codeposited, and the codeposited OC drives As release from the sediments.

Co-deposition of organic carbon and arsenic in Bengal Delta aquifers.

We present data showing that arsenic (As) was codeposited with organic carbon (OC) in Bengal Delta sediments as As and OC concentrations are highly (p < 0.001) positively correlated in core profiles collected from widely dispersed geographical sites with different sedimentary depositional histories. Analysis of modern day depositional environments revealed that the As/OC correlations observed in cores are due to As retention and high OC inputs in vegetated zones of the deltaic environment. We hypothesize that elevated concentrations of As occur in vegetated wetland sediments due to concentration and retention of arsenate in aerated root zones and animal burrows where copious iron(III) oxides are deposited. On burial of the sediment, degradation of organic carbon from plant and animal biomass detritus provides the reducing conditions to dissolve iron(III) oxides and release arsenite into the porewater. As tubewell abstracted aquifer water is an invaluable resource on which much of Southeast Asia is now dependent, this increased understanding of the processes responsible for As buildup and release will identify, through knowledge of the palaeosedimentary environment, which sediments are at most risk of having high arsenic concentrations in porewater. Our data allow the development of a new unifying hypothesis of how As is mobilized into groundwaters in river flood plains and deltas of Southeast Asia, namely that in these highly biologically productive environments, As and OC are codeposited, and the codeposited OC drives As release from the sediments. We present data showing that arsenic (As) was codeposited with organic carbon (OC) in Bengal Delta sediments as As and OC concentrations are highly (p < 0.001) positively correlated in core profiles collected from widely dispersed geographical sites with different sedimentary depositional histories. Analysis of modern day depositional environments revealed that the As?OC correlations observed in cores are due to As retention and high OC inputs in vegetated zones of the deltaic environment. We hypothesize that elevated concentrations of As occur in vegetated wetland sediments due to concentration and retention of arsenate in aerated root zones and animal burrows where copious iron(III) oxides are deposited. On burial of the sediment, degradation of organic carbon from plant and animal biomass detritus provides the reducing conditions to dissolve iron(III) oxides and release arsenite into the porewater. As tubewell abstracted aquifer water is an invaluable resource on which much of Southeast Asia is now dependent, this increased understanding of the processes responsible for As buildup and release will identify, through knowledge of the palaeosedimentary environment, which sediments are at most risk of having high arsenic concentrations in porewater. Our data allow the development of a new unifying hypothesis of how As is mobilized into groundwaters in river flood plains and deltas of Southeast Asia, namely that in these highly biologically productive environments, As and OC are codeposited, and the codeposited OC drives As release from the sediments.

Alternate wetting and drying irrigation for rice in Bangladesh: Is it sustainable and has plant breeding something to offer?

The crop management practice of alternate wetting and drying (AWD) is being promoted by IRRI and the national research and extension program in Bangladesh and other parts of the world as a water-saving irrigation practice that reduces the environmental impact of dry season rice production through decreased water usage, and potentially increases yield. Evidence is growing that AWD will dramatically reduce the concentration of arsenic in harvested rice grains conferring a third major advantage over permanently flooded dry season rice production. AWD may also increase the concentration of essential dietary micronutrients in the grain. However, three crucial aspects of AWD irrigation require further investigation. First, why is yield generally altered in AWD? Second, is AWD sustainable economically (viability of farmers' livelihoods) and environmentally (aquifer water table heights) over long-term use? Third, are current cultivars optimized for this irrigation system? This paper describes a multidisciplinary research project that could be conceived which would answer these questions by combining advanced soil biogeochemistry with crop physiology, genomics, and systems biology. The description attempts to show how the breakthroughs in next generation sequencing could be exploited to better utilize local collections of germplasm and identify the molecular mechanisms underlying biological adaptation to the environment within the context of soil chemistry and plant physiology.

Numerical assessment of potential impacts of hydraulically fractured Bowland Shale on overlying aquifers

Natural gas extracted from hydraulically fractured shale formations potentially has a big impact on the global energy landscape. However, there are concerns of potential environmental impacts of hydraulic fracturing of the shale formations, particularly those related to water quality. To evaluate the potential impact of hydraulically fractured shale on overlying aquifers, we conduct realizations of numerical modeling simulations to assess fluid flow and chloride transport from a synthetic Bowland Shale over a period of 11,000 years. The synthetic fractured shale was represented by a three-dimensional discrete fracture model that was developed by using the data from a Bowland Shale gas exploration in Lancashire, UK. Chloride mass exchange between fractures and the rock matrix was fully accounted for in the model. The assessment was carried out to investigate fluid and chloride mass fluxes before, during, and after hydraulic fracturing of the Bowland Shale. Impacts of the upward fracture height and aperture, as well as hydraulic conductivity of the multilayered bedrock system, are also included this assessment. This modeling revealed that the hydraulically fractured Bowland Shale is unlikely to pose a risk to its overlying groundwater quality when the induced fracture aperture is ≤200 µm. With the fracture aperture ≥1000 µm, the upward chloride flux becomes very sensitive to the upward fracture height growth and hydraulic conductivity of the multilayered bedrock system. In the extremely unlikely event of the upward fracture growth directly connecting the shale formation to the overlying Sherwood Sandstone aquifer with the fracture aperture ≥1000 µm, the upward chloride mass flux could potentially pose risks to the overlying aquifer in 100 years. The model study also revealed that the upward mass flux is significantly intercepted by the horizontal mass flux within a high permeable layer between the Bowland Shale and its overlying aquifers, reducing further upward flux toward the overlying aquifers.

Geophysical and hydrogeological characterisation of the impacts of on-site wastewater treatment discharge to groundwater in a poorly productive bedrock acquifer

Contaminants discharging from on-site wastewater treatment systems (OSWTSs) can impact groundwater quality, threatening human health and surface water ecosystems. Risk of negative impacts becomes elevated in areas of extreme vulnerability with high water tables, where thin unsaturated intervals limit vadose zone attenuation. A combined geophysical/hydrogeological investigation into the effects of an OSWTS, located over a poorly productive aquifer (PPA) with thin subsoil cover, aimed to characterise effluent impacts on groundwater. Groundwater, sampled from piezometers down-gradient of the OSWTS percolation area displayed spatially erratic, yet temporally consistent, contaminant distributions. Electrical resistivity tomography identified an area of gross groundwater contamination close to the percolation area and, when combined with seismic refraction and water quality data, indicated that infiltrating effluent reaching the water table discharged to a deeper more permeable zone of weathered shale resting on more competent bedrock. Subsurface structure, defined by geophysics, indicated that elevated chemical and microbiological contaminant levels encountered in groundwater samples collected from piezometers, down-gradient of sampling points with lower contaminant levels, corresponded to those locations where piezometers were screened close to the weathered shale/competent rock interface; those immediately up-gradient were too shallow to intercept this interval, and thus the more impacted zone of the contaminant plume. Intermittent occurrence of faecal indicator bacteria more than 100 m down gradient of the percolation area suggested relatively short travel times. Study findings highlight the utility of geophysics as part of multidisciplinary investigations for OSWTS contaminant plume characterisation, while also demonstrating the capacity of effluent discharging to PPAs to impact groundwater quality at distance. Comparable geophysical responses observed in similar settings across Ireland suggest the phenomena observed in this study are more widespread than previously suspected.