Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Enviornmental Monitoring and Assessment”

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

Subduction zone earthquake catalog separation tool: Implementation in the USGS 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

The U.S. Geological Survey (USGS) periodically releases updates to National Seismic Hazard Model (NSHM) for the United States and its territories leveraging current scientific knowledge and methodologies to guide public policy, building codes, and risk assessments regarding potential ground shaking due to earthquakes that may result in infrastructure damage. In subduction zones, there is a need to separate the earthquake catalog into tectonic regimes to create specific seismicity models for which the most appropriate ground‐motion models are then applied. Here, we describe newly developed methods and software, called CatSep, that classifies subduction zone events into three primary tectonic regimes: crustal, interface, and intraslab. This method incorporates information about the location of the earthquake relative to the subducting slab, the depth of the Mohorovičić discontinuity, and the earthquake’s moment tensor. Applying this method is a first step in the NSHM workflow for regions covering U.S. subduction zones. Results using this subduction zone earthquake catalog separation tool for the 2025 Puerto Rico and U.S. Virgin Islands NSHM earthquake catalog are presented and analyzed.

Puerto Rico, U.S. Virgin Islands↗
USGS2026

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

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

JGR Machine Learning and Computation↗
USGS2026

Computational electromagnetic geophysics for groundwater system studies: A review on established practices and recent advances

Identifying effective solutions for locating groundwater resources and ensuring the quality of drinking water is increasingly urgent, given the challenges posed by climate change and population growth. This review investigates electromagnetic geophysical imaging techniques, in both time- and frequency-domain, that can provide valuable insights for groundwater assessment. We explore computational electromagnetic methods used to evaluate electromagnetic data and several recent hydrogeophysical case studies. As open-source frameworks for modeling electromagnetic geophysical problems become available, a broader range of researchers can interpret their data with computationally advanced software. We provide an overview of documented open-source codes for evaluating electromagnetic data and analyze various hydrological targets in relation to their electromagnetic surveying technique and the computational method applied. Furthermore, we evaluate the potential of advanced computational techniques, including three-dimensional modeling, non-deterministic inversion and machine learning, to couple geophysical with numerical groundwater modeling and apply it in groundwater system studies. Despite obstacles such as complexity and resource demands, our findings indicate that the quantification and integration of predictive uncertainties from both electromagnetic and hydrological data and simulations would significantly improve the reliability of hydrogeophysical models. This can lead to a deeper understanding of groundwater systems and improved management practices.

Journal of Hydrology↗
USGS2023

Paleomagnetism of the Harrat Rahat volcanic field, Kingdom of Saudi Arabia—Geologic unit correlations and geomagnetic cryptochron identifications

Paleomagnetic rock samples were collected from 173 drill sites in the Quaternary alkali basaltic volcanic field of northern Harrat Rahat, Kingdom of Saudi Arabia. Laboratory measurements on these samples established that lava flows and vent complexes—identified and mapped from field characteristics, rock types, and compositions as products of single or temporally close eruptions—typically record single, or very similar, directions of remanent magnetization. Correlations defined through geologic mapping, spatial association, geochemistry, geochronology, and identical mean remanent directions indicate at least 16 brief episodes of temporally clustered eruptions. These episodes had durations of a few centuries or less. Anomalous remanent magnetic directions were found for at least 13 mapped lavas of northern Harrat Rahat, which demonstrate that they were acquired during brief geomagnetic cryptochrons during the Brunhes Normal Polarity Chron. These uncharacteristic directions enhance the opportunity to identify common eruptive episodes, and to better understand and evaluate assessments of eruption ages based on 40 Ar/ 39 Ar geochronology. Combining paleomagnetic and regional archaeomagnetic results for the youngest eruptions allows us to evaluate their historical age assignments and, in one case, refute a previously assigned provisional age.

Harrat Rahat volcanic field↗
USGS2026

Production of mineral commodities and geospatial map of the mineral industries and related infrastructure of China

As part of the U.S. Geological Survey’s (USGS) mission to distribute global mineral information and analyze supply chains, this study provides a comprehensive review of the global significance of China’s mineral production and capacity in 2023. Of 77 mineral commodities in the USGS dataset, China produced 74 and was the world’s first-ranked producer for 39 of the 74. Compared to the high share of global mineral production, including up to 98 percent of global gallium production, the country’s share of global mineral reserves was relatively small, ranging from 20 percent (zinc ore) to 52 percent (tungsten ore). China’s imports of metal ores, slag, and ash accounted for 64 percent of global imports of such commodities by value. The country’s exports of base metals and articles of base metal accounted for 17 percent of the global exports. To help nongeographic information system users assess the spatial distribution of mineral mines, processing facilities, and ports for trades in China, this study created a geospatial (also called “georeferenced”) portable document format (GeoPDF) map. In addition, the GeoPDF contains mineral resource tracts (such as antimony, copper, potash, coal, and oil and gas), exploration sites, and energy infrastructure based on the preexisting USGS data.

Open-File Report↗
USGS2026

Coastal barrier resilience and resistance: Analysis and metrics for characterizing coastal state

Barrier islands are shaped by a variety of short- and long-term environmental processes such as storms and relative sea-level rise. These islands, found along the estuarine-marine interface, provide ecosystem services including storm surge and wave attenuation, erosion protection to inland marshes, habitat for fish and wildlife, and recreation. Natural resource managers require actionable information on how barrier island resilience and resistance changes over time to understand how an island’s current state relates to past conditions and to inform restoration prioritization and implementation. The U.S. Geological Survey and The Water Institute collaborated on a study to develop indicators of resilience and resistance for barrier islands in Louisiana. Here, resilience captures island persistence on yearly to decadal time scales, and resistance captures persistence on event time scales of days to weeks. The indicators fall in two categories: Tier 1 Screening Metrics, that can be readily calculated from available data, are easily interpretable as an evaluation of barrier condition, and provide a high-level snapshot of overall barrier resilience and resistance; and Tier 2 Analysis Metrics, which are detailed metrics that required specialized analysis and interpretation and are more applicable to answering specific questions managers may have about barrier state. The research team derived Tier 1 resilience indicators from subaerial land and vegetation cover calculated from publicly available maps and products based on satellite imagery. By benchmarking the total land and vegetation extent against their respective historical maxima, this metric provides a snapshot of an island’s current state in the context of its long-term trajectory. The research team developed Tier 1 resistance indicators based on subaerial island configuration and water level recurrence as a proxy for evaluating island resistance to storms, which are the primary driver of short-term change. These Tier 1 metrics can be analyzed over time to provide a high-level assessment of how an island’s resistance decreases because of elevation loss or sea-level rise or increases due to restoration or natural recovery. The research team developed Tier 2 resilience and resistance indicators and associated analyses to provide detailed information for specific time periods or applications (e.g., wildlife management). These metrics include habitat coverage from high-resolution maps, which show composition changes over time to capture the evolving resilience of specific habitat types; high tide flooding analysis, which evaluate island area relative to specified flooding thresholds to characterize resistance in the short-term or, if analyzed over time, indicate changes in resilience; and hypsometric curve analysis, which allows managers to evaluate island area changes above their own elevation benchmarks of interest and similarly characterize resistance in the short-term or indicate changes in resilience if assessed over time. The research team calculated Tier 1 metrics of the barrier islands and headlands along the coast of Louisiana for the period of 1984 through 2021 and Tier 2 metrics for select times during that period depending on data available and quality. The results were captured in a report card for each barrier, which also includes an overview of the metrics and their interpretation; a restoration and storm history; and Tier 1 and Tier 2 metric analysis, including benchmarking against coastwide and regional values as well as to an island’s pre-restoration trajectory. These report cards provide a readily digestible synthesis of barrier condition and trajectory that coastal managers can use to support restoration prioritization and other decisions.

Louisiana↗
USGS2023

3D Dynamic rupture modeling of the 6 February 2023, Kahramanmaraş, Turkey Mw 7.8 and 7.7 earthquake doublet using early observations

The 2023 Turkey earthquake sequence involved unexpected ruptures across numerous fault segments. We present 3D dynamic rupture simulations to illuminate the complex dynamics of the earthquake doublet. Our models are constrained by observations available within days of the sequence and deliver timely, mechanically consistent explanations of the unforeseen rupture paths, diverse rupture speeds, multiple slip episodes, heterogeneous fault offsets, locally strong shaking, and fault system interactions. Our simulations link both earthquakes, matching geodetic and seismic observations and reconciling regional seismotectonics, rupture dynamics, and ground motions of a fault system represented by 10 curved dipping segments and embedded in a heterogeneous stress field. The M w 7.8 earthquake features delayed backward branching from a steeply branching splay fault, not requiring supershear speeds. The asymmetrical dynamics of the distinct, bilateral M w 7.7 earthquake are explained by heterogeneous fault strength, prestress orientation, fracture energy, and static stress changes from the previous earthquake. Our models explain the northward deviation of its eastern rupture and the minimal slip observed on the Sürgü fault. 3D dynamic rupture scenarios can elucidate unexpected observations shortly after major earthquakes, providing timely insights for data‐driven analysis and hazard assessment toward a comprehensive, physically consistent understanding of the mechanics of multifault systems.

The Seismic Record↗
USGS2025

ARCHI: A new R package for automated imputation of regionally correlated hydrologic records

Missing data in hydrological records can limit resource assessment, process understanding, and predictive modeling. Here, we present ARCHI (Automated Regional Correlation Analysis for Hydrologic Record Imputation), a new, open-source software package in R designed to aggregate, impute, cluster, and visualize regionally correlated hydrologic records. ARCHI imputes missing data in “target” records by linear regression using more complete “reference” records as predictors. Automated imputation is implemented using a novel, iterative algorithm that allows each site to be considered a target or reference for regression, growing the pool of complete references with each imputed record until viable gap-filling ceases. Users can limit artifacts from spurious correlations by specifying model-acceptance criteria and applying geospatial, correlation, and group-based filters to control reference selection. ARCHI provides additional functions for visualizing results, clustering records with similar correlation structures, evaluating holdout data, and interactive parameterization with an accessible and intuitive graphical user interface (GUI). This methods brief provides an overview of the ARCHI package, modeling guidelines, and benchmarking on two regional groundwater-level datasets from the Central Valley, CA and Long Island, NY. We evaluate ARCHI alongside widely used multivariate imputation software to highlight and contextualize its computational efficiency, imputation accuracy, and model transparency when applied to large, groundwater-level datasets.

California, New York↗
USGS2025

Three-dimensional seismic velocity model for the Cascadia Subduction Zone with shallow soils and topography, version 1.7

The U.S. Geological Survey’s seismic velocity model for the Cascadia Subduction Zone provides P- and S-wave velocity ( V P and V S , respectively) information from 40.2° to 50.0° N. latitude and −129.0° to −121.0° W. longitude, and is used to support a variety of research topics, including three-dimensional (3D) earthquake simulations and seismic hazard assessment in the Pacific Northwest. This report describes an update to the previous version (v) 1.6 of the 3D seismic velocity model for the Cascadia Subduction Zone. This new model (herein referred to as v1.7) contains more detailed near-surface structure for improved earthquake ground motion modeling. Updated features include the addition of a new shallow soil velocity model in the top few hundred meters and the option of adding user-specified topography. Although v1.6 of the Cascadia seismic velocity model has a minimum V S of 600 meters per second (m/s), the new model (v1.7) has a minimum V S of approximately 40 m/s. Overall, this update will allow for more accurate ground motion estimates from 3D simulations of scenario earthquakes in the Cascadia Subduction Zone region.

British Columbia, California, Oregon, Washington↗
USGS2024

Flood-inundation maps for the Cuyahoga River in and near Independence, Ohio, 2024

Digital flood-inundation maps for a 9.9-mile reach of the Cuyahoga River in and near Independence, Ohio, were created by the U.S. Geological Survey (USGS) in cooperation with the Northeast Ohio Regional Sewer District Board of Trustees. Water-surface profiles were computed for the stream reach by using a one-dimensional steady-state step-backwater model. The model was calibrated to the current (2024) stage-streamflow relation (rating curve 43.0) for the USGS streamgage 04208000, Cuyahoga River at Independence, Ohio. The resulting hydraulic model was then used to compute 13 water-surface profiles for water levels (flood stages) ranging from 14.00 to 26.00 feet. The flood stages range from “action stage” to above “major flood stage” as reported by the National Weather Service. The simulated water-surface profiles were then used in combination with a digital elevation model derived from light detection and ranging data to map the inundated areas associated with each flood profile. The flood-inundation maps and the supporting hydraulic model produced by this study can be used by emergency managers and local officials to assess flood mitigation strategies and to define flood hazard areas to protect life and property, to coordinate flood response activities such as evacuations and road closures, and to aid postflood recovery efforts.

Ohio↗
USGS2025

Global maps of critical mineral production in 2023

Introduction The global production of many mineral commodities, especially critical minerals, is concentrated in a few countries that have mineral resources and the infrastructure necessary to mine and process those resources. For this reason, the type and amount of mineral production differ by country. For example, many countries produce such metallic ores as gold and silver, whereas only a few countries produce magnesium, niobium, platinum-group metals, and rare earths. The concentration of mining and processing in certain countries necessitates the existence of a global supply chain. A mineral supply chain is the sequence of mining and processing of minerals and manufacturing of products. Mineral supply chains are global in scale, complex, and dynamic. Supply chain data can be used to understand how a country’s mineral resources and various economic, technical, and environmental factors affect the complexity of global supply chains. This fact sheet summarizes the world’s leading countries (those accounting for 5 percent or more of a commodity’s global production in 2023) for production of select mineral commodities (mainly critical minerals) in the mining and processing stages. These countries and the minerals they produce are synthesized on global maps to communicate the status of, and potential risk to, mineral commodity supply chains from geographic production concentration. Trade data from United Nations Statistics Division (2025) is used to support assessments of the observed production data.

Fact Sheet↗
USGS2026

The damability function: A probabilistic approach to regional landslide dam susceptibility analysis applied to the Oregon Coast Range, USA

Landslides can dam rivers and require rapid response to mitigate catastrophic outburst floods. Here, we present a workflow to map landslide dam formation susceptibility at a regional scale. We define a probabilistic function that combines river valley width and landslide volume to efficiently determine the likelihood of a landslide dam or “damability”. We combine damability values with landslide susceptibility to estimate landslide dam susceptibility. The valley width measurements are automated using a new elevation threshold-based algorithm. Landslide volume is represented as a statistical distribution from mapped landslides. We validate and apply our approach to the Oregon Coast Range, USA and find that 36 % of river stretches exceed a dam potential threshold; these are in river headwaters and steeper terrain, which in this case correlate with more resistant lithologies. We also estimate volumes of the potential dammed lakes and find that most rivers with high dam susceptibility are less likely to impound large lakes because they have low drainage areas. However, widespread susceptibility, and the potential impacts from exceptionally large landslides, suggest that this hazard should be considered in the Pacific Northwest. The damability function workflow can ingest new data and be applied more broadly to assess future landslide dam hazards.

Oregon↗
USGS2025

3D viscoelastic models of slip-deficit rate along the Cascadia subduction zone

Interseismic deformation in the Pacific Northwest is constrained by the horizontal crustal velocity field derived from the Global Positioning System (GPS) in addition to vertical rates derived from GPS, leveling, and tide gauge measurements. Such measurements were folded in to deformation models of fault slip rates as part of the 2023 National Seismic Hazard Model (NSHM) update. Here I build upon one of the contributing models, the viscoelastic earthquake-cycle model of Pollitz [2022]. This model permits inclusion of effects of time-dependent viscoelastic relaxation within earthquake cycles (i.e., ‘ghost transients’) and laterally variable elastic and/or ductile material properties. I lever-age these capabilities to incorporate the Cascadia megathrust into Western U.S.-wide deformation models in which crustal fault slip rates are estimated simultaneously with slip deficit rates along the interplate boundary between the descending Juan de Fuca plate and North American plate. This effort includes construction of a margin-wide model of viscoelastic structure founded on the Slab 2.0 model and probes different models of the ductile properties of the surrounding oceanic asthenosphere, continental lower crust, and mantle asthenosphere. This results in new estimates of the distribution of slip deficit rate along the ∼ 1000 km long margin, highlights the importance of correcting for glacial-isostatic adjustment effects, and permits assessment of sensitivity of results to assumed ductile properties.

California, Oregon, Washington↗
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↗
USGS2023

Stream sediment geochemistry in mineral exploration: A review of fine-fraction, clay-fraction, bulk leach gold, heavy mineral concentrate and indicator mineral chemistry

Stream sediment surveys support early-stage reconnaissance mineral exploration and regional assessment programmes, enhanced by recent improvements in analytical method detection limits, continuously improving mineral chemistry, and new approaches to the interpretation of geochemical data. Sediment surveys may be used to predict catchment basin lithology, mineralization type based on pathfinder geochemistry, and geological features based on indicator mineral chemistry. Sediment surveys that target a finer-fraction sediment sample led to the discovery of the La Colosa gold deposit, Colombia. The Batu Hijau porphyry Cu–Au deposit in Indonesia was discovered based on an anomalous clay-sized fraction sample 12 km downstream. In an arid region with poorly developed drainages and minor topographic relief, the Ag-base-metal Navidad District in Argentina was discovered with clay-fraction sediment geochemistry. Heavy mineral concentrate (HMC) sediment surveys that include mineral chemistry determinations have led to global diamond discoveries. HMC surveys contributed to discovery of the Ring of Fire Ni–Cu–PGE and chromite district, Ontario, Canada. Discoveries and geochemical mapping can assist advancement of the application of stream sediment geochemistry in those global areas for which lithologies and deposits are exposed. Stream sediment surveys continue to be one of the most cost-effective geochemical methods for covering large areas for mineral exploration.

Geochemistry: Exploration, Environment, Analysis↗
USGS2026

USGS Geochron—A database of geochronological and thermochronological dates and data—Technical documentation

Geochronological and thermochronological data are essential for constraining the timing and rates of geological processes, supporting geologic mapping, natural hazard assessment, and resource exploration. The U.S. Geological Survey (USGS) Geochron Database is a centralized, relational database that integrates USGS and State geological survey data in accordance with the National Geologic Mapping Act of 1992 and its reauthorizations. This report documents the structure and development of the database, including its standardized schema, controlled vocabularies, and protocols for compiling legacy and newly published data. The database is designed to adhere to FAIR (Findable, Accessible, Interoperable, and Reusable) data principles and currently (2026) includes data from a wide range of geochronological and thermochronological methods. Public access is provided through the USGS Geochron Database Explorer, a browser-based geographic information system (GIS) interface, and through versioned data releases available on ScienceBase, which provide sample-level summary data and detailed analytical information. The USGS Geochron Database is structured to be extensible, supporting the inclusion of additional methods and data types as they are compiled. This format ensures long-term utility and aligns with the Mapping Act’s directive to create a national archive of geochronological information that adds interpretive value to geologic map data. Its design reflects a commitment to data transparency, scientific reproducibility, and national geoscience priorities.

Data Report↗
USGS2024

Revision of ModelMuse to support the use of PEST software with MODFLOW and SUTRA models

Executive Summary ModelMuse is a graphical user interface for several groundwater modeling programs. ModelMuse was updated to generate the input files for the parameter estimation software suite PEST. The software is used with MODFLOW or SUTRA models to run PEST-based parameter estimation and display the updated model inputs after parameter estimation. The PEST input files can also be used with the PEST++ version 5 software suite. Parameter estimation typically requires defining the parameters being adjusted during calibration and observations for assessing calibration quality. After a parameter is defined in ModelMuse, it can be applied to all or part of a model dataset. Pilot points—a parameterization device that facilitates higher levels of parameterization—can be used to assign spatially variable distributions of model inputs. Parameters can be applied to temporally varying features, such as boundary conditions, by either applying them to all the values in a series in one step or by applying separate parameters to individual members of a series. ModelMuse allows the definition of many observation types from various model output files. For MODFLOW 6 and SUTRA models, new options were added to ModelMuse to allow it to display the changed input after parameter estimation is complete. For MODFLOW–2005 and MODFLOW–NWT models, ModelMuse can import an entire model for visualization. An example illustrates the use of PEST with a MODFLOW 6 model in ModelMuse.

Techniques and Methods↗
USGS2025

A 700-year rupture sequence of great eastern Aleutian earthquakes from tsunami modeling of stratigraphic records

Great Aleutian underthrusting earthquakes produced destructive tsunamis impacting Hawaiʻi in 1946 and 1957. Prior modeling of the 1957 tsunami deposit and runup records on eastern Aleutian and Hawaiian Islands jointly with tide-gauge observations across the Pacific Ocean constrained a rupture model with shallow slip up to 26 m along 600 km of the plate boundary. Here we implement this modeling approach to older deposits and show alternating deep and shallow megathrust slip up to 26, 32, and 22 m for great earthquakes along the same segment in the 18 th , 15 th , and 14 th centuries. All three modeled prehistoric Aleutian earthquakes produce tsunami inundation in Hawaiʻi with the most severe, 14 th century event having impacts exceeding the 1957 event. The along-dip variability of these four ruptures spanning seven centuries provides insights on earthquake cycles for engineering design and hazard assessment. The 15 th century and 1957 rupture models provide evidence for recurrence of tsunami earthquakes, which can produce disproportionately large tsunamis for a given moment magnitude due to reduced rigidity in the shallow megathrust. The 14 th and 18 th century events likely ruptured deeper regions that did not slip in 1957, suggesting potential for corresponding deeper failure in the next great eastern Aleutian earthquake.

Alaska, Hawaii↗