SEARCH · USGS Science

Results for “Journal of Environmental Quality”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,662 records · Page 69Linked to original sources

USGS2025

Airborne geophysics for geologic mapping of critical mineral systems in the United States southern midcontinent

The increased demand for clean energy technology and a significant reliance on foreign supply chains have given impetus to understanding critical mineral systems and locating potential resources within the United States. At least thirteen critical mineral-bearing systems have been identified throughout the U.S. southern Midcontinent (Hofstra and Kreiner, 2020) but much of the region’s geologic framework is concealed by vegetation and sedimentary cover that hinder traditional geologic mapping efforts. Airborne geophysical data provide an effective way to overcome these obstacles and to provide additional insight into the deeper structures that underlie shallow mineralization. However, legacy airborne magnetic and radiometric data were collected using now-outdated instruments and methods, inconsistent survey parameters, and large flight-line spacings resulting in low-resolution data that present challenges to regional-scale study and interpretation. Over the last decade, the U.S. Geological Survey Earth Mapping Resources Initiative (EMRI) and National Cooperative Geologic Mapping Program have conducted a series of high-resolution airborne magnetic and radiometric surveys across the southern Midcontinent (Fig. 1) as part of an effort to improve understanding of the geophysical framework and natural resource potential in the region. These surveys are designed using modern survey methods and instruments with consistent parameters for flight-line spacing and flight height relative to magnetic sources. The EMRI airborne surveys are planned in collaboration with State geological surveys based on focus areas (Dicken et al., 2022) according to the presence of or potential for critical mineral deposits. High-resolution airborne magnetic and radiometric data cover focus areas such as the southeast Missouri iron metallogenic province and South-Central iron-oxide-apatite (IOA) – iron-oxide-copper-gold (IOCG) province, the Magnet Cove alkaline-carbonatite complex, the Midwest Permian ultramafic dike district, the Illinois-Kentucky fluorspar district, and several Mississippi Valley-type lead-zinc deposits and districts (Fig. 1). These focus areas represent known deposits or prospective host systems of critical minerals including rare earth elements (REEs), platinum-group elements (PGEs), cobalt, lithium, fluorspar, niobium, titanium, vanadium, lead, zinc, gallium, germanium, and many more. Other significant geologic and geophysical features covered include the Reelfoot rift, the New Madrid seismic zone, the Illinois basin, the Arkoma basin, the South-Central magnetic lineament, and the Kentucky-Tennessee magnetic anomaly (Fig. 1). This presentation focuses on new airborne magnetic and radiometric data with continuous coverage across parts of six states, preliminary interpretations, examples of geologic mapping applications, and discussion of newly discovered magnetic anomalies and follow-up investigations.

Conference Paper↗

USGS2024

A benchmark dataset and workflow for landslide susceptibility zonation

Landslide susceptibility shows the spatial likelihood of landslide occurrence in a specific geographical area and is a relevant tool for mitigating the impact of landslides worldwide. As such, it is the subject of countless scientific studies. Many methods exist for generating a susceptibility map, mostly falling under the definition of statistical or machine learning. These models try to solve a classification problem: given a collection of spatial variables, and their combination associated with landslide presence or absence, a model should be trained, tested to reproduce the target outcome, and eventually applied to unseen data. Contrary to many fields of science that use machine learning for specific tasks, no reference data exist to assess the performance of a given method for landslide susceptibility. Here, we propose a benchmark dataset consisting of 7360 slope units encompassing an area of about 4,100 km 2 "> 4,100 km 2 in Central Italy. Using the dataset, we tried to answer two open questions in landslide research: (1) what effect does the human variability have in creating susceptibility models; (2) how can we develop a reproducible workflow for allowing meaningful model comparisons within the landslide susceptibility research community. With these questions in mind, we released a preliminary version of the dataset, along with a “call for collaboration,” aimed at collecting different calculations using the proposed data, and leaving the freedom of implementation to the respondents. Contributions were different in many respects, including classification methods, use of predictors, implementation of training/validation, and performance assessment. That feedback suggested refining the initial dataset, and constraining the implementation workflow. This resulted in a final benchmark dataset and landslide susceptibility maps obtained with many classification methods. Values of area under the receiver operating characteristic curve obtained with the final benchmark dataset were rather similar, as an effect of constraints on training, cross–validation, and use of data. Brier score results show larger variability, instead, ascribed to different model predictive abilities. Correlation plots show similarities between results of different methods applied by the same group, ascribed to a residual implementation dependence. We stress that the experiment did not intend to select the “best” method but only to establish a first benchmark dataset and workflow, that may be useful as a standard reference for calculations by other scholars. The experiment, to our knowledge, is the first of its kind for landslide susceptibility modeling. The data and workflow presented here comparatively assess the performance of independent methods for landslide susceptibility and we suggest the benchmark approach as a best practice for quantitative research in geosciences.

Earth-Science Reviews↗

USGS2024

An orbital comparison of a late mantling unit on Aeolis Mons with other erosion-resistant strata explored by MSL in Gale crater, Mars

Previous investigations along Curiosity 's traverse in Gale crater have explored the relationship between orbital and in situ observations. This work aims to better understand the geologic environment of units only observable from orbit and compare them to the properties of units examined by Curiosity . Here, we map an erosion-resistant dark-toned mantling unit that overlies the modern topography of Aeolis Mons (informally known as Mt. Sharp) and compare this unit to two other previously mapped dark-toned resistant units, the marker band and the mound skirting unit (MSU), that have been inferred to represent different geologic environments (lacustrine and aeolian, respectively). Visible to short wave infrared spectra from the Compact Reconnaissance Imaging Spectrometer for Mars and visual images from the High Resolution Imaging Science Experiment and Context Cameras aboard the Mars Reconnaissance Orbiter are used for this comparison. Spectral data suggest a mafic composition with minor alteration, although the composition varies more with location around Mt. Sharp rather than between units. Morphologically, the mantling unit has strong similarities to the marker band based on their consistent low-albedo, erosion-resistance, and smooth appearance, contrasting with the highly variable surface texture of the MSU. We hypothesize that all three units had a similar sediment source but experienced aqueous alteration at different times: early ubiquitous cementation in a surface aqueous environment in the mantling unit and marker band versus patchy late diagenesis in the MSU. If true, these results suggest that water activity continued within the Gale crater long after the erosion of Mt. Sharp.

JGR Planets↗

USGS2026

Wavelet Inversion for SliP (WISP): Open-source earthquake slip modeling software

Models of the spatiotemporal evolution of earthquake slip, termed finite-fault models, are a critical component of rapid earthquake and tsunami response, earthquake forecasting, seismic ground-motion estimates, and studies of earthquake kinematics. Here, we detail a newly released finite-fault modeling software, Wavelet Inversion for SliP (WISP), in use at the U.S. Geological Survey’s National Earthquake Information Center (NEIC) and available to the public. WISP version 1.1.0 allows inversion of teleseismic body and surface waves, as well as local strong-motion, static and dynamic Global Navigation Satellite System, and satellite imagery (e.g., Interferometric Synthetic Aperture Radar) observations on single or multiple planar fault segments. The software is used in NEIC rapid response of earthquakes M w ≥ 7, generally resulting in a published model within the first few hours after the event origin time. The rupture location and dimensions are then used as inputs to downstream products to estimate earthquake shaking, predict loss, and model the likelihood of secondary hazards, namely landslides and liquefaction. WISP is also used in research studies to evaluate the characteristics of complex ruptures including multifault ruptures and earthquake doublets, among others. The WISP version 1.1.0 software release is composed of Python-wrapped FORTRAN code to accomplish the inversion procedure. A simple command line interface facilitates ease of use even for those with only a cursory knowledge of Python scripting. WISP version 1.1.0 includes a Jupyter Notebook tutorial demonstrating use of the software for modeling the 2015 M w 8.3 Illapel, Chile, earthquake. In parallel with the tutorial, we demonstrate the typical usage of the WISP software using the M w 8.3 Illapel earthquake example here.

Seismological Research Letters↗

USGS2025

Mineral commodity summaries 2025

Introduction Each mineral commodity chapter of the 2025 edition of the U.S. Geological Survey (USGS) Mineral Commodity Summaries (MCS) includes information on events, trends, and issues for each mineral commodity as well as discussions and tabular presentations on domestic industry structure, Government programs, tariffs, 5-year salient statistics, and world production, reserves, and resources. The MCS is the earliest comprehensive source of 2024 mineral production data for the world. More than 90 individual minerals and materials are covered by two-page synopses. Abbreviations and units of measure and definitions of selected terms used in the report are in Appendix A and Appendix B, respectively. Reserves and resources information is in Appendix C, which includes “Part A—Resource and Reserve Classification for Minerals” and “Part B—Sources of Reserves Data.” A directory of USGS minerals information country specialists and their responsibilities is in Appendix D. The USGS continually strives to improve the value of its publications to users. Constructive comments and suggestions by readers of the MCS 2025 are welcomed.

Mineral Commodity Summaries↗

USGS2026

Mineral commodity summaries 2026

Introduction Each mineral commodity chapter of the 2026 edition of the U.S. Geological Survey (USGS) Mineral Commodity Summaries (MCS) includes information on events, trends, and issues for each mineral commodity as well as discussions and tabular presentations on domestic industry structure, Government programs, tariffs, 5-year salient statistics, and world production, reserves, and resources. The MCS is the earliest comprehensive source of 2025 mineral production data for the world. More than 90 individual minerals and materials are covered by two-page synopses. Abbreviations and units of measure and definitions of selected terms used in the report are in Appendix A and Appendix B, respectively. Reserves and resources information is in Appendix C, which includes “Part A—Resource and Reserve Classification for Minerals” and “Part B—Sources of Reserves Data.” A directory of USGS minerals information country specialists and their responsibilities is in Appendix D. The USGS continually strives to improve the value of its publications to users. Constructive comments and suggestions by readers of the 2026 MCS are welcomed.

Mineral Commodity Summaries↗

USGS2025

Relationship between peak and cumulative ground motions from 49 Mw 3-6 earthquakes in the San Francisco Bay Area

We examine the relationship between peak ground velocity (PGV) and cumulative absolute displacement (CAD) for 49 M w 3 – 6 earthquakes in the San Francisco Bay Area (SFBA) and gain insight into the spatiotemporal partitioning of seismic energy in ground motion records with respect to source, path, and site effects. PGV and CAD are positively correlated, but there can be large deviations from the average trend. For example, ground motion records with either very long duration resonance or short pulse-like motions will have higher or lower CAD, respectively, but could have very similar PGV. We perform principal component analysis (PCA) on PGV-CAD for >7000 records in the SFBA with the goal of investigating what influences positive or negative anomalies in cumulative motions. PCA rotates the PGV-CAD datapoints into two principal components, where the one with the larger variance, which we call the “primary intensity” represents mostly the distance-dependence of ground motion amplitudes. The other principal component, which we call the “excess motion”, represents the deviation from cumulative motions that would be typical for a ground motion record with a given PGV. The excess motion will be positive in the case of records with long duration ringing and will be negative for short duration pulse-like ground motions. We find that excess motion is generally positive at sites in sedimentary basins and in soft sediments around the SF Bay. Excess motion is generally negative in the very near field, as well as at sites on hard bedrock. We discuss the findings here in terms of implications for seismic hazard applications and other wave propagation phenomena.

California↗

USGS2025

U.S. Geological Survey science strategy to address white-nose syndrome and bat health in 2025–2029

Since its discovery in 2006, the fungal disease known as white-nose syndrome (WNS) has killed millions of bats. Of the 47 bat species native to the conterminous United States, Alaska, Hawaii, and Canada, 12 have been affected by WNS, including 3 endangered species and 1 proposed endangered species. WNS has also been detected in 40 States and 9 Canadian Provinces. U.S. Geological Survey (USGS) scientists have been critical in identifying the causal fungus for WNS ( Pseudogymnoascus destructans [Pd]), characterizing the effects of WNS, and tracking the spread of Pd in many bat populations in North America. The mission of the USGS WNS and Bat Health Science Team is to deliver integrated science in order to build resiliency into free-ranging bat populations through more effective WNS management, build capacity for bat health science, and enhance bat health information sharing across USGS science centers and cooperative research units as well as with stakeholders. The USGS can play an important role in supporting regional and national capacity building by providing resources and guidance to local, State, and Tribal management entities and by providing tools to enhance disease management. The USGS Ecosystems Mission Area’s Biological Threats and Invasive Species Research Program is the lead Federal program for free-ranging wildlife disease research and surveillance. As of 2024, guided by the science priorities set by the WNS Steering Committee, USGS scientists are engaged in a nationwide response to WNS. This work is done in close coordination with our partners at the U.S. Fish and Wildlife Service, National Park Service, Bureau of Land Management, U.S. Forest Service of the U.S. Department of Agriculture, U.S. Department of Defense, as well as State and Tribal agencies. In addition to conducting WNS research, the USGS is mapping the spread of WNS and coordinating the North American Bat Monitoring Program (NABat) to understand how WNS and other stressors affect the status and trends of native bats across their range. The USGS is supporting the national WNS response through four science goals: (1) provide situational awareness on the health of bat populations; (2) conduct ecological studies of bats along the gradient of disease vulnerability; (3) contribute actionable science to enhance the resiliency of bat populations; and (4) implement an adaptive, holistic approach to bat health.

Circular↗

USGS2025

Lithium from magma to mine in an early Yellowstone hotspot caldera

Renewable energy technologies rely on the extraction of metals not historically in high demand, such as lithium (Li), for which ore deposit models are incompletely understood. One of the world’s largest Li deposits is hosted in lake sediments of the 16.4 Ma McDermitt caldera, which formed during the early stages of Yellowstone hotspot volcanism in the western United States. Eruptive and posteruptive mobility of Li are major challenges in elucidating deposit formation. Melt inclusions preserved in quartz crystals provide a means to assess pre-eruptive magmatic Li contents. Concentrations of Li determined by ion microprobe for melt inclusions in a McDermitt rhyolite lava are 400−1350 ppm, compared to 20−70 ppm Li in matrix rhyolite glasses. Synthesis with melt inclusion data for eight additional calderas demonstrates a recurrence of Li-rich rhyolitic magmas (200−2000 ppm Li) in the western part of the Yellowstone hotspot track. However, unlike the multicyclic caldera complexes with overlapping fault networks that may have compromised Li retention, the McDermitt caldera remained a closed hydrologic system throughout its evolution. Modeling indicates 100 km 3 of resurgent magma could yield 25−150 Mt Li in a magmatic fluid and supports accumulation of Li-rich magmatic fluid in a closed intracaldera lake, followed by evaporative concentration and sequestration of Li within clay minerals to generate the McDermitt deposit.

Idaho, Nevada, Oregon, Wyoming↗

USGS2026

Postglacial ecosystem development of a hydrothermal landscape in Yellowstone National Park

The Yellowstone geo-ecosystem has been the subject of much research, but the ecological history of the Yellowstone Plateau volcanic field and its iconic geyser basins is less known. In this investigation, paleoenvironmental analyses of sediment cores from lakes in Lower Geyser Basin were compared with regional records and paleoclimate model simulations to reconstruct the vegetation, wildfire, limnology, hydrothermal dynamics, and climate drivers since deglaciation, 15,000 to 14,000 years ago. Pollen data from Lower Geyser Basin lakes reveal the strong influence of infertile rhyolitic soils on vegetation history: an initial late-glacial steppe was replaced by lodgepole pine forest from 12.8 to 11.0 ka, with little change in forest composition or cover thereafter despite changing climate. This stability contrasts with the more dynamic vegetation response on nonrhyolite substrates in the Yellowstone region where nutrient and moisture availability is greater. Highest wildfire activity and low lake nutrient levels in Lower Geyser Basin occurred from 12 to 4 ka, when summers were substantially warmer and drier and fire-inducing vapor pressure deficits were 29 to 56% higher. The hydrothermal history, inferred from sedimentary arsenic and cesium abundances, was spatially and temporally variable but lake-forming hydrothermal events align with periods of abundant moisture. Thus, long-term changes in wildfire, limnology, and, to some extent, hydrothermal activity were governed by insolation-driven climate variations, whereas the vegetation response was muted and constrained by geologic processes. These findings suggest that warmer, drier conditions in the future could result in less hydrothermal activity yet little change in forest cover across the Yellowstone Plateau volcanic field despite more wildfires.

Idaho, Montana, Wyoming↗

USGS2025

Tracing fluid evolution during batholith and deposit formation: In-situ and bulk rock Mo isotopes from the Questa porphyry Mo deposit, New Mexico, USA

The Oligocene Questa porphyry Mo deposit, located in New Mexico, USA, is a well-characterized system that formed over 400 kyr, punctuating the 8 Myr formation of the Questa-Latir batholith. Detailed sampling of the intrusions and minerals associated with magmatichydrothermal mineralization, as well as bookending barren intrusions, allows for examination of the Mo isotope record of the deposit and its host batholith. Through this established chronologic framework, such data can yield insights into the metallogenesis of deposit formation and long-term magma genesis evolution for upper crustal magmatism. The Mo isotopic compositions of barren intrusions yield a younging trend towards lighter δ 98 Mo. However, intrusions associated with Mo deposit formation yield heavier δ 98 Mo relative to barren magmatism, by up to 2.0‰. This indicates that they were affected by an isotopically heavy and Mo-rich fluid. The in-situ data from molybdenite covers the range observed in mineralized bulk rocks. However, individual samples consistently yield Mo isotopic variations of 0.4‰ to 0.5‰. The consistency of the δ 98 Mo variation in in-situ analyses of molybdenite within hand samples and individual mineral grains is likely the result of MoS 2 crystallization in a closed system and reflects an almost complete precipitation of Mo.

New Mexico↗

USGS2026

Global pegmatite-hosted lithium, cesium, and rubidium resources: A dataset for grade and tonnage modeling

Quantitative mineral resource assessments of potential undiscovered deposits can inform future mineral supply scenarios, but their accuracy is conditional on building robust grade and tonnage models of known deposits. This study presents an up-to-date global compilation and analysis of recently discovered and original, in-situ pegmatite-hosted Li, Cs, and Rb resources prior to historic production. Our analysis yields a median tonnage of 21.2 million tons (Mt) and grade of 1.12% Li 2 O, respectively, for global Li pegmatite deposits ( n = 73). The grades and tonnages of Li pegmatite resources vary depending on the age of the bedrock host domain, pegmatite crystallization age, and primary ore mineralogy. Lithium pegmatite resources hosted in Archean to transitional Archean-Paleoproterozoic domains have the largest median tonnage (29.8 Mt; n = 38), and those hosted in Paleoproterozoic to Mesoproterozoic domains have smaller median tonnages (6.5 Mt; n = 16). Cesium deposits where pollucite is the primary ore mineral have a bimodal grade distribution, with modes of 2.40 and 0.035 wt% Cs 2 O for high- and low-grade deposits, respectively, while Rb deposits are more unimodal with a median grade of 0.247 wt% Rb 2 O. Pegmatite-hosted Cs and Rb resources have median tonnages of 7.6 and 6.3 Mt, respectively. Covariation between ore mineralogy and the degree of crustal enrichment in pegmatite-hosted deposits is diagnostic of petrogenetic differences, including melt source characteristics, magma evolution, or variable degrees of volatile solubility. The Li pegmatite compilation is suitable for fitting robust numerical models to support quantitative assessments. More well-defined Rb and Cs pegmatite resources are required for quantitative assessments, but these data provide useful information about original in-place resources for framing supply discussions.

Ore Geology Reviews↗

USGS2025

Improved prediction of postfire debris flows through rainfall anomaly maps

Predicting where runoff-generated debris flows might occur during rainfall on steep, recently burned terrain is challenging. Studies of mass-movement processes in unburned areas indicate that event locations are well-predicted by rainfall anomaly, R* , in which peak observed rainfall is normalized by local rainfall climatology. Here, we use remote and field methods to map debris flows triggered within the 2020 Dolan Fire burn area in coastal California, demonstrate that a short-duration R* metric predicts debris-flow occurrence more effectively than absolute peak intensity or longer-duration rainfall metrics, and show that incorporating an R* criterion into an existing debris-flow likelihood model can reduce false positive predictions and improve accuracy. We test R * at three other climatically distinct fires in California, demonstrating its utility for mapping likely debris-flow locations in different climates. We also consider how R* can benefit postfire debris-flow prediction given recent increases in climatological variability within individual burn perimeters.

Callifornia↗

USGS2025

Strong shaking from past Cascadia Subduction Zone earthquakes encoded in coastal landforms

Strong earthquakes along subduction zones are often devastating events, but sparse records along some tectonic margins limit our understanding of seismic hazards. Constraining shaking intensities is critical, especially in subduction zones with infrequent but large-magnitude earthquakes like the Cascadia Subduction Zone (CSZ), where the lack of recorded ground motions has led to uncertainty in the severity and potential impacts of future earthquakes. Here we fill this observational gap with a novel inventory of quantitative estimates of past shaking intensities from geotechnical modeling of coastal landforms. One hundred fifty-four deep-seated landslides and 65 fragile geologic features constrain minimum and maximum peak ground accelerations, respectively. These estimates are broadly consistent with model predictions of M9 ruptures, suggesting strong shaking of 0.4–0.8 g during past CSZ earthquakes. Local discrepancies between our geologic shaking constraints and earthquake simulations may inform past rupture behavior, leading to better predictions of shaking intensity for future earthquakes.

California, Oregon, Washington↗

USGS1977

Leachate plumes in ground water from Babylon and Islip landfills, Long Island, New York

Landfills operated by the Towns of Babylon and Islip in southwest and central Suffolk County contain urban refuse, incinerated garbage, and scavenger (cesspool) waste; some industrial refuse is deposited at the Babylon site. The Islip landfill was started in 1933, the Babylon landfill in 1947. The landfills are in contact with and discharge leachate into the highly permeable upper glacial aquifer (hydraulic conducitvity 190 to 500 ft/d [feet per day]). The aquifer is 74 feet thick at the Babylon landfill and 170 feet thick at the Islip landfill. The leachate-enriched water occupies the entire thickness of the aquifer beneath both landfills, but hydrologic boundaries retard downward migration of the plumes to deeper aquifers. The Babylon plume is 1,900 feet wide at the landfill and narrows to about 700 feet near its terminus 10,000 feet from the landfill. The Islip plume is 1,400 feet wide at the landfill and narrows to 500 feet near its terminus 5,000 feet from the landfill. Hydrochemical maps and sections show the distribution of the major chemical constituents of the plumes. The most highly leachate-enriched ground water obtained was from the Babylon site; it contained 860 mg/L (milligrams per liter) sodium, 110 mg/L potassium, 565 mg/L calcium, 100 mg/L magnesium, 2,700 mg/L bicarbonate, and 1,300 mg/L chloride. Sulfate was notably absent or in low concentration in most parts of both plumes. Nitrogen in plume water was mostly in the form NH4+, and concentrations as high as 90 mg/L were found; concentrations of nitrogen as N in the plume were less than 10 mg/L. As much as 440 mg/L iron and 190 mg/L manganese were found in the leachate-enriched water. Samples were also tested for arsenic, boron, cadmium, cobalt, chromium, copper, mercury, nickel, lead, selenium, strontium, and zinc. Boron was more or less ubiquitous and was found in concentrations as high as 2 mg/L. Organic carbon was found in concentrations as high as 2,250 mg/L in the most highly leachate-enriched water but attenuated rapidly to less than 20 mg/L. Dissolved-solids concentrations near the landfills were between 400 and 3,000 mg/L at Babylon and between 500 and 1,500 mg/L at Islip. Ground-water temperatures near the landfills exceed those in ambient water by as much as 7°C at Babylon and 16°C at Islip. Heat contributed by the landfills was mostly dissipated within 0.4 mi of the landfill, butat Islip, the warm leachate-enriched water extended 0.5 mi downgradient. The entrance of leachate into the less dense ground water as pulsations after rainfall may explain the presence of high leachate enrichment near the bottom of the aquifer. A comparison of the physical characteristics of leachate-enriched ground water with those of ambient water suggests that the downward movement of leachate results from its greater density. Simulation of the movement and dispersion of the Babylon plume with a mathematical dispersion model indicated the coefficient of longitudinal dispersion to be about 60 ft 2 /d (feet squared per day) and the ground-water velocity to be 1 ft/d. However, the velocity determined from the hydraulic gradient and public-supply wells in the area was 4 ft/d, which would cause a plume four times as long as that predicted by the mathmatical dispersion model. At the Islip site, the plume was one-third the length calculated on the basis of the age of the landfill. The shortness of the plumes has not been explained; it may be a result of the leachate's having been too dilute to form a plume during the early years of the landfills.

New York↗

USGS2024

U.S. Geological Survey Earthquake Hazards Program decadal science strategy, 2024–33

Executive Summary Earthquakes represent one of our Nation’s most significant and costly natural hazards, with estimated annual loses from earthquakes close to $15 billion in 2023. Over the past two centuries, 37 U.S. States have experienced an earthquake exceeding a magnitude of 5, and 50 percent of States have a significant potential for future damaging shaking; these statistics speak to the need for nationwide interest and investment in earthquake hazard characterization and risk reduction. Authorized under the Earthquake Hazards Reduction Authorization Act, the U.S. Geological Survey (USGS) Earthquake Hazards Program (EHP) provides the scientific information, situational awareness, and knowledge necessary to reduce deaths, injuries, and economic losses from earthquakes and earthquake-induced tsunamis, landslides, and soil liquefaction. The EHP supports activities in three focused topical areas: (1) earthquake monitoring, (2) hazard assessment, and (3) applied research, using the results of each—and the coordination among them—to further support risk translation and communication in regions at risk nationwide. For earthquake monitoring, the Advanced National Seismic System (ANSS), a cooperative effort of USGS networks, university partner regional seismic networks, and real-time geodetic networks, collects and analyzes data on earthquakes; issues timely, reliable notifications of their occurrence and impacts; and provides data for earthquake research, hazard, and risk assessment as a foundation for building an earthquake-resilient Nation. The USGS-operated ShakeAlert Earthquake Early Warning system is a recent addition to EHP’s ANSS infrastructure. In the realm of earthquake hazard assessment, the EHP contributes to earthquake risk mitigation strategies by developing the National Seismic Hazard Model and maps, and other related products, that describe the likelihood and potential effects of earthquakes nationwide, especially in the urban areas of highest risk. The EHP also conducts research on the causes, characteristics, and effects of earthquakes and prioritizes work that directly increases the accuracy and precision of earthquake hazards assessments, earthquake forecasts, and earthquake monitoring and situational-awareness products and that supports the Nation’s earthquake mitigation practices. Bridging the EHP’s efforts across research, hazard assessments, and earthquake monitoring is a broad and comprehensive collection of earthquake information products, including the National Seismic Hazard Model, ShakeAlert, and other products describing impact, such as ShakeMap and PAGER (Prompt Assessment of Global Earthquakes for Response), which have been developed and integrated into EHP’s real-time monitoring systems. EHP funds external partners to carry out many important collaborative activities through an active external grants program—one of the largest in the USGS—and through cooperative agreements with other partners such as the university-operated regional seismic networks, funded as part of the ANSS. To continue its support of earthquake hazard characterization and risk reduction, the EHP aims to strengthen its foundational products and practices while positioning itself to respond to the evolving needs of the Nation and follow best practices of the scientific community. This document describes a strategy for the program to ensure it can meet these demands. The foundational priorities outlined in this strategy represent those activities that remain critical to the core functionality of the program and those that can be supported under current fiscal year 2024-level appropriations. Priorities described as aspirational are important for future growth, and to maintain the program’s position as a leading global resource in earthquake science, but would require increases in appropriated funding to be fully realized. Across the program’s portfolio of activities, several major themes have been identified as the most critical activities to advance EHP science over the coming decade. Together, these activities provide the framework necessary to integrate critical hazard characterization and risk reduction activities across the program. They provide the structure for research to advance the understanding of where, when, and why earthquakes occur and how we can use improved knowledge to drive short-term and actionable forecasts of seismic activity. They expand the usefulness of critical earthquake products and advance the sophistication of those products to keep pace with the rapidly evolving needs of an ever-expanding user base while maintaining the position of the USGS as a global leader in earthquake science. Focus on system-level science. Establish an automated earthquake-processing pipeline. Enhance the accuracy and reliability of the ShakeAlert Earthquake Early Warning system and plan for extension to other regions. Implement time-dependent earthquake forecasting. Develop physically realistic models. Expand computational capacity. This science strategy is organized into three primary sections. The first section provides an overview of the EHP and its budget, governance, and program council. Readers familiar with the program may wish to focus on the second section, which describes the core of the science strategy, including priorities across each of the EHP’s major program activities in monitoring, hazard assessment, and targeted research. The third section outlines science priorities that cut across program activities, including those involving collaborations external to the EHP.

American Samoa, Guam, Northern Marianas Islands, U↗

USGS2025

Igneous differentiation in arcs

This chapter assesses crystallization-differentiation of convergent-margin magmas both from the standpoints of natural rock suites and from laboratory phase-equilibrium experiments. Crystallization-differentiation of dry subalkaline and alkaline basalts are summarized at the beginning for reference and context. Experiments crystallizing arc basalts through andesites at middle to upper crustal pressures (100–400 MPa), saturated either with H 2 O or mixed H 2 O-CO 2 , produce liquid lines of descent similar to arc igneous suites, with higher f O 2 s giving results closer to calc-alkaline magma series. Early enrichment of melts in SiO 2 with little or no enrichment in FeO* (total Fe as FeO) comes about by the combined effects of H 2 O suppressing crystallization of plagioclase relative to olivine, shifting plagioclase to more calcic – therefore lower SiO 2 – compositions, and suppressing silicate minerals to greater extents than Fe-rich spinel. Also contributing to SiO 2 enrichment is crystallization near the liquidus of higher proportions of olivine relative to calcic clinopyroxene, compared with MORB, but it is unclear if this is brought about by H 2 O or elevated pressure. Exhumed arc terranes, however, generally fail to expose evidence for primitive liquids having reached and differentiated at middle to upper crustal depths, at least in any abundance. Instead, most differentiation in arcs appears to take place in the lower crust or near the crust-mantle transition, and although hydrous experiments at such pressures produce liquids broadly similar to arc suites in their SiO 2 enrichment, the experimental liquids differ in becoming peraluminous at early stages of differentiation while melts have basaltic, basaltic andesitic, or andesitic SiO 2 concentrations. Crystallization of abundant calcic clinopyroxene close to the liquidus causes the high-pressure experimental liquids to become peraluminous with only small increases in SiO 2 concentration. In contrast, natural arc volcanic and plutonic suites become gradually less metaluminous with increasing SiO 2 , crossing to peraluminous at dacitic to rhyolitic SiO 2 concentrations. Hybridism, encompassing magma-mixing, assimilation of old crust, and incorporation of near-solidus melts from only shortly older deep, antecedent arc intrusions may accompany natural crystallization-differentiation, preventing melts from attaining peraluminous compositions at low-SiO 2 concentrations. If so, then nearly all arc suites are to substantial extents composites, rather than preserving simple, progressive liquid lines of descent. Some comments and assessments are also made on the participation of garnet in arc differentiation and on the production of arc batholiths.

Book chapter↗

USGS2025

Critical minerals in mine waste

Introduction Critical minerals are commodities with vulnerable supply chains that play a vital role in supporting the United States’ economy, national defense and security, emerging technologies, and energy independence. The prosperity of our Nation depends on generating a resilient supply of domestic critical minerals; mine waste may be an untapped source of these commodities. Mine waste from centuries of legacy mining persist on the landscape and may contain critical minerals and other valuable commodities previously deemed uneconomic to recover. At modern mines, the financial viability of recovering byproduct critical minerals, which are not the primary target, may be marginal and can ultimately destine them to mine waste. Further, mine waste can be a liability for the mining company or, at legacy mines, the taxpayer because of its effect on the landscape. The U.S. Geological Survey (USGS) has several initiatives to evaluate critical mineral resources in various types of waste. This factsheet highlights studies of mine waste carried out by USGS scientists at the Geology, Energy & Minerals Science Center in collaboration with other science centers funded through the USGS Mineral Resources Program. Recovery of critical minerals from mine waste can aid in remediation efforts and increase domestic supply of vital mineral resources.

Fact Sheet↗