Groundwater flow velocities in a fractured carbonate aquifer-type: Implications for contaminant transport.

Contaminants that are highly soluble in groundwater are rapidly transported via fractures in mechanically resistant sedimentary rock aquifers. Hence, a rigorous methodology is needed to estimate groundwater flow velocities in such fractured aquifers. Here, we propose an approach using borehole hydraulic testing to compute flow velocities in an un-faulted area of a fractured carbonate aquifer by applying the cubic law to a parallel plate model. The Cadeby Formation (Yorkshire, NE England) - a Permian dolostone aquifer present beneath the University of Leeds Farm - is the fractured aquifer selected for this hydraulic experiment. The bedding plane fractures of this dolostone aquifer, which are sub-horizontal, sub-parallel and laterally persistent, largely dominate the flow at shallow (<~40 mBGL) depths. These flowing bedding plane discontinuities are separated by a rock matrix which is relatively impermeable (Kwell-test/Kcore-plug~104) as is common in fractured carbonate aquifers. In the workflow reported here, the number of flowing fractures - mainly bedding plane fractures - intersecting three open monitoring wells are found from temperature/fluid conductivity and acoustic/optical televiewer logging. Following well installation, average fracture hydraulic apertures for screened intervals are found from analysis of slug tests. For the case study aquifer, this workflow predicts hydraulic apertures ranging from 0.10 up to 0.54 mm. However, groundwater flow velocities range within two order of magnitude from 13 up to 242 m/day. Notably, fracture apertures and flow velocities rapidly reduce with increasing depth below the water table; the upper ~10 m shows relatively high values of hydraulic conductivity (0.30-2.85 m/day) and corresponding flow velocity (33-242 m/day). Permeability development around the water table in carbonate aquifer-types is common, and arises where high pCO2 recharge water from the soil zone causes calcite/dolomite dissolution. Hence, agricultural contaminants entering the aquifer with recharge water are laterally transported rapidly within this upper part. Computation of groundwater flow velocities allows determination of the Reynolds number. Values of up ~1, indicating the lower limit of the transition from laminar to turbulent flow, are found at the studied site, which is situated away from major fault traces. Hence, turbulent flow is likely to arise in proximity to tectonic structures, such as normal faults, which localize flow and enhance karstification. The occurrence of turbulent flow in correspondence of such tectonic structures should be represented in regional groundwater flow simulations.

Oxalate secretion by ectomycorrhizal Paxillus involutus is mineral-specific and controls calcium weathering from minerals.

Trees and their associated rhizosphere organisms play a major role in mineral weathering driving calcium fluxes from the continents to the oceans that ultimately control long-term atmospheric CO2 and climate through the geochemical carbon cycle. Photosynthate allocation to tree roots and their mycorrhizal fungi is hypothesized to fuel the active secretion of protons and organic chelators that enhance calcium dissolution at fungal-mineral interfaces. This was tested using (14)CO2 supplied to shoots of Pinus sylvestris ectomycorrhizal with the widespread fungus Paxillus involutus in monoxenic microcosms, revealing preferential allocation by the fungus of plant photoassimilate to weather grains of limestone and silicates each with a combined calcium and magnesium content of over 10 wt.%. Hyphae had acidic surfaces and linear accumulation of weathered calcium with secreted oxalate, increasing significantly in sequence: quartz, granite < basalt, olivine, limestone < gabbro. These findings confirmed the role of mineral-specific oxalate exudation in ectomycorrhizal weathering to dissolve calcium bearing minerals, thus contributing to the geochemical carbon cycle.

Plant-driven weathering of apatite--the role of an ectomycorrhizal fungus.

Ectomycorrhizal (EcM) fungi are increasingly recognized as important agents of mineral weathering and soil development, with far-reaching impacts on biogeochemical cycles. Because EcM fungi live in a symbiotic relationship with trees and in close contact with bacteria and archaea, it is difficult to distinguish between the weathering effects of the fungus, host tree and other micro-organisms. Here, we quantified mineral weathering by the fungus Paxillus involutus, growing in symbiosis with Pinus sylvestris under sterile conditions. The mycorrhizal trees were grown in specially designed sterile microcosms in which the supply of soluble phosphorus (P) in the bulk media was varied and grains of the calcium phosphate mineral apatite mixed with quartz, or quartz alone, were provided in plastic wells that were only accessed by their fungal partner. Under P limitation, pulse labelling of plants with (14)CO(2) revealed plant-to-fungus allocation of photosynthates, with 17 times more (14)C transferred into the apatite wells compared with wells with only quartz. Fungal colonization increased the release of P from apatite by almost a factor of three, from 7.5 (±1.1) × 10(-10) mol m(-2) s(-1) to 2.2 (±0.52) × 10(-9) mol m(-2) s(-1). On increasing the P supply in the microcosms from no added P, through apatite alone, to both apatite and orthophosphate, the proportion of biomass in roots progressively increased at the expense of the fungus. These three observations, (i) proportionately more plant energy investment in the fungal partner under P limitation, (ii) preferential fungal transport of photosynthate-derived carbon towards patches of apatite grains and (iii) fungal enhancement of weathering rate, reveal the tightly coupled plant-fungal interactions underpinning enhanced EcM weathering of apatite and its utilization as P source.

Biological weathering and the long-term carbon cycle: integrating mycorrhizal evolution and function into the current paradigm.

The dramatic decline in atmospheric CO2 evidenced by proxy data during the Devonian (416.0-359.2 Ma) and the gradual decline from the Cretaceous (145.5-65.5 Ma) onwards have been linked to the spread of deeply rooted trees and the rise of angiosperms, respectively. But this paradigm overlooks the coevolution of roots with the major groups of symbiotic fungal partners that have dominated terrestrial ecosystems throughout Earth history. The colonization of land by plants was coincident with the rise of arbuscular mycorrhizal fungi (AMF),while the Cenozoic (c. 65.5-0 Ma) witnessed the rise of ectomycorrhizal fungi (EMF) that associate with both gymnosperm and angiosperm tree roots. Here, we critically review evidence for the influence of AMF and EMF on mineral weathering processes. We show that the key weathering processes underpinning the current paradigm and ascribed to plants are actually driven by the combined activities of roots and mycorrhizal fungi. Fuelled by substantial amounts of recent photosynthate transported from shoots to roots, these fungi form extensive mycelial networks which extend into soil actively foraging for nutrients by altering minerals through the acidification of the immediate root environment. EMF aggressively weather minerals through the additional mechanism of releasing low molecular weight organic chelators. Rates of biotic weathering might therefore be more usefully conceptualized as being fundamentally controlled by the biomass, surface area of contact, and capacity of roots and their mycorrhizal fungal partners to interact physically and chemically with minerals. All of these activities are ultimately controlled by rates of carbon-energy supply from photosynthetic organisms. The weathering functions in leading carbon cycle models require experiments and field studies of evolutionary grades of plants with appropriate mycorrhizal associations. Representation of the coevolution of roots and fungi in geochemical carbon cycle models is required to further our understanding of the role of the biota in Earth's CO2 and climate history.

Assessing the natural attenuation of organic contaminants in aquifers using plume-scale electron and carbon balances: model development with analysis of uncertainty and parameter sensitivity.

A quantitative methodology is described for the field-scale performance assessment of natural attenuation using plume-scale electron and carbon balances. This provides a practical framework for the calculation of global mass balances for contaminant plumes, using mass inputs from the plume source, background groundwater and plume residuals in a simplified box model. Biodegradation processes and reactions included in the analysis are identified from electron acceptors, electron donors and degradation products present in these inputs. Parameter values used in the model are obtained from data acquired during typical site investigation and groundwater monitoring studies for natural attenuation schemes. The approach is evaluated for a UK Permo-Triassic Sandstone aquifer contaminated with a plume of phenolic compounds. Uncertainty in the model predictions and sensitivity to parameter values was assessed by probabilistic modelling using Monte Carlo methods. Sensitivity analyses were compared for different input parameter probability distributions and a base case using fixed parameter values, using an identical conceptual model and data set. Results show that consumption of oxidants by biodegradation is approximately balanced by the production of CH4 and total dissolved inorganic carbon (TDIC) which is conserved in the plume. Under this condition, either the plume electron or carbon balance can be used to determine contaminant mass loss, which is equivalent to only 4% of the estimated source term. This corresponds to a first order, plume-averaged, half-life of > 800 years. The electron balance is particularly sensitive to uncertainty in the source term and dispersive inputs. Reliable historical information on contaminant spillages and detailed site investigation are necessary to accurately characterise the source term. The dispersive influx is sensitive to variability in the plume mixing zone width. Consumption of aqueous oxidants greatly exceeds that of mineral oxidants in the plume, but electron acceptor supply is insufficient to meet the electron donor demand and the plume will grow. The aquifer potential for degradation of these contaminants is limited by high contaminant concentrations and the supply of bioavailable electron acceptors. Natural attenuation will increase only after increased transport and dilution.

Processes controlling the distribution and natural attenuation of dissolved phenolic compounds in a deep sandstone aquifer.

Processes controlling the distribution and natural attenuation (NA) of phenol, cresols and xylenols released from a former coal-tar distillation plant in a deep Triassic sandstone aquifer are evaluated from vertical profiles along the plume centerline at 130 and 350 m from the site. Up to four groups of contaminants (phenols, mineral acids, NaOH, NaCl) form discrete and overlapping plumes in the aquifer. Their distribution reflects changing source history with releases of contaminants from different locations. Organic contaminant distribution in the aquifer is determined more by site source history than degradation. Contaminant degradation at total organic carbon (TOC) concentrations up to 6500 mg l(-1) (7500 mg l(-1) total phenolics) is occurring by aerobic respiration NO3-reduction, Mn(IV)-/Fe(III)-reduction, SO4-reduction, methanogenesis and fermentation, with the accumulation of inorganic carbon, organic metabolites (4-hydroxybenzaldehyde, 4-hydroxybenzoic acid), acetate, Mn(II), Fe(II), S(-II), CH4 and H2 in the plume. Aerobic and NO3-reducing processes are restricted to a 2-m-thick plume fringe but Mn(IV)-/Fe(II)-reduction, SO4-reduction, methanogenesis and fermentation occur concomitantly in the plume. Dissolved H2 concentrations in the plume vary from 0.7 to 110 nM and acetate concentrations reach 200 mg l(-1). The occurrence of a mixed redox system and concomitant terminal electron accepting processes (TEAPs) could be explained with a partial equilibrium model based on the potential in situ free energy (deltaGr) yield for oxidation of H2 by specific TEAPs. Respiratory processes rather than fermentation are rate limiting in determining the distribution of H2 and TEAPs and H2 dynamics in this system. Most (min. 90%) contaminant degradation has occurred by aerobic and NO3-reducing processes at the plume fringe. This potential is determined by the supply of aqueous O2 and NO3 from uncontaminated groundwater, as controlled by transverse mixing, which is limited in this aquifer by low dispersion. Consumption to date of mineral oxides and SO4 is, respectively, <0.15% and 0.4% of the available aquifer capacity, and degradation using these oxidants is <10%. Fermentation is a significant process in contaminant turnover, accounting for 21% of degradation products present in the plume, and indicating that microbial respiration rates are slow in comparison with fermentation. Under present conditions, the potential for degradation in the plume is very low due to inhibitory effects of the contaminant matrix. Degradation products correspond to <22% mass loss over the life of the plume, providing a first-order plume scale half-life >140 years. The phenolic compounds are biodegradable under the range of redox conditions in the aquifer and the aquifer is not oxidant limited, but the plume is likely to be long-lived and to expand. Degradation is likely to increase only after contaminant concentrations are reduced and aqueous oxidant inputs are increased by dispersion of the plume. The results imply that transport processes may exert a greater control on the natural attenuation of this plume than aquifer oxidant availability.