Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research: Solid Earth”

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,665 records · Page 5Linked to original sources
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↗
USGS2026

Emergent headwater types for watershed monitoring, research, and planning

Headwaters—that is headwater streams and the watersheds they drain—strongly influence water quality, ecosystem services, and hydrologic connectivity, yet remain poorly mapped, monitored, and understood at large spatial scales. Using available spatial data to classify headwater “types” could therefore enhance the utility of sparse monitoring data to advance research and watershed planning—and help update literature-based perceptual models of headwaters as steep, forested, and minimally disturbed watersheds. Here, we analyze headwaters across the contiguous United States to quantify how their attributes (climate, topography, land cover, soils) vary regionally and differ from their downstream watersheds. Our analysis (a) identified seven dominant headwater types with distinct combinations of biophysical and climatic attributes and (b) revealed systematic contrasts with downstream watersheds, with strongest differences in attributes representing potential flowpath depth. Our findings offer a transferable framework for classifying, modeling, and managing these critical yet understudied components of the hydrologic network.

conterminous United States↗
USGS2025

Reaction kinetics and accelerant effects of sulfides in early mature hydrocarbon generation using hydrous pyrolysis

Hydrocarbon generation in organic-rich sediments is influenced by the molecular organic composition and relative abundance of associated minerals. Certain mineral-derived elements act as catalysts and reaction intermediaries, facilitating early-stage hydrocarbon formation in potential source rocks. This study investigated the role of sulfur contributed from pyrite as an accelerant in thermal reaction, focusing on its effects on early maturation and consequent hydrocarbon generation from gilsonite (low-sulfur solid petroleum). Hydrous pyrolysis (HP) experiments were conducted on mixtures of gilsonite and pyrite in varying ratios (1:0.1, 1:0.5, 1:1, 1:2, and 1:10 w/w gilsonite:pyrite) at 320, 350, and 370 °C for 72 h. Untreated and thermally altered residues were analyzed using solid bitumen reflectance (BR o , %), total organic carbon (TOC) content, programmed temperature pyrolysis, scanning electron microscopy with energy-dispersive spectroscopy (SEM-EDS), and X-ray diffraction (XRD) to evaluate the potential accelerant role of pyritic sulfur in hydrocarbon formation. The results show HP residues at 320 and 350 °C with greater pyrite concentrations had higher BR o , while reflectance values were similar in the 370 °C residues, regardless of pyrite concentration, suggesting enhanced reaction at lower thermal conditions. Increasing pyrite content systematically decreased hydrogen index (HI) values while increasing the transformation ratio (TR) and production index (PI), indicating enhanced conversion of organic matter to hydrocarbons with increasing pyrite concentrations. Gas yields increased with pyrite addition, particularly at 350 °C, confirming secondary cracking effects. However, gas production stabilized or declined at higher pyrite loadings (1:10), suggesting alternative reaction pathways such as coke formation. Our data indicate the presence of pyrite lowers the activation energy for thermal cracking, shifting peak experimental hydrocarbon generation temperatures downward by 20–30 °C, with the most pronounced accelerant effects observed at moderate pyrite concentrations (1:0.5 and 1:1). The thermodynamic framework reveals that pyrite stability is influenced by experimental conditions, with pyrrhotite formation favored in the presence of gilsonite due to reduced oxygen fugacity. Pyrite transformation to pyrrhotite, as observed through XRD, SEM-EDS, and predicted by thermodynamic data, further supports the accelerant role of S, as pyrrhotite exhibits a higher hydrogen transfer potential, promoting early oil generation. These findings highlight the importance of pyrite in modulating hydrocarbon generation pathways in organic-rich systems.

Journal of Analytical and Applied Pyrolysis↗
USGS2024

Groundwater-Surface water interactions research: Past trends and future directions

Interactions between groundwater and surface water sustain groundwater-dependent ecosystems and regulate river temperature and biogeochemical cycles, amongst many other processes. These interactions occur in freshwater environments including rivers, springs, lakes, and wetlands, and in coastal environments via tidal pumping, submarine groundwater discharge, and seawater intrusion. Here, we explore groundwater-surface water interactions research using bibliometric analyses of titles, abstracts, and keywords from 20,275 journal papers published between 1970 and 2023 extracted from Scopus. Analyses show that research into groundwater-surface water interactions is highly multi-disciplinary, with growing contributions from the social and biological sciences. The number of groundwater-surface water interactions papers is rapidly increasing with over 1200 papers published per year since 2020. Drawing on our data-driven approach and expert knowledge, we synthesise current research trends and identify critical future research directions. Despite the thousands of papers on groundwater-surface water interactions, important processes are still difficult to quantify or predict at meaningful spatial scales to inform water-resources management. We see benefits in future groundwater-surface water interactions research focusing on: (1) using new technologies including internet-of-things-based sensors, uncrewed vehicles, and remote-sensing approaches for data collection to inform groundwater-surface water interactions at large scales, (2) seeking approaches to upscale site-specific findings to better inform management, and (3) continuing the movement towards multi-disciplinary investigations to better inform the understanding of groundwater-surface water interactions and processes that will enable better management outcomes.

Journal of Hydrology↗
USGS2025

Mechanics and statistics of postseismic shaking

Analysis of two weeks of continuous post-seismic shaking after the 2019 M7.1 Ridgecrest, CA earthquake sequence using 4 nearby borehole seismometers reveals that continuous ground motions decay as Omori’s law in time and follow the Gutenberg-Richter distribution in logarithmic amplitude. The measured temporal decay in amplitudes agrees with predictions of the rate-and-state framework and indicates shaking amplitudes are proportional to the velocity of afterslip. Our ground motion-based statistical framework provides a basis to forecast shaking intensity in the minutes to hours after a large earthquake.

California↗
USGS2025

New constraints on location and timing of the Great Lakes tectonic zone, central Upper Peninsula, Michigan, USA

The Great Lakes tectonic zone (GLTZ) forms the boundary between the Wawa–Abitibi and Minnesota River Valley subprovinces within the Archean Superior Province. The GLTZ is concealed for all of its 1100 km length, except for a segment in the central Upper Peninsula of Michigan. There, it is exposed as a northwest-striking mylonite zone along a 11 km segment, extending to the onlap of Paleozoic rocks to the east. Farther east, its location has been unknown. Here, we use aeromagnetic and gravity data to develop interpretations of the expression of the GLTZ and to define its extent under cover. Aeromagnetic gradients over the mylonite zone are interpreted to be produced by structurally juxtaposed rocks with varying magnetizations. Gravity data show a regional gradient along the GLTZ, produced by the juxtaposition of a dense greenstone belt on the north against lower-density gneisses and granites on the south. The GLTZ is interpreted to extend ∼55 km under cover to the east. The GLTZ is terminated on the east by the buried eastern arm of the ca. 1100 Ma Midcontinent Rift. An undeformed granitic dike that cuts the mylonitic foliation produces a U–Pb apatite age of 2523 ± 33 Ma, implying no major post-Archean shearing occurred, and is at odds with previous interpretations of major Proterozoic reactivation. A granite intrusion in the Minnesota River Valley subprovince produces a Pb–Pb zircon age of 2606.9 ± 3.6/7.4 Ma. This suggests that magmatism related to the Sacred Heart orogeny, previously known in Minnesota, extended to Michigan.

Michigan↗
USGS2025

S/P amplitude ratios with Distributed Acoustic Sensing and application to earthquake focal mechanisms

Distributed acoustic sensing (DAS), which transforms a fiber optic cable into an array of high frequency strainmeters, has the potential to help us characterize earthquakes with a dense sampling of measurements. While earthquake focal mechanisms are frequently determined using P -wave polarities and S / P amplitude ratios with inertial seismometers, the dense sampling of DAS over potentially large portions of the focal sphere can aid our solutions. Here, we consider ∼200 regional earthquakes and thousands of S / P measurements on a DAS cable colocated with a network of inertial seismometers near Arcata, California. We demonstrate the S / P ratio measurements made on DAS are similar to those made on the inertial seismometers and can be used to constrain focal mechanism solutions.

California↗
USGS2024

GRAPES: Earthquake early warning by passing seismic vectors through the grapevine

Estimating an earthquake's magnitude and location may not be necessary to predict shaking in real time; instead, wavefield-based approaches predict shaking with few assumptions about the seismic source. Here, we introduce GRAph Prediction of Earthquake Shaking (GRAPES), a deep learning model trained to characterize and propagate earthquake shaking across a seismic network. We show that GRAPES’ internal activations, which we call “seismic vectors”, correspond to the arrival of distinct seismic phases. GRAPES builds upon recent deep learning models applied to earthquake early warning by allowing for continuous ground motion prediction with seismic networks of all sizes. While trained on earthquakes recorded in Japan, we show that GRAPES, without modification, outperforms the ShakeAlert earthquake early warning system on the 2019 M7.1 Ridgecrest, CA earthquake.

Shimane/HiroshimaPrefectures↗
USGS2025

Airborne geophysical analysis to decipher salinization for coastal Louisiana

Coastal Louisiana is known for saltwater intrusion that threatens wetlands, aquifers, and rivers. However, the extent of saltwater intrusion is not well understood. This study develops an innovative framework with airborne electromagnetic (AEM) data to map chloride concentration distributions for wetlands in the Mississippi River deltaic plain and Chenier plain as well as for the Mississippi River Valley alluvial aquifer (MRVA) and Chicot aquifer. Moreover, the framework maps chloride concentrations along the Mississippi River and Atchafalaya River. Key components in the framework include the establishment of resistivity-to-chloride concentration transformation, 3D resistivity architecture building through geostatistics, and the employment of a lithologic model. The transformation functions correlate AEM resistivity data with porewater salinity measurements and groundwater and river chloride samples. The results show that AEM data reliably infers soil water chloride concentrations and correlates well with the distribution of various marsh types. AEM data reveals extensive saltwater presence at depth and near the coast, originating from salt domes and the Gulf of Mexico, respectively. The saltwater upconing pattern in the Chicot aquifer is likely due to excessive groundwater withdrawals. The AEM data also confirms a distinct tongue of saltwater intruding into the Atchafalaya Basin from the Gulf. The AEM data helps to identify faults that are obscured or eroded at the surface, which appear as leaky barriers in the subsurface where dramatic changes in chloride concentration are apparent. Finally, this study uses the AEM data to infer the presence of an extensive seawater wedge in the Mississippi River and Atchafalaya River.

Louisiana↗
USGS2025

The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought

Quantitative understanding is lacking on how the depth of active groundwater circulation in bedrock affects mountain streamflow response to a multi‐year drought. We use an integrated hydrological model to explore the sensitivity of a variety of streamflow metrics to bedrock circulation depth and porosity under a plausible extreme drought scenario lasting up to 5 years. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. With drought, a deeper circulation system with higher drainable porosity more effectively buffers minimum flow and significantly limits perennial stream loss in comparison to a shallow circulation system. Streamflow buffering is accomplished through extensive groundwater storage loss. However, deeper circulation systems experience prolonged recovery from drought in comparison to storage‐limited shallow systems. Research highlights the importance of characterizing the deeper bedrock hydrogeology in mountainous watersheds to better understand and predict drought impacts on stream ecosystem health and water resource sustainability.

Colorado↗
USGS2024

A robust quantitative method to distinguish runoff-generated debris flows from floods

Debris flows and floods generated by rainfall runoff occur in rocky mountainous landscapes and burned steeplands. Flow type is commonly identified post-event through interpretation of depositional structures, but these may be poorly preserved or misinterpreted. Prior research indicates that discharge magnitude is commonly amplified in debris flows relative to floods due to volumetric bulking and increased frictional resistance. Here, we use this flow amplification to develop a metric ( Q* ) to separate debris flows from floods based on the ratio of observed peak discharge to the theoretical maximum water discharge from rainfall runoff. We compile 642 observations of floods and debris flows and demonstrate that Q* distinguishes flow type to ∼92% accuracy. Q* allows for accurate identification of debris flows through simple channel cross-section surveys rather than through qualitative interpretation of deposits, and therefore should increase the performance of models and engineered structures that require accurate flow-type observations.

Geophysical Research Letters↗
USGS2024

Detrital zircons and the magmatic history of Viti Levu, Fiji

We integrate the existing detrital zircon data from multiple modern river sediment samples on Viti Levu, Fiji, with the most current available geological and topographic mapping of the respective river drainage basins to compare detrital populations with potential bedrock sources. The temporal and spatial variations in zircon geochemistry supplement what is known from igneous rocks and confirm the petrological differences between plutonic and volcanic rocks from the Eocene to early Oligocene (Yavuna age, >30 Ma), middle Oligocene to middle Miocene (Wainimala age, 30–12.5 Ma), late Miocene (Colo age, 12.5–6.5 Ma) and latest Miocene (Namosi age, 6.5–5 Ma). The >30 Ma Yavuna-age zircons are restricted to areas that drain the previously mapped Yavuna Group. The 30–12.5 Ma zircons are found across central Viti Levu from west to east, and the 30–15 Ma zircons have distinctively low U/Yb and high Dy/Yb ratios. They are the best radiometric evidence of widespread early to middle Miocene arc magmatism in Fiji that was relatively U-poor. Peak deconvolution of the Colo age zircons from individual basins suggests the following ages for undated or poorly dated plutons from central Viti Levu. The large Mavuvu pluton is probably composed of multiple intrusions in the 12–10 Ma range, the Waiqa pluton is probably ca 10 Ma, and the Noikoro pluton is probably ca 9 Ma. There are zircons from unknown plutonic or volcanic sources between 8 and 7 Ma in western Viti Levu that have distinct Eu/Eu* ratios. We attribute the highest U/Yb ratios in some Colo age zircons to crustal anatexis. Namosi-age zircons are abundant in the Medrausucu Group and can be found scattered across Viti Levu.

Viti Levu↗
USGS2026

Decadal shifts in groundwater age detected by environmental tracers across California, USA

Groundwater age offers important insight into recharge, storage, and contamination risk. Although models predict age changes can be driven by pumping and climate variability, direct observational evidence remains limited. Here, we analyzed paired environmental tracer suites (tritium, carbon-14, and tritiogenic helium-3) collected a decade apart from 268 wells across California to assess the prevalence of groundwater age transience. Travel-time distribution models and statistical tests indicated age transience at 29% of sites, occurring most often in agricultural regions, such as the San Joaquin Valley and Southern Coast Ranges, where large carbon-14 changes coincided with substantial nitrate and chloride shifts. Sites with tritiogenic helium-3 data showed more frequent age transience, underscoring the value of multi-tracer data sets. These results provide the first regional evidence of widespread groundwater age change and a method for detecting changing water balances with implications for groundwater sustainability and water quality.

California↗
USGS2025

Anomalous shear stress variation in wet granular medium: Implications for landslide lateral faults

Landslide assessments typically focus on the mechanical properties of the basal shear zone, but lateral faults are frequently overlooked, possibly due to their lower normal stresses and variably saturated conditions. Using double-cylinder shear experiments on wet granular systems as analogs for landslide lateral faults, we observe anomalous shear stress variations with fluid volume fractions, defying an expected unimodal relationship associated with capillary cohesion. At low fluid volume fractions, shear strength weakens as the wet grain assembly experiences reduced lateral pressure and increased boundary slip. This boundary slip subsequently vanishes, with an abrupt strengthening due to the dilation of the grain assembly against fluid surface tension as saturation approaches. Strike-slip motion and confinement in this system explain the strength anomaly, highlighting a critical role of lateral faults in landslide stability, particularly in cases where dynamics cannot be adequately explained by monitored pore-water pressure or basal friction.

Geophysical Research Letters↗
USGS2024

Time varying crustal anisotropy at Whakaari/White Island volcano

Whakaari/White Island has been the most active New Zealand volcano in the 21st century, producing small phreatic and phreatomagmatic eruptions, which are hard to predict. The most recent eruption occurred in 2019, tragically claiming the lives of 22 individuals and causing numerous injuries. We employed shear-wave splitting analyses to investigate variations in anisotropy during quiescence, unrest, and an eruption between 2018 and 2020. We examined spatial and temporal variations in 7972 shear-wave splitting measurements and derived 1402 V p /IV s ratio measurements. Comparing shear-wave splitting parameters from similar earthquake paths across different times indicates that the observed temporal changes are unlikely to result from variations in earthquake paths through media with spatial variability. Instead, these changes may stem from variations in anisotropy over time, likely caused by changes in crack alignment due to stress or varying fluid content.

Whakaari/White Island volcano↗
USGS2025

Surface-wave relocation and characterization of the October 2023 tsunamigenic seismic unrest near Sofugan volcano, Izu Islands, Japan

A moderate-magnitude earthquake swarm occurred in the remote Izu Islands region of Japan between October 1 and 8, 2023. The swarm included 151 shallow earthquakes cataloged by the U.S. Geological Survey, which notably included a roughly 2.5-hr episode of 15 successive magnitude (M) < 5.5 earthquakes. Origin times were coincident with regionally recorded tsunami waves, but tsunamigenesis for moderate-magnitude earthquakes is uncommon, indicating that volcanic activity generated the ocean displacements. Leveraging a surface-wave relative relocation approach, we estimate precise epicentroid locations for the remote swarm. Final epicentroids and caldera analogs indicate a three-stage model to explain swarm activity: (a) caldera pressurization due to magma intrusion, (b) depressurization via dike propagation away from the caldera, and (c) eruption corresponding with caldera reactivation either by collapse or additional intrusion.

Sofugan volcano, Izu Islands↗
USGS2026

Asynchronous landslide seasonality across the United States

Mid-range landslide outlooks can facilitate weather-related landslide preparedness and disaster response planning, but seasonal landslide activity remains poorly quantified at continental scales. Leveraging >55,000 reported landslides from across the United States (U.S.), we used circular statistics to quantify landslide seasonality in 67 National Weather Service County Warning Areas (CWAs). We found regional differences in landslide season timing and duration, with transitions between domains variably corresponding to climate class or river basin. We assessed differences in seasonality by movement type for slides, flows, and falls, detecting apparent, but uncertain, differences between slide and fall seasonalities in 27 of 35 (77%) of CWAs with both types reported. In the Pacific Northwest, where long records exist, we found a credible shift toward a later mean landslide season in western Washington from 1990 to 2020, but no trend in western Oregon. Our results can provide emergency planners a resource to assess seasonal landslide probability nationwide.

Geophysical Research Letters↗
USGS2025

Multi-scale predictors of Northern Long-eared Bat (Myotis septentrionalis) occupancy in the United States

Historically, Myotis septentrionalis (Northern Long eared Bat) was among the most common forest-interior species in North America. Largely due to high mortality from white-nose syndrome, this species has experienced severe population declines across its range. To create an updated species distribution map representing summer occupancy probabilities from 2017 to 2022, we integrated stationary acoustic data with live-capture data from the database of the North American Bat Monitoring Program into a multi-scale, multi-method occupancy modeling framework. Our results provide data-driven predictions with quantified uncertainty for summer occupancy probabilities for Northern Long-eared Bats at 2 spatial scales across the range of the species, while also accounting for inherent observation biases (e.g., imperfect detection).

Journal of North American Bat Research↗