SEARCH · USGS Science

Results for “Journal of the Air & Waste Management Association”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,689 records · Page 45Linked to original sources

USGS2024

Global variability of the composition and temperature at the 410-km discontinuity from receiver function analysis of dense arrays

Seismic boundaries caused by phase transitions between olivine polymorphs in Earth's mantle provide thermal and compositional markers that inform mantle dynamics. Seismic studies of the mantle transition zone often use either global averaging with sparse arrays or regional sampling from a single dense array. The intermediate approach of this study utilizes many densely spaced seismic arrays distributed around the globe. We systematically compute teleseismic P-to-S receiver functions for each seismic array and invert for the 1-D seismic velocity structure of the mantle transition zone beneath each array to facilitate a comparison between densely sampled regions. We stack 3,600 receiver functions on average at 67 arrays in total. The stack is used in a probabilistic inversion to estimate the mantle transition zone interface depths and velocities beneath each array. We focus on the 410-km discontinuity (410) because it is a prominent seismic interface that is clearly linked to a single mineral phase transition between olivine and wadsleyite. The depths and velocity contrasts of the 410 are mapped to temperatures and compositions using mineral physics constraints. The depth of the 410 ranges from ∼405–440 km, which is consistent with a ∼360 K temperature range in a dry mantle and a ∼260 K temperature range in a wet mantle (2 wt. % water). The Vs contrast across the 410 ranges from ∼2.5–8 %, which is consistent with ∼20–70 vol. % olivine composition in a dry mantle and ∼25–80 vol. % in a wet mantle. The bulk composition of the upper mantle near the 410-km discontinuity is typically considered to be well-mixed because there is no thermodynamic impediment to convection at the olivine to wadsleyite phase transition. However, the wide range of inferred olivine content from our study suggests that there are large lateral variations in the bulk composition of the upper mantle near the 410-km discontinuity.

Earth and Planetary Science Letters↗

USGS2024

Bedrock geologic map of the Woodstock quadrangle, Grafton County, New Hampshire

The bedrock geology of the Woodstock 7.5-minute quadrangle consists of highly deformed metasedimentary rocks of the Central Maine trough, including the Silurian Rangeley and Perry Mountain Formations and the Devonian Littleton Formation. The central, northern, and eastern parts of the quadrangle are underlain by the oldest rocks in the area, the Rangeley Formation. In the southwest and south-central part of the quadrangle, metaturbidites of the Perry Mountain Formation and subsequent Littleton Formation overly the Rangeley Formation in a deformed F 1 synform, herein informally called the Bagley Brook basin. The metasedimentary rocks were intruded by widespread syn- to post-tectonic granitoids of the Devonian New Hampshire Plutonic Suite and minor post-metamorphic Jurassic-Cretaceous mafic dikes of the White Mountain Plutonic-Volcanic Suite. The metasedimentary rocks were affected by at least two episodes of deformation in the Devonian Acadian orogeny. The dominant regional foliation is second-generation (S 2 /D 2 ) and formed during the development of sillimanite-muscovite mineral assemblages. Large bodies of the Early Devonian Kinsman Granodiorite intruded the metasedimentary rocks semi-concordantly during D 2 deformation. Dikes of the Late Devonian Concord Granite cut the Kinsman Granodiorite and the metasedimentary rocks and were emplaced either syn- or post-D 2 . The map pattern in the Rangeley Formation is dominated by northeast to northwest trending, moderately to steeply north-dipping F 2 and F 3 folds. Map-scale F 1 folds are defined by the Bagley Brook basin. Previous division of Rangeley Formation stratigraphy in this region into “upper” and “lower” parts was not corroborated by 1:24,000-scale mapping of lithodemic units, and rocks previously mapped as part of the Smalls Falls and Madrid Formations are here reassigned to the Rangeley Formation. Some rocks previously mapped as the lower part of the Littleton Formation are now assigned to the Perry Mountain Formation. The Littleton Formation on this map is approximately equivalent to rocks previously mapped as the upper part of the same formation. Steeply dipping fractures in the quadrangle show a preferred northeast orientation, consistent with subsurface fracture orientations in the well fields near Mirror Lake. Jurassic-Cretaceous mafic dikes and normal faults show preferred northeast orientations, similar to the fractures, suggesting that the extensional stress field that controlled dike orientation during the Mesozoic also produced the dominant brittle fabrics in the area.

New Hampshire↗

USGS2025

The geology of Canadian potash: A critical mineral for feeding the world

Potash, potassium-bearing water-soluble salt, is the primary global economic source of potassium. Potash is recognized as a critical mineral in Canada as it is the largest source of potassium used in fertilizers. It is essential for global agricultural productivity and food security. Canada is the world’s largest potash exporter with vast deposits in the widely mined Prairie Evaporite of Saskatchewan, which formed in the epicontinental Elk Point Basin during the Middle Devonian. Potash is also found in the Windsor Group of Atlantic Canada where it formed in a series of tectonically active basins during the Mississippian that have undergone substantial post-depositional subsurface alteration and deformation. Potash deposits were mined in New Brunswick up until 2016. Both deposits are salt giants, recording times in the geologic record of extensive and long-lasting evaporite genesis under arid conditions in restricted seas. This paper reviews the geological and economic significance of Canadian potash, including (1) the genesis of each deposit, (2) diagenetic, erosional, and tectonic modification, and (3) exploration and mining in each basin. Underdeveloped regions, possible undiscovered resources, environmental considerations, and the importance of sustainable practices in light of climate change and socioeconomic risks are also addressed.

Facets↗

USGS2026

Earthquake catalog for the Fairbanks region of central Alaska, 2014–2024, based on waveform cross-correlation

The Fairbanks region of central Alaska is part of a broad zone of intraplate crustal deformation, situated north of the Denali fault and north of the ongoing collision and flat‐slab subduction of the Yakutat oceanic plateau. Seismicity in the Fairbanks region occurs both in diffuse areas as well as in well‐defined lineaments, such as the left‐lateral Salcha fault, which hosted the 1937 M 8 7.3 earthquake. Starting with the regional seismicity catalog, we perform waveform cross‐correlation, network‐matched filtering, and relative relocation to obtain an enhanced seismicity catalog over the time period 2014–2024. Based on the relocated catalog, we interpret a set of 15 fault segments, including two conjugate faults and two new faults east of the previously documented fault system. Considering the combined seismicity in the Minto and Fairbanks regions, the median depth of seismicity decreases from east (6 km) to west (20 km). Our interpreted faults provide guidance for future tectonic modeling and assessment of seismic hazards in this region.

Alaska↗

USGS2026

Constraining source and path effects of large magnitude earthquakes using ground motion simulations

The purpose of this study is to use ground‐motion simulations to investigate ways in which source and path effects for large‐magnitude earthquakes can be represented in nonergodic ground‐motion models (GMMs). To achieve this, we designed a ground‐motion study in the San Francisco Bay Area that includes earthquakes with a broad range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of kinematic simulations (magnitude 4–7), we then develop a nonergodic GMM with the simulated data. We find that trends in the within‐site residuals are affected significantly by the earthquake radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within‐site residuals among multiple source realizations. Finally, after removing the source effects from the within‐site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site as the large events, whereas the second approach considers all small events on the fault plane. The results indicate that it is difficult to satisfactorily approximate the path effects of large events with those of small events using either approach, at least in the case of simulations.

California↗

USGS2025

The effects of line simplification on planform geometry

Data on maps should retain accuracy regardless of scale. Yet, as cartographic lines are generalized, there can be impacts on properties such as topology, density, and planform geometry. Here, we investigate the use of the Scale Specific Sinuosity (S3) metric (Stanislawski et al., 2023) to evaluate the effects of line simplification on planform geometry, which is the bends of streams in map view. We employ an open-source Python S3 workflow to characterize the geometry of five diverse stream channels in the United States. The original data are extracted from the U.S. Geological Survey National Hydrography Dataset 1:24,000-scale vector data (U.S. Geological Survey, 2000) (Table 1), and the simplification is done using the Visvalingam and Whyatt method (2017) with a simplification tolerance of 0.5, 1.0, 1.5, and 2.0 km. The S3 analysis is calculated at each level of simplification and S3 derivatives are generated. Derivatives include measures of sinuosity, fractal dimension, and the dominant bend wavelength. The findings show that the change in planform geometry is scale-dependent, though simplification will have little effect on straighter lines. The change becomes more apparent in complex lines as the degree of simplification aligns with the scale of the dominant bend geometries. These logical conclusions are evidence that the S3 is a useful metric for automated characterization of bend geometry regardless of line complexity.

Conference Paper↗

USGS2024

Seismic response comparison of a historical masonry church subject to real and simulated ground motions

In recent years, advanced numerical models and high-performance computing have facilitated the utilization of ground motion time series in the assessment of the non-linear dynamic behavior of historic masonry structures. Since recorded accelerograms can be sparse for specific analysis conditions, stochastic ground motion simulations have become a viable alternative to overcome this limitation. This study simulates the recorded acceleration time series of the Central Italy 2016 earthquake event at the closest station to the town of Macerata using a site-based stochastic approach. The simulated motions are seismologically evaluated using a goodness-of-fit method in terms of various intensity measures. The simulated records, in conjunction with real records, are used to study the non-linear dynamic behavior of San Filippo Neri church located in Macerata. The church of San Filippo represents an important example of Baroque religious architecture in central Italy, which was damaged and closed off to the public after the 2016 earthquake events. The construction was investigated with a vast diagnostic campaign which included on-site testing and dynamic identification tests. The collected data is used to calibrate the dynamic response of a three-dimensional finite element model of the church. The model is finally used to compare the non-linear seismic responses under real and simulated ground motions with the site recorded damage. The results of structural responses demonstrate a strong agreement between the real and simulated records, providing evidence to support the validation of the site-based stochastic simulation.

Macerata↗

USGS2025

Apatite geo-thermochronology and geochemistry constrain Oligocene-Miocene growth and geodynamics of the northeastern Tibetan Plateau

Understanding the geodynamics of plateau evolution requires examining the spatial and temporal aspects of mountain building in the northeastern Tibetan Plateau, which are still under debate. Here we integrate apatite geo-thermochronological and geochemical data from the Oligocene-Miocene succession of the Xunhua Basin to elucidate the evolution of the regional topography. The results suggest sediment provenance changes, at ca. 28, 20 and 12 Ma, indicating topographic growth of the West Qinling, Laji Shan, and Jishi Shan source areas, respectively. Our findings reveal Oligocene-Miocene stepwise deformation and middle Miocene stress reorganization within the northeastern Tibetan Plateau. We attribute this to Oligocene lithospheric foundering beneath the south-central Tibetan Plateau and the middle Miocene initiation of sinistral strike-slip faults due to the underthrusting of India and North China. This study highlights the roles of lithospheric removal and continental underthrusting in plateau growth, and the value of multi-proxy apatite analyses in provenance studies.

Northeastern Tibetan Plateau↗

USGS2025

To heal or not to heal?: 2. The moment-recurrence time behavior of repeating earthquakes in the 2011 Prague, Oklahoma aftershock sequence is consistent with laboratory healing rates

The timing and failure conditions of an earthquake are governed by the interplay between fault reloading and restrengthening. The moment-recurrence time behavior of repeating earthquakes can give observational estimates of fault healing rates; however, it is difficult to link these observed healing rates to laboratory studies of frictional healing in part because of uncertainty in lithology. Here, we study the 2011 Prague earthquake sequence, which includes repeating earthquakes in the Arbuckle group and the granitic basement, and compare them to laboratory experiments on samples of the Arbuckle and Troy granite (representative of the basement rock) (Okamoto et al., 2025, https://doi.org/10.1029/2024JB030573 ). We find three spatially distinct groups of repeating earthquakes with different moment-recurrence behavior: (a) constant moment-recurrence time in the Arbuckle group, (b) scattered moment-recurrence time at the intersection of the foreshock-mainshock fault in the granitic basement, and (c) moment-predictable behavior outside of the foreshock-mainshock fault intersection also in the granitic basement. Our observation of stagnant healing for repeating sequences in the Arbuckle group is consistent with laboratory observations of low healing rates for moderately high pore fluid pressures in Arbuckle samples. For the moment-predictable group, the source radius that is required in order to match healing rates is consistent with source radius estimations when taking into account reasonable attenuation of the P- pulse width. Overall, we observe diverse healing behaviors in the seismic families that are consistent with laboratory healing rates, providing seismic evidence that contact-scale frictional mechanisms are relevant to large-scale earthquake dynamics.

Oklahoma↗

USGS2025

Rebuilding a volcano one lava flow at a time—Visualizing the lava dome-building eruption in the crater of Mount St. Helens, 1982–1986

Between 1980 and 1986, the U.S. Geological Survey made a series of 1:2,000-scale topographic contour maps from aerial photographic surveys to monitor the eruption. These maps were made for operational purposes and were not intended for publication. Since then, advances in technology made it possible to digitize the original, highly detailed hardcopy maps and derive new digital data elevation models of the surface of the lava dome. These digital elevation models allow for the visualization of the progression of the eruption and reveal the rubbly, chaotic surface of the lava flows and dome. Additionally, these new data help fill gaps in the long-term record of topographic changes that have occurred at the volcano since the May 18, 1980, eruption.

Washington↗

USGS2024

Temporal, spatial, and chemical evolution of Quaternary high-silica rhyolites in the Mineral Mountains, Utah

The Mineral Mountains in southwestern Utah are a structurally controlled core complex at the confluence of the Colorado Plateau and the Basin and Range physiographic provinces. Aside from hosting Utah’s largest batholith, the Mineral Mountains host some of the State’s youngest high-silica rhyolites, which have been linked to a magma source that is presently being utilized as an enhanced geothermal system. The high-silica rhyolites take the form of effusive lavas and domes, and explosive products are rare. Previous K-Ar dating of these Pleistocene rhyolites placed eruptions between about 790 and 500 kilo-annum (ka) with contemporaneous basalts erupting in the valley to the east of the Mineral Mountains. Large uncertainties on these ages obscured the tempo of eruptions and thus hindered attempts to constrain the timescales of the petrogenetic processes that produced the rhyolites. In this study, we build on previous studies conducted in the 1970s and 1980s by using new geochronologic and geochemical data to investigate the temporal and spatial evolution of the youngest phase of volcanism in the Mineral Mountains. We identify two major eruptive periods, from approximately 850 to 750 ka and from approximately 590 to 480 ka. The older phase is characterized by the eruption of several basaltic lavas, two obsidian flows, and a series of coalescing porphyritic rhyolite domes. The younger phase included the eruption of six evolved high-silica rhyolite domes and one pyroclastic deposit, followed by the eruption of trachyandesite in the adjacent valley to the east. Whole-rock geochemical data indicate that the rhyolites can be divided into three chemical groups, with more evolved compositions erupting through time. The youngest rhyolites along the range crest have the lowest total iron and TiO 2 concentrations and the highest incompatible element concentrations, indicative of increasing differentiation with time and elevation. Improved precision on the eruption ages indicates a recurrence interval of approximately 20 thousand years. The eruptive flux for both periods of rhyolitic volcanism is about 0.01 cubic kilometers per thousand years, which is less than the magma resurgence flux rates for syn-caldera and post-caldera eruptions of the Valles Caldera and Yellowstone Caldera volcanic systems. Collectively, these geochemical, geochronological, and volumetric data may facilitate a better understanding of heat flux and the longevity of magmatic sources related to geothermal resources in similar small-volume, silicic systems.

Utah↗

USGS2025

Introduction to the special section on improving measurements of earthquake source parameters

Earthquake source parameters such as magnitude, seismic moment, source dimension, stress drop, and radiated energy are fundamental to understanding earthquake physics, and are also key ingredients in earthquake ground‐motion modeling, rupture simulation, and statistical seismology. However, the uncertainties in these parameters estimated from the radiated seismic wavefield are large due to variability in approaches, including site and attenuation characterizations, and so estimates for an individual earthquake made by different studies can vary greatly. Estimating spectral source parameters remains a popular topic, due to a combination of their intrinsic importance and their apparent ease of measurement, but also a controversial one, due to many sources of variability and large uncertainty. The available methods coupled with necessary parameter choices and assumptions in the analysis make it challenging for researchers to apply methods or understand the reliability in results or reported source parameters. This Special Section on Improving Measurements of Earthquake Source Parameters showcasing comparisons between methods and studies seeks to alleviate some of these difficulties to help the community identify the important components and trade‐offs of decomposing recorded seismograms into their source, path, and site components.

Bulletin of the Seismological Society of America↗

USGS2025

Formation of the Mount Weld rare earth element deposit, Western Australia: A carbonatite-derived laterite

Carbonatite-hosted rare earth element (REE) deposits are the primary source of the world’s light REEs. The Mount Weld REE deposit in Western Australia is hosted in a lateritic sequence that reflects supergene enrichment of the underlying carbonatite. Water-rock interaction is a key to the formation of this world-class deposit. REE enrichment in the laterite is controlled by the breakdown of primary minerals, the release and transport of REEs, and the formation of secondary minerals. Secondary REE-bearing phosphate minerals are the primary REE-host phases in the laterite ore with monazite as the dominant phase; other REE-bearing phases include rhabdophane, cerianite, churchite, florencite, and crandallite subgroup minerals. Profiles through the laterite show that in the REE-rich zone, apatite and primary calcite and dolomite have broken down such that the loss of Ca and Mg, as well as Si and K, leads to a relative increase in the REEs. Sequestering of REEs in secondary mineral phases formed by groundwater further enhances the REE concentration.

Mount Weld Mine↗

USGS2026

Altitude of the potentiometric surface and depth to water in the Mississippi River Valley alluvial aquifer, spring 2024

Potentiometric-surface and depth-to-water maps for spring 2024 were created for the Mississippi River Valley alluvial aquifer using groundwater-altitude data from 1,151 wells completed in the Mississippi River Valley alluvial aquifer and from the altitude of the top of the water surface in area rivers from 160 U.S. Geological Survey and U.S. Army Corps of Engineers streamgages. The potentiometric-surface and depth-to-water maps for 2024 were created to support investigations to characterize the Mississippi River Valley alluvial aquifer as part of the U.S. Geological Survey Water Availability and Use Science Program. Sufficient data were available to map the potentiometric surface and depth to water of the Mississippi River Valley alluvial aquifer for spring 2024 for about 83 percent of the aquifer area. The potentiometric contours ranged from 0 to 330 feet (ft) above the North American Vertical Datum of 1988. The regional direction of the groundwater gradient was generally to the south-southwest, except in areas of groundwater-altitude depressions, where the groundwater gradient direction was into the depression, and near rivers, where the groundwater gradient direction was either from the aquifer to the river or from the river into the aquifer. Groundwater depressions in the potentiometric-surface map are in the lower one-half of the Cache region and in most of the Grand Prairie and Delta regions. Depth to water by well in the Mississippi River Valley alluvial aquifer in spring 2024 ranged from 0.27 ft above land surface to 145.65 ft below land surface.

Arkansas, Illinois, Kentucky, Louisiana, Mississip↗

USGS2024

A decade of shaking in the Garden City: The dynamics of preparedness, perceptions, and beliefs in Canterbury, New Zealand, and implications for earthquake information

This study explored earthquake preparedness over time - before, during, and 10 years after the Canterbury Earthquake Sequence (CES) in Aotearoa New Zealand (NZ). Surveys of Canterbury residents were conducted in 2009, 2013 and again in 2021, using variables derived from Community Engagement Theory (CET). The surveys measured earthquake perceptions and beliefs, participation and engagement, and preparedness actions. Results were compared across the three samples. Findings indicate that perceptions and beliefs (e.g. risk perception, outcome expectancy beliefs), and types of preparedness actions taken (e.g. collection of survival items, structural preparedness, community and agency relationships), differed over time, depending on people’s experiences before, during, and after the CES. For example, during and after the CES people were more likely to believe that preparing provided a benefit to daily life, but less likely to think it could reduce property damage, perhaps due to people’s experiences of disruption and damage during the earthquakes. An understanding of such dynamics can assist with the provision and timing of risk and preparedness information. This study highlights the importance of providing applicable and actionable preparedness information, that is relevant to people’s experiences, throughout an earthquake sequence. Such information might evolve and change in focus over time depending on risks and needs. Focus could also be given to information that builds peoples beliefs and capacities to undertake preparedness in evolving situations. Understanding preparedness in the context of different experiences and timeframes is useful in helping update models such as the CET, where the dynamics of time might be better incorporated.

Canterbury↗

USGS2025

Comparative properties of saponitic fault gouge and serpentinite muds cored from mud volcanoes of the Mariana subduction zone

We obtained 12 core samples for physical and chemical characterization from three serpentinite mud volcanoes (Yinazao, Asùt Tesoru, and Fantangisña) located on the forearc of the Mariana subduction system, that were drilled during International Ocean Discovery Program Expedition 366. Two samples from the Fantangisña mud volcano are interpreted to be clay-rich fault gouges derived from the subduction channel. Their bulk compositions are intermediate between the serpentinites and oceanic basalts. The oceanic crustal materials in the gouges have been thoroughly metasomatized and the serpentinites extensively altered to the trioctahedral, Mg-rich smectite clays saponite and corrensite. The only relict phases in clasts of crustal rock are accessory Ti- and P-bearing minerals. The two fault gouge samples have lower frictional strengths ( μ < 0.2) than the serpentinites ( μ = 0.2–0.4), and their measured permeabilities are also somewhat lower. Their physical and compositional properties correspond to saponitic gouges from other faults that juxtapose serpentinite against crustal rocks, in particular gouges from the two creeping traces of the San Andreas Fault recovered in the core from the San Andreas Fault Observatory at Depth. The décollement beneath Fantangisña mud volcano is thus expected to be very weak and likely characterized by stable slip.

Geophysics, Geochemistry, Geosystems↗

USGS2026

Importance of fish in the diet of Adélie penguins across multiple life stages

Over recent decades, the Adélie penguin ( Pygoscelis adeliae ) population has grown across most of its range, the exception being the northern coast of the western Antarctic Peninsula (WAP). In the Ross Sea, the very large Cape Crozier colony grew to reach ~9% of the global population. Demographic factors driving the Cape Crozier trend have yet to be identified. However, low chick fledging mass - a negative influence in the northern WAP where the diet is mostly less-energy-dense krill - is not showing an effect in the Cape Crozier colony. Previously, we hypothesized that Cape Crozier fledglings, and post-breeding adults, must be finding sufficient higher-quality prey once free of the prey-depleted colony foraging area, with subsequent higher survival. Here we explore penguin diet within and outside the colony foraging area using stable isotope analysis (SIA) of feathers: those grown by post-breeding adults outside the foraging area, prior to moult, compared to those of near-fledged chicks raised on meals obtained within the foraging area. Second, we explore whether a major difference occurs between breeding and post-breeding diet of adults from a small, nearby colony (Cape Royds) that exploits a smaller, little-depleted preyscape. Finally, we compare the results to SIA analyses published by others to explore whether the pre-moult diet of Ross Island adults, mostly foraging within the north-eastern Ross Sea, compares with that of adults nesting in the southern WAP (Bellingshausen Sea coast). At the latter, populations are not decreasing, with diet containing an appreciable contribution of fish. The results confirmed that fish contributed greatly to chick diets at both Ross Island colonies (δ 15 N 10.6–11.4). Once adults were released from central-place foraging, SIA levels (δ 15 N 9.7–12.9) were similar to those of the southern WAP. Foraging on energy-dense fish in coastal Antarctic waters could be one reason Adélie penguins are among the most abundant and therefore most ecologically successful penguin species.

Antarctic Science↗

USGS2025

A geospatial analysis of water-quality threats from orphan wells in principal and secondary aquifers of the United States

Throughout the history of oil and gas production in the United States, millions of wells have been drilled for exploration and energy production. Hundreds of thousands of unplugged wells are no longer actively producing and are currently under orphan status, with no responsible party obligated for plugging. Orphan wells can pose threats to water resources by providing pathways for contaminants such as hydrocarbons and brines to migrate into water-supply aquifers. In this study, we investigate the potential threats to groundwater resources posed by orphan wells at the national scale. Water-quality data is extremely sparse in relation to orphan wells nationally and may not be suitable for identifying contamination from oil and gas development. We used geospatial and statistical methods to evaluate which principal and secondary aquifer systems may be most susceptible to contamination from orphan wells. Analysis involved three sets of susceptibility factors including: 1) factors related to the number and density of orphan wells; 2) factors that can threaten well integrity and contribute to transport of contaminants; and 3) factors related to groundwater withdrawal rates and the affected populations/communities in the event of water quality disturbances. From a dataset of 117,672 documented orphan wells, 64,203 fall within a principal aquifer system, while the remainder fall within a secondary aquifer system. By assessing the combination of well integrity and hydrogeologic factors within these aquifer systems, five groupings of principal aquifers were identified, where groups ranged from aquifer systems with high numbers of orphan wells, multiple well integrity threats and high withdrawals, to aquifers with a relatively low number of orphan wells, limited well integrity threats and minimal water use. Three regions of the country emerge containing aquifers with higher susceptibility to contamination from orphan oil and gas wells. These regions include 1) The Appalachian Basin (including the Pennsylvanian Aquifer System), 2) The Gulf Coast Aquifers (including the Coastal Lowlands Aquifer system) and 3) The California Aquifers (including the California Coastal Basin Aquifer system). This work is the first multivariate geospatial investigation of orphan wells and groundwater resources on a national scale, and sheds light on which aquifers are most susceptible to groundwater contamination from orphan wells.

Science of the Total Environment↗