USGS Science⌕ Search

SEARCH · USGS Science

Results for “PLoS ONE”

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

Critical Minerals in Ores (CMiO) database

Critical minerals are commodities essential to modern industrial and strategic technologies and are highly vulnerable to supply chain disruption. The Critical Minerals Mapping Initiative (CMMI) is a collaboration among the U.S. Geological Survey (USGS), the Geological Survey of Canada, and Geoscience Australia that aims to deepen global understanding of where critical minerals are located. A key output of this initiative is the Critical Minerals in Ores (CMiO) database that is advancing our collective understanding of critical minerals distributions. For instance, publicly available data on the concentrations of many critical minerals are sparse because these commodities can only be produced in small, yet essential, quantities compared to the primary commodities like copper and zinc. The CMiO database helps bridge this gap by offering high-quality, multielement geochemical data from a wide variety of critical mineral-bearing deposits around the world. Importantly, it uses a novel consensus deposit environment, group, and type classification scheme developed by the agencies that allows comparisons among ore deposits from different regions. The CMiO database contains geochemical data for more than 20,000 samples from more than 100 deposit types comprising 10 deposit environments.

Fact Sheet↗
USGS1993

Hydraulic properties of the Madison aquifer system in the western Rapid City area, South Dakota

Available information on hydrogeology, data from borehole geophysical logs, and aquifer tests were used to determine the hydraulic properties of the Madison aquifer. From aquifer-test analysis, transmissivity and storage coefficient were determined for the Minnelusa and Madison aquifers, and vertical hydraulic conductivity (K v ') along with specific storage (Ss') for the Minnelusa confining bed. Borehole geophysical well logs were used to determine the thickness and location of the Minnelusa aquifer, the lower Minnelusa confining bed, and the Madison aquifer within the Madison Limestone. Porosity values determined from quantitative analysis of borehole geophysical well logs were used in analyzing the aquifer-test data. The average porosity at the two aquifer-test sites is about 10 percent in the Minnelusa aquifer, 5 percent in the lower Minnelusa confining bed, and 35 percent in the Madison aquifer. The first aquifer test, which was conducted at Rapid City production well #6, produced measured drawdown in the Minnelusa and Madison aquifers. Neuman and Witherspoon's method of determining the hydraulic properties of leaky two-aquifer systems was used to evaluate the aquifer-test data by assuming the fracture and solution-opening network is equivalent to a porous media. Analysis of the aquifer test for the Minnelusa aquifer yielded a transmissivity value of 12,000 feet squared per day and a storage coefficient of 3 x 10 -3 . The specific storage of the Minnelusa confining bed was 2 x 10 -7 per foot, and its vertical hydraulic conductivity was 0.3 foot per day. The transmissivity of the Madison aquifer at this site was 17,000 feet squared per day, and the storage coefficient was 2 x 10 -3 . The second aquifer test, which was conducted at Rapid City production well #5 (RC-5) produced measured drawdown only in the Madison aquifer. Hantush and Jacob's method of determining the hydraulic properties of leaky confined aquifers with no storage in the confining bed was used to evaluate the aquifer-test data by assuming the fracture and solution-opening network is equivalent to a porous media. The analysis of data from the RC-5 aquifer test showed that transmissivity was not equal in all directions. Hantush's method was used to determine the direction of radial anisotropy and magnitude of the major and minor axes of transmissivity. The major axis of transmissivity is at an angle of 42° east of north, and the transmissivity along this axis is about 56,000 feet squared per day. The minor axis of transmissivity is at an angle of 48° west of north, and the transmissivity along this axis is about 1,300 feet squared per day. The major axis of transmissivity intersects Cleghorn Springs, a large resurgent spring on the west edge of Rapid City. The shape of the potentiometric contours of the Madison aquifer near RC-5 agree with the orientation of the transmissivity ellipse. The average value of the storage coefficient from the isotropic analysis of the aquifer-test data was 3.5 x 10 -4 , and the average vertical hydraulic conductivity of the lower Minnelusa confining bed was 9.6 x 10 -3 foot per day.

South Dakota↗
USGS2025

REDPy: A Python tool for automated repeating earthquake detection and visualization

Detecting and cataloging seismic events are among the most fundamental tasks in seismology. Many standardized tools for these tasks exist, including the open‐source package repeating earthquake detector in Python (REDPy). REDPy generates an organized catalog of seismic events from continuous waveform data, in which events are automatically separated into groups (“families”) by their waveform similarity through cross‐correlation. REDPy also automatically generates various outputs that allow a user to visualize important trends in the catalog, which may be used in real time or in retrospective analyses to allow rapid identification of interesting features. The code was designed for near‐real‐time volcano monitoring but is applicable across a broad range of use cases in seismology and seismoacoustics. In this article, the utility and performance of REDPy are demonstrated on two highly seismogenic volcanic eruption sequences: the onset of the dome‐building eruption of Mount St. Helens, Washington, from 2004 to 2005, and the entirety of the summit caldera collapse sequence of Kīlauea, Hawai‘i, in 2018. This article is meant to be a companion to the documentation of the code; in addition to detailing the basic required inputs, script functionality, and resulting outputs, the reasonings behind several important design decisions are also discussed.

Seismological Research Letters↗
USGS2024

Quantitative risk of earthquake disruption to global copper and rhenium supply

Earthquakes have the potential to substantially affect mining operations, potentially leading to supply chain disruptions and adversely affecting the global economy. This study explores the quantification of earthquake risk to copper and rhenium commodity supply by examining the spatial concentration of high earthquake hazard areas and the commodity-specific mining, smelting, and refining operations across the globe. Because many of the largest facilities are concentrated geographically near the highly seismic regions of South America, East Asia, and the Pacific, there is a potential for cascading effects on the entire supply chain. The analysis indicates that the expected annual disruption of global production is 0.3–1.1 percent for copper mines, 1.8–4.0 percent for smelters, and 1.5–3.3 percent for refineries. Expected annual disruption of global rhenium production capacity is 0.32–1.32 percent. The research highlights that the potential lost revenue from earthquake disruptions is from $315 million to $1.29 billion for copper mining, $1.92 billion to $4.33 billion for copper smelting, $2.06 billion to $4.52 billion for copper refining, and $337,000 to $1.40 million for rhenium production capacity.

Open-File Report↗
USGS1999

Estimation of magnitude and frequency of floods for streams in Puerto Rico: New empirical models

Flood-peak discharges and frequencies are presented for 57 gaged sites in Puerto Rico for recurrence intervals ranging from 2 to 500 years. The log-Pearson Type III distribution, the methodology recommended by the United States Interagency Committee on Water Data, was used to determine the magnitude and frequency of floods at the gaged sites having 10 to 43 years of record. A technique is presented for estimating flood-peak discharges at recurrence intervals ranging from 2 to 500 years for unregulated streams in Puerto Rico with contributing drainage areas ranging from 0.83 to 208 square miles. Loglinear multiple regression analyses, using climatic and basin characteristics and peak-discharge data from the 57 gaged sites, were used to construct regression equations to transfer the magnitude and frequency information from gaged to ungaged sites. The equations have contributing drainage area, depth-to-rock, and mean annual rainfall as the basin and climatic characteristics in estimating flood peak discharges. Examples are given to show a step-by-step procedure in calculating a 100-year flood at a gaged site, an ungaged site, a site near a gaged location, and a site between two gaged sites.

Puerto Rico↗
USGS2024

Diel temperature signals track seasonal shifts in localized groundwater contributions to headwater streamflow generation at network scale

Groundwater contributions to streamflow sustain aquatic ecosystem resilience; streams without significant groundwater inputs often have well-coupled air and water temperatures that degrade cold-water habitat during warm low flow periods. Widespread uncertainty in stream-groundwater connectivity across space and time has created disparate predictions of energy and nutrient fluxes across headwater networks, hindering predictions of cold-water habitat resilience under climate change scenarios. Recently, annual paired air and water temperature signals have been harnessed to indicate stream water thermal sensitivity and the dominance of deep versus shallow groundwater influence, although the utility of diel air–water temperature signal metrics for hydrologic inference has remained unexplored. Here we analyzed two consecutive years of locally paired, air–water temperature data from 47 headwater stream sites in the Catskill Mountains, New York, USA, and discovered characteristic seasonal patterns in diel temperature signal sinusoid metrics (amplitude ratio, phase lag, and mean ratio) driven by shifts in streamflow generation mechanisms and stream network position. Hydrologic interpretations of observed patterns were supported by stream heat budget model scenarios and additional analysis of paired air–water temperature data from two streams in Shenandoah National Park, Virginia, USA, with well characterized stream-groundwater connectivity. We found that within smaller tributaries, streamflow generation transitions from runoff to groundwater dominance were driven by hillslope drying during seasonal periods of lower precipitation. This was evidenced by significant correlations (p < 0.01) between daily water:air temperature signal amplitudes (non-linear decreases of ∼ 50 %) and derived base-flow index at 22 of the 28 sites, indicating enhanced local groundwater influence on streamflow promotes decoupling of diel air–water temperature signals. Additionally, ratios between daily water:air temperature signal means were lower in tributaries (∼0.68) when compared to main-stem (∼0.8) sites, increasing linearly throughout the observational period. In conceptual stream heat budget models, groundwater inflow had minimal effects on daily phase lags (∼0.2 hr), but increases in fractional groundwater discharge (0–50 %) depressed daily amplitude (∼20 % to 50 %) and mean ratios (∼15 %), supporting the sensitivity of daily metrics to interpreted changes in seasonal groundwater contributions to streamflow. During observational periods (i.e., April through October 2021 and 2022), significant differences (p < 0.01) between tributary and main-stem air–water metrics occurred when base-flow contributions were highest (∼0.93 vs. ∼ 0.68), as sites lower in the network had daily temperature metrics dominated by stream channel thermal inertia, rather than local groundwater connectivity, showing enhanced air–water diel signal coupling during warmer, drier periods. Divergent air temperature coupling across the network was interpreted as being driven by distance from local groundwater source zones, additional lateral groundwater inflows do not contribute a meaningful fraction to channel discharge lower in the network. Given the growing footprint of stream temperature observations, diel air–water temperature signals can provide distributed metrics sensitive to upstream groundwater discharge. Consequently, these metrics can support ongoing efforts by resource managers and researchers seeking to forecast the resilience of cold-water habitat to climate warming and changing precipitation regimes in mountain headwater streams.

New York, Virginia↗
USGS2025

Cancer risk and estimated lithium exposure in drinking groundwater in the US

Importance Lithium is a naturally occurring element in drinking water and is commonly used as a mood-stabilizing medication. Although clinical studies have reported associations between receiving lithium treatment and reduced cancer risk among patients with bipolar disorder, to our knowledge, the association between environmental lithium exposure and cancer risk has never been studied in the general population. Objectives To evaluate the association between exposure to lithium in drinking groundwater and cancer risk in the general population. Design, Setting, and Participants This cohort study included participants with electronic health record and residential address information but without cancer history at baseline from the All of Us Research Program between May 31, 2017, and June 30, 2022. Participants were followed up until February 15, 2023. Statistical analysis was performed from September 2023 through October 2024. Exposure Lithium concentration in groundwater, based on kriging interpolation of publicly available US Geological Survey data on lithium concentration for 4700 wells across the contiguous US between May 12, 1999, and November 6, 2018. Main Outcome and Measures The main outcome was cancer diagnosis or condition, obtained from electronic health records. Stratified Cox proportional hazards regression models were used to estimate the hazard ratios (HRs) and 95% CIs for risk of cancer overall and individual cancer types for increasing quintiles of the estimated lithium exposure in drinking groundwater, adjusting for socioeconomic, behavioral, and neighborhood-level variables. The analysis was further conducted in the western and eastern halves of the US and restricted to long-term residents living at their current address for at least 3 years. Results A total of 252 178 participants were included (median age, 52 years [IQR, 36-64 years]; 60.1% female). The median follow-up time was 3.6 years (IQR, 3.0-4.3 years), and 7573 incident cancer cases were identified. Higher estimated lithium exposure was consistently associated with reduced cancer risk. Compared with the first (lowest) quintile of lithium exposure, the HR for all cancers was 0.49 (95% CI, 0.31-0.78) for the fourth quintile and 0.29 (95% CI, 0.15-0.55) for the fifth quintile. These associations were found for all cancer types investigated in both females and males, among long-term residents, and in both western and eastern states. For example, for the fifth vs first quintile of lithium exposure for all cancers, the HR was 0.17 (95% CI, 0.07-0.42) in females and 0.13 (95% CI, 0.04-0.38) in males; for long-term residents, the HR was 0.32 (95% CI, 0.15-0.66) in females and 0.24 (95% CI, 0.11-0.52) in males; and the HR was 0.01 (95% CI, 0.00-0.09) in western states and 0.34 (95% CI, 0.21-0.57) in eastern states. Conclusions and Relevance In this cohort study of 252 178 participants, estimated lithium exposure in drinking groundwater was associated with reduced cancer risk. Given the sparse evidence and unknown mechanisms of this association, follow-up investigation is warranted.

contiguous United States↗
USGS2025

Disturbance is the primary determinant of food chain length when the top predator is constant

Food chain length (FCL) is a primary determinant of food web structure and is hypothesized to be influenced by habitat size, productivity, and disturbance. Understanding the environmental characteristics that determine food chain length can assist in understanding how food webs may be impacted due to changes in habitats and environmental characteristics. This study examines the impact of hydrologic disturbance on stream food webs when the top predator is constant. We analyzed FCL in less disturbed groundwater flashy streams and more disturbed runoff flashy streams using stable isotopes. Despite no difference in species richness or fish density, food chains in more disturbed streams had a lower FCL compared to food chains in more stable streams. Assemblage analysis showed that flow regime and drainage area significantly impacted individual species abundances. The more disturbed runoff flashy streams had higher proportions of primary consumer fish, such as the algivorous Campostoma sp. (Stonerollers), which likely drives the reduced FCL. Drainage area and land cover had non-significant relationships with FCL. Shifting community structure due to hydrologic variability likely leads to differences in diet of Micropterus dolomieu (Smallmouth Bass), and thus a difference in FCL.

Arkansas, Missouri, Oklahoma↗
USGS2026

CRESCENT earthquake dynamic rupture, earthquake cycle, and tsunami code verification platform

Physics-based simulations are critical for understanding natural hazards. The increasing complexity of numerical codes requires benchmark exercises to verify that different computational methods yield consistent results when solving the same governing equations. Here, we present an open-access web platform designed for the verification of earthquake dynamic rupture, seismic cycle, and tsunami simulations. The platform architecture utilizes a modular, serverless backend on Amazon Web Services (AWS) to provide scalable file processing and visualization. A lightweight static web application provides a secure interface for uploading and managing results, while the browser-based data visualization enables interactive analysis of time series and surface grid data. By using structured JavaScript Object Notation (JSON) text files to define benchmark structures, the system remains fully extensible, allowing the addition of new scenarios without modifying the underlying software logic. The platform hosts the "The Tsunami Problem Versions" (TTPV) 1 & 2, two benchmarks for 3D fully coupled earthquake dynamic rupture and tsunami generation, and provides a framework for earthquake cycle models. This community resource aims to build trust in numerical simulations and facilitate long-term collaborative code verification as modeling software continues to evolve.

Seismica↗
USGS2025

Origin of the high Pd/Pt ratio of the J-M Reef, Stillwater Complex Montana USA

The J-M Reef of the Stillwater Complex exhibits a high and consistent Pd/Pt ratio (~3.8). This ratio results from the equilibration of an immiscible sulfide liquid with a relatively high Pd/Pt silicate melt rather than an unusually Pd- and Pt-enriched parental melt. Numerical modeling suggests that the original silicate melt contained typical mantle-derived concentrations of Pd and Pt (~10–20 ppb Pd/Pt ~1). The partitioning of Pt and Pd between sulfide liquid and silicate melt alone cannot explain the consistently high Pd/Pt ratio across variable melt-to-sulfide mass ratios (R factors). Instead Pt-depletion caused by the early fractionation of Pt-alloy from S-undersaturated silicate magma likely established the high Pd/Pt signature. High Pd/Pt ratios can form through batch equilibration of sulfide liquid with silicate melt if partition coefficients are extremely high (>10⁶). Alternatively Pd enrichment may result from sulfide upgrading within the resident footwall mush under smaller partition coefficients (10⁴–10⁶) in this model the instantaneous R factors remain low (R ≈ 100–700). This limits the impact of Pt and Pd partitioning on sulfide composition and helps explain the Pd-enriched character of the J-M Reef.

Montana↗
USGS2025

Deformation mechanisms in quartz veins and shear zones elucidate the origin of gold mineralization at Pogo, Alaska

Pogo is a quartz vein hosted, ca. 8 Moz gold deposit. Although it has similarities to orogenic and magmatic-hydrothermal deposits, its origin remains enigmatic. Observations from surface exposures, underground workings, and drill core provide new constraints on quartz vein origins with implications for mineralization. Abundant, largely barren metamorphic segregation quartz veins are found throughout interior Alaska and share characteristics with mineralized quartz veins at Pogo. Pogo quartz veins show crystal plastic deformation fabrics consistent with greenschist to amphibolite facies metamorphism and a lack of internal shear. In contrast, gold in Pogo quartz veins occurs in brittle microfaults and fractures with sulphides such as arsenopyrite and pyrrhotite. Major and minor cataclastic shear zones also exist at Pogo. Cataclastic shear zones commonly cut the mineralized quartz veins and porphyroclasts are dominated by quartz with the same plastic deformation fabrics found in the major Pogo quartz veins. The porphyroclasts, and the carbonaceous clay-rich matrix they sit in, also contain gold indicating that the shear zones postdate quartz and earliest gold deposition. These observations suggest that competency contrasts between the quartz veins and their phyllosilicate-rich host rocks controlled transient permeability formed by late, preferential brittle deformation events localizing sulphide ± gold deposition.

Alaska↗
USGS2026

A practical framework for identifying genetic subpopulations and ESUs: Insights for IUCN assessments and broader management

Species conservation assessments evaluate extinction risk, and recovery potential, advancing species persistence through guiding resource prioritization and planning. Assessment frameworks, including the International Union for Conservation of Nature Red List and Green Status of Species, typically focus on species as a whole. Importantly, they do not routinely account for genetically distinct units or do not have standardized methods of unit delineation. This limits the representation of genetically distinct components, including adaptive genetic diversity that underpins long-term resilience and recovery. Incorporating standardized within-species units like subpopulations and Evolutionarily Significant Units (ESUs) into species assessments could help address this oversight. However, identifying and delineating such units remain challenging, particularly when molecular data are limited. Here, we propose a flexible framework that integrates molecular and non-molecular evidence to identify both subpopulations and ESUs across taxa, providing a practical tool to incorporate within-species diversity into conservation assessments.

BioScience↗
USGS2008

Geology and geochemistry of Jurassic plutonic rocks, Baboquivari Mountains, south-central Arizona

Among the plutons of the Jurassic magmatic arc segment in southern Arizona, northern Sonora, and southern California, two distinctive rock types dominate, occurring together in numerous mountain ranges: porphyritic, titanite-bearing hornblende-biotite granodiorite; and biotite leucogranite. These characteristic rocks are particularly well preserved and well exposed around Kitt Peak, in the Baboquivari Mountains, southern Arizona. The Middle to early Late Jurassic Kitt Peak Plutonic Suite (KPPS) comprises three units: Aguirre Peak Quartz Diorite (APQD;<170 Ma, >=165 Ma, U-Pb), chiefly hornblende mesodiorite; Kitt Peak Granodiorite (KPGD; 165 Ma), hornblende-biotite granodiorite and monzogranite; and Pavo Kug Granite (PKG; 159 Ma), equigranular biotite leucogranite. This range also hosts a fourth, regionally unusual, Jurassic plutonic unit: the Baboquivari Peak Perthite Granite, related to post-arc, Late Jurassic crustal extension. Binary variation diagrams, incompatible element patterns, and REE spectra indicate APQD and KPGD are consanguineous. APQD basaltic mesodiorite evidently approximates the parental magma for the APQD–KPGD series. More evolved rocks of this series evidently formed by combinations of fractional crystallization and mixing with or assimilation of crustal material. PKG is less directly related to APQD and KPGD than these two are to one another.

Arizona, Sonora↗
USGS2025

Comparing magmatism and hydrothermal alteration using magnetic modelling and stable isotopes at the Clementine porphyry copper prospect, Montana, USA

Recent mapping of the underexplored Clementine prospect in southwestern Montana has revealed evidence of hydrothermal alteration and mineralized breccia vein gossan interpreted to represent the upper expression of a Cretaceous, sediment-hosted copper porphyry system. The prospect is at the nexus of several Cretaceous granites, including the pre-mineralization Butte Granite and Big Hole Canyon plutons and local pre- to post-mineralization granite intrusions. Here, 3D magnetic inversions and stable isotope data are used to evaluate the spatial and genetic relationship of mineralization and Cretaceous magmatism. Magnetic inversions reveal a zone of high magnetic susceptibility beneath the prospect that is related to an exposure of unaltered, likely post-mineralization granite. The granite appears to connect to the Big Hole Canyon pluton at depth, suggesting it is a late-stage differentiate of the cooling, pre-mineralization pluton that may have been one of several potential heat sources for the mineral system at Clementine. This is supported by carbon isotope values of graphite that indicate organic material in shale of the Cretaceous Kootenai Formation was subject to metamorphism. Preliminary light hydrogen isotope values of actinolite hydrothermal orbs suggest formation following magmatic outgassing during thermal advection from an underlying intrusion.

Montana↗
USGS2024

The importance of nodal plane orientation diversity for earthquake focal mechanism stress inversions

Inversions of earthquake focal mechanisms are among the most accessible and reliable methods for determining crustal stress. However, the use of this method varies widely, and assumptions that underpin it are often violated, potentially compromising stress estimates. We investigate the consequences of violating the little-studied assumption that the focal mechanisms have diverse orientations. Our approach is to employ data-informed synthetic mechanisms, with nodal plane orientations defined by recent earthquake lineaments in the Midland Basin, western Texas, and rakes consistent with slip in the mapped stress field. Using both the traditional stress inversion method that assumes constant shear stress magnitudes on the causative faults as well as a recently published variable shear stress method, we show that low fault plane diversity can cause maximum horizontal stress ( S Hmax ) orientation and relative principal stress magnitude (faulting regime) estimates to differ markedly from the true values. This problem is compounded for catalogs with even modest amounts of noise (≤15°) or few (e.g., 20) mechanisms. Significantly, traditional approaches for quantifying uncertainty such as the bootstrap can severely underestimate the true uncertainty under these circumstances. To remedy this, we provide simple tools to quantify nodal plane orientation diversity and stress inversion reliability.

Geological Society of London Special Publications↗
USGS2024

Reduced injection rates and shallower depths mitigated induced seismicity in Oklahoma

The proximity of wastewater disposal to the Precambrian basement is a critical factor influencing induced earthquake rates in the Central United States, but the impact of reducing injection depths has not been widely demonstrated. Beginning in 2015, state regulatory efforts in Oklahoma and Kansas mandated that wells injecting into the lower Arbuckle Group, a basal sedimentary unit, be backfilled with cement (i.e. “plugged back”) so that they inject into shallower formations. This plug back activity gives us a unique opportunity to investigate the relationship between injection depth and induced seismicity rate. To evaluate the impact that decreased injection rates and plug backs had on the seismicity rates, we create a suite of rate-state earthquake models. Observed seismicity rates are best fit when only lower Arbuckle volumes are considered, suggesting the lower Arbuckle injectors were primarily responsible for the seismicity and that plug backs were effective at isolating the injected volumes to shallower formations. Our models demonstrate that if these wells had not been plugged back, seismicity rates would be multiple times larger than they are today. We find that the combination of well plug backs and injection volume decreases can be effective strategies for reducing induced seismicity rates.

Oklahoma↗
USGS2024

Uncertainty in ground-motion-to-intensity conversions significantly affects earthquake early warning alert regions

We examine how the choice of ground‐motion‐to‐intensity conversion equations (GMICEs) in earthquake early warning (EEW) systems affects resulting alert regions. We find that existing GMICEs can underestimate observed shaking at short rupture distances or overestimate the extent of low‐intensity shaking. Updated GMICEs that remove these biases would improve the accuracy of alert regions for the ShakeAlert EEW system for the West Coast of the United States. ShakeAlert uses ground‐motion prediction equations (GMPEs), which calculate spatial distributions of peak ground acceleration (PGA) and peak ground velocity (PGV) from earthquake source estimates, combined with GMICEs to translate GMPE output into modified Mercalli intensity (MMI). We find significant epistemic uncertainty in alert distances; near‐source MMI estimates from different GMICEs can differ by over 1 MMI unit, and MMI extents used for public EEW alerts can differ by hundreds of kilometers for larger magnitude earthquakes ( M ∼6.5+). We use a catalog of “Did You Feel It?” shaking reports to evaluate how well GMICEs predict observed shaking. Our preferred GMICE is the one that computes MMI using PGV for high intensities and transitions to using PGA for nondamaging intensities. These results motivate updating GMICE relationships more generally, including in ShakeMap applications.

The Seismic Record↗
USGS2025

Fingerprinting magmatic REE deposit sources with zircon petrochronology

Carbonatites and associated alkaline silicate rocks are of considerable economic interest due to their enrichments in rare earth elements. The petrogenesis and source(s) of these complexes, however, are poorly understood. Models propose either mantle plume-derived carbon-rich melts or a mantle source enriched by subduction-related metasomatism. We use zircon trace elements to fingerprint carbonatite-alkaline complex sources, focusing on the economically significant 1.4 Ga Mountain Pass intrusive suite (MPIS). Autocrystic zircon from MPIS alkaline silicate rocks are enriched in Th and U relative to Nb suggesting a subduction influenced, oxidized source region. Lower Sc/Yb, higher Ti concentrations, and the absence of Eu anomalies in autocrystic zircon suggest derivation from a less hydrous, hotter, and deeper mantle source relative to arc-related 1.8–1.6 Ga inherited zircon. These data are best explained by syn- to post-collisional Mesoproterozoic reactivation of a subduction-metasomatized mantle source. The source and tectonic setting of the MPIS contrasts with Mesoproterozoic Gifford Creek and Bayan Obo carbonatites which exhibit plume-like compositions and are associated with rifts, suggesting different sources and tectonic settings for these economically significant deposits. We demonstrate that zircon petrochronology is a robust method for distinguishing carbonatite sources and can inform more targeted exploration strategies for critical mineral resources.

Conference Paper↗