USGS Science⌕ Search

SEARCH · USGS Science

Results for “Data Series”

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

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

887 records · Page 46Linked to original sources
USGS2024

Application of non-stationary shear-wave velocity randomization approach to predict 1D seismic site response and its variability at two downhole array recordings

Accounting for uncertainties in seismic site response is crucial to improving the performance of one-dimensional (1D) ground response analyses (GRAs) at downhole array recording sites. In addition to site effects, uncertainties in 1D-GRAs can also be contributed from the seismic source and/or path. Though often representing not more than one percent of the distance (path) from the source, site conditions are known to have an enormous influence on ground shaking. In this study, we focus on the site shear-wave velocity ( V S ) structure, which is the main ingredient for estimating the variability of site response. As such, V S can manifest aleatory uncertainties related to the effects of small-scale spatial heterogeneities within the near surface, thus V S can substantially modify ground shaking during earthquakes. We apply a novel V S randomization approach to propagate the small-scale heterogeneities of V S to estimate seismic site response within a non-stationary probabilistic framework. The randomization approach generates samples of V S profiles that are used to perform several 1D-GRAs and obtain an averaged site response and related variability. The proposed method is implemented on data recorded at two downhole array sites with different subsurface soil conditions: a soft soil site on Treasure Island (California, United States of America) and a rock outcrop site in Cadarache (South-East France). We show that synthetic surface-to-borehole transfer functions from 1D-GRAs provide an acceptable fit to the empirical transfer functions from low-motion earthquake records and succeed in reproducing most of the site-specific seismic response variability. The remaining mismatch between transfer functions is likely due to insufficient precision on the seismic bedrock and the impedance contrast. The variability in site response is discussed with emphasis on the role of V S small-scale heterogeneities, attenuation, and input motion incidence angle in ground motion variability for the site and soil conditions at both locations.

California↗
USGS2024

Indications of preferential groundwater seepage feeding northern peatland pools

Groundwater seepage from underlying permeable glacial sedimentary structures, such as eskers, has been hypothesized to directly feed pools in northern peat bogs. These hypotheses directly contradict classical peat bog models for ombrogenous systems, wherein meteoric water is the sole water input to these systems. Variations in the underlying mineral sediment in contact with the peat imply that unrecognized hydrogeologic connectivity may exist with pools in northern peat bogs, particularly where high permeability materials are in contact with the peat. Seepage dynamics originating from these structural variations were investigated using a suite of thermal and hydrogeophysical methods deployed around pools in a peat bog of northeastern Maine, USA. Thermal characterization methods mapped anomalies that were confirmed as matrix seepage or preferential flow pathways (PFPs). Geochemical methods were employed at identified thermal anomalies to confirm upwelling of solute-rich groundwater. Conduits around pools were associated with surficial terminations of suspected peat pipes, based on the inference of pathways extending down into the peat, that focus flow through PFPs in the peat matrix. Discharge also occurred through the peat matrix adjacent to suspected pipe structures and matrix seepage rates were quantified using analysis of diurnal temperature signals recorded at multiple depths. Seepage rates, with a maximum of nearly 0.4 m/d, were measured at localized points around pools. Periods of synchronized temperatures paired with highly muted diurnal temperature signals, recorded in diurnal temperature with depth data, were interpreted qualitatively as activation of strong upward discharge rates through suspected peat pipes. These time periods correlated strongly with local precipitation events around the peatland. Ground-penetrating radar surveys revealed discontinuities in the low permeability glacio-marine clay at the mineral sediment-peat interface, interpreted to be regional glacial esker deposits, which were located beneath and around pools. Heat tracing, specific conductance contrasts, seepage rates, and trace metal concentrations all imply groundwater seepage originating from underlying permeable glacial esker deposits and directly sourcing pools. Preferential groundwater inputs into northern peat bogs may play a key role in developing and maintaining pool systems, with enhanced solute transport impacting peatland ecology, water resources, and carbon cycling.

Maine↗
USGS2024

Distribution of ancient carbon in groundwater and soil gas from degradation of petroleum near the Red Hill Bulk Fuel Storage Facility, O‘ahu, Hawai‘i

The groundwater below the Red Hill Bulk Fuel Storage Facility (the facility) in Oʻahu, Hawaiʻi, contains fuel compounds from past spills. This study used carbon-14 analyses to distinguish fuel-derived carbon from background carbon, along with other biodegradation indicators, to address two goals: (1) determine the extent and migration direction of groundwater affected by residual fuel below the facility and (2) determine if residual fuel locations in the subsurface could be identified by analyzing soil gas at the surface above the facility. Groundwater from 19 wells was sampled between September 2022 and April 2023. Nonvolatile dissolved organic carbon (NVDOC) from a well presumed to be unaffected by past spills contained 38 percent ancient carbon indicating a natural source of ancient carbon in the subsurface. The NVDOC concentrations and ancient carbon percentages indicate fuel biodegradation products are likely present on the north and south of Red Hill with the greatest effects at well RHMW02 near the 2014 spill site. The NVDOC concentrations are almost three times higher than diesel range organic (DRO) concentrations in groundwater from the same sites. Major ion data indicate that iron reduction is an important biodegradation process. Soil probe samples and soil carbon traps were used to determine the carbon-14 content of soil carbon dioxide. Ancient carbon from fuel biodegradation was not detected at any soil probe or carbon trap site in contrast to a 2017 study which reported ancient carbon detections. A reanalysis of the 2017 results using a range of local values for background carbon-14 indicates that ancient carbon from fuel biodegradation was probably only detected in lower tunnel exhaust system samples and not in any soil carbon trap samples. Measurements of carbon dioxide efflux with a dynamic closed chamber were highly variable. The soil gas results indicate that soil gas measurements at land surface were not useful for detecting residual fuel at the facility.

Hawaii↗
USGS2025

Stratigraphy, structure, and geomorphology of the central Appalachians across the North Mountain fault zone near Harrisonburg, Virginia, USA

This field trip focuses on the geology of the central Appalachian Valley and Ridge province near Harrisonburg, Virginia, USA. Recent geologic mapping utilizing 1-m resolution lidar data has revealed new insights into the Paleozoic stratigraphy, structural geology, and Neogene landscape evolution of the region. The detailed mapping reveals the presence of the Big Spring Station Member and multiple thrombolite zones in the Cambrian Conococheague Formation extending as far south as the Briery Branch 7.5 min quadrangle, providing insights into Late Cambrian sea-level fluctuations. Multiple outcrop exposures in the study area of this guidebook confirm recent work in Pennsylvania, USA, showing that the Ordovician Reedsville Shale overlies the Martinsburg Formation and that the two are distinct and mappable as separate formations rather than laterally equivalent units as previously interpreted. Our work extends the Silurian Williamsport Sandstone into Shenandoah County, Virginia, and describes its facies relationships with the Bloomsburg Formation along strike and across the Adams Run anticline. Mapping within the thick Devonian siliciclastic sequence reveals the presence of the Mahantango Formation on the western limb of Supin Lick syncline and illustrates its complex facies relationship with the Millboro Shale. In addition, we highlight new mapping criteria for the Brallier and Foreknobs Formations and demonstrate how the specific changes to the placement of the contact between them addresses previous challenges in their differentiation. We present cosmogenic burial ages of broad alluvial fan sediments in the Shenandoah Valley near Timberville and Briery Branch, Virginia, and erosion rates estimated for the Briery Branch stream basin. Both analyses provide new constraints on the timing of landscape evolution and karst development since the middle Pliocene. This field guide also highlights some significant structural features within the North Mountain fault zone, such as evidence of imbricated thrust sheets cut by cross-strike faults that have been exploited by Eocene igneous intrusions. Map-scale horses of Silurian and Ordovician rocks hold up ridges that are oblique to the regional strike. Deformation internal to one of these horse blocks is shown to be non-coaxial with respect to the main regional northwest directed transport.

Virginia↗
USGS2026

Mountain goat declines in a protected, interior, native population

A shifting climate poses threats to alpine-adapted species including mountain goats. We used long-term (12 years) citizen science monitoring data and Bayesian N-mixture modeling to estimate population trends and drivers of population metrics among mountain goats in Glacier National Park (GNP). Median goats per site ( n = 37 sites) declined by 45% (95% credible interval [CRI] = 32%, 57%) from 77.8 (95% CRI = 64.4, 95.1) in 2008 to 42.3 (95% CRI = 34.3, 52.2) in 2019, with consistent declines from 2008 until 2015, when the number of estimated goats stabilized. The decline exceeds IUCN criteria for classifying a population as vulnerable, >30% declines over only two generations. Across years, relatively few goats occupied northwestern GNP. Goat numbers declined the most at northeastern sites, trended toward decline in most southern sites, and increased at only two west-central sites. The proportion of permanent snow and glaciers, the presence of natural mineral licks, and habituation strongly increased the initial abundance of goats in the area. Weather variables had the greatest influence on population growth rates, particularly precipitation between May 15 and June 15 of the previous summer, the neonatal period. Lower growth occurred with less snow water equivalent and lower mean winter temperature, early summer temperature, and early summer precipitation. Projected reductions of permanent snow, increasing spring and summer temperatures, and insufficient and variable spring precipitation raise concerns for the future of native goats in this region. Our analyses reveal ways to improve detection rates of goats during surveys, which is important for optimizing the precision of estimates and the power to detect future trends. Detection increased with goat habituation, retention of observers with experience, use of binoculars, and conducting surveys at lower temperatures and earlier dates. Improving detection will be particularly important given the lower number of goats currently observed in the park. Research to estimate park-wide population size, evaluate genetic structure and diversity, assess changing habitat, human recreation levels and forage, and forward-project climate effects on persistence will be crucial to understanding the context of these results and conserving this iconic, metapopulation at the southern edge of the distribution of native mountain goats.

Montana↗
USGS2024

U.S. Geological Survey Earthquake Hazards Program decadal science strategy, 2024–33

Executive Summary Earthquakes represent one of our Nation’s most significant and costly natural hazards, with estimated annual loses from earthquakes close to $15 billion in 2023. Over the past two centuries, 37 U.S. States have experienced an earthquake exceeding a magnitude of 5, and 50 percent of States have a significant potential for future damaging shaking; these statistics speak to the need for nationwide interest and investment in earthquake hazard characterization and risk reduction. Authorized under the Earthquake Hazards Reduction Authorization Act, the U.S. Geological Survey (USGS) Earthquake Hazards Program (EHP) provides the scientific information, situational awareness, and knowledge necessary to reduce deaths, injuries, and economic losses from earthquakes and earthquake-induced tsunamis, landslides, and soil liquefaction. The EHP supports activities in three focused topical areas: (1) earthquake monitoring, (2) hazard assessment, and (3) applied research, using the results of each—and the coordination among them—to further support risk translation and communication in regions at risk nationwide. For earthquake monitoring, the Advanced National Seismic System (ANSS), a cooperative effort of USGS networks, university partner regional seismic networks, and real-time geodetic networks, collects and analyzes data on earthquakes; issues timely, reliable notifications of their occurrence and impacts; and provides data for earthquake research, hazard, and risk assessment as a foundation for building an earthquake-resilient Nation. The USGS-operated ShakeAlert Earthquake Early Warning system is a recent addition to EHP’s ANSS infrastructure. In the realm of earthquake hazard assessment, the EHP contributes to earthquake risk mitigation strategies by developing the National Seismic Hazard Model and maps, and other related products, that describe the likelihood and potential effects of earthquakes nationwide, especially in the urban areas of highest risk. The EHP also conducts research on the causes, characteristics, and effects of earthquakes and prioritizes work that directly increases the accuracy and precision of earthquake hazards assessments, earthquake forecasts, and earthquake monitoring and situational-awareness products and that supports the Nation’s earthquake mitigation practices. Bridging the EHP’s efforts across research, hazard assessments, and earthquake monitoring is a broad and comprehensive collection of earthquake information products, including the National Seismic Hazard Model, ShakeAlert, and other products describing impact, such as ShakeMap and PAGER (Prompt Assessment of Global Earthquakes for Response), which have been developed and integrated into EHP’s real-time monitoring systems. EHP funds external partners to carry out many important collaborative activities through an active external grants program—one of the largest in the USGS—and through cooperative agreements with other partners such as the university-operated regional seismic networks, funded as part of the ANSS. To continue its support of earthquake hazard characterization and risk reduction, the EHP aims to strengthen its foundational products and practices while positioning itself to respond to the evolving needs of the Nation and follow best practices of the scientific community. This document describes a strategy for the program to ensure it can meet these demands. The foundational priorities outlined in this strategy represent those activities that remain critical to the core functionality of the program and those that can be supported under current fiscal year 2024-level appropriations. Priorities described as aspirational are important for future growth, and to maintain the program’s position as a leading global resource in earthquake science, but would require increases in appropriated funding to be fully realized. Across the program’s portfolio of activities, several major themes have been identified as the most critical activities to advance EHP science over the coming decade. Together, these activities provide the framework necessary to integrate critical hazard characterization and risk reduction activities across the program. They provide the structure for research to advance the understanding of where, when, and why earthquakes occur and how we can use improved knowledge to drive short-term and actionable forecasts of seismic activity. They expand the usefulness of critical earthquake products and advance the sophistication of those products to keep pace with the rapidly evolving needs of an ever-expanding user base while maintaining the position of the USGS as a global leader in earthquake science. Focus on system-level science. Establish an automated earthquake-processing pipeline. Enhance the accuracy and reliability of the ShakeAlert Earthquake Early Warning system and plan for extension to other regions. Implement time-dependent earthquake forecasting. Develop physically realistic models. Expand computational capacity. This science strategy is organized into three primary sections. The first section provides an overview of the EHP and its budget, governance, and program council. Readers familiar with the program may wish to focus on the second section, which describes the core of the science strategy, including priorities across each of the EHP’s major program activities in monitoring, hazard assessment, and targeted research. The third section outlines science priorities that cut across program activities, including those involving collaborations external to the EHP.

American Samoa, Guam, Northern Marianas Islands, U↗
USGS2024

Geophysical constraints on continental rejuvenation in central China: Implications for outward growth of the Tibetan Plateau

Continental rejuvenation results from the tectonic reactivation of crustal structures and lithospheric reworking by mantle flow. Geochemical observations and field mapping have traditionally provided the primary evidence for the secular evolution of crustal composition and tectonic processes during continental rejuvenation. Nonetheless, the impact of continental rejuvenation on the observed present-day strain rate and orogenic-scale lithospheric structure has not been well constrained. The pre-existing E-W–trending Central China Orogenic Belt has been overprinted by the N-S–trending Central Longitudinal Seismic Belt and constitutes the intracontinental West Qinling Syntaxis in central China, where the tectonic setting changes eastward from contraction to extension. Combining updated global positioning system data and high-resolution crustal seismic tomography, we reveal a modern continental rejuvenation process within the West Qinling Syntaxis in central China. The northward extrusion of the Tibetan Plateau's weak lithospheric layer (middle-lower crust and lithospheric mantle) of southwestern China relative to the rigid Sichuan Basin/Ordos Block of the eastern West Qinling Syntaxis results in regional dextral shearing that shapes the Central Longitudinal Seismic Belt and defines the eastern Tibetan Plateau margin. The pre-existing E-W–trending Central China Orogenic Belt has been preserved above the brittle-ductile transition zone, and the northward movement of the deep lithospheric layer drives the deformation of the upper crust in the West Qinling Syntaxis. Our results, along with previous studies, suggest the presence of an intracontinental lithospheric interchange structure in central China. The continental rejuvenation of the West Qinling Syntaxis results from a combination of fault reactivation in the upper crust (Stage I, Eocene–Oligocene) and reworking of the deep lithosphere (Stage II, middle–late Miocene) related to the plateau-wide shift in stress accommodation ultimately driven by the redistribution of mass outward from the central Tibetan Plateau. At present, the transition zone between the high- and low-velocity anomalies along the Central Longitudinal Seismic Belt not only shapes the landscape boundary but controls the size and recurrence interval of earthquakes within the West Qinling Syntaxis in central China.

GSA Bulletin↗
USGS2024

Status and performance of the ShakeAlert® earthquake early warning system: 2019-2023

The U.S. Geological Survey (USGS)‐operated ShakeAlert® system is the United States West Coast earthquake early warning system ( Given et al. , 2018 ). In this study we detail ShakeAlert’s performance during some of the largest events seen by the system thus far. Statewide public alerting using ShakeAlert messages was authorized in California in October 2019. Over the next few years, public alerts were expanded into Oregon and then into Washington ( U.S. Geological Survey, 2024 ). ShakeAlert source results are routinely compared to the USGS Comprehensive Catalog (ComCat; Guy et al. , 2015 ; U.S. Geological Survey, Earthquake Hazards Program, 2017 ), which contains the earthquake location and magnitude determined using complete waveform data. M 4.5 and larger is the threshold used for public alerting and was deliberately set below the level where damage is likely to compensate for cases where the system underestimates the magnitude. Between 17 October 2019 and 1 September 2023, the ShakeAlert system created 95 events with maximum magnitude estimates of M ≥4.5, the public alerting threshold. 94 of the 95 events were due to real earthquakes. Seven were categorized “false” per ShakeAlert’s internal definition that there was no matching catalog event within 100 km and 30 s of origin time; however, all but one of these were real earthquakes that were poorly located, primarily because they were at the edges of the seismic network. Three detected events were labeled “missed” because they were very poorly located (>100 km location error). In addition, the system did not produce solutions for four ComCat events M ≥4.5 ( U.S. Geological Survey, Earthquake Hazards Program, 2017 ), which were all at the edge of the alerting and network boundaries. The ShakeAlert system has accurately detected the majority of earthquakes that have occurred within the operational region since completing the public rollout, and alerts from the system have been delivered to millions of cell phone users throughout the West Coast.

California, Oregon, Washington↗
USGS2025

Overview of The SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest earthquake sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress‐drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude‐spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress‐drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two‐week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. In this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite‐fault modeling, ground‐motion‐based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between‐method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress‐drop values.

California↗
USGS2025

Reconstruction of Holocene and Last Interglacial vegetation dynamics and wildfire activity in Southern Siberia

Wildfires are a rapidly increasing threat to boreal forests. While our understanding of the drivers behind wildfires and their environmental impact is growing, it is mostly limited to the observational period. Here we focus on the boreal forests of southern Siberia and exploit a U–Th-dated stalagmite from Botovskaya Cave, located in the upper Lena region of southern Siberia, to document wildfire activity and vegetation dynamics during parts of two warm periods: the Last Interglacial (LIG; specifically part of the Last Interglacial maximum between 124.1 and 118.8 ka) and the Holocene (10–0 ka). Our record is based on levoglucosan (Lev), a biomarker sensitive to biomass burning, and on lignin oxidation products (LOPs) that discriminate between open and closed forest and hard- or softwood vegetation. In addition, we used carbonate carbon stable isotope ratios ( δ 13 C ), which reflect a dominant control of the host rock, to evaluate soil respiration and local infiltration changes. Our LOP data suggest that, during the Last Interglacial, the region around Botovskaya Cave was characterised by open forest, which by ca. 121.5 ka underwent a transition from fire-resistant hardwood to fire-prone softwood. The Lev record indicates that fire activity was high and increased towards the end of Last Interglacial just before 119 ka. In contrast, the Holocene was characterised by a closed-forest environment with mixed hard- and softwood vegetation. Holocene fire activity varied but at a much lower level than during the Last Interglacial. We attribute the changes in wildfire activity during the intervals of interest to the interplay between vegetation and climate. The open forests of the Last Interglacial were more likely to ignite than their closed Holocene equivalents, and their flammability was aided by warmer and drier summers and a stronger seasonal temperature contrast due to the increase in seasonal insolation difference compared to the Holocene. Our comparison of the last two interglacial intervals suggests that, with increasing global temperatures, the boreal forest of southern Siberia may become progressively more vulnerable to higher wildfire activity.

Botovskaya Cave, Siberia↗
USGS2026

Linking distribution and return-on-investment models to optimize woody management for prairie grouse in Nebraska

Grasslands in Nebraska, USA, face threats from agricultural conversion, urban development, and woody encroachment, all of which negatively affect prairie grouse ( Tympanuchus spp.) populations. To optimize conservation planning, Nebraska wildlife agencies developed probabilistic area-based surveys for greater prairie-chicken ( T. cupido ) and sharp-tailed grouse ( T. phasianellus ) to sample landscapes across a range of environmental conditions. This design improves historical surveys and enables the development of distribution models that quantitatively define habitat associations and support scenario-based conservation planning. Using survey data collected during 2020–2022, we modeled prairie grouse occurrence and abundance as functions of land cover, topography, and climate using Bayesian logistic and zero-inflated negative binomial models with regularized horseshoe priors. We then conducted a maximum potential return-on-investment analysis of woody cover treatments, assuming sustained treatment success, relative to projected impacts of woody encroachment on prairie grouse populations by 2050. Among modeled associations were a positive association with grasslands having low woody cover and a negative association with grasslands having high woody cover. Across the 3-year period, median estimated annual populations were 142,380 for greater prairie-chicken (range of 95% CIs across years = 68,821–277,615) and 64,154 for sharp-tailed grouse (range of 95% CIs across years = 27,550–146,559). Under projected woody encroachment, mean predicted population declines were 10% for greater prairie-chicken (range of 95% CIs = 7–14%) and 6% for sharp-tailed grouse (range of 95% CIs = 5–7%). Areas with high prairie grouse density and low treatment costs in 2021, and high projected woody encroachment and population loss by 2050, offered the greatest return on investment for woody management. Return on investment was greatest in the northwestern Shortgrass Prairie ecoregion (northwestern Nebraska) for sharp-tailed grouse and the eastern Sandhills ecoregion (central Nebraska) for both species. These models underscore the value of evidence-based, quantitative approaches for prioritizing conservation actions on working lands. Scenario-based modeling could be extended to guide other treatments, such as optimizing restoration (e.g., Conservation Reserve Program) or incentivizing grassland persistence in areas with predicted climate resilience.

Nebraska↗
USGS2025

Fluid inclusion constraints on the geometry of the magmatic plumbing system beneath Mauna Loa – Part 2: Xenoliths

Mauna Loa volcano erupts crystal-poor material at its summit and more crystal-rich material on its rift zones. Some of the more olivine-rich lava flows contain xenoliths with diverse mineralogy, including cumulate harzburgites with high-Mg# orthopyroxenes and high-Fo olivines (both > 84). Previous experimental work and thermodynamic modelling has proposed that high-Mg# orthopyroxenes only crystallize from Mauna Loa melts at high pressures (> 6 kbar, > 20 km), leading to suggestions that there is a region of sub-Moho magma storage at Mauna Loa in addition to the geophysically imaged magma reservoir at 2–5 km depth below the summit. We use melt and fluid inclusion barometry combined with thermodynamic models to further investigate this suggestion. Fluid inclusion data from harzburgites and dunitic xenoliths yield storage depths remarkably similar to those found in non-xenolithic crystals from lavas and tephras, with a clear peak at ~ 2–3 km (below the summit). Depths from melt inclusions in these xenoliths overlap with fluid inclusion pressures, ruling out the possibility of fluid inclusion re-equilibration during a period of stalling in a shallower reservoir. We examine five different thermodynamic models and find that the minimum pressure of olivine-orthopyroxene co-saturation varies by ~ 4 kbar (~ 12 km). These models also fail to predict that orthopyroxene is stable in ~ 15–80% of compositionally relevant experimental charges which grew orthopyroxene. Overall, this shows that phase stability modelling is an unreliable method of determining magma storage depth at Mauna Loa. We suggest that model discrepancies reflect a lack of experimental constraints on orthopyroxene stability at > 1200 ℃ and 0.01–5 kbar. Based on the presence of large oikocrystic orthopyroxenes completely enclosing rounded olivine chadacrysts, we suggest that these harzburgitic xenoliths formed through the reaction of intruding melts with olivine mush piles within the Mauna Loa edifice at ~ 3 km depth below the summit, with no need for a deeper storage reservoir. The predominance of pre-eruptive shallow storage means that there is more chance of detecting reservoir destabilization with geophysical monitoring techniques compared to a scenario where melts are supplied from sub-Moho reservoirs.

Hawaii↗
USGS2025

Groundwater structures fish growth and production across a riverscape

Landscapes are composed of habitat patches and conditions that vary across space and time. While habitat variability and complexity can support important ecological processes and ecosystem services, the dynamic nature of habitats can also constrain organismal growth and production as optimal conditions are fleeting. In riverine ecosystems, groundwater discharge to streams stabilises water temperature and flow regimes, thus mediating how habitat complexity is expressed. Yet, how stable habitats structure growth and production within the broader landscape matrix is not well understood. In this study, we explored the effects of groundwater on spatiotemporal variation in growth and production for juvenile Yellowstone cutthroat trout ( Oncorhynchus virginalis bouvieri ) across the upper Snake River catchment, Wyoming, USA. We combined machine learning techniques and remotely sensed landscape data to estimate groundwater availability across the river network, which we linked to stream temperature regimes and conspecific density. We then used Bayesian hierarchical models to quantify the effects of temperature, density and groundwater on spatiotemporal variation in fish growth and production in 52 focal reaches. Finally, we predicted body size trajectories and trends in total production continuously over both space and time to understand the effect of groundwater at the riverscape scale. Groundwater discharged to streams where topography changes abruptly in valley-bottom areas underlain by coarse glacial deposits. Groundwater stabilised temperature regimes and was associated with high trout densities. Temperature and density, in turn, interacted to influence growth rates: growth increased strongly with temperature, but this effect was reduced when density was high. Accordingly, variation in groundwater availability among stream reaches diversified growth and production regimes. In reaches with low groundwater availability, growth and production declined over time from summer maxima. In contrast, in reaches with high groundwater availability, temporal trends in growth and production were hump-shaped—peaking in autumn—and mean production was greater. At the riverscape scale, temporal asynchrony in growth rates generated convergent spatial variation in growth capacity, but—when combined with density—led to the formation of distinct hotspots of production. Our results demonstrate how groundwater, an important driver of aquatic ecosystem heterogeneity, structures trout growth and production across space and time. Importantly, rare, but stable habitats may disproportionately affect ecological processes and serve as key sources of population diversity at larger spatial scales.

Wyoming↗
USGS2025

Metal fingerprints of Eocene rhyolite magmas coincident with Carlin-type gold deposition in Nevada USA

Eocene magmatic systems contemporaneous with world-class Carlin-type Au deposits in Nevada (USA) have been proposed by some researchers as a key ingredient for Au mineralization, though evidence conclusively demonstrating their genetic relationship remains tenuous. This study provides the first direct evidence of the pre-eruptive metal budget of volatile- and metal-charged silicic magmas coincident in time (~41 to 34 Ma) and space (within 5 km) with Carlin-type Au deposits. We characterize the pre-eruptive metal fingerprints of these diverse magmatic systems to assess their potential as sources of metals for Carlin-type Au mineralization. Metal abundances from quartz-hosted melt inclusions (Au, Te, Ag, Sb, Tl, Mo, W, Sn, As, Pb, Co, Cu, Ni, and Zn) characterized in situ by SHRIMP-RG and LA-ICP-MS represent our best (and only) estimates for the pre-eruptive metal budget in these systems. Median metal concentrations are generally within one order of magnitude of average upper crust and average continental rhyolite values. But there are two notable exceptions, with median Au contents extending >1 order of magnitude higher than average upper crust and median Cu contents ranging >1 order of magnitude lower than upper crust. Despite this, melts contain lower Au/Cu (<0.1), Au/Ag (<5), and Au/Tl (<0.3) than most ore-grade Carlin-type rock samples and quartz-hosted fluid inclusions, regardless of their age and timing relative to nearby Carlin-type Au mineralization. The metal fingerprints of these magmatic systems, de-fined both by traditional and multivariate compositional data analysis techniques, are distinct from one another. Yet none are particularly specialized, e.g., high Au/Cu, in terms of being ideal ingredients as postulated by magmatic models for Carlin-type Au mineralization. Magmatic Au contents do not appear to be correlated with rhyolite “flavors” in the way that Cu, Sn, and Nb contents are. Fluid/melt partitioning modeling and magma volume estimates support the idea that a diverse array of non-specialized silicic magmas could feasibly contribute some or potentially all of the Au, Ag, and Cu in Carlin-type systems. The compositional diversity among contemporaneous magmatic systems could possibly contribute to some of the diversity observed across Carlin-type Au districts in Nevada.

Nevada↗
USGS2024

The addition of 144Nd atomic mass to routine ICP-MS analysis as a Quick Screening Tool for Approximating Rare Earth Elements (Q-STAR) in natural waters

Rare earth elements (REEs) are a class of critical minerals, all of which can have supply chain vulnerability that impacts economic security. These elements are widely measured in environmental matrices via inductively coupled plasma mass spectrometry (ICP-MS); however, successful quantification can require time-consuming, sample-specific optimization. While a sample-by-sample approach is appropriate for targeted quantification studies, this approach is not suitable for mineral exploration efforts where rapidly screening thousands of samples for the presence of REEs is desired. Here, we demonstrated the use of a Quick Screening Tool for Approximating REEs (Q-STAR) to detect REEs in surface water and groundwater matrices, collected as part of existing environmental studies. A mass-to-charge ratio of 144 ( m / z = 144) was added to an ICP-MS method to screen for REEs in filtered water samples submitted for metals analyses to the U.S. Geological Survey (USGS) National Water Quality Laboratory. We detected the presence of REEs above a reference threshold of 1200 counts per second in 18 % of pre-selected 6626 samples. Using this screened dataset, we mapped estimated dissolved REE concentrations across the United States in relation to ecoregions and underlying geology. Data are constrained to where sample collection took place but nevertheless show estimated aqueous dissolved REE concentrations on a geographic scale that has not yet been studied. To validate Q-STAR, REEs were measured in a USGS standard reference sample, a subset of 88 archived filtered water samples, and in fresh filtered surface water samples. Our targeted analyses demonstrated a strong linear relationship between Q-STAR predicted and measured values in all archived samples for Nd (r 2 = 0.94), and light REEs (LREEs) such as lanthanum (La) (r 2 = 0.93), praseodymium (Pr) (r 2 = 0.94) and samarium (Sm) (r 2 = 0.94). Using Q-STAR screen values, nine field sites were identified and surface water samples recollected to confirm the continued presence of Nd and LREEs. Q-STAR can be used to screen an unlimited number of water samples for the presence of REEs prior to time-intensive and costly quantitative analyses and to generate large REE datasets for further investigation.

Journal of Geochemical Exploration↗
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↗
USGS2025

Reconnaissance of potential alternate water supply sources for the City of Gary, West Virginia

Seven potential sources of water, consisting of free-flowing discharge from abandoned coal mines at six locations and one abandoned flooded underground coal mine air shaft, were sampled for chemical analysis to assess the quality of the groundwater emanating from the seven mine sources. The six free-flowing mine discharge sources were also assessed for discharge by current-meter measurements on two separate occasions. The U.S. Geological Survey assessed these seven sources to provide information to the City of Gary, West Virginia (W. Va.), and the City of Gary’s consulting engineer with groundwater-quality and flow data to allow them to assess the seven sites as potential alternate sources of water for the City of Gary to augment its existing supply. For the six sites where discharge could be measured, discharge ranged from a minimum of 0.082 cubic feet per second (ft 3 /s) to a maximum of 3.685 ft 3 /s. Of the six sites measured, only two, Harmon Branch at Thorpe, W. Va. (USGS site 372201081303501) and the abandoned public-supply water wells near Havaco, W. Va. (USGS site 372358081344601), had discharge in excess of 1.00 ft 3 /s. Discharge from the abandoned public supply wells was 3.685 ft 3 /s on September 20, 2023, and 2.888 ft 3 /s on October 16, 2023, and discharge from Harmon Branch at Thorpe, W. Va., was 1.049 ft 3 /s on September 22, 2023, and 1.038 ft 3 /s on October 17, 2023. Discharge in the abandoned underground mine air shaft (USGS site 372224081340901) could not be assessed, but the air shaft drains an abandoned mine that likely contains water stored in approximately 1.7 square miles (mi 2 ) of abandoned underground coal mines in the Pocahontas No. 3 coal seam, and possibly an additional 0.9 mi 2 of leakage from the overlying Pocahontas No. 4 coal seam. Discharge for the six sites measured for the study was measured during a period between September 20 and October 18, 2023, and corresponded to the 12th to the 15th percentile of flow-duration statistics for the Tug Fork downstream of Elkhorn Creek at Welch, W. Va. streamgage (USGS site 03212750). Water-quality data for the seven sites sampled overall were acceptable with respect to drinking water standards. Of the 203 constituents analyzed, only a few failed to meet applicable U.S. Environmental Protection Agency (EPA) drinking water standards. Iron exceeded the 300 micrograms per liter (μg/L) secondary maximum contaminant level (SMCL) at only 1 of the 7 sites (14.3 percent) sampled. Iron concentrations ranged from a minimum of less than (<) 5.00 μg/L to a maximum of 724 μg/L with a median concentration of 7.62 μg/L. Manganese exceeded the 50.0 μg/L SMCL at 2 of the 7 sites (28.6 percent) sampled. Manganese concentrations ranged from a minimum of 1.93 μg/L to a maximum of 271 μg/L with a median concentration of 4.03 μg/L. No sites sampled exceeded the arsenic maximum contaminant level (MCL) of 10 μg/L. Arsenic concentrations ranged from a minimum of <0.100 μg/L to a maximum of 2.35 μg/L with a median arsenic concentration of 0.200 μg/L. None of the seven sites sampled for selenium for this study exceeded the EPA MCL of 50.0 μg/L. Selenium concentrations ranged from a minimum of <0.050 μg/L to a maximum of 5.26 μg/L with a median concentration of 3.21 μg/L. All seven sites were sampled for volatile organic compounds (VOCs), semivolatile organic compounds (SVOCs), and polychlorinated biphenyls (PCBs), but most had concentrations below the detection limit. Of the 10 PCB compounds analyzed for the seven sites sampled, none contained detectable concentrations of PCBs or Aroclor compounds. Of the 44 SVOCs analyzed at each of the seven sites sampled, only 1 SVOC, acenaphthene, was detected, at a concentration of 0.02 μg/L. Of the 96 VOCs analyzed, from each of the seven sites sampled, only two were found at detectable concentrations. Trichloromethane was detected only at 1 of the 7 (14.3 percent) sites sampled at a concentration of 0.027 μg/L, and benzene was detected at the same site and 3 additional sites (4 of the 7 sites or 57.1 percent of the sites sampled) at concentrations of 0.028, 0.029, 0.021, and 0.035 μg/L, but none exceeded the EPA MCL for benzene of 5.00 μg/L. Total coliform bacteria are ubiquitous in the environment, and their presence only suggests the potential for contamination by near-surface processes. Escherichia coli ( E. coli ) bacteria are derived from either human or animal fecal material and can be an indicator of potential contamination by pathogenic bacteria or viruses. Total coliform bacteria were detected at all 7 sites sampled at concentrations ranging from 17.5 to greater than (>) 2,420 most probable number per 100 mL (MPN/100 mL) of sample, with a median total coliform concentration of 1,553 MPN/100 mL. Escherichia coli bacteria were detected at 4 of the 7 sites sampled at concentrations ranging from <1 to 11.9 MPN/100 mL, with a median E. coli concentration of 5.1 MPN/100 mL.

West Virginia↗
USGS2025

Water supply in the conterminous United States, Alaska, Hawaii, and Puerto Rico, water years 2010–20

We present an assessment of water supply across the conterminous United States (CONUS), Alaska, Hawaii, and Puerto Rico covering water years 2010–20. Our analysis drew on two national hydrologic models, the National Hydrologic Model Precipitation-Runoff Modeling System and the Weather Research and Forecasting model hydrologic modeling system. Both models produced estimates of streamflow, evapotranspiration, soil moisture, snow water equivalent, and other hydrologic states and fluxes. The models were driven by the bias-adjusted 4-kilometer-resolution, long-term regional hydroclimate simulation over the conterminous United States dataset (CONUS404). We assessed spatial and temporal error distributions by comparing monthly simulations at the 12-digit hydrologic unit code and regional scale from both models against external benchmarking datasets. Results showed that average annual rainfall across the CONUS was 857 millimeters per year for the period of analysis, with water year 2012 the driest year (729 millimeters) and water year 2019 the wettest year (995 millimeters). Key interannual variability results included the following: (1) the California–Nevada hydrologic region had the highest variability in precipitation and snow accumulation, and (2) the Texas hydrologic region was among hydrologic regions with the highest variability in precipitation. We related interannual variability in precipitation to storage volumes in soil moisture, snow water equivalent, and lakes and reservoirs to highlight areas with little storage and large year-to-year variability in precipitation. These areas included the Southern High Plains, Central High Plains, Texas, Souris–Red–Rainy, Mississippi Embayment, and Midwest regions. Our analysis of groundwater-level data showed that several of these areas overlap aquifers where groundwater levels were considerably lower than historical averages, including the Colorado Plateaus aquifers, the Rio Grande aquifer system, and the Central and Southern regions of the High Plains aquifer. Many of these lowered groundwater levels are continuations of decades-long declines from overpumping that started well before the assessment period. The resulting water budgets and their analyses provide a high-resolution foundational assessment of the mean state and variability of the terrestrial hydrologic cycle across the CONUS and Alaska, Hawaii, and Puerto Rico to support a wide range of water resource management applications.

Professional Paper↗