USGS Science⌕ Search

SEARCH · USGS Science

Results for “Resources, Conservation, and Recycling”

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,663 records · Page 70Linked to original sources
USGS2026

Efficient physics‐informed ground‐motion simulations with reduced‐order models: CyberShake implications and high‐resolution site terms for southern San Andreas fault earthquakes

Recent advances in Probabilistic Seismic Hazard Analysis (PSHA) leverage physics‐based ground‐motion simulations to estimate seismic hazard, such as the CyberShake project. However, computational costs quickly escalate when performing PSHA for numerous faults or sites and can become prohibitively expensive. To reduce computational demands, CyberShake uses reciprocity and interpolates physics‐informed corrections from simulations conducted at fewer locations, but the accuracy of these interpolations remains poorly quantified. To quantify the interpolation accuracy, we derive high‐resolution, frequency‐dependent site terms for southern California and compare them with interpolated site terms using the CyberShake approach. We accomplish this by performing a set of earthquake point‐source simulations distributed along the nonplanar fault geometry for the southern San Andreas fault (SSAF) extending from Bombay Beach to Lake Hughes. Using SeisSol, we simulate three minutes of viscoelastic seismic wave propagation for these sources and store the horizontal‐component Green’s functions for 480,000 sites. We then use a scientific machine learning approach based on interpolated proper orthogonal decomposition to construct an accurate reduced‐order model of the Green’s functions to efficiently predict effective amplitude spectra (EAS) for finite‐source rupture models of SSAF earthquakes. Using minimum curvature interpolation with tension, as used in CyberShake, we compare the interpolated site terms against our high‐resolution site terms. We identify local discrepancies with EAS differing by up to a factor of approximately three. Furthermore, we identify locations where unexpectedly high or low ground motions are missed when using the interpolated dataset for these earthquakes. We estimate that our approach may be used within CyberShake to reduce the time‐to‐solution by a factor of 336 for the entire earthquake rupture forecast. Our analysis of physics‐based site terms provides more insight into the seismic hazard due to SSAF ruptures and guides future developments by combining high‐performance computing and reduced‐order modeling techniques for PSHA.

California↗
USGS2025

Paleoproterozoic vein graphite mineralization caused by decarbonation in the Ruby Range, Montana, USA

Hydrothermal graphite veins are a possible source for modern battery materials and require better understanding of their carbon source(s) and absolute timing to develop mapable criteria for exploration models. We present new observations of graphite vein and alteration paragenesis and U-Pb LA-ICP-MS titanite age data from the Ruby prospect, Montana, USA, that constrain mineralization timing and source. The graphite veins cut high-temperature metamorphic rocks of the lower Christensen Range suite and are associated with intense diopside (Di0.69Hd0.27Jo0.04) alteration of marble. The oldest titanite ages in fresh marble and partially altered calc-silicate gneiss are ca. 2,500 – 2,450 Ma, show elevated REE values, and Eu/Eu* <1.5, consistent with growth during early regional metamorphism. Titanite in diopside-graphite alteration, interpreted as part of the hydrothermal vein-forming episode, cluster around 1,750 Ma, are characterized by lower REE values, and show Eu/Eu* >2; suggesting a low fO 2 fluid generated from metamorphism during the Big Sky orogeny. Our paragenetic observations and titanite ages indicate graphite vein formation via skarnoid decarbonation reactions in marble late in regional orogenesis. Granulite-facies, carbonate-bearing, supracrustal rock terranes are thus favorable for hydrothermal graphite vein deposits.

Montana↗
USGS2026

Increasing artificial light at night enhances salmonid predator foraging efficiency in an urbanized lake

Artificial light at night (ALAN) poses a threat to ecosystems globally. It includes both direct and indirect light, or skyglow, which occurs when ALAN scatters in the atmosphere, extending beyond its original source. We analyzed ALAN trends in Lake Washington, WA, from 2014 to 2023 using Visible Infrared Imaging Radiometer Suite (VIIRS) nighttime light measurements, evaluated the relationship between in situ and satellite measurements, and modeled predator search volumes for a juvenile salmon predator, cutthroat trout ( Oncorhynchus clarki ), as a function of ambient light and turbidity conditions. Open water regions experienced significant increases in ALAN, while nearshore areas primarily showed no or negative trends, revealing the role of skyglow in shaping open water light environments. Using a visual foraging model for cutthroat trout, we found that juvenile salmon at the shallow southern pelagic site experienced light 28 times brighter, resulting in a 168% greater predation vulnerability than those at the northern site. In the Ship Canal, a narrow corridor for outmigrating salmon, predator search volumes were 249% higher than at the southern site. These contrasts in predation vulnerability demonstrate how local conditions influence predator–prey dynamics and provide critical insight for targeting mitigation of both nearshore and distant light sources.

Washington↗
USGS2025

Refining the earthquake history of south-central Alaska through lake records

The Alaska–Aleutian subduction zone (AASZ) is one of the world's most seismically active plate boundaries and the source of the 1964 M w 9.2 Great Alaska earthquake–the second largest instrumentally recorded earthquake in the world. Understanding the nature and frequency of such earthquakes is necessary for seismic and tsunami hazard assessment, but instrumental and historical records that span less than 150 years are too short to allow a statistically reliable analysis of earthquake recurrence times. This calls for studies of evidence of past earthquakes, extending the earthquake catalog further back in time. Subduction-zone paleoseismology in south-central Alaska is predominantly based on coastal evidence of land-level changes and tsunamis generated by megathrust earthquakes and preserved in the geological record. A complementary approach is lacustrine paleoseismology, which is still a relatively young discipline in Alaska. However, globally, lake basins are well-established high-resolution and continuous recorders of paleoseismic activity along subduction zones, relying on the identification of underwater landslide deposits and turbidites generated by seismic shaking. As a result, lake basins not only register ground shaking from megathrust earthquakes, but also from intraslab and crustal earthquakes, which are typically not accompanied by significant land-level changes. In this review paper, we combine coastal and lacustrine paleoseismology approaches to refine the south-central Alaskan earthquake history by comparing the paleoseismic records from two lakes (i.e., Eklutna Lake, located in the Chugach Mountain Range, and Skilak Lake, situated on the Kenai Peninsula) with the coastal and crustal earthquake catalog in Alaska. The resulting age ranges of all known megathrust earthquakes involving the Alaskan megathrust between the Kodiak and Prince William Sound (PWS) sections are more precise and accurate for the last 1.3 kyrs BP than the previously published age ranges from coastal records. As a result, this study supports the following key conclusions: (1) The 1964 CE earthquake was an exceptionally strong and unique event in the last 2000 years, rupturing the PWS, Kenai, Barren Islands, and Kodiak sections simultaneously. (2) The high-resolution and seasonal markings of the varved lake records now disentangle for the first time closely timed earthquakes, which was not possible based on the coastal evidence alone. (3) No persistent megathrust rupture boundaries exist. So, the possibility of a full rupture of the entire eastern AASZ, from PWS to Semidi cannot be excluded. (4) The rupture pattern in the eastern AASZ reveals superimposed cycles of multi-asperity ruptures (1964 earthquake) and clustered complementary partial ruptures, or rupture cascades. (5) The PWS section hosts the largest asperity in the eastern AASZ. (6) The shaking record of megathrust earthquakes indicates a time-dependent (quasiperiodic) behavior for the study area, but the observation of complementary clusters means that the hazards will not drop to zero but instead may even increase for a neighboring section. (7) The time-independent behavior of intraplate earthquakes implies that the intraslab hazard did not decrease following the 2016 and 2018 earthquakes. This study utilizes an integrated approach for subduction zone paleoseismology as a solution for unraveling recurrence and rupture patterns in Alaska, which can be applied worldwide.

Alaska↗
USGS2025

Global patterns of coseismic landslide runout mobility differ from aseismic landslide trends

Coseismic landslides significantly contribute to human and economic losses during and immediately following earthquakes, yet very little data on the runout of such landslides exist. While well-established behavior of aseismic (e.g., hydrologically triggered) landslide runout mobility suggests strong correlation between landslide size and mobility, limited studies of coseismic landslide runout find conflicting mobility trends. We present a global dataset of runout lengths produced from a new automated method for estimating landslide runout, developed and validated using 1726 manually mapped landslides from five unique earthquakes. We then apply the automated runout tool to 23 global earthquake-induced landslide inventories, producing a compiled database of 73,665 measured and estimated runout lengths of coseismic landslides to assess mobility trends. We find a significant divergence between well-established aseismic mobility trends and that of coseismic landslides, with far greater scatter and more complex mobility patterns in earthquake-triggered landslides. As a function of landslide size, we observe global coseismic landslide mobility patterns are bilinear, becoming increasingly less mobile with increasing size above some threshold. This discordance between aseismic and coseismic landslide mobility may be a function of landslide type, kinematics, hydrology, and or setting that systematically differ between triggering mechanisms and should be explored in more depth to develop predictive models of these unique runout patterns. These results suggest hazard and risk models for coseismic landslides may significantly under-predict or over-predict impacts, depending on the size of triggered landslides.

Engineering Geology↗
USGS2026

Latest Pleistocene to 19th-century earthquakes on bending-moment reverse faults of the Seattle fault zone, Washington

Fault-related folds and their associated secondary faults play a critical yet often underrecognized role in accommodating strain and generating earthquakes in active fold-and-thrust belts. In the Seattle fault zone (SFZ), Washington, USA, we present new paleoseismic, geomorphic, and geophysical evidence for late Pleistocene and Holocene earthquakes on shallow, south-dipping secondary faults—the Lytle Beach and Vasa Park faults—that lie within the hanging wall of the greater SFZ and are on trend with the primary, blind Blakely Harbor fault. Our data show that these structures have ruptured independently, producing localized uplift and deformation at the surface, with the most recent event (RH2) likely occurring in the early nineteenth century. While a temporal overlap between the late Pleistocene RH1 and VP1 earthquakes raises the possibility of a ≥35 km rupture along the Blakely Harbor fault, structural and temporal evidence instead supports independent rupture on individual faults related to folding. We interpret these faults as bending-moment reverse faults that formed within a synclinal hinge zone of the main fault, reflecting mechanical and kinematic influences of the broader fault system. Combined with prior studies, our findings indicate that faulting related to folding dominates the mode of strain release within the SFZ since the late Pleistocene with more frequent earthquake recurrence (∼350 yr) over the past ∼2500 yr.

Washington↗
USGS2026

Long-term intermittent connection between the western Snake River Plain and Columbia basin: A two-phased incision history of Hells Canyon

For more than a century, researchers have debated the route of the Snake River across the northern Cordillera and U.S. Pacific Northwest, including the associated incision history of Hells Canyon. Here, we use detrital zircon U-Pb provenance analysis of Miocene strata upstream and downstream of Hells Canyon to constrain the evolution in fluvial pathways over time. Downstream of Hells Canyon, we show that the ca. 10−6 Ma Clarkston Heights gravel was dominantly sourced from the nearby Salmon and Clearwater rivers but also with substantial input from drainages that now flow into the western Snake River Plain. These same western Snake River Plain drainages were contributing to time-equivalent Lake Idaho strata, upstream of Hells Canyon, indicating a Hells Canyon fluvial connection between the western Snake River Plain and Columbia basin before 6 Ma. Results from Lake Idaho strata suggest a two-phased lacustrine history in the western Snake River Plain: Before ca. 6.0 Ma, lake strata were derived from local tributaries, suggesting that the western Snake River Plain was isolated from the modern upper Snake River system. Between ca. 4.3 Ma and 2.2 Ma, the source area expanded to include eastern tributaries in conjunction with migration of high-standing topography of the Yellowstone hotspot. Our results challenge the long-held hypothesis of “capture” of the Snake River through Hells Canyon. Instead, we infer that Hells Canyon was a long-established route for outflow of lakes occupying the western Snake River Plain, possibly intermittently, during a ca. 10−2 Ma phase of slow regional incision, followed by rapid incision in Hells Canyon starting ca. 2 Ma in conjunction with erosion of the topographic barrier impounding Lake Idaho.

Idaho, Oregon↗
USGS2025

Dendroseismological investigation of redwood trees along the North Coast section of the San Andreas Fault

Sequoia sempervirens (coast redwood) tree rings have the potential to annually resolve late-Holocene earthquakes on the northern San Andreas Fault based on direct (e.g., physical damage) and indirect (e.g., co-seismic environmental change) impacts, but scarcity of suitable samples and challenges crossdating this long-lived species have limited progress. More precise dating of the pre-1906 (penultimate) earthquake can improve hazard assessment and understanding of rupture segmentation. We target old trees (maximum >815 yr) along the North Coast section of the fault (increment cores via rope-climbing, 11 living trees; plunge cuts, 23 stumps) and employ complementary disturbance detection methods including radial-growth averaging (tree- and series-level), cataloging anatomical indicators (e.g., traumatic resin ducts, TRD), and dating structural components (e.g., reiterated trunks, leans). Multi-centennial ring-width chronologies at Fort Ross (1569−2023) and Gualala (1397−2023) promote continued study with incomplete crossdating limiting utilization of some series. Growth pulses (reductions, releases) and TRD dispersed across the record reflect dynamic environments that obfuscate detection of earthquake signals. The 1906 earthquake did not leave strong signatures on most trees, and when it did, within-tree response varied from normal presentation to discoloration, TRD, and missing rings. Synchrony of indicators at both locations identified 1678−1680 (6 of 15 trees) and 1698−1700 (8 of 16 trees) as the strongest disturbances among dated rings in the time range of the penultimate earthquake, peaking at 1698 (15.7 % of possible growth and anatomical indicators), but the triggering mechanisms for these events are unknown.

California↗
USGS2025

Seamless Geologic Map Database for the Intermountain West, United States: A foundational dataset for mineral systems analysis

The Intermountain West has a complex geologic history, resulting in the formation of a diverse array of mineral deposits. Effective mineral exploration requires understanding the spatial and temporal relationships among geologic processes and events, a key focus of the mineral systems approach to exploration. This paper presents the Intermountain West Seamless Geologic Map Database, a unified dataset designed to support mineral exploration. Integrating geologic provinces, structural settings, and hydrothermal alteration, the database leverages the Seamless Integrated Geologic Mapping (SIGMa) extension to the USGS Geologic Map Schema (GeMS) to standardize geologic data from varied sources. SIGMa's hierarchical stratigraphic organization and feature-level metadata enhance data interoperability and reusability, enabling seamless query, analysis, and visualization of lithology, structural features, mineral deposits, geochronology, hydrothermal alteration. volcanic activity, and By providing a regionally consistent and dynamically evolving geologic map, this database provides a foundational framework for mineral exploration and geologic research. It also allows for an efficient workflow that expedites the publication of integrated geologic map databases.

Conference Paper↗
USGS2025

Mineralization processes at Escanaba Trough

Escanaba Trough off northern California is the southernmost and only sediment covered segment of the Gorda Ridge. Hydrothermal mineralization was first discovered in Escanaba Trough in 1985, and subsequent investigations included the delineation of the NESCA (northern Escanaba) and SESCA (southern Escanaba) sections of Escanaba Trough, as well as drilling by Ocean Drilling Program (ODP) Leg 169 in 1996. Preliminary research conducted in the 1980s distinguished two classes of sulfide mineral assemblages within the Escanaba Trough: one dominated by pyrrhotite (Fe 1-x S) and one containing ‘mixed sulfides’ including chalcopyrite (CuFeS 2 ) and sphalerite (ZnS). In 2022, we conducted a research expedition to the Escanaba Trough to more thoroughly characterize mineralization processes and consider preservation of precipitated minerals throughout hydrothermally active and inactive sections of Escanaba Trough, with the specific objective of characterizing the host minerals, and processes of enrichment for specific elements which are currently listed as critical elements and can also be toxic under certain conditions (Sb, As, Ba, Co, Ga, Ge, Te). The following will include initial results from that work.

California↗
USGS2026

Transit time modeling framework for predicting freshwater salinization in urban catchments

The salinity of inland freshwaters is rising globally, particularly in urban watersheds where winter road deicers are widely applied. Attributing stream salinity dynamics to specific sources and transport pathways remains challenging due to episodic salt inputs, engineered drainage, and strong coupling between hydrology and subsurface storage. We present a modeling framework that couples climate-driven deicer build-up and wash-off with transient transit time distribution theory to simulate salt transport through drainage, interflow, and groundwater pathways. Applied to an urban watershed in Northern Virginia (USA), the model reproduces ten years of high-frequency stream salinity measurements across daily-to-decadal timescales. The calibrated model implies an average deicer application of 206 tonnes Cl −1 yr −1 , or roughly one 20 kg bag of rock salt person −1 yr −1 when normalized by the 20,000 people living in the watershed. In winter months, higher infiltration routes a large fraction of snowmelt and deicers into shallow subsurface pathways, enhancing vadose-zone and interflow contributions to stream salinity. Limited subsurface storage capacity and seasonal hydrologic turnover flush excess chloride from the vadose zone and groundwater during subsequent summer storms. By linking climate-driven deicer inputs, hydrologic connectivity, and stream water age, the framework provides a transferable basis for diagnosing and managing freshwater salinization in urban watersheds.

Virginia↗
USGS2024

Potential corrosivity of untreated groundwater in Louisiana

Corrosive groundwater can cause lead, copper, and other metals to leach from pipes and plumbing fixtures in water distribution systems. Metals, if ingested, could lead to serious health implications to the nearly 2.9 million people in Louisiana who obtain their drinking water from groundwater sources. Four indices—the Langelier Saturation Index (LSI), Ryznar Stability Index (RSI), Puckorius Scaling Index (PSI), and the Potential to Promote Galvanic Corrosion (PPGC)—in addition to an analysis which normalized the results from the existing indices, the Combined Index (CI), were used to assess the corrosivity of groundwater in Louisiana and identify areas within eight major aquifers and aquifer systems with moderate to high corrosivity potential. The purpose of this study is to provide State and local governments, public water system managers, and the nearly 500,000 private well owners in Louisiana with information needed to manage drinking-water supplies and mitigate potential health risks related to leaching of metals from water pipes and fixtures. The average scores of untreated groundwater samples from approximately 375 wells by index are as follows: LSI, −1.28; RSI, 9.78; PSI, 9.34; and CI, 4.14. The PPGC does not produce a numerical score, but the total percentage of class counts can be used to assign a classification; overall, samples in Louisiana were classified as significant concern. The percentages of groundwater samples from wells classified as potentially corrosive, by index, are as follows: LSI, 53 percent; RSI, 94 percent; PSI, 81 percent; PPGC, 98 percent; and CI, 81 percent. The percentages of samples classified as indeterminate, by index, are as follows: LSI, 46 percent; RSI, 5 percent; PSI, 12 percent; PPGC, 0 percent; and CI, 18 percent.

Louisiana↗
USGS2026

Site response models based on geometric parameters for southern California sedimentary basins

Site response in sedimentary basins is influenced by complex three-dimensional (3D) features, including trapping of seismic waves, focusing of seismic energy and basin resonance. Current ground motion models (GMMs) incorporate basin effects using one-dimensional parameters like V S30 and shear wave velocity isosurface depths, which are limited in capturing lateral and 3D effects. To address these limitations, we develop seismic site response models based on novel parameters that represent multi-dimensional properties of the Los Angeles Basin (LAB) geometry and shear wave velocity. We define a basin shape for the LAB using depth to subsurface geologic interfaces associated with the oldest sedimentary deposits (depth to a particular shear wave velocity horizon, i.e., 1.5 km/s - z 1.5 ) and the depth to the crystalline basement ( z cb ) which are determined using geologic cross sections and community seismic velocity model profiles. We explore a suite of geometric descriptors computed for the LAB and southern California, from which three parameters with the greatest predictive potential are selected and evaluated using empirical ground motion residual analyses in combination with the Boore et al. GMM. The results demonstrate that the zonal heterogeneity index ( ), standard deviation of the absolute difference between z 1.5 and z cb ( ) and standard deviation of z cb ( ) each provide a reduction in site-to-site variability ( ϕ S2S ) of empirical GMMs. The reduction in ϕ S2S is period-dependent, with average decreases of 3%, 26% and 6% for , , and , respectively. Although these reductions are modest from an engineering application perspective, they are statistically significant, underscoring the inherent difficulty in fully characterising complex basin effects. Collectively, these findings indicate that the inclusion of basin-specific geometric parameters yields measurable, albeit incremental, improvements in site response prediction and establishes a framework for the progressive refinement of seismic hazard characterisation within sedimentary basins.

California↗
USGS2025

Revisiting an enigma on California's north coast: The Mw6.5 Fickle Hill earthquake of 21 December 1954

Many earthquakes occur along the North Coast of California in the vicinity of the Mendocino Triple Junction (MTJ), where the Pacific, Gorda, and North American (NA) plates meet, and on the adjacent plate boundaries. The MTJ marks the nexus of the Mendocino and San Andreas faults with the Cascadia subduction zone (CSZ). Historically, most large earthquakes around the MTJ have been within the offshore Gorda plate and its subducted portion beneath the NA plate. North of the MTJ, active faults mapped in the NA plate are part of the CSZ fold‐and‐thrust belt. Although some events have been detected in the NA plate, no large historic events have been associated with mapped surface faults. The 21 December 1954 M w 6.5 earthquake in Humboldt County is one possible exception. Using published data from catalogs and articles, unpublished data from Berkeley’s archives, and S ‐ P times interpreted from two U.S. Coast and Geodetic Survey (USCGS) accelerometers, we determine a probability cloud for the earthquake’s hypocenter using NonLinLoc. The highest probability location lies beneath Fickle Hill just east of the city of Arcata, California, at 40.87° N, 124.03° W, and ∼11 km depth. Using P ‐wave polarities from Berkeley stations and the digitized waveforms from the accelerometers, we find that the focal mechanism most consistent with the data indicates thrust movement with strike, dip, and rake of 350°, 10°, and 90°, respectively, at a depth of 14 km. Given the depth uncertainties of both this event and the megathrust, this implies that the earthquake most likely took place on the subduction interface rather than on the mapped faults in the Mad River fault zone that trend 322° and dip to the northeast. The revisited intensity in the epicentral region also supports a location beneath Fickle Hill to the east of the city of Arcata, California.

California, Idaho, Nevada, Oregon, Washington↗
USGS2024

Recent applications of the USGS National Crustal Model for Seismic Hazard Studies

The U.S. Geological Survey is developing the National Crustal Model (NCM) for seismic hazard studies to facilitate modeling site, path, and source components of seismic hazard across the conterminous United States. The NCM is composed of a 1km grid of geophysical profiles, extending from the Earth’s surface into the upper mantle. It is constructed from a threedimensional (3D) geologic framework and geophysical rules that use (1) a petrologic and mineral physics database; (2) a 3D temperature model; and (3) a calibrated rock type- and age-dependent porosity model. Parameters needed to estimate site response for existing ground motion models (GMMs), including the time-averaged velocity in the upper 30 meters (VS30), the depths to 1.0 and 2.5 km/s shear-wave velocity (Z1.0 and Z2.5), and sediment thickness, can be computed from the NCM. As GMMs continue to improve in the future, other metrics could also be extracted or derived from the NCM, such as fundamental period, site attenuation (ko), a fully frequency-dependent site response function, or 3D geophysical volumes for wavefield simulations. Application of the NCM may also benefit other aspects of seismic hazard analysis, including better accounting for path-dependent attenuation and geometric spreading, more accurate estimation of earthquake source properties such as hypocentral location and stress drop, and calculation of crustal strength profiles that inform estimates of the base of seismicity.

conterminous United States↗
USGS2026

Surface variable‐based machine learning for scalable arsenic prediction in undersampled areas

In the United States, private wells are not federally regulated, and many households do not test for Arsenic (As). Chronic exposure is linked with multiple health outcomes, and risk can change sharply over short distances and with well depth. Coarse maps or sparse sampling often miss exceedances. Most existing models operate at ∼1 km resolution and use groundwater chemistry or detailed geologic logs, which limits their use in undersampled areas where improved guidance is most needed. We overcome these limitations by developing a machine learning model for Minnesota, USA, that predicts As exposure risk using only surficial variables from remote sensing and global data sets. Variables related to surface water hydrology and geomorphology are selected based on mechanistic links that control redox conditions and As mobilization. Local training was essential, and surficial geology variables that are more sensitive to local conditions were needed to maximize model accuracy. The resulting complete model was sufficiently sensitive to generate accurate and detailed risk maps and depth profiles of As concentrations above the 10 μg/L maximum contaminant level. Accuracy depended on local training data density. We identified a training data density of 0.07 wells/km 2 as a practical target for stable county-level performance. Maps of exceedance probabilities highlight priority areas for testing that are particularly important in rural communities that have received less sampling. These results support public health action by guiding where to install wells and where to test them, how much new sampling is needed, and where treatment outreach is most urgent.

Minnesota↗
USGS2025

Earthquake recurrence estimates for northern Caribbean faults from combinatorial optimization

We use combinatorial optimization to find the optimal spatial distribution of random samples of earthquakes (≥6.5) that minimize the misfit in target slip rates for all faults in the northeast Caribbean, and we derive magnitude-frequency relationships with uncertainties for these faults. Slip rates for many faults are derived from GPS block models, not direct measurements, because of their underwater locations. Predicted recurrence rates for eastern Hispaniola and Puerto Rico Trench faults are 220-450 yr for M7 and 3-5 kyr for M8, with maximum feasible magnitude of M8.2. The most frequent earthquakes with magnitudes ≥7.0 are predicted on the large upper plate strike-slip faults, Enriquillo (EF) and Septentrional Fault, commensurate with the historical record. Calais et al. (2023) suggested that shortening in western Hispaniola is accommodated on the offshore Jérémie and onshore Malpasse faults north and south of EF, instead of on terrestrial faults in western Hispaniola and EF. Because of our system-modeling approach, such a configuration predicts less frequent earthquakes on EF and on western Hispaniola and Muertos convergent zones. Recurrence times of a few 100s yr for M6.7 earthquakes is predicted on the submerged faults in Mona Passage, and infrequent M>7 earthquakes are predicted on the Virgin Islands faults.

Lesser Antilles↗
USGS2025

Dietary bioavailability of uranium to a model freshwater invertebrate

Uranium (U) mining increases environmental exposures. Understanding how U is taken up by organisms can aid in evaluating the potential for bioaccumulation and toxicity. Although the importance of aqueous geochemical speciation is well recognized for U bioavailability after dissolved exposures, far less is known about the processes controlling U bioavailability after dietary exposures. This study characterizes the biogeochemical drivers of dietary U uptake in the freshwater snail Lymnaea stagnalis in laboratory experiments. Solids tested included benthic diatoms pre-exposed to dissolved U(VI), soils from contaminated U mine sites, and colloidal hydrous ferric oxide (HFO) synthesized in the presence of dissolved U(VI) or with U complexed by natural organic matter (NOM). Results showed that U was bioavailable from all solids. Uranium assimilation efficiency (AE), a proxy for dietary U bioavailability, varied among solids. AE was lowest for the U-contaminated soils (25 ± 17%) and highest for the U-laden diatoms (71 ± 13%). AE varied slightly among HFO preparations, suggesting modest influences of NOM and iron on U bioavailability. Increases in dietary U exposures reduced feeding rates, and the extent of feeding inhibition appeared inversely related to U bioavailability. The high U assimilation and range of bioavailability have implications for toxicity risks inferred without considering dietary uptake.

Environmental Science and Technology↗