Susan Brantley

While unconventional oil and gas development (UOGD) is changing the world economy, processes that are used during UOGD such as high-volume hydraulic fracturing (“fracking”) have been linked with water contamination. Water quality risks include leaks of gas and salty fluids (brines) that are co-produced at wellpads. Identifying the cause of contamination is difficult, however, because UOG wells are often co-located with other contaminant sources. We investigated the world’s largest shale gas play with publicly accessible groundwater data (~29,000 analyses from the Marcellus Shale in Pennsylvania, U.S.A.) and discovered that concentrations of brine-associated species barium ([Ba]) and strontium ([Sr]) show small regional increases within 1km of UOGD. Higher concentrations in groundwaters are associated with greater proximity to and density of UOG wells. Concentration increases are even larger when considering the locations of i) spill-related violations and ii) some wastewater impoundments. These statistically significant relationships persist even after correcting for other natural and anthropogenic sources of salts. The most likely explanation is that UOGD slightly increases salt concentrations in regional groundwaters not because of fracking but because of the ubiquity of wastewater management issues. The high frequency of spills and leaks across shale gas basins suggests other plays could show similar effects.

Earth’s climate may be stabilized over millennia by solubilization of atmospheric carbon dioxide (CO2) as minerals weather, but the temperature sensitivity of this thermostat is poorly understood. We discovered that the temperature dependence of weathering expressed as an activation energy increases from laboratory to watershed as transport, clay precipitation, disaggregation, and fracturing increasingly couple to dissolution. A simple upscaling to the global system indicates that the temperature dependence decreases to ~22 kilojoules per mole because (i) the lack of runoff limits weathering and retains base metal cations on half the land surface and (ii) other landscapes are regolith-shielded and show little weathering response to temperature. By comparing weathering from laboratory to globe, we reconcile some aspects of kinetic and thermodynamic controls on CO2 drawdown by natural or enhanced weathering.

Although metal redox reactions in soils can strongly affect carbon mineralization and other important soil processes, little is known about temporal variations in this redox cycling. Recently, potentiostatically poised electrodes (fixed-potential electrodes) have shown promise for measuring the rate of oxidation and reduction at a specific reduction potential in situ in riparian soils. Here for the first time, we used these electrodes in unsaturated soils to explore the fine-scale temporal redox fluctuations of both iron and manganese in response to environmental conditions. We used three-electrode systems with working electrodes fixed at 100 mV (vs. SHE) and 400 mV at 50 cm and 70 cm in the valley floor soil of a headwater watershed. Electrodes fixed at 100 mV to mimic iron oxides and at 400 mV to mimic manganese oxides allowed real-time reduction and oxidation rates to be calculated from temporal variations in the …

As drinking‐water scarcity grows worldwide, we need to improve predictions of the quantity and quality of our water resources. An overarching problem for model improvement is that we do not know the geological structure of aquifers in sufficient detail. In this work, we demonstrate that mineral‐water reactions imprint structure in the subsurface that impacts the flow and transport of some chemical species. Specifically, pyrite, a ubiquitous mineral, commonly oxidizes and depletes in the upper layers of the weathering profile in most humid watersheds, only remaining at depths of meters. We hypothesize that variations in concentrations (C) of pyrite‐derived sulfate released into rivers as a function of discharge (q) reflect the rate‐limiting step and depth of pyrite‐oxidizing layers. We found that log C − log q behaviors thus differ in small and large watersheds in the Susquehanna River Basin as well as in selected …

The oxidation of petrogenic organic carbon (OCpetro) is a source of carbon dioxide to the atmosphere over geological timescales. The rates of OCpetro oxidation in locations that experience low rates of denudation remain poorly constrained, despite these landscapes dominating Earth's continental surface area. Here, we track OCpetro oxidation using radiocarbon and the trace element rhenium (Re) in the deep weathering profiles, soils and stream waters of the Susquehanna Shale Hills Critical Zone Observatory (PA, USA). In a ridge‐top borehole, radiocarbon measurements reveal the presence of a broad OCpetro weathering front, with a first‐order assessment of ~40% loss occurring over ~6 m. However, the low OCpetro concentration (< 0.05 wt%) and inputs of radiocarbon throughout the deepest parts of the profile complicate the assessment of OCpetro loss. The OCpetro weathering front coincides with a …

Despite its importance, only a few researchers have incorporated the effects of fracturing into models of reactive transport for rock weathering. Here we explore 2D simulations that describe weathering under conditions of diffusive and advective transport within heterogeneous media consisting of rocky blocks and fractures. In our simulations, the Darcy velocities vary in space and time and depend on weathering processes within the rock matrix. We explore simulations with saturated and unsaturated flow for weathering bedrock that consists of blocks separated by inert or weathered material. The simulations show that a simplified homogenized model can approximate exact solutions for some of the simulated columns and hills and can allow exploration of coupling between flow and reaction in fractured rock. These hillslope simulations document that, even in the presence of 2D water flow, i) an increase in fracture density results in faster weathering advance rates; and ii) the water table locates deeper for a rock system that is weathered and fractured rather than weathered and unfractured. Some of these patterns have also been observed for natural systems. But these simulations also highlight how simplified models that do not use appropriate averaging of heterogeneities can be inaccurate in predicting weathering rate for natural systems. For example, if water flows both vertically and laterally through the vadose zone of a hill, then a prediction of the depth of regolith that is based on modeling strictly unidirectional downward infiltration will be unrealistically large. Likewise, if the fracture density observed near the land surface is used in a model to …

To understand geological systems that are important to our national well-being -- such as nuclear waste repositories or hydraulically fractured rocks -- requires the ability to make quantitative projections of the evolution of rock-water systems forward in time. A fundamental goal of much of our geochemical research was focused on development of numerical models of reactive transport to simulate the evolution of natural water-rock systems. The geochemical community has provided such models and they are rapidly being improved with field-model cross-testing. This project was especially impactful in field-model testing for well-constrained weathering systems to provide new conceptual understanding of how these systems function.

Understanding emissions of methane from legacy and ongoing shale gas development requires both regional studies that assess the frequency of emissions and case studies that assess causation. We present the first direct measurements of emissions in a case study of a putatively leaking gas well in the largest shale gas play in the United States. We quantify atmospheric methane emissions in farmland >2 km from the nearest shale gas well cited for casing and cementing issues. We find that emissions are highly heterogeneous as they travel long distances in the subsurface. Emissions were measured near observed patches of dead vegetation and methane bubbling from a stream. An eddy covariance flux tower, chamber flux measurements, and a survey of enhancements of the near-surface methane mole fraction were used to quantify emissions and evaluate the spatial and temporal variability. We combined …