USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Geophysical Research Solid Earth”

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

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

1,663 records · Page 33Linked to original sources

Stable isotope composition and geochemistry of calcite and dolomite in the Mountain Pass carbonatite: A lens into petrogenesis

Carbonatites host most of the global rare earth element (REE) deposits. The petrogenesis of these rocks, including magmatic and post-magmatic processes, are poorly understood but critical in forming and upgrading these deposits. The Mountain Pass carbonatite, which hosts the only active REE mine in North America, is lithologically variable but consistently contains >50% calcite and/or dolomite. Multiple generations of carbonates are evident, suggesting bulk analysis of carbonate minerals would generate a mixed isotopic and geochemical signature. In-situ stable carbon and oxygen isotopes, combined with geochemical analysis and petrography, distinguish carbonate generations, trace the source region of the carbonatite magma, and link carbonate genesis to REE mineralization throughout the intrusion. Primary carbon and oxygen isotope values outside the accepted mantle range may suggest a subduction-related carbon and oxygen source in the mantle. Deformation textures in primary calcite indicate syn- or post-crystallization ductile flow. Carbonates with secondary textures and heavy δ 13 C and δ 18 O values are associated with ore mineralization, suggesting fluid mobilization was important in concentrating REE and upgrading of the Mountain Pass deposit.

California

Measuring stress In high pressure deformation experiments with high speed fiber-optics

High‐pressure, high‐temperature rock deformation experiments are essential for understanding deep Earth processes, but accurately measuring stress is challenging due to the inaccuracy introduced by seal friction within the apparatus and large inertia, which affects earthquake process measurements made far from the experimental fault. To overcome these limitations, we developed and implemented a simple, compact fiber‐optic sensor based on an External Cavity Fabry–Perot Interferometer for in situ load measurement on a piston inside the pressure vessel of a Griggs‐type high‐pressure apparatus. The sensor can be used at sample temperatures up to 800°C for both slow creep and fast rupture testing with bandwidth capability from DC to 6 MHz. Two important experimental results are described: (1) the first direct measurements of seal friction in this apparatus type; these measurements record seal friction approximately one‐tenth of the confining pressure under both low (150 MPa) and high (1 GPa) pressure conditions; and (2) high‐temporal resolution capture of dynamic stress drops during stick‐slip faulting events at high‐confining pressure, revealing near‐fault wave propagation details, high‐frequency oscillations, and implying high‐transient slip rates (6–12 m/s) previously inaccessible with external sensors.

The Seismic Record

The influence of pre-existing structures on geothermal springs: Inferences from potential field mapping in Surprise Valley, CA and other sites In the northwestern Great Basin

Surprise Valley, located in the northwestern Great Basin, is an asymmetric extensional basin that marks a major tectonic transition between the relatively un-extended volcanic Modoc Plateau to the west, and the Basin and Range to the east that has undergone 10-15% extension. In addition, it sits just north of the Walker Lane which accommodates up to 20% of dextral slip associated with Pacific-North American plate interactions. Thermal springs issue from eight areas within Surprise Valley. Most of these occur within the basin and are not situated on the main basin forming range-front faults. As a result, efforts to resolve the structural setting of the valley’s hydrothermal system have relied on geophysics to characterize basin structure and geology. Extensive efforts to map the basin with ground and airborne magnetics have revealed a >35 km-long linear, intra-basin magnetic high, interpreted as a buried dike swarm. Geothermal springs on the eastern side of the valley, including Seifert hot springs, Leonards hot springs, and Surprise Valley hot springs (SVHS), are all situated along the magnetic high and occur at local breaks and bends in the anomaly, suggesting that fracture permeability is enhanced along the feature and particularly at these discontinuities. Recent studies, including drilling over the anomaly near SVHS that likely intersected dike material, as well as subsequent mapping and sampling of dikes outcropping along the anomaly on the playa surface south of SVHS, confirm (as previously inferred) that mafic intrusives are the principal source of the anomaly. Similar interpretations made in two other valleys (in southern Oregon and northwestern Nevada), where inferred intra-basin dikes appear to be spatially correlated with hot springs or prospective geothermal resource areas, suggest that the impact of pre-existing basement structure on hydrothermal activity may pertain more generally to other hydrothermal settings throughout the Great Basin. If so, efforts to map basement may enhance understanding structural controls on some geothermal systems. Furthermore, similarities across these disparate sites suggest that magmatism may play a much larger role in accommodating extension and influencing basin evolution across the western Great Basin than previously recognized.

Conference Paper

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

Hydrogeologic conceptual model of groundwater occurrence and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River is a major concern in the southwestern United States where the river provides water to about 40 million people for municipal and industrial use and is used to irrigate about 5.5 million acres of land. Much of the salinity in the Colorado River Basin is derived from natural interactions of surface water and groundwater with various geologic materials (rocks, soils, and alluvial deposits). The Dolores River in southwest Colorado is a major tributary of the Colorado River that historically accounts for about 6 percent of the salinity load to the Upper Colorado River Basin with the Paradox Valley being the primary source of salinity to the Dolores River. The Paradox Valley, one of several salt-anticline valleys in the region, is a fault-bounded topographic basin aligned with and exposing an underlying salt-anticline core. Salt deposits in the Pennsylvanian Paradox Formation of the Hermosa Group form an elongated salt diapir oriented northwest to southeast that is up to 12,000 feet (ft) thick beneath the present valley floor. Surface erosion, groundwater circulation, and weathering during Tertiary and Quaternary valley formation contributed to development of a cap rock, collapse features, breccia, and brine at the top of the exposed salt diapir. Today (2023), brine occurring in the brecciated cap rock and underlying salt deposits is in hydraulic connection with an overlying freshwater alluvial aquifer, and depending on seasonal river stage and hydrologic conditions, the brine discharges to the Dolores River causing the observed increase in salinity as the river crosses the Paradox Valley. To reduce salinity concentrations in the Dolores River, the Bureau of Reclamation (Reclamation) operates the Paradox Valley Unit (PVU). The PVU project consists of nine shallow brine pumping wells near the Dolores River and one deep disposal well where the brine is injected for disposal. When operational, the PVU pumping wells extract brine from the base of the alluvial aquifer that is piped and injected into a deep disposal well about 3 miles southwest of the PVU. The PVU became fully operational July 1, 1996, and by 2015, operation of the PVU had reduced salinity concentrations in the Dolores River by as much as 70 percent compared to pre-PVU conditions. In response to a 4.5 magnitude earthquake, injection operations, and thus PVU pumping, were ceased from March 2019 to June 2022. A trial period of PVU operation began in June 2022 with a reduced injection rate, and thus PVU pumping rate, of about two-thirds capacity to gather additional information and guide future operational decisions. In cooperation with Reclamation, the U.S. Geological Survey (USGS) developed this report to present the current (2023) understanding of groundwater and brine occurrence and discharge to the Dolores River in the Paradox Valley. Results from the compilation of spatial datasets, groundwater sampling and age dating, and aquifer tests are presented to provide improved understanding of the Paradox Valley hydrogeology, to supply datasets for a numerical groundwater-flow and brine-transport model, and to support future operations of the PVU. The hydrogeologic data provided herein, along with the most recent loading analysis for the Dolores River in the Paradox Valley, and a previous conceptual model for brine discharge to the river are used to present a conceptual understanding of groundwater occurrence in the Paradox Valley.

Colorado

A diatom-based quantitative sea-ice proxy for the Bering and Chukchi seas

Sea ice affects Earth's climate system on both regional and global scales. Its incorporation into climate can be used to achieve more accurate predictions of future climate. However, instrumental records of sea-ice concentration do not extend earlier than 1978. In an effort to extend this record, we constructed a proxy using the generalized additive model based on relative abundances of five easy-to-identify diatom species found in sediment samples across the Bering and Chukchi seas. Here we present the first quantitative diatom-based sea-ice proxy developed for Beringia. The developed proxy has been applied to two sediment cores in the Bering Sea ranging from 0 to 25.7 ka (HLY0204 51JPC) and 369 to 430 ka (IODP Exp 323 Site U1345) and one in the Chukchi Sea ranging from 2.7 to 10 ka (HLY0204 24JPC). The obtained reconstructions of sea-ice concentrations are similar, but not identical to previously published qualitative and nearby records based on other proxies. Because our results are quantitative, they can be incorporated into regional climate models. The proxy is publicly available as an R Shiny application (app) and can be applied to any diatom count from marine sediments in the region.

Bering Sea, Chukchi Sea

Assessment of extreme subsurface hydrologic conditions captured during atmospheric river storms in the San Francisco Bay area (California, USA) with applications to shallow landslide early warning

An increase in soil pore water pressure is the typical trigger for the majority of landslides caused by rainfall. Atmospheric river storms, common to the west coast of North America during the winter season, can deliver landslide triggering rainfall resulting in severe impacts to coastal communities. Using a network of hydrologic monitoring stations situated within landslide-prone terrain in the San Francisco Bay area of California (USA), we assess the meteorologic conditions and resulting hydrologic and landslide response resulting from eight consecutive storm events that caused thousands of shallow landslides during the winter of 2022–2023. We find disparate hydrological responses and resultant degrees of observed landsliding ranging from < 1 landslide/km2 to 18 landslides/km2 at the monitoring sites that reflect the interplay and differences between rainfall delivery, subsurface hydrological characteristics, and geotechnical properties at each site. Antecedent soil moisture from both early season rainfall and the first storm in the sequence played a critical role in setting up some hillslopes for failure. Subsequent storms then generated elevated pore water pressures for several hours with associated landsliding. However, we find that the occurrence of widespread landsliding required not only sufficient pore water pressure magnitude in susceptible hillslopes, but also full and prolonged effective soil saturation throughout hillslope profiles. Landslides may still occur at lower values and durations of effective saturation but are likely to be less extensive regionally. We present these findings within the context of research directions and improvements to landslide early warning systems first suggested by researchers 40 years ago.

California

The feasibility of using lidar-derived digital elevation models for gravity data reduction

Gravity data require submeter elevation accuracy for data processing, and differential global navigation satellite system (dGNSS) equipment is commonly used to acquire three-dimensional positional data to achieve such accuracy. However, lidar (light detection and ranging) data are commonly used to develop digital elevation models (DEMs) of Earth’s surface. Therefore, using elevations from lidar-derived DEMs for gravity-data acquisition and reduction may improve field efficiency and reduce cost. This study examines the feasibility of using DEMs for gravity-data reduction by comparing dGNSS elevation data from 435 gravity stations in Michigan, Wyoming, and Colorado with their respective DEM elevations. The results show that the average difference between DEM and dGNSS elevations is 13 centimeters (cm) and that 93 percent of those differences are less than 50 cm, even in areas with steep terrain. Because an elevation discrepancy of 50 cm corresponds to an error of roughly 0.1 milligals (mGal) in the simple Bouguer gravity anomaly, the results suggest that lidar-derived DEMs are a viable source for acquiring the elevation data needed to process gravity data, thus improving both the cost and efficiency of data collection for regional surveys where an accuracy of less than 1.0 mGal is desired.

Open-File Report

The 3D National Topography Model Call for Action—Part 2: The Next Generation 3D Elevation Program

The three-dimensional (3D) National Topography Model initiative to integrate elevation and hydrography data includes the next generation of hydrography data from the 3D Hydrography Program and the next generation of elevation data from the 3D Elevation Program (3DEP). The first-ever collection of light detection and ranging (lidar) data for the nation (IfSAR for Alaska) provides a critical baseline reference, and the addition of multiple repeat elevation mapping projects as part of the next generation of 3DEP would substantially expand analysis capabilities. As the U.S. Geological Survey (USGS) is closing in on our goal of complete coverage with 3DEP data available or in progress for 98.3 percent of the Nation at the end of fiscal year 2024, the USGS is already transitioning to the next generation of 3DEP. Based on the 3D Nation Study results and input from a broad range of stakeholders, the USGS National Geospatial Program has finalized a new design for 3DEP that provides increased lidar quality levels and refresh rates. The new program is designed with more flexibility to meet changing user needs and take advantage of improvements in mapping technologies. The program will aim to expand the level of interagency coordination for topobathymetric lidar acquisition for inland rivers. The next generation of 3DEP will also aim to emphasize research, including advancing program design, products, and services and engaging and leveraging the evolving 3D industry. Research goals also include becoming more flexible in meeting user needs and taking advantage of evolving remote-sensing technologies. The program also plans to move from focusing on producing standard products to producing a concept of a 3D Nation Ecosystem with a variety of inputs, products, and services.

Circular

Four-band image mosaic of the Colorado River Corridor downstream of Glen Canyon Dam in Arizona, derived from the May 2021 airborne image acquisition

In May 2021, the U.S. Geological Survey’s Grand Canyon Monitoring and Research Center acquired airborne multispectral high-resolution data for the Colorado River in the Grand Canyon, Arizona. The image data, which consist of four spectral bands (red, band 1; green, band 2; blue, band 3; and near infrared, band 4) with a ground resolution of 20 centimeters, are available as 16-bit unsigned-integer GeoTIFF files in Sankey and others (2024) (available online at https://doi.org/10.5066/P9BBGN6G ). The image files are projected in the State Plane Coordinate System, using the central Arizona zone (202) with the North American Datum of 1983 National Adjustment of 2011. The assessed spatial accuracy for these data is based on 47 ground-control points that were independent from the ground-control points used by the contractor for aerotriangulation and is reported at the 95-percent confidence level as 0.514 meter (m) and a root mean square error of 0.297 m. The intended uses of this dataset are primarily in support of scientific research and monitoring applications. Examples of these applications include high-resolution spatial and temporal change detection of the river channel, geomorphic landforms, riparian vegetation, and backwater and nearshore habitat, as well as other ecosystem-wide mapping. These imagery data also serve as reference material for field science mission planning, as base data for field data collection including community science activities, and as a highly detailed guide for technical boat operation during science activities such as reconnaissance for nighttime missions and navigating rapids during low flows.

Arizona, Nevada, Utah

Woods Hole Coastal and Marine Science Center—2023 annual report

The 2023 annual report of the U.S. Geological Survey Woods Hole Coastal and Marine Science Center highlights accomplishments of 2023, includes a list of 2023 publications, and summarizes the work of the center, as well as the work of each of its science groups. This product allows readers to gain a general understanding of the focus areas of the center’s scientific research and learn more about specific projects and progress made throughout 2023, all while enjoying photographs taken in various environments and laboratories, and applicable maps and figures.

Circular

Woods Hole Coastal and Marine Science Center—2024 annual report

The 2024 annual report of the U.S. Geological Survey Woods Hole Coastal and Marine Science Center highlights accomplishments of 2024, includes a list of 2024 publications, and summarizes the work of the center, as well as the work of each of its science groups. This product allows readers to gain a general understanding of the focus areas of the center’s scientific research and learn more about specific projects and progress made throughout 2024, all while enjoying photographs taken in various environments and laboratories, and applicable maps and figures.

Circular

The 3D Elevation Program—Supporting New Hampshire’s economy

Introduction The topography of New Hampshire ranges from the Coastal Lowlands to the Eastern New England Upland to the White Mountains region. High-quality statewide elevation data are useful in managing this very diverse landscape. For example, the short coastline, including the Great Bay estuary and the Hampton-Seabrook marshes, is of disproportionately high value to New Hampshire’s tourist economy. The vulnerability of the coast to the effects of sea-level rise underscores the need for accurate, high-quality nearshore topographic elevation data and offshore bathymetric data to effectively manage the coast’s valuable resources, which include important fisheries, habitat, and infrastructure. Another important use for accurate elevation data in New Hampshire is in the evaluation of flood hazards and their potential environmental and infrastructure effects. This evaluation includes mapping of inundation and sediment transport, and assessing the associated costs of flooding. Addressing this challenge requires detailed knowledge of both surface topography and inland bathymetry. Other important activities having a substantial economic element and needing accurate elevation data include geologic resource assessment and hazard mitigation, urban and regional planning, infrastructure and construction management, and cultural resources preservation and management. Critical applications that meet the State’s management needs depend on light detection and ranging (lidar) data that provide a highly detailed three-dimensional model of the Earth’s surface and aboveground features.

New Hampshire

New developments at the Center for Engineering Strong-Motion Data (CESMD)

The Center for Engineering Strong-Motion Data (CESMD), an internationally utilized joint center of the U.S. Geological Survey (USGS) and the California Geological Survey (CGS), provides a single access point for earthquake strong-motion records and station metadata from the CGS California Strong-Motion Instrumentation Program (CSMIP), the USGS National Strong-Motion Project (NSMP), the USGS Advanced National Seismic System, and other affiliates. The CESMD has been continuously improving its webtools to facilitate the access of strong-motion data and metadata for use in post-earthquake response and for scientific and engineering research applications. The Center provides raw and processed strong-motion data via the Engineering Data Center (EDC) and the Virtual Data Center (VDC) web portals. This paper focuses on the strong-motion products provided by the EDC where more than 48,000 records with peak ground accelerations greater than 0.1% g from over 2400 earthquakes are currently hosted. and on the ongoing efforts to develop data access tools and applications. The new developments and ongoing efforts in the EDC include: 1) enhancements to the CESMD webservices to facilitate access to station metadata, earthquake information, and strong motion records 2) new features to the interactive map interface, improving the visualization and access to earthquake, station, and record information, 3) efforts to develop a new web application tool for data format conversion from a number of data formats, 4) efforts to unify varying waveform data formats into a consistent format, 5) ongoing efforts to compile seismic station site geology, measured or inferred Vs30 values, shear-wave profiles, NEHRP site class, and available structural instrument deployment schematics, and 6) a special studies pages for research topic-specific ground motion datasets that offer uniform processing of records from a variety of sources.

Conference Paper

Geochemical processes related to mined, milled, or natural metal deposits in a rapidly changing global environment

The demand for metals and raw materials, such as nickel and copper, has been projected to expand in the coming decades, driven by the global energy transition, the need for green technologies, and expanding infrastructure. Consequently, the increasing extraction and production of mining waste can have adverse impacts on surrounding environments and human health. The aim of this thematic collection is to fill critical knowledge gaps in the present-day cycles of metal(loid)s from source to larger sinks, and the effect of environmental management, anthropogenic development, and climate change. Altogether, the studies have been conducted in different natural settings around the world and comprise investigations in laterites, a soil-medicinal plant system, watersheds, and banded iron formations, among others. The geochemical applications in tracing mineralization, its secondary products, and/or potential impact on the immediate environment are highly diverse with applied tools ranging from isotope tracers to major and trace element systematics. Particularly the use of rare earth elements, their patterns and anomalies are methods employed by several studies in this collection. We summarize the findings to offer a potential future direction for the use of geochemical tracing techniques in resource exploration in the context of climate change and environmental challenges.

Geochemistry: Exploration, Environment, Analysis

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

2025 USGS National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands: Overview of model and hazard results

The U.S. Geological Survey recently updated the National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The first version of the PRVI NSHM was released in 2003, and therefore this 2025 update includes over 20 years of new geologic, geophysical, and engineering data, methods, and models. Updates follow similar efforts performed in the recent 2023 50-state NSHM. However, this is the first NSHM in which we: (1) apply an inversion methodology to subduction interface fault sources in the earthquake rupture forecast (ERF) model; (2) develop scaled backbone median ground-motion models and independent aleatory variability models that are applied in the ground-motion characterization (GMC) model; and (3) calculate epistemic uncertainty related to alternative scenarios in the ERF and GMC models for all grid points in the study region. Long-term time-independent mean hazard calculations were performed for peak ground acceleration and 5%-damped pseudospectral acceleration at 21 spectral periods from 0.01- to 10.0-s, for eight National Earthquake Hazards Reduction Program site conditions ranging from V S30 = 150 to 1500 m/s, and for 2%, 5%, and 10% in 50-year probabilities of exceedance (return periods of 2475, 975, and 475 years, respectively). Epistemic uncertainty, in the form of selected percentiles, is also provided for a suite of test sites and all grid points in the study region for limited periods, site conditions, and probabilities of exceedance. Selected results, including comparisons with the 2003 PRVI NSHM, are shown and discussed for selected periods, site conditions, and probabilities of exceedance. When comparing the 2025 PRVI NSHM with the 2003 PRVI NSHM, hazard is generally higher at shorter periods and lower at longer periods, as a result of updates in both ERF and GMC models. The 2025 PRVI NSHM is applicable for return periods greater than ∼475 or less than ∼10,000 years.

Puerto Rico, U.S. Virgin Islands

Oblique contraction along the fastest ocean-continent transform plate boundary focuses rock uplift west of the Fairweather fault, southeast Alaska

Contraction along the Yakutat–North America plate boundary drives 4.6–9.0 mm/year Holocene rock uplift rates along Earth's fastest slipping (≥49 mm/year) ocean–continent transform fault, the Fairweather Fault. Between Icy Point and Lituya Bay, the near-vertical Fairweather fault focuses rock uplift and rapid right-lateral slip by accommodating both vertical and fault-parallel strain during oblique-slip and separate, predominantly strike-slip ruptures. Unusually high uplift rates, indicated by radiocarbon and luminescence dating, result from a 10-km-wide, asymmetric, positive flower structure along a 20°, ∼30-km-long restraining double bend in the Fairweather fault. The principal reverse fault in the flower structure, the offshore, blind Icy Point–Lituya Bay fault, ruptures no more than every 460–1040 years evidenced by uplifted Holocene shorelines. Maximum 3–5 m coseismic uplifts imply 3.1–10 m dip slip per event and earthquake magnitudes of M w 7.0–7.5. The Yakutat block collides obliquely into North America, and our model entails oblique slip on the Fairweather fault with and without corupture on the reverse fault. Oblique slip is evident by vertically offset (&gt;25 m) fluvial and marine terraces and by the primary Fairweather fault strand that strikes &gt;20° to the west of plate-boundary motion.

Alaska