USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geological Society of America Bulletin”

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.

1,666 records · Page 75Linked to original sources
USGS2026

Detecting snow avalanche activity using infrasound: Hooker Valley, New Zealand

Snow avalanches pose considerable hazards to people and infrastructure in alpine environments. Traditional avalanche monitoring relies on meteorological data and visual observations, which can be limited in scope and timeliness. Infrasound offers a promising complementary monitoring tool by detecting the low-frequency sound waves generated by avalanches. Here, we present infrasound and camera observations during a 50-day field campaign in the Hooker Valley of Aoraki/Mount Cook National Park, New Zealand. Our study detected seven avalanches with the cameras, whereas the infrasound system identified only one of these events, which was the largest and occurred under conditions that likely favoured infrasound propagation. The infrasound system recorded numerous other events not captured by the cameras, indicating the benefit of further investigation to determine their sources. These findings highlight the potential of infrasound technology for detecting avalanches and providing broad spatial coverage, capturing events in areas not monitored by cameras, while also showcasing limitations in infrasound capabilities. The limited detection of smaller avalanches underscores the opportunity for further research to enhance detection capabilities and understand environmental influences such as snow cover and wind noise. Overall, this study emphasises the utility of multidisciplinary monitoring techniques to improve avalanche detection in alpine environments.

Hooker Valley↗
USGS2025

Evidence for nonlocal sediment transport on hillslopes from fault scarp morphology

Hillslope sediment transport processes such as bioturbation, rainsplash, and granular mechanics occur across the entire planet. Yet, it remains uncertain how these small-scale processes act together to shape landscapes. Longstanding hillslope diffusion theory posits that hillslope processes are spatially limited, whereas new concepts of nonlocal sediment transport argue otherwise. However, each theory produces subtly different, but distinct, predictions for the evolution of fault scarps. We use the topographic change of fault scarps to demonstrate that hillslope processes produce nonlocal sediment transport. Analysis of a global compilation of 340 dated single-earthquake scarp profiles reveals a statistically significant ( p < 0.05 ) relationship between scarp age and scarp asymmetry, here defined as the ratio of imaginary to real components of the Fourier transform of absolute slope. Numerical simulations show that nonlocal models predict this relationship, whereas hillslope diffusion models do not. To further investigate this result, we examined the depositional geometry of a well-exposed colluvial wedge along the Wasatch fault in central Utah, United States. Our quantitative comparison between the exposure and numerical simulations reveals better agreement with the nonlocal model. Nonlocal sediment transport theory appears to best capture the physics of how hillslope processes shape fault scarps, yet hillslope diffusion provides a useful approximation in many cases. As the processes that act on fault scarps are nearly identical to those acting on hillslopes, our results provide evidence supporting nonlocality as a generalized model of hillslope sediment transport.

Utah↗
USGS2024

Examining 22 years of ambient seismic wavefield at Mount St. Helens

An increase in seismic activity precedes most volcanic eruptions. Whereas event-based forecasting approaches have been successful, some eruptions remain unanticipated, resulting in casualties and damage. Our study leverages the recent advancements in ambient field seismology. We explore features extracted from continuous ambient fields using traditional methods, for example, peak ground velocity, peak ground acceleration, root mean square, root median square, real-time seismic amplitude measurement, and novel methods (displacement seismic amplitude ratio and spectral width). In addition, we explore unsupervised learning of higher order wavelet features using scattering networks. We find that combining all the methods was necessary to disentangle the effects of seismic sources from structural changes at Mount St. Helens. Although the ambient wavefield-based approach does not yield additional or more significant precursory signals than event-based methods at Mount St. Helens, our study demonstrates that the ambient wavefield provides supplementary information, mainly about structural changes and complements traditional methods. The ambient seismic wavefield offers additional insights into long-lasting processes. We find enhanced wave attenuation correlating with geochemical measurements. We interpret this as ongoing structural changes, such as dome growth or the evolution of the volcanic conduit system. On annual and decadal timescales, we interpret seasonal seismic attenuation in the shallow subsurface as groundwater fluctuations, corroborated by observations at the nearby Spirit Lake level. This multimethod approach at Mount St. Helens sheds light on a volcanic system’s underlying dynamics and structure.

Washington↗
USGS2024

Methods for quantifying interactions between groundwater and surface water

Driven by the need for integrated management of groundwater (GW) and surface water (SW), quantification of GW–SW interactions and associated contaminant transport has become increasingly important. This is due to their substantial impact on water quantity and quality. In this review, we provide an overview of the methods developed over the past several decades to investigate GW–SW interactions. These methods include geophysical, hydrometric, and tracer techniques, as well as various modeling approaches. Different methods reveal valuable information on GW–SW interactions at different scales with their respective advantages and limitations. Interpreting data from these techniques can be challenging due to factors like scale effects, heterogeneous hydrogeological conditions, sediment variability, and complex spatiotemporal connections between GW and SW. To facilitate the selection of appropriate methods for specific sites, we discuss the strengths, weaknesses, and challenges of each technique, and we offer perspectives on knowledge gaps in the current science.

Annual Review of Environment and Resources↗
USGS2025

A spatial analysis of the groundwater emergence flood hazard in Long Island, New York and near coastal areas surrounding Long Island Sound in New York, Connecticut, and Rhode Island

Long Island, New York and near coastal areas surrounding Long Island Sound are densely populated and, like other coastal areas, are susceptible to flooding from several potential sources, including stormwater from precipitation events, tidal flooding and storm surge, and groundwater inundation or groundwater emergence flooding. The latter refers to the intersection of a rising water table with land surface or critical infrastructure. Many studies of flood drivers either neglect or only briefly discuss how shallow groundwater conditions may contribute to or exacerbate flood conditions. As part of a comprehensive study of compound flood hazards in the near coastal areas surrounding Long Island and Long Island Sound, a spatial analysis was completed, in cooperation with the Environmental Protection Agency’s Long Island Sound Study, using available regional datasets to characterize the potential hazard for groundwater emergence flooding. The approximately 3,100 square mile study area was subdivided into 11,407 900-meter by 900-meter (approximately 3,000-feet by 3,000-feet) grid cells, for the purposes of integrating the spatial datasets to calculate and map the groundwater emergence flood hazard. The depth to the water table, hydrologic soil groups, and National Land Cover Database were harmonized to the common grid. A groundwater emergence flood hazard rank was calculated for each grid cell for current average conditions following a set of rules accounting for the depth to the water table and the percent of area within each cell with slow infiltrating soils. A higher sea level position scenario was also calculated for the Long Island part of the study area. The calculated groundwater emergence flood hazard rank was reviewed in concert with the National Land Cover Data Base to identify developed areas and associated infrastructure that may be at risk to groundwater emergence flooding. Study results indicate that the groundwater emergence flood hazard is highest in coastal areas and near surface water where the water table is close to ground surface. Inland areas away from surface water bodies are not likely to be exposed to groundwater emergence flooding. For Long Island, under a scenario with higher sea level position, a greater groundwater emergence flood hazard is calculated in some locations closer to the coast and where land is submerged. Away from the coast and surface-water drainage, the groundwater emergence flood hazard is similar between the current average sea level condition and a higher sea level position scenario.

EarthArXiv↗
USGS2025

Modeling seawater intrusion along the Alabama coastline using physical and machine learning models to evaluate the effects of multiscale natural and anthropogenic stresses

Seawater intrusion threatens groundwater resources in coastal regions, including southern Baldwin County, Alabama, where the freshwater-saltwater interface dynamics remain poorly understood. To address this gap, this study uses combined physics-based and machine-learning models to quantify seawater intrusion caused by natural (storm surges) and anthropogenic (human activities) perturbations. The long short-term memory network and wavelet analysis were used to assess vertical aquifer vulnerabilities, revealing that the shallow part of the Coastal lowlands aquifer system (CL1) in the southern Baldwin County region is more susceptible to sea level rise and groundwater extraction than deeper aquifers. Based on these findings, a cross-sectional numerical model (physics approach) for the CL1 aquifer was developed to evaluate tidal and storm surge effects, using Tropical Storm Claudette (June 2021) as a case study. Results showed that tidal fluctuations had a minimal impact on the saltwater-freshwater interface location, whereas storm surges caused substantial inland movement, with effects lasting for nine months. The steady-state version of the three-dimensional (3D) physical model predicted seawater intrusion across the entire area, and convolutional neural network-based modeling further validated the model results. The 3D physical model was also applied to a smaller area to assess human impact on the saltwater interface due to two groundwater pumping scenarios (± 50% of the baseline pumping rate). Results revealed that a 50% increase in groundwater withdrawals caused seawater to advance ~ 320 m inland, whereas a 50% reduction led to a ~ 270-meter retreat. This study highlights the vulnerability of Alabama’s shallow coastal aquifers to seawater intrusion due to storm surges and human activities, and demonstrates that combining physics-based models with machine learning approaches can improve groundwater predictions, though its accuracy depends on the availability of site-specific data.

Alabama↗
USGS2025

Bayesian mapping of regionally grouped, sparse, univariate earth science data

Some earth science data are naturally grouped by region, and it is often desirable to map these data by region. However, if there are only a few samples within each region, then the map should be smoothed in an appropriate way to mitigate the problems that arise from having only a few samples. A smoothing algorithm based on a Bayesian hierarchical model is developed and presented in this report. This algorithm has several features that make it especially suitable for mapping earth science data: it can account for measurements that are censored, it can process multiple datasets with different measurement errors and different censoring thresholds, and it can calculate the uncertainty in any statistic that is mapped. The algorithm is demonstrated by mapping gold concentrations that are measured in streambed sediments in the Taylor Mountains quadrangle in southwestern Alaska.

Alaska↗
USGS2025

Speleothem evidence for Late Miocene extreme Arctic amplification – An analogue for near-future anthropogenic climate change?

The Miocene provides an excellent climatic analogue for near-future runaway anthropogenic warming, with atmospheric CO 2 concentrations and global average temperatures similar to those projected for the coming century under extreme-emissions scenarios. However, the magnitude of Miocene Arctic warming remains unclear due to the scarcity of reliable proxy data. Here we use stable oxygen isotope and trace element analyses, alongside clumped isotope and fluid inclusion palaeothermometry of speleothems to reconstruct palaeo-environmental conditions near the Siberian Arctic coast during the Tortonian (8.68 ± 0.09 Ma ). Stable oxygen isotope records suggest warmer-than-present temperatures. This is supported by temperature estimates based on clumped isotopes and fluid inclusions giving mean annual air temperatures between + 6.6 and + 11.1 °C , compared with − 12.3 °C today. Trace elements records reveal a highly seasonal hydrological environment. Our estimate of > 18 °C of Arctic warming supports the wider consensus of a warmer-than-present Miocene and provides a rare palaeo-analogue for future Arctic amplification under high-emissions scenarios. The reconstructed increase in mean surface temperature far exceeds temperatures projected in fully coupled global climate models, even under extreme-emissions scenarios. Given that climate models have consistently underestimated the extent of recent Arctic amplification, our proxy data suggest Arctic warming may exceed current projections.

Lena River delta↗
USGS2025

Did they feel it? Legacy maroseismic data illuminates an engimatic 20th century earthquake

The challenges and the importance of preserving legacy instrumental records of earthquakes are now well-recognized (e.g., Richards & Hellweg, 2020, https://doi.org/10.1785/0220200053 ). Seismologists may not be aware of parallel challenges and opportunities with legacy macroseismic data for earthquakes in the United States. For much of the 20th century, macroseismic data were collected by a series of U.S. government agencies using a standard questionnaire distributed on postcards. Published summaries of postcards provide macroseismic data akin to modern Did You Feel It? questionnaire responses. In this paper we focus on the M 6.5 Fickle Hill, California earthquake, on 21 December 1954 (Hellweg et al., 2025) as a proof-of-concept, illustrating the potential of what we dub Did They Feel It? (DTFI) data to improve our understanding of significant 20th century U.S. earthquakes for which instrumental data are sparse. Legacy macroseismic data interpreted following modern conventions can potentially constrain traditional ShakeMaps at a level of detail and accuracy that in some respects rival maps for modern earthquakes. The updated ShakeMap for the 1954 Fickle Hill earthquake, also drawing from recently published media and first-person accounts, supports the location, depth, and stress drop value estimated from available instrumental data (Hellweg et al., 2025).

California↗
USGS2024

Taking heat (downstream): Simulating groundwater and thermal equilibrium controls on annual paired air–water temperature signal transport in headwater streams

Headwater stream temperature often exhibits spatial variation at the kilometer-scale, but the relative importance of the underlying hydrogeological processes and riverine perturbations remains poorly understood. In this study, we investigated the relative importance of groundwater (GW) and other processes on downstream annual stream temperature signal characteristics using deterministic heat budget model (HFLUX) scenarios within an idealized stream reach representative of mountainous forested conditions. We summarized annual stream thermal regimes from the relationship of paired sinusoidal air and water temperature signals (amplitude ratio, phase lag, and mean ratio). Results showed that downstream changes in annual temperature depended on the thermal gradient between water and the hypothetical equilibrium temperature (where all heat fluxes sum to zero). GW inflow, riparian shading, and the boundary input signal were the most significant factors affecting downstream annual water temperature signals, while flow volume and channel dimensions impacted how quickly annual temperature signals changed. Effects of GW were dominated by advective rather than conductive heat exchange processes, but conduction played a larger role when GW input was more spatially diffuse. Our results indicated several mechanisms by which local processes may affect stream thermal resilience to disturbances and can help guide management of wildfire and climate change.

Journal of Hydrology↗
USGS2026

Suitability of ground-motion models for seismic hazard assessment in Puerto Rico and the U.S. Virgin Islands

We perform linear mixed-effects analyses with a ground-motion dataset to evaluate how well ground-motion models (GMMs) fit active crustal, subduction interface, and subduction intraslab earthquakes in Puerto Rico and the US Virgin Islands. Most of the GMMs reproduce the spatial variation in peak ground motions with earthquake magnitude and rupture distance but predict ground motions 0.3–1.0 natural log units (35%–270%) greater than observations. Two GMMs developed for Puerto Rico that are based on ground-motion records from mostly magnitude 4–5 earthquakes do not perform as well. We attribute the period-dependent overprediction in observed ground motions to differences between the observed site response and the linear site response in the GMMs. Consequently, we developed region-specific GMMs by adjusting the period-dependent linear site response coefficients and period-dependent constant coefficients to remove most of the bias between predicted ground motions and observations. For the analysis, we compile ground-motion records and process waveforms to build a dataset with 10,127 records at 72 stations from 849 magnitude 4.0–6.4 earthquakes between 1 January 2006 and 31 March 2024. The earthquakes include active crustal, subduction interface, and subduction intraslab events. We evaluate the GMMs using the time-averaged shear wave speed in the top 30 m ( ), which we compile from site surveys and proxy values computed from horizontal to vertical spectral ratios. Site terms exhibit strong consistency across GMMs and crustal and subduction earthquakes, indicating that the linear mixed-effects analysis successfully isolates the effects of local site response. The event terms show little spatial correlation and more substantial variability than in other regions, which we attribute to uncertainties in the earthquake magnitudes. This analysis guides the selection of the GMMs for the 2025 update of the National Seismic Hazard Model for Puerto Rico and the US Virgin Islands.

Puerto Rico, U.S. Virgin Islands↗
USGS2024

The 17 January 1994 Northridge, California, earthquake: A retrospective analysis

The 17 January 1994 Northridge, California, earthquake was a watershed event, with far-reaching societal and scientific impacts. The earthquake, which occurred in the early days of both broadband seismic networks and the Internet, spurred advances in seismic monitoring, real-time systems, and development of data products. Motivated by the 30th anniversary of the earthquake, we present a brief retrospective of the earthquake and its impact, and reconsider both ground motions and the aftershock distribution using modern tools and best-available data. With improvements in instrumentation and analysis methodology, recent earthquakes continue to reveal increasing complexity of ground motions, fault systems, and earthquake ruptures. Even in the absence of data from state-of-the art instrumentation, a retrospective consideration of ground motion data from the Northridge earthquake reveals complexities beyond what could be characterized (and modeled) thirty years ago. Aftershock relocations for both the 1971 Sylmar and 1994 Northridge earthquakes also reveal an updated view of fault complexity. Our study does provide a cautionary tale regarding legacy data sets and research results that are not easily accessible, which can result in discrepancies between catalog data and products from best-available science. We also briefly describe outreach products produced as part of the anniversary commemoration.

California↗
USGS2026

Cultural benefits knowledge and social learning in Elwha River dam removal and ecosystem decision-making

In a retrospective case study of the seminal Elwha River dam removal in Washington State, USA, we explore pathways through which ecosystem services (ES) knowledge can inform environmental decision-making. We specifically investigate how knowledge of cultural benefits linked to the Elwha River ecosystem was made available. Through interviews and historical document review we identify diverse forms of cultural benefits knowledge (i.e., how individuals or groups know the contributions of an ecosystem to their identities, experiences, and capabilities), including both knowledge in the form of practices – embodied or enacted forms of knowledge, and knowledge in the form of products – quantitative or written documentation. Recognizing and learning from the diverse forms in which benefits knowledges are conveyed can support more accurate and equitable consideration of the diverse values, well-beings, and human-nature relationships linked to ecosystems. With this case study, we illuminate pathways for meaningful consideration of multiple knowledge forms, and the plural ecosystem values they convey, with relevance for environmental decision processes in the Puget Sound and beyond.

Washington↗
USGS2026

GST-1: A high-resolution global sediment thickness model

Global Sediment Thickness 1 (GST-1) is a high-resolution sedimentary thickness model calculated on a 0.125° x 0.125° grid. It modifies the sediment thickness of the 1° x 1° Earth Crustal Model 1 (ECM1) by means of 3D inversions of free air gravity anomalies. GST-1 is calculated by performing structural inversions on high-density contrasts across two crustal boundaries: the sediment – basement interface and the crystalline crust – upper mantle interface. The inversions are calculated in each of ten overlapping 3D models that span the globe. These ten models are merged to obtain the GST-1 global model, providing an eight-fold increase in lateral spatial resolution in comparison with ECM1 and CRUST 1.0. Our sediment thickness model exploits the nearly continuous sampling of gravity data when compared to the irregular, sparse sampling of seismic refraction data. Sediment thickness values in GST-1 are in excellent agreement with independently derived cross sections from well-studied sedimentary basins, and within expected resolution limits of seismic refraction data. GST-1 offers a robust, high resolution global model of sedimentary thickness to support studies of sedimentary basins.

Tectonophysics↗
USGS2025

Haunted Summerville: Ghostly lights or earthquake lights?

Among the colorful local lore in the Charleston, South Carolina, area, are a number of ghost stories, shared not only over campfires but also in published books. Among the most well-known of the stories is the tale of the Summerville Light. Local lore holds that a strange light sometimes seen in a remote area is a lantern carried by the ghost of a local woman who once waited hours for her husband, who turned out to have been decapitated earlier that day in a train accident (DePoppe, 2023). Extant sources suggest the ghost stories began to circulate in the 1950s to 1960s. So pervasive was the lore that (Old) Sheep Island Road became known among local residents as Light Road, with a local stretch of road known today as Old Light Road. Reviewing the location where the lights appear as well as the nature of accounts, I suggest that many if not all of the anecdotal observations can be most readily attributed to natural phenomena, including earthquake lights from earthquakes that were too small to be felt. Accounts of lights near Summerville cluster in proximity to the generally accepted epicenter of the 1886 Charleston, South Carolina, earthquake, where foreshocks to the 1886 mainshock were apparently concentrated, and within a few kilometers of three M3.5 – 4.4 earthquakes in 1959 and 1960.

South Carolina↗
USGS2026

Hydrological control on soil redox condition and carbon loss of coastal wetland under sea-kevel rise

Coastal wetlands are critical carbon sinks with their biogeochemical and ecological functioning shaped by dynamic hydrological conditions that are increasingly influenced by climate change. A key unresolved question is how hydrologic flow, vegetation response, and rising sea levels interact to regulate soil redox condition and carbon loss in coastal wetlands. Using a field-tested hydrological–biogeochemical–ecological modeling framework, we reveal how the interplay between terrestrial groundwater discharge and tidal fluctuations generates complex groundwater flow patterns at the terrestrial–aquatic interface, and how these patterns modulate soil redox conditions, in turn influencing soil organic matter decomposition and carbon loss. Notably, rising sea levels suppress soil CO 2 emissions while reducing lateral dissolved carbon losses, thereby enhancing litter carbon sequestration under anoxic conditions. As vegetation responds to sea-level rise and carbon inputs diminish, litter carbon subsequently declines. These findings underscore a critical hydrological control on carbon cycling, advancing our understanding of coastal ecosystem resilience in a warming world.

Louisiana↗
USGS2025

ARCHI: A new R package for automated imputation of regionally correlated hydrologic records

Missing data in hydrological records can limit resource assessment, process understanding, and predictive modeling. Here, we present ARCHI (Automated Regional Correlation Analysis for Hydrologic Record Imputation), a new, open-source software package in R designed to aggregate, impute, cluster, and visualize regionally correlated hydrologic records. ARCHI imputes missing data in “target” records by linear regression using more complete “reference” records as predictors. Automated imputation is implemented using a novel, iterative algorithm that allows each site to be considered a target or reference for regression, growing the pool of complete references with each imputed record until viable gap-filling ceases. Users can limit artifacts from spurious correlations by specifying model-acceptance criteria and applying geospatial, correlation, and group-based filters to control reference selection. ARCHI provides additional functions for visualizing results, clustering records with similar correlation structures, evaluating holdout data, and interactive parameterization with an accessible and intuitive graphical user interface (GUI). This methods brief provides an overview of the ARCHI package, modeling guidelines, and benchmarking on two regional groundwater-level datasets from the Central Valley, CA and Long Island, NY. We evaluate ARCHI alongside widely used multivariate imputation software to highlight and contextualize its computational efficiency, imputation accuracy, and model transparency when applied to large, groundwater-level datasets.

California, New York↗