The changing nature of groundwater in the global water cycle.

In recent decades, climate change and other anthropogenic activities have substantially affected groundwater systems worldwide. These impacts include changes in groundwater recharge, discharge, flow, storage, and distribution. Climate-induced shifts are evident in altered recharge rates, greater groundwater contribution to streamflow in glacierized catchments, and enhanced groundwater flow in permafrost areas. Direct anthropogenic changes include groundwater withdrawal and injection, regional flow regime modification, water table and storage alterations, and redistribution of embedded groundwater in foods globally. Notably, groundwater extraction contributes to sea level rise, increasing the risk of groundwater inundation in coastal areas. The role of groundwater in the global water cycle is becoming more dynamic and complex. Quantifying these changes is essential to ensure sustainable supply of fresh groundwater resources for people and ecosystems.

Development of complex patterns of anthropogenic uplift and subsidence in the Delaware Basin of West Texas and southeast New Mexico, USA.

The Delaware Basin in west Texas and southeast New Mexico is now the largest global oil producing basin, averaging ∼400,000 m3 (∼2,500,000 barrels) per day in 2022. The shale-dominated strata targeted for production can co-produce 4-5 times more water than oil, necessitating disposal by injection of ∼1,400,000 m3 (∼8,700,000 barrels) of water per day in 2022. Through a comprehensive assimilation of regional Sentinel-1 satellite radar data and analysis of production and injection, we show how petroleum operations have caused the development of complex and accelerating patterns of surface deformation from 2015 through 2021. We observe uplift from reservoir swelling, subsidence from reservoir contraction, and the development of linear features that are indicative of faulting. Subsidence is predominantly caused by production, and an important finding of this study is that the magnitude of subsidence is linearly proportional to total production. Uplift is caused by pressurization from wastewater injection of shallow permeable strata. The patterns of uplift are complex and extend laterally well-beyond areas where injection was performed. Linear surface deformation features are observed throughout the Delaware Basin, and they are lengthening and densifying as uplift and subsidence accelerate. Many of the lineations can be linked to known strata-bounded faults and shallow seismicity in the southern Delaware Basin where they serve as permeable and anisotropic conduits for pore pressure migration. In the Northern Delaware Basin, co-seismic rupture is hosted along basement-rooted faults that may link to the linear surface features. Understanding these dynamic changes in Delaware Basin is a pressing concern for management of subsurface reservoirs and safeguarding the surface environment. Concerns include ongoing induced seismicity, hazard of drilling through over-pressured strata, maintenance of integrity for newer wellbores, mitigation of flows of brine and petroleum fluids at the surface of old wellbores, and management of the pore space resource for wastewater injection.

Mobilization of Arsenic and Other Naturally Occurring Contaminants during Managed Aquifer Recharge: A Critical Review.

Population growth and climate variability highlight the need to enhance freshwater security and diversify water supplies. Subsurface storage of water in depleted aquifers is increasingly used globally to alleviate disparities in water supply and demand often caused by climate extremes including floods and droughts. Managed aquifer recharge (MAR) stores excess water supplies during wet periods via infiltration into shallow underlying aquifers or direct injection into deep aquifers for recovery during dry seasons. Additionally, MAR can be designed to improve recharge water quality, particularly in the case of soil aquifer treatment and riverbank filtration. While there are many potential benefits to MAR, introduction of recharge water can alter the native geochemical and hydrological conditions in the receiving aquifer, potentially mobilizing toxic, naturally occurring (geogenic) contaminants from sediments into groundwater where they pose a much larger threat to human and ecosystem health. On the basis of the present literature, arsenic poses the most widespread challenge at MAR sites due to its ubiquity in subsurface sediments and toxicity at trace concentrations. Other geogenic contaminants of concern include fluoride, molybdenum, manganese, and iron. Water quality degradation threatens the viability of some MAR projects with several sites abandoning operations due to arsenic or other contaminant mobilization. Here, we provide a critical review of studies that have uncovered the geochemical and hydrological mechanisms controlling mobilization of arsenic and other geogenic contaminants at MAR sites worldwide, including both infiltration and injection sites. These mechanisms were evaluated based on site-specific characteristics, including hydrological setting, native aquifer geochemistry, and operational site parameters (e.g., source of recharge water and recharge/recovery cycling). Observed mechanisms of geogenic contaminant mobilization during MAR via injection include shifting redox conditions and, to a lesser extent, pH-promoted desorption, mineral solubility, and competitive ligand exchange. The relative importance of these mechanisms depends on various site-specific, operational parameters, including pretreatment of injection water and duration of injection, storage, and recovery phases. This critical review synthesizes findings across case studies in various geochemical, hydrological, and operational settings to better understand controls on arsenic and other geogenic contaminant mobilization and inform the planning and design of future MAR projects to protect groundwater quality. This critical review concludes with an evaluation of proposed management strategies for geogenic contaminants and identification of knowledge gaps regarding fate and transport of geogenic contaminants during MAR.

A screening approach to improve water management practices in undeveloped shale plays, with application to the transboundary Eagle Ford Formation in northeast Mexico.

Hydraulic fracturing (HF) operations have transformed the unconventional energy industry, leading to a global increase in hydrocarbon production. Despite this, only the US, China, Canada and Argentina currently dominate production of unconventional resources, with the majority of shale basins globally remaining unprofitable to develop. An important gap in current water-energy nexus research, which this study addresses, is the assessment of potential water use to satisfy HF procedures in emergent plays. This work presents a screening tool for assessing first-order estimates of water impacts in undeveloped shale plays, testing the approach in the transboundary Eagle Ford (EF) play in northeast Mexico. We couple surface water and groundwater stress indicators derived from global hydrological variables to depict a baseline water stress index. Relative water stress is mapped for proposed blocks to be leased by the Mexican government in the future. We simulate four HF scenarios to assess new total water stress indicators for each block, considering shale production schemes using representative well drilling density (well lateral length(s) per unit area) and HF water intensity (HF water volume per unit lateral length) from existing EF development in Texas. Results suggest that the most feasible management scenario would consider the drilling of ∼1360 new unconventional wells/yr with projected HF water use of ∼57 Mm3/yr (0.7% of the total water withdrawals). The remaining scenarios will largely affect groundwater resources. Though applied to the EF in Mexico, this screening tool can assess water use constraints in emerging unconventional plays globally.

Water Stress from High-Volume Hydraulic Fracturing Potentially Threatens Aquatic Biodiversity and Ecosystem Services in Arkansas, United States.

Demand for high-volume, short duration water withdrawals could create water stress to aquatic organisms in Fayetteville Shale streams sourced for hydraulic fracturing fluids. We estimated potential water stress using permitted water withdrawal volumes and actual water withdrawals compared to monthly median, low, and high streamflows. Risk for biological stress was considered at 20% of long-term median and 10% of high- and low-flow thresholds. Future well build-out projections estimated potential for continued stress. Most water was permitted from small, free-flowing streams and "frack" ponds (dammed streams). Permitted 12-h pumping volumes exceeded median streamflow at 50% of withdrawal sites in June, when flows were low. Daily water usage, from operator disclosures, compared to median streamflow showed possible water stress in 7-51% of catchments from June-November, respectively. If 100% of produced water was recycled, per-well water use declined by 25%, reducing threshold exceedance by 10%. Future water stress was predicted to occur in fewer catchments important for drinking water and species of conservation concern due to the decline in new well installations and increased use of recycled water. Accessible and precise withdrawal and streamflow data are critical moving forward to assess and mitigate water stress in streams that experience high-volume withdrawals.

Trace Element Behavior in Methane-Rich and Methane-Free Groundwater in North and East Texas.

There is concern about adverse impacts of natural gas (primarily methane) production on groundwater quality; however, data on trace element concentrations are limited. The objective of this study was to compare the distribution of trace elements in groundwater samples with and without dissolved methane in aquifers overlying the Barnett Shale (Hood and Parker counties, 207 samples) and the Haynesville Shale (Panola County, 42 samples). Both shales have been subjected to intensive hydraulic fracturing for gas production. Well clusters with high dissolved methane were previously found in these counties and are thought to be of natural origin. Overall, groundwater in these counties is of excellent quality with typically low elemental concentrations. Several statistical analyses strongly suggest that most trace element concentrations, generally at low background levels, are no higher and even reduced when dissolved methane is present. In addition, trace element concentrations are not correlated with distance to gas wells. The reduction in trace element concentrations is attributed to anaerobic microbial degradation of methane, is associated with a higher pH (>8.5), and, likely, with precipitation of carbonates and pyrite and formation of clays. Trace and other elements are likely incorporated within the precipitating mineral crystalline network or sorbed. High pH values are found throughout these high-methane clusters (e.g., Parker-Hood cluster), even in subregions where methane is not present, which is consistent with a pervasive natural origin of dissolved methane rather than a limited gas well source.

Characterizing the Noble Gas Isotopic Composition of the Barnett Shale and Strawn Group and Constraining the Source of Stray Gas in the Trinity Aquifer, North-Central Texas.

This study presents the complete set of stable noble gases for Barnett Shale and Strawn Group production gas together with stray flowing gas in the Trinity Aquifer, Texas. It places new constraints on the source of this stray gas and further shows that Barnett and Strawn gas have distinct crustal and atmospheric noble gas signatures, allowing clear identification of these two sources. Like stray gas, Strawn gas is significantly more enriched in crustal 4He*, 21Ne*, and 40Ar* than Barnett gas. The similarity of Strawn and stray gas crustal noble gas signatures suggests that the Strawn is the source of stray gas in the Trinity Aquifer. Atmospheric 22Ne/36Ar ratios of stray gas mimic also that of Strawn, further reinforcing the notion that the source of stray gas in this aquifer is the Strawn. While noble gas signatures of Strawn and stray gas are consistent with a single-stage water degassing model, a two-stage oil modified groundwater exsolution fractionation model is required to explain the light atmospheric noble gas signature of Barnett Shale production gas. These distinct Strawn and Barnett noble gas signatures are likely the reflection of distinct evolution histories with Strawn gas being possibly older than that of Barnett Shale.

Methane Occurrences in Aquifers Overlying the Barnett Shale Play with a Focus on Parker County, Texas.

Clusters of elevated methane concentrations in aquifers overlying the Barnett Shale play have been the focus of recent national attention as they relate to impacts of hydraulic fracturing. The objective of this study was to assess the spatial extent of high dissolved methane previously observed on the western edge of the play (Parker County) and to evaluate its most likely source. A total of 509 well water samples from 12 counties (14,500 km2 ) were analyzed for methane, major ions, and carbon isotopes. Most samples were collected from the regional Trinity Aquifer and show only low levels of dissolved methane (85% of 457 unique locations <0.1 mg/L). Methane, when present is primarily thermogenic (δ13 C 10th and 90th percentiles of -57.54 and -39.00‰ and C1/C2+C3 ratio 10th, 50th, and 90th percentiles of 5, 15, and 42). High methane concentrations (>20 mg/L) are limited to a few spatial clusters. The Parker County cluster area includes historical vertical oil and gas wells producing from relatively shallow formations and recent horizontal wells producing from the Barnett Shale (depth of ∼1500 m). Lack of correlation with distance to Barnett Shale horizontal wells, with distance to conventional wells, and with well density suggests a natural origin of the dissolved methane. Known commercial very shallow gas accumulations (<200 m in places) and historical instances of water wells reaching gas pockets point to the underlying Strawn Group of Paleozoic age as the main natural source of the dissolved gas.

Controls on Methane Occurrences in Aquifers Overlying the Eagle Ford Shale Play, South Texas.

Assessing natural vs. anthropogenic sources of methane in drinking water aquifers is a critical issue in areas of shale oil and gas production. The objective of this study was to determine controls on methane occurrences in aquifers in the Eagle Ford Shale play footprint. A total of 110 water wells were tested for dissolved light alkanes, isotopes of methane, and major ions, mostly in the eastern section of the play. Multiple aquifers were sampled with approximately 47 samples from the Carrizo-Wilcox Aquifer (250-1200 m depth range) and Queen City-Sparta Aquifer (150-900 m depth range) and 63 samples from other shallow aquifers but mostly from the Catahoula Formation (depth <150 m). Besides three shallow wells with unambiguously microbial methane, only deeper wells show significant dissolved methane (22 samples >1 mg/L, 10 samples >10 mg/L). No dissolved methane samples exhibit thermogenic characteristics that would link them unequivocally to oil and gas sourced from the Eagle Ford Shale. In particular, the well water samples contain very little or no ethane and propane (C1/C2+C3 molar ratio >453), unlike what would be expected in an oil province, but they also display relatively heavier δ13 Cmethane (>-55‰) and δDmethane (>-180‰). Samples from the deeper Carrizo and Queen City aquifers are consistent with microbial methane sourced from syndepositional organic matter mixed with thermogenic methane input, most likely originating from deeper oil reservoirs and migrating through fault zones. Active oxidation of methane pushes δ13 Cmethane and δDmethane toward heavier values, whereas the thermogenic gas component is enriched with methane owing to a long migration path resulting in a higher C1/C2+C3 ratio than in the local reservoirs.

Unconventional Oil and Gas Spills: Risks, Mitigation Priorities, and State Reporting Requirements.

Rapid growth in unconventional oil and gas (UOG) has produced jobs, revenue, and energy, but also concerns over spills and environmental risks. We assessed spill data from 2005 to 2014 at 31 481 UOG wells in Colorado, New Mexico, North Dakota, and Pennsylvania. We found 2-16% of wells reported a spill each year. Median spill volumes ranged from 0.5 m3 in Pennsylvania to 4.9 m3 in New Mexico; the largest spills exceeded 100 m3. Seventy-five to 94% of spills occurred within the first three years of well life when wells were drilled, completed, and had their largest production volumes. Across all four states, 50% of spills were related to storage and moving fluids via flowlines. Reporting rates varied by state, affecting spill rates and requiring extensive time and effort getting data into a usable format. Enhanced and standardized regulatory requirements for reporting spills could improve the accuracy and speed of analyses to identify and prevent spill risks and mitigate potential environmental damage. Transparency for data sharing and analysis will be increasingly important as UOG development expands. We designed an interactive spills data visualization tool ( http://snappartnership.net/groups/hydraulic-fracturing/webapp/spills.html ) to illustrate the value of having standardized, public data.

Controls on Methane Occurrences in Shallow Aquifers Overlying the Haynesville Shale Gas Field, East Texas.

Understanding the source of dissolved methane in drinking-water aquifers is critical for assessing potential contributions from hydraulic fracturing in shale plays. Shallow groundwater in the Texas portion of the Haynesville Shale area (13,000 km2 ) was sampled (70 samples) for methane and other dissolved light alkanes. Most samples were derived from the fresh water bearing Wilcox formations and show little methane except in a localized cluster of 12 water wells (17% of total) in a approximately 30 × 30 km2 area in Southern Panola County with dissolved methane concentrations less than 10 mg/L. This zone of elevated methane is spatially associated with the termination of an active fault system affecting the entire sedimentary section, including the Haynesville Shale at a depth more than 3.5 km, and with shallow lignite seams of Lower Wilcox age at a depth of 100 to 230 m. The lignite spatial extension overlaps with the cluster. Gas wetness and methane isotope compositions suggest a mixed microbial and thermogenic origin with contribution from lignite beds and from deep thermogenic reservoirs that produce condensate in most of the cluster area. The pathway for methane from the lignite and deeper reservoirs is then provided by the fault system.

Unconventional oil and gas spills: Materials, volumes, and risks to surface waters in four states of the U.S.

Extraction of oil and gas from unconventional sources, such as shale, has dramatically increased over the past ten years, raising the potential for spills or releases of chemicals, waste materials, and oil and gas. We analyzed spill data associated with unconventional wells from Colorado, New Mexico, North Dakota and Pennsylvania from 2005 to 2014, where we defined unconventional wells as horizontally drilled into an unconventional formation. We identified materials spilled by state and for each material we summarized frequency, volumes and spill rates. We evaluated the environmental risk of spills by calculating distance to the nearest stream and compared these distances to existing setback regulations. Finally, we summarized relative importance to drinking water in watersheds where spills occurred. Across all four states, we identified 21,300 unconventional wells and 6622 reported spills. The number of horizontal well bores increased sharply beginning in the late 2000s; spill rates also increased for all states except PA where the rate initially increased, reached a maximum in 2009 and then decreased. Wastewater, crude oil, drilling waste, and hydraulic fracturing fluid were the materials most often spilled; spilled volumes of these materials largely ranged from 100 to 10,000L. Across all states, the average distance of spills to a stream was highest in New Mexico (1379m), followed by Colorado (747m), North Dakota (598m) and then Pennsylvania (268m), and 7.0, 13.3, and 20.4% of spills occurred within existing surface water setback regulations of 30.5, 61.0, and 91.4m, respectively. Pennsylvania spills occurred in watersheds with a higher relative importance to drinking water than the other three states. Results from this study can inform risk assessments by providing improved input parameters on volume and rates of materials spilled, and guide regulations and the management policy of spills.

Methane Sources and Migration Mechanisms in Shallow Groundwaters in Parker and Hood Counties, Texas-A Heavy Noble Gas Analysis.

This study places constraints on the source and transport mechanisms of methane found in groundwater within the Barnett Shale footprint in Texas using dissolved noble gases, with particular emphasis on 84Kr and 132Xe. Dissolved methane concentrations are positively correlated with crustal 4He, 21Ne, and 40Ar and suggest that noble gases and methane originate from common sedimentary strata, likely the Strawn Group. In contrast to most samples, four water wells with the highest dissolved methane concentrations unequivocally show strong depletion of all atmospheric noble gases (20Ne, 36Ar, 84Kr, 132Xe) with respect to air-saturated water (ASW). This is consistent with predicted noble gas concentrations in a water phase in contact with a gas phase with initial ASW composition at 18 °C-25 °C and it suggests an in situ, highly localized gas source. All of these four water wells tap into the Strawn Group and it is likely that small gas accumulations known to be present in the shallow subsurface were reached. Additionally, lack of correlation of 84Kr/36Ar and 132Xe/36Ar fractionation levels along with 4He/20Ne with distance to the nearest gas production wells does not support the notion that methane present in these groundwaters migrated from nearby production wells either conventional or using hydraulic fracturing techniques.

Source and fate of hydraulic fracturing water in the Barnett Shale: a historical perspective.

Considerable controversy continues about water availability for and potential impacts of hydraulic fracturing (HF) of hydrocarbon assets on water resources. Our objective was to quantify HF water volume in terms of source, reuse, and disposal, using the Barnett Shale in Texas as a case study. Data were obtained from commercial and state databases, river authorities, groundwater conservation districts, and operators. Cumulative water use from ∼ 18,000 (mostly horizontal) wells since 1981 through 2012 totaled ∼ 170,000 AF (210 Mm(3)); ∼ 26 000 AF (32 Mm(3)) in 2011, representing 32% of Texas HF water use and ∼ 0.2% of 2011 state water consumption. Increase in water use per well by 60% (from 3 to 5 Mgal/well; 0.011-0.019 Mm(3)) since the mid-2000s reflects the near-doubling of horizontal-well lengths (2000-3800 ft), offset by a reduction in water-use intensity by 40% (2000-1200 gal/ft; 2.5-1.5 m(3)/m). Water sources include fresh surface water and groundwater in approximately equal amounts. Produced water amount is inversely related to gas production, exceeds HF water volume, and is mostly disposed in injection wells. Understanding the historical evolution of water use in the longest-producing shale play is invaluable for assessing its water footprint for energy production.

Water use for Shale-gas production in Texas, U.S.

Shale-gas production using hydraulic fracturing of mostly horizontal wells has led to considerable controversy over water-resource and environmental impacts. The study objective was to quantify net water use for shale-gas production using data from Texas, which is the dominant producer of shale gas in the U.S. with a focus on three major plays: the Barnett Shale (~15,000 wells, mid-2011), Texas-Haynesville Shale (390 wells), and Eagle Ford Shale (1040 wells). Past water use was estimated from well-completion data, and future water use was extrapolated from past water use constrained by shale-gas resources. Cumulative water use in the Barnett totaled 145 Mm(3) (2000-mid-2011). Annual water use represents ~9% of water use in Dallas (population 1.3 million). Water use in younger (2008-mid-2011) plays, although less (6.5 Mm(3) Texas-Haynesville, 18 Mm(3) Eagle Ford), is increasing rapidly. Water use for shale gas is <1% of statewide water withdrawals; however, local impacts vary with water availability and competing demands. Projections of cumulative net water use during the next 50 years in all shale plays total ~4350 Mm(3), peaking at 145 Mm(3) in the mid-2020s and decreasing to 23 Mm(3) in 2060. Current freshwater use may shift to brackish water to reduce competition with other users.

Unsaturated zone arsenic distribution and implications for groundwater contamination.

Arsenic compounds have been applied at the land surface as pesticides in agricultural areas globally. The purpose of this study was to evaluate the fate of anthropogenic arsenic applications related to agriculture, using arsenic applications on cotton in the southern High Plains (SHP), Texas, as a case study and examining possible linkages with contamination of the underlying Ogallala aquifer in this region, where 36% of wells exceed the new EPA 10 microg/L standard. Unsaturated zone soil samples were collected from boreholes beneath natural ecosystems (grassland/ shrubland) to provide a control (no arsenic application) (5 profiles) and cotton cropland (20 profiles) for analyses of water-extractable arsenic, vanadium, phosphate, chloride, and nitrate. Natural ecosystem profiles have high arsenic concentrations at depth (maximum of 7.2-69.6 microg As/ kg dry soil at 5.9-21.4 m depth) that are attributed to a geologic source. Most profiles beneath cotton cropland have high arsenic concentrations within the upper meter (profile means 1.7 to 31.6 microg/kg) that correlate with phosphate (r = 0.70, p < 0.01) and are attributed to anthropogenic arsenic application associated with phosphate fertilizer application. High arsenic concentrations at >1 m depth (profile means < or =36.3 microg/kg) found in cropland profiles are attributed to a geologic source because of similarity with profiles beneath natural ecosystems, lack of correlation with phosphate, and pore-water ages that predate anthropogenic arsenic application in many profiles. GIS analyses showed poor correlations between groundwater arsenic and percent cultivated land (r = -0.15, p < 0.01), groundwater nitrate (r = 0.30, p < 0.01), and water table depth (r= -0.31, p < 0.01), further supporting the idea that anthropogenic-derived arsenic in the shallow subsurface is not linked to groundwater arsenic contamination in this region.

Scoping analyses of geochemical sealing of early cracks in a waste container and associated drip shield, Yucca Mountain, Nevada.

Early after final emplacement of the nuclear waste containers at the proposed Yucca Mountain, Nevada, high-level-waste repository, tiny cracks (less than 200 microm wide, 1 to 2 cm deep, and a few centimeters long at most) could appear in the containers and in the drip shield protecting them. Modeling calculations were performed to understand how fast those cracks could be sealed. Under dripping conditions, they are expected to be bridged with water. If cracks are located in the drip shield, any further dripping on the waste containers located underneath will be limited. If cracks are located in a container, potentially harmful radionuclides could only travel by diffusion. In addition, water-bridged cracks will be sealed through at least two processes: precipitation of calcite with minor silica following evaporative concentration of the water residing in the cracks and continuous corrosion of the crack walls. The sealing rate is calculated as the intersection of the time of emergence of the cracks, the water dripping rate, and the decreasing evaporation rate. The evaporative driving force declines as short-lived radioactive elements, having given up much of the heat affecting the repository, are progressively depleted from the waste. Depending on the crack initiation time and environmental conditions, crack sealing varies from a few tens of years to a few thousand years. Because environmental conditions in the vicinity of the cracks and at the crack scale have not been produced, a parametric method scaling drift scale conditions is used.