USGS Science⌕ Search

SEARCH · USGS Science

Results for “Mining, Metallurgy and Exploration”

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,691 records · Page 48Linked to original sources
USGS2025

Groundwater salinity: Applying the specific conductance and water type proxy

Groundwater is increasingly needed for water supplies but may have limited utility in some locations because of its salinity. Salinity, often expressed as total dissolved solid (TDS), is frequently estimated using specific conductance (SC) measurements. However, the commonly used proxy (0.65 multiplied by SC to indicate TDS, common in many handheld meters) can result in inaccurate TDS estimates. First, the TDS–SC relationship is not linear over the entire concentration range of groundwater. Furthermore, the TDS (and salinity)–SC relationships vary substantially depending on the major-ion composition. Here we develop a proxy method utilizing SC and major-ion water type to estimate TDS and salinity specifically for groundwaters. Compared to most surface waters, groundwater tends to have a wider range of salinity (fresh to highly saline) and higher concentrations of bedrock-derived solutes such as carbonate ions, silica, and many other ions. The dataset used to develop the proxies includes water chemistry data from 149,059 discrete groundwater samples. The groundwater proxies, which employ nonlinear log–log relations, utilize five water types (HCO 3 , Cl, Ca-Mg-SO 4 , Na-K-SO 4 , and mixed waters), are accurate (median percent difference between TDS and salinity determined using the proxy compared to discrete measurements was <±0.8%) over a wide range of SC (up to 200 mS/cm), rapid, cost-effective, and can be measured on-site.

Groundwater↗
USGS2024

Seismic techniques and suggested instrumentation to monitor volcanoes

Introduction Changes in the pressure or location of magma can stress or break surrounding rocks and trigger flow of nearby waters and gases, causing seismic signals, such as discrete earthquakes and tremor. These phenomena are types of seismic unrest that commonly precede eruption and can be used to forecast volcanic activity. Mass movements at the surface, including avalanches, debris flows, and lahars, may also generate seismic signals that are specifically addressed in chapter H , this volume (Thelen and others, 2024). Our focus in this chapter is to determine the levels of instrumentation recommended to produce high-quality, well-constrained seismic observations important for early warning of impending eruptions, detecting changes in ongoing eruptions, and characterizing other hazardous volcanic events. There are emerging techniques and new types of instrumentation, such as distributed acoustic sensing or rotational seismometers, that we do not consider here. These types of instrumentation show promise for monitoring but still require maturation before being considered more generally in volcano monitoring. Most of the capabilities mentioned below are universal for all types of volcanic systems, although some are best applied to stratovolcanoes with an apical single vent. In some settings, such as calderas or shield volcanoes, we must broaden coverage to include multiple possible storage regions or vent locations. As an example, Thelen (2014) discretized the long rift zones of shield volcanoes in Hawaiʻi as a set of evenly spaced “vents.” In this construct, each vent comes with recommendations, and several thousand network configurations were simulated to assess the effect on network quality levels and to determine the most efficient network design. The same process could be applied in a caldera setting or a volcanic field, where an evenly spaced grid of potential vents is considered. Localized recommendations for each unique system are beyond the scope of this report and left up to local experts to assess based on the conditions, restrictions, and requirements of each volcano.

Scientific Investigations Report↗
USGS2026

Status and understanding of groundwater quality in the San Joaquin Valley Kern County subbasin domestic-supply aquifer study unit, 2022—California GAMA Priority Basin Project

The quality of water accessed by domestic wells (here referred to as domestic groundwater resources) in the San Joaquin Valley Kern County subbasin (basin number 5-022.14) was assessed as part of the California Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project (GAMA-PBP), in cooperation with the California State Water Resources Control Board. Kern County is at the southern end of the San Joaquin Valley in California, and about 30,000 residents are estimated to use privately owned domestic wells for drinking water. Domestic wells typically draw from shallower parts of the aquifer system than public-supply wells and can be more vulnerable to effects from surface activities. Kern County is host to a highly productive agricultural industry, with Bakersfield as the main urban center. The Kern River runs through Bakersfield from the southern Sierra Nevada and intersects the Kern Water Bank, one of the largest groundwater banking operations in California, at the Kern River Intertie. The section of the Kern River running through the Kern Water Bank is dry most years. Kern County also encompasses some of the most productive oil and gas basins in California, with extensive underground and surface disposal of oil-field wastewater. This study was based on data collected from 33 sites sampled by the U.S. Geological Survey for the GAMA-PBP in 2022. To provide context for the water quality assessment, measured concentrations were compared to regulatory and non-regulatory health-based and aesthetic benchmarks. A grid-based method was used to estimate the proportions of the groundwater resources used for domestic-supply wells that have water-quality constituents below (low relative concentration), approaching (moderate relative concentration), or above (high relative concentration) benchmark concentrations. At least one measured constituent with a regulatory benchmark was categorized as having a high relative concentration in 72 percent of the aquifer area used for domestic groundwater resources. Inorganic constituents were detected at high concentrations in 45 percent of the domestic groundwater resources, and the constituents detected above regulatory benchmarks were arsenic, nitrate, and uranium. At least one organic constituent was detected at high concentrations in 41 percent of the domestic groundwater resources, and the constituents exceeding regulatory benchmarks were the fumigants 1,2,3-trichloropropane (1,2,3-TCP), 1,2-dibromo-3-chloropropane (dibromochloropropane [DBCP]), 1,2-dibromoethane (EDB), and the per-and polyfluoroalkyl substance (PFAS) perfluorooctanesulfonate. The disinfection by-product chloroform, the fumigant 1,2-dichloropropane, the herbicides atrazine and hexazinone, and the herbicide degradates 2-chloro-6-ethylamino-4-amino-s-triazine, 2-chloro-4,6-diamino-s-triazine, 4-hydroxychlorothalonil, and metolachlor sulfonic acid were detected in more than 10 percent of domestic groundwater resources, but concentrations did not exceed regulatory benchmarks. Land use, groundwater age (fraction of modern water and mean age), and geochemical environment (oxic or anoxic conditions, pH, alkalinity) were associated with the distribution of high relative concentrations of inorganic and organic constituents. Young, oxygenated water is recharged along the Kern River and adjacent recharge ponds, or as irrigation water in the agricultural areas. High concentrations of nitrate and volatile organic compounds occurred in the oxic water in urban and agricultural areas. The fumigants 1,2,3-TCP, DBCP, and EDB were reported throughout the agricultural areas, whereas chloroform, tetrachloroethene, and PFAS were associated with urban land use. High uranium concentrations were associated with young, modern groundwater in agricultural areas with low pH and high bicarbonate. Total dissolved solids increased with distance from the Kern River, as the contributions of fresh, oxic water decreased. High concentrations of arsenic were present in older anoxic or alkaline groundwater away from areas of recharge. Overall, groundwater age, redox conditions, and the source of recharge as a result of different land uses contribute to large aquifer-scale portions of domestic groundwater resources that exceed health-based benchmarks for nitrate, uranium, and fumigant concentrations.

California↗
USGS2026

SlideDetect: Spatio-temporal landslide detection using a three-dimensional convolutional neural network

Landslides pose a serious and ongoing threat to both human lives and infrastructure worldwide; therefore, it is of interest to predict where and when landslides are likely to occur. Advances in machine learning techniques have spurred numerous studies aimed at estimating relative landslide propensity, but are limited to spatial (as opposed to temporal) prediction due to the sparsity of landslide timing data. We address this data gap by training SlideDetect, a 3-dimensional convolutional neural network (3D CNN), to identify landslides based on their spatial and temporal occurrence within multitemporal image stacks. We use an inventory of landsides triggered by the 2018 Hokkaido earthquake and two years of monthly composite optical imagery spanning this event. The model can identify not only landslide location but also landslide date with an area under the precision-recall curve (PR-AUC) of 0.84. We further present a new standard for presenting PR curve results that explicitly compares model performance at different confidence thresholds, allowing for clearer model evaluation and comparison. Our new approach to constraining landslide timing paired with this more consistent and objective method for evaluating model performance shows considerable promise, and with further application and testing, SlideDetect could enhance the data availability and tools needed to advance landslide hazard and risk assessments.

JGR Machine Learning and Computation↗
USGS2025

Flood-inundation maps for Río de la Plata in and near Comerío, Puerto Rico, 2025

Digital flood-inundation maps for a 3.1-mile reach of Río de la Plata in and near Comerío, Puerto Rico, were created by the U.S. Geological Survey (USGS). Water-surface profiles were computed for the stream reach by using a one-dimensional steady-state step-backwater model. The model was calibrated to the current (2025) stage-streamflow relation (rating curve 11.0) for the USGS streamgage 50043800, Río de la Plata at Comerío, Puerto Rico. The resulting hydraulic model was then used to compute 16 water-surface profiles for water levels (flood stages) ranging from 10.00 to 40.00 feet at the streamgage and ranging from “action stage” to above “major flood stage” as reported by the National Weather Service. The 40.00-foot stage was selected because it exceeds the peak stage of 34.86 ft recorded during Hurricane Maria at the USGS streamgage 50043800, Río de Plata at Comerío, Puerto Rico. The simulated water-surface profiles were then used in combination with a digital elevation model derived from light detection and ranging data to map the inundated areas associated with each flood profile. The flood-inundation maps and the supporting hydraulic model produced by this study can be used by emergency managers and local officials to assess flood mitigation strategies and to define flood hazard areas to help protect life and property, to coordinate flood response activities such as evacuations and road closures, and to aid post-flood recovery efforts.

Comerio↗
USGS2026

Upper Triassic igneous rocks of the southern Kenai Peninsula, Alaska – Prelude to Early Jurassic subduction along the western Wrangellia composite terrane margin

New U–Pb zircon geochronology identifies a latest Triassic (ca 214–201 Ma) igneous suite of tuff, hypabyssal dikes, and a pluton on the southern Kenai Peninsula, Alaska. The igneous suite was emplaced within Upper Triassic sedimentary rocks along the southern margin of Western Wrangellia, the western-most fragment of the Wrangellia composite terrane. The igneous rocks range from mafic (50.6% SiO 2 ) to felsic (78.3% SiO 2 ), characteristically have less than 1.55% K 2 O, and generally have low trace element abundances. The tonalitic and trondhjemitic magmas were largely sourced in mafic-rich lower crust and incompletely assimilated quartz and other mineral xenocrysts are common. Fractionation involving plagioclase and amphibole is indicated for some magmas and composite intrusions and igneous xenoliths indicate magma mixing was possible. Paleozoic and Precambrian inherited zircons and initial 87 Sr/ 86 Sr (0.704103–0.705609) and 143 Nd/ 144 Nd (0.512396–0.512777) ratios indicate that the Western Wrangellia crustal sources are heterogeneous and contain sialic components. The latest Triassic magmatism reflects processes that preceded Early Jurassic subduction along the Wrangellia composite terrane and Pacific Ocean plate boundary. These processes involved heating and melting of mantle lithosphere and lower crust as mantle instabilities accompanied the breaking of the plate boundary linkages. The Late Triassic transition to subduction along the Wrangellia composite terrane margin coincided with the transition to subduction cessation in the Late Triassic arcs of the western Intermontane terranes of Canada. The shift to subduction along the outboard Wrangellia composite terrane margin marks the beginning of the Pacific Ocean–Cordillera plate interactions that came to dominate the tectonic evolution of the northern Cordillera from the Early Jurassic to today.

Alaska↗
USGS1998

Transmissivity and water quality of water-producing zones in the intermediate aquifer system, Sarasota County, Florida

The intermediate aquifer system is an important water source in Sarasota County, Florida, because the quality of water in it is usually better than that in the underlying Upper Floridan aquifer. The intermediate aquifer system consists of a group of up to three water-producing zones separated by less-permeable units that restrict the vertical movement of ground water between zones. The diverse lithology, that makes up the intermediate aquifer system, reflects the variety of depositional environments that occurred during the late Oligocene and Miocene epochs. Slight changes in the depositional environment resulted in aquifer heterogeneity, creating both localized connection between water-producing zones and abrupt culmination of water-producing zones that are not well documented. Aquifer heterogeneity results in vertical and areal variability in hydraulic and water-quality properties. The uppermost water-producing zone is designated producing zone 1 but is not extensively used because of its limited production capability and limited areal extent. The second water-producing zone is designated producing zone 2, and most of the domestic- and irrigation-supply wells in the area are open to this zone. Additionally, producing zone 2 is utilized for public supply in southern coastal areas of Sarasota County. Producing zone 3 is the lowermost and most productive water-producing zone in the intermediate aquifer system. Public-supply well fields serving the cities of Sarasota and Venice, as well as the Plantation and Mabry Carlton Reserve well fields, utilize producing zone 3. Heads within the intermediate aquifer system generally increase with aquifer depth. However, localized head-gradient reversals occur in the study area, coinciding with sites of intense ground-water withdrawals. Heads in producing zones 1, 2, and 3 range from 1 to 23, 0.2 to 34, and 7 to 42 feet above sea level, respectively. Generally, an upward head gradient exists between producing zones 3 and 2. However, an upward head gradient between producing zones 2 and 1 does not consistently occur throughout Sarasota County, probably the result of greater ground-water withdrawals from producing zone 2 than from producing zone 1. The transmissivity of the intermediate aquifer system is spatially variable. Specific-capacity data from selected wells penetrating producing zones 2 and 3, were used to estimate transmissivity. Estimated transmissivity values for producing zones 2 and 3 range from about 100 to 26,000 feet squared per day and from about 1,300 to 6,200 feet squared per day, respectively. Because the capacity of specific water-producing zones is highly variable from site to site, estimating the performance of a specific water-producing zone as a water resource is difficult. Water samples collected during the study were analyzed for major-ion concentrations. Generally, bicarbonate type water from rock interaction occurs in northern Sarasota County; enriched calcium-magnesium-sulfate type water from deeper aquifers occurs in central Sarasota County; and sodium-chloride type water from saltwater mixing occurs in southern Sarasota County. In some areas of northern Sarasota County, the major-ion concentrations in water are lower in producing zone 2 than in producing zone 1. Major-ion concentrations in water are higher in producing zone 3 throughout the study area. A major objective of the study was to evaluate hydraulic and water-quality data to determine distinctions that could be used to characterize a particular producing zone. However, data indicate that both hydraulic and water-quality properties are highly variable within and between zones, and are more related to the degree of connection between and areal extent of water-producing zones than to aquifer depth and distance from the coast.

Florida↗
USGS2026

Moment magnitude for small earthquakes in the Delaware basin of west Texas and southeast New Mexico, USA

The Delaware Basin region of west Texas and southeast New Mexico has become one of the most prolific regions of seismic activity in the continental United States due to widespread hydraulic fracturing and wastewater disposal injection. In response to the increased number of earthquakes in this region, rapid and accurate characterization of earthquake sources is necessary to understand the evolution of seismic activity and level of seismic hazard associated with these earthquakes. This study re-evaluates earthquake magnitudes, estimating moment magnitude (MW) for small earthquakes in the Delaware Basin using 1) moment-rate spectra derived from S-wave coda envelopes, and 2) a relative magnitude method that relies exclusively on the ratio of waveform amplitudes between highly correlated waveform pairs. The coda-envelope method produces accurate M W estimates for small earthquakes ( M 1.5 – 3) that are consistent with independent, waveform modeled moment magnitudes for events with M W > 3 . Using the relative amplitudes method to extend these M W magnitudes to many other events, we successfully provide relative moment magnitude ( M W,rel ) values for 81% of the Texas Seismological Network catalog in the Delaware Basin region, and 45% of the USGS Induced Seismicity Project’s catalog of events in southeast New Mexico. The adoption and integration of the calibrated M W,rel method with current magnitude estimation methods offers valuable insights into the relationships between local and moment magnitude and will contribute to improved characterization of widespread induced seismicity.

New Mexico, Texas↗
USGS2026

Polystyrene microplastics alter the accumulation and elimination dynamics of silver nanoparticles in Daphnia magna

Rationale Microplastics (MPs) can interact with engineered nanomaterials and alter their environmental fate and bioavailability. However, their influence on the bioaccumulation dynamics of silver nanoparticles (AgNPs) in aquatic filter feeders under environmentally relevant conditions remains poorly understood. We hypothesized that polystyrene microplastics (PS-MPs) alter the uptake, elimination and overall bioaccumulation dynamics of AgNPs in Daphnia magna . Methodology Adult D. magna were exposed to isotopically labeled citrate-coated 109 AgNPs in the presence and absence of PS-MPs. Waterborne uptake, dietary uptake, assimilation efficiency, food ingestion and elimination were quantified experimentally and incorporated into a biodynamic model to predict steady-state silver (Ag) concentrations under environmentally relevant exposure scenarios. Results PS-MPs increased waterborne Ag accumulation by approximately eightfold compared with AgNPs exposure alone. However, Ag elimination was substantially faster, with 98% of accumulated Ag eliminated after five days of depuration compared with 66% without PS-MPs. PS-MPs also reduced food ingestion rates (IRs), while assimilation efficiency remained largely unchanged. Biodynamic modeling predicted that steady-state Ag accumulation was approximately fivefold greater in the presence of PS-MPs, with waterborne exposure becoming the dominant accumulation pathway. Discussion These findings suggest that PS-MPs alter bioaccumulation dynamics of AgNPs in filter-feeding organisms by enhancing organism-associated Ag during waterborne exposure while accelerating Ag elimination, likely through particle-associated transport and gut egestion. Collectively, these results underscore the complexity of Ag bioaccumulation processes when MPs and AgNPs interact and highlight the importance of studies conducted under environmentally relevant conditions.

Environmental Chemistry↗
USGS2024

Fish beta diversity associated with hydrologic and anthropogenic disturbance gradients in contrasting stream flow regimes

Understanding the role of hydrologic variation in structuring aquatic communities is crucial for successful conservation and sustainable management of native freshwater biodiversity. Partitioning beta diversity into the additive components of spatial turnover and nestedness can provide insight into the forces driving variability in fish assemblages across stream flow regimes. We examined stream fish beta diversity across hydrologic and anthropogenic disturbance gradients using long-term (1916–2016) site occurrence records ( n = 17,375) encompassing 252 species. We assessed total beta diversity (Sørensen dissimilarity), spatial turnover, and nestedness of fish assemblages in contrasting stream flow regimes across a gradient of decreasing flow stability: groundwater stable ( n = 77), groundwater ( n = 67), groundwater flashy ( n = 175), perennial runoff ( n = 141), runoff flashy ( n = 255), and intermittent ( n = 63) streams. Differences in total beta diversity among the stream flow regimes were driven predominantly (>86 %) by spatial turnover (i.e. species replacement) as opposed to nestedness (i.e. species loss or gain). Total fish beta diversity and spatial turnover were highest in streams with intermediate flow stability (groundwater flashy), while more flow-stable streams (groundwater stable and groundwater) had lower turnover and higher nestedness. Species turnover was also strongly associated with seasonal variation in hydrology across all flow regimes, but these relationships were most evident for assemblages in intermittent streams. Distance-based statistical comparisons showed significant correlations between beta diversity and anthropogenic disturbance variables, including dam density, dam storage volume and water withdrawals in catchments of groundwater stable streams, while hydrologic variables were more strongly correlated with beta diversity in streams with runoff-dominated and flashy flow regimes. The high spatial turnover of species implies that fish conservation actions would benefit from watershed-focused approaches targeting multiple streams with wide spatial distribution, as opposed to simply focusing on preserving sites with the greatest number of species.

Science of the Total Envionrment↗
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↗
USGS2026

Mineral chemistry perspective on remobilization of stored magma at Kamakai'a Hills, Southwest Rift Zone of Kilauea, Island of Hawai'i, USA

Differentiated magmas stored in the rift zones of Kīlauea have received more attention in recent years following eruption of andesite during the early phase of 2018 lower East Rift Zone activity. Despite this growing interest, some of the most voluminous eruptions of differentiated rift zone magmas remain poorly studied. One such eruption, and the most voluminous exposed differentiated flow field at Kīlauea, is the Kamakaiʻa Hills. This eruption took place in the Southwest Rift Zone of Kīlauea, a region that is hypothesized to contain a long-lived rift zone reservoir. The Kamakaiʻa Hills flow field encompasses >250 × 10 6 m 3 of basaltic andesite and basalt compositions with a mineral assemblage of orthopyroxene + clinopyroxene + plagioclase during its early ʻaʻā phase and clinopyroxene + plagioclase + olivine during its late pāhoehoe phase. To better understand storage conditions and magma accumulation, this study focuses on major, minor, and trace elements from the mineral assemblage present within the early ʻaʻā and late pāhoehoe phases. The diversity of clinopyroxene and plagioclase compositions within the early ʻaʻā and late pāhoehoe phases, as well as diverse compositions of plagioclase and orthopyroxene within the early ʻaʻā phase, suggest multiple magma bodies and limited pre-eruption magma mixing within the broader Kamakaiʻa Hills reservoir. Oscillatory zoning patterns (particularly in clinopyroxene) imply processes such as recharge events, magma mixing or mingling, or convection within a differentially cooling, chemically stratified reservoir over protracted time intervals, whereas only limited resorbed mineral textures indicate incomplete mixing of heat and chemically distinct magmas during the dike intrusion that triggered the eruption. Mineral-mineral and mineral-melt thermobarometry indicate predominantly shallow (≤2.5 km depth) crustal storage conditions of the cooled, differentiated magma (∼1100 °C and cooler for the basaltic andesites) to hotter temperatures for the basalts (all >1100 °C). Despite the known large standard errors estimated for mineral-melt and mineral-mineral barometry (10s to >100 MPa), the calculated pressures and depths broadly correspond with earthquake swarm depths beneath the Kamakaiʻa Hills, and drill core and fluid inclusion barometry storage depths of differentiated magmas within the lower East Rift Zone. The Kamakaiʻa Hills differentiated magmas have H 2 O contents (∼0.5 wt%, using plagioclase-melt hygrometry) equivalent to typical Kīlauea basalts. Our data and interpretations demonstrate a complex, long-lived rift zone storage system that consisted of multiple magma bodies and was mobilized into eruption through intrusion of a hotter and more primitive summit-derived (uprift) magma.

Hawaii↗
USGS2024

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

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

Earth and Planetary Science Letters↗
USGS2026

Top Elevation of Glacial Till and Thickness of the Big Sioux Aquifer Delineated From Electrical Resistivity Tomography Surveys Near Sioux Falls, South Dakota, 2022 and 2025

The City of Sioux Falls, South Dakota, requested the U.S. Geological Survey perform electrical resistivity surveys on three parcels of land north of the city. Electrical resistivity data were collected along a total of 22 transects during March 14–18, 2022, and November 17–21, 2025. Results from electrical resistivity surveys were used to delineate the top of glacial till deposits for the purpose of characterizing the Big Sioux aquifer near the city. Delineating geologic contacts provides important information on groundwater storage, flow dynamics, well design and placement, contaminant transport, groundwater–surface-water interactions, and regional water modeling. The top elevation of glacial till and the thickness of the Big Sioux aquifer varied among the three survey areas. The interpreted top elevation of glacial till in the North survey area decreases from east to west toward a slough, with elevations ranging from 1,403 to 1,418 feet (ft). The estimated thickness of the Big Sioux aquifer in the North survey area increased from east to west, with thicknesses ranging from 23 to 38 ft. The top elevation of glacial till in the Well 72 survey area generally decreases from northwest to southeast. Top elevations of the glacial till in the Well 72 survey area ranged from 1,400 to 1,409 ft along the southern end of transect W72_2. The estimated thickness of the Big Sioux aquifer in the Well 72 survey area was greatest along a southeast to northwest trending channel, with thicknesses ranging from 28 to 40 ft. The top elevation of glacial till in the Nose survey area generally decreases west toward the Big Sioux River. Top elevations of the glacial till in the Nose survey area ranged from 1,362 to 1,395 ft. The estimated thickness of the Big Sioux aquifer in the Nose survey area ranged from 33 to 70 ft.

South Dakota↗
USGS2025

Assessment of coastal and fluvial morphodynamic changes using Structure-for-Motion: A case study of the Sfȃntu Gheorghe Mouth (Danube Delta, Romania)

The ability to accurately map erosion, flooding, and habitat loss in coastal environments is crucial for formulating national strategies aimed at preventing and mitigating the impacts of natural disasters. A fundamental component of this process is the implementation of coastal morphodynamics monitoring through Structure-from-Motion (SfM) techniques, utilizing high-resolution 2D/3D data obtained from aerial photogrammetry. To assess morphodynamic changes over a three-year period (2022 – 2024), several SfM-based photogrammetric studies were conducted, each year, in the Romanian sector of the Danube-Black Sea coastal zone, specifically at the mouth of one of the Danube River distributaries (Sf Gheorghe branch) into the Black Sea, and along the left bank, near Sf Gheorghe locality, located within the Danube Delta Biosphere Reserve (DDBR). The essential equipment for aerial photogrammetry comprises Unmanned Aerial Vehicles (UAVs) and Global Navigation Satellite Systems (GNSS). In this study, the UAV used was a DJI Mavic 3T (Enterprise/Thermal) drone, complemented by two Trimble R12i and R4 GNSS systems, as well as approximately 10 Ground Control Points (GCPs). Data acquisition and processing were carried out using specialized photogrammetric software (Agisoft Metashape) along with various GIS tools (e.g., Blue Marble Geographics Global Mapper and ESRI ArcMap). The photogrammetric products generated for the study, as detailed in this paper, include Digital Elevation Models (DEMs), Digital Terrain Models (DTMs), orthomosaics (orthophotos), and others. At Sfântu Gheorghe beach, a comparison between 2023 and 2024 photogrammetric surveys revealed that the left bank of the Sf. Gheorghe Arm, at the river mouth into the Black Sea, suffered from a twist (erosion) of up to 64 metres. Additionally, on the selected perimetre (total area of 31,910 square meters ) from the beach and dune zone of Sf. Gheorghe, an area of up to 16,202 square meters was eroded between 2023 and 2024. This contrasts with the period between 2022 and 2023, during which deposition predominated. Erosion at the Danube mouths and the adjacent Black Sea coastline is driven by a complex interaction of natural and anthropogenic factors. Natural processes, including subsidence, sea-level rise, and episodic extreme storm events, contribute significantly to coastal dynamics. Meanwhile, human-induced factors, such as upstream hydrotechnical works that limits sediment transport, cutting of navigation canals, as well as the exacerbating effects of climate change, further accelerate erosion. The recent Structure-from-Motion (SfM) surveys provide essential quantitative data, enabling a detailed analysis of both short-term and long-term morphodynamic changes influenced by seasonal variations and extreme hydrometeorological events in this highly dynamic coastal system.

Sfȃntu Gheorghe Mouth (Danube Delta)↗
USGS2024

Influence of inherited structure on flexural extension in foreland basin systems: Evidence from the northern Arkoma basin and southern Ozark dome, USA

Extensional faults are key components of foreland basin systems. They form within the upper crust in response to flexure of the lithosphere and accommodate subsidence within the foredeep and forebulge depozones. Such faults are excellent proxies for orogenic system evolution and control the distribution of natural resources and hazards. However, the spatiotemporal evolution of flexural extension has not been documented previously at a regional scale, thereby limiting our understanding of underlying geodynamic controls. Here, we resolve late Paleozoic flexural extension in the northern Arkoma basin and southern Ozark dome, USA. We synthesize a large database of previous mapping, existing research, subsurface data, and geophysical data into 3D geologic and 2D kinematic models. Mesh surfaces representing several key horizons from the Carboniferous Period (ca. 335-306 Ma) were constructed. These surfaces were built from oil and gas well tops (n = ∼10,000) and surface geologic map contacts using an advanced kriging method. The mesh surfaces are offset by a complex 3D fault network, allowing detailed analysis of along-strike and down-dip variations in fault displacement. Analysis of the 3D model reveals a regular and repeated fault segmentation pattern wherein E -W striking, left- and foreland-stepping en échelon normal faults are segmented by inherited NE striking basement faults. Maximum vertical separation along the E -W normal faults is generally focused between the inherited NE-trending faults. This suggests that the inherited basement faults delocalized extensional strain during late Paleozoic normal faulting. Maximum vertical separation and fault localization may correlate to areas with high-amplitude positive magnetic anomalies interpreted as Mesoproterozoic granitic rocks. Speculative covariance of magnetic anomalies and fault displacements implies that the relatively strong basement granite concentrated stress, leading to localized faulting within the relatively thin sedimentary cover. Lastly, we show that flexural extension migrated southeast to northwest from the Chesterian-Morrowan (ca. 335-319 Ma) to the Desmoinesian (ca. 306 Ma). The migratory flexural extension may be explained by diachronous loading during Pangean assembly, or by synchronous loading but variable load compensation due to inherent factors.

Arkansas, Oklahoma↗
USGS2026

The U.S. Geological Survey 2025 Puerto Rico and U.S. Virgin Islands time-independent earthquake rupture forecast

We present the 2025 U.S. Geological Survey Puerto Rico and U.S. Virgin Islands (PRVI) time‐independent earthquake rupture forecast (ERF), developed for the 2025 update to the National Seismic Hazard Model (NSHM) for PRVI. The updated ERF improves upon a prior model from 2003, including an expanded fault inventory with slip‐rate estimates, updated seismicity catalogs, and refined subduction zone geometries and deformation models. It applies the fault‐system inversion methodology to solve for rates of ruptures on modeled faults, adapted from the 2023 NSHM (NSHM23) for the western United States, including the first application of the inversion to model rates on a U.S. subduction interface. Off‐fault and intraslab seismicity are constrained by observed seismicity and use updated methods developed for NSHM23. Uncertainties in model components are substantial, and the ERF represents epistemic uncertainties through a comprehensive logic tree consisting of 1.7 billion logic‐tree branches combined across all sources.

Puerto Rico, U.S. Virgin Islands↗