Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Data Report”

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,150 records · Page 53Linked to original sources
USGS2024

Methods to evaluate and improve the modeling of rupture directivity in assessment of seismic hazard

In recent years, there have been several advancements related to the modelling of near-source effects of earthquake rupture on strong ground shaking, leading to an improved characterization of ground motions and resulting seismic hazard. Some of these modifications have stemmed from physics-based numerical modelling of the earthquake rupture process, using physics-based dynamic rupture simulations. These contributions have led to a better understanding of how fault rupture characteristics, geometry, and the style of faulting can interact with the hypocenter-dependence on the path from source to site that may ultimately guide the development of seismic directivity models. Moving forward, the application of modern techniques can be used to incorporate these source characteristics and near-fault ground motion behavior that contribute to the azimuthally varying effects that result in rupture directivity. One example is the application of machine learning methods to support more automated integration of new predictor variables in model development and open more evaluation opportunities to access residuals. Here, we utilize several techniques to take advantage of the plethora of synthetic data and its ability to supplement preexisting trends observed in data. We showcase two examples of how models can be either developed, expanded upon, or constrained using artificial neural network model (ANNs). We evaluate the performance of the ANN with existing methods, comparing misfit, potential limitations, and ability to continue to improve upon these methods in the future. One approach uses a set of simulations with corresponding synthetic ground motions from the Southern California Earthquake Center (SCEC) CyberShake study to develop a ground motion model adapted to incorporate seismic directivity information using an ANN. This large database (TBs) enables us to train the model to capture magnitude, period, and distance variations and how these parameters relate to amplification from hypocenters located along finite-faults. In some cases, there is reduced misfit from better representing source features that aren’t included in base ground motion models that neglect hypocenter location (e.g. azimuthal variation, source-to-site terms). Another ANN method uses a shallow-layered neural network model to better fit a hypocenter-independent model. This method adjusts the median and aleatory variability to account for the averaged impact of various hypocenter distributions to fit the underlying directivity adjustment model. This method serves as a template to apply to other directivity models, improving computational efficiency and more readily enabling integration in hazard codes.

California↗
USGS2024

A history of cryohydrogeology modeling and recent advancements through the integration of solute transport

Groundwater flow systems and permafrost are interrelated because permafrost thaw enhances permeability, while groundwater flow can advect heat and accelerate permafrost thaw (McKenzie et al. 2021). Given amplified climate change in cold regions, there is renewed interest in ‘cryohydrogeology’, the study of groundwater in cold regions. Many data-driven studies have shown that permafrost thaw is leading to activated aquifers and increased baseflow across the pan-Arctic region (e.g. Walvoord and Striegl 2007, Evans et al. 2020). Empirical evidence of a subsurface ‘replumbing’ (Walvoord and Kurylyk 2016) in permafrost regions raises questions about the fate of sequestered contaminants in the North (Langer et al. 2023). We will discuss the history of and emerging opportunities in cryohydrogeological modeling, with a focus on recent contaminant transport modeling.

Conference Paper↗
USGS2025

Reaction kinetics and accelerant effects of sulfides in early mature hydrocarbon generation using hydrous pyrolysis

Hydrocarbon generation in organic-rich sediments is influenced by the molecular organic composition and relative abundance of associated minerals. Certain mineral-derived elements act as catalysts and reaction intermediaries, facilitating early-stage hydrocarbon formation in potential source rocks. This study investigated the role of sulfur contributed from pyrite as an accelerant in thermal reaction, focusing on its effects on early maturation and consequent hydrocarbon generation from gilsonite (low-sulfur solid petroleum). Hydrous pyrolysis (HP) experiments were conducted on mixtures of gilsonite and pyrite in varying ratios (1:0.1, 1:0.5, 1:1, 1:2, and 1:10 w/w gilsonite:pyrite) at 320, 350, and 370 °C for 72 h. Untreated and thermally altered residues were analyzed using solid bitumen reflectance (BR o , %), total organic carbon (TOC) content, programmed temperature pyrolysis, scanning electron microscopy with energy-dispersive spectroscopy (SEM-EDS), and X-ray diffraction (XRD) to evaluate the potential accelerant role of pyritic sulfur in hydrocarbon formation. The results show HP residues at 320 and 350 °C with greater pyrite concentrations had higher BR o , while reflectance values were similar in the 370 °C residues, regardless of pyrite concentration, suggesting enhanced reaction at lower thermal conditions. Increasing pyrite content systematically decreased hydrogen index (HI) values while increasing the transformation ratio (TR) and production index (PI), indicating enhanced conversion of organic matter to hydrocarbons with increasing pyrite concentrations. Gas yields increased with pyrite addition, particularly at 350 °C, confirming secondary cracking effects. However, gas production stabilized or declined at higher pyrite loadings (1:10), suggesting alternative reaction pathways such as coke formation. Our data indicate the presence of pyrite lowers the activation energy for thermal cracking, shifting peak experimental hydrocarbon generation temperatures downward by 20–30 °C, with the most pronounced accelerant effects observed at moderate pyrite concentrations (1:0.5 and 1:1). The thermodynamic framework reveals that pyrite stability is influenced by experimental conditions, with pyrrhotite formation favored in the presence of gilsonite due to reduced oxygen fugacity. Pyrite transformation to pyrrhotite, as observed through XRD, SEM-EDS, and predicted by thermodynamic data, further supports the accelerant role of S, as pyrrhotite exhibits a higher hydrogen transfer potential, promoting early oil generation. These findings highlight the importance of pyrite in modulating hydrocarbon generation pathways in organic-rich systems.

Journal of Analytical and Applied Pyrolysis↗
USGS2024

Land-use interactions, Oil-Field infrastructure, and natural processes control hydrocarbon and arsenic concentrations in groundwater, Poso Creek Oil Field, California, USA

Like many hydrocarbon production areas in the U.S., the Poso Creek Oil Field in California includes and is adjacent to other land uses (agricultural and other developed lands) that affect the hydrology and geochemistry of the aquifer overlying and adjacent to oil development. We hypothesize that the distributions of hydrocarbons and arsenic in groundwater in such areas will be controlled by complex interactions between mixed land uses, oil-field infrastructure, and natural processes. In 2020–2021, samples of groundwater and surface water were collected and analyzed for a large suite of inorganic and organic chemicals and isotope and gas tracers to test this hypothesis. Those data are supplemented with ancillary data on historical geochemistry, hydrology, geology, and oil-field infrastructure. Hydrocarbons in groundwater (e.g., methane through pentane gases and benzene) are associated with natural processes (e.g., fault offsets or transition in sediment depositional environment) and oil-field infrastructure (e.g., fluid-migration pathways associated with uncemented annulus in oil wells or unlined pits). Arsenic concentrations >10 μg per liter (μg/L; maximum concentration 12.9 μg/L) are associated with natural processes in old, high-pH groundwater, and more recent recharge of water from natural and/or engineered recharge processes. Along the southwest margin of the oil field, pumping for drinking-water and irrigation supplies in combination with engineered groundwater recharge produce a depression in groundwater elevations where groundwater with elevated sulfate concentrations from agricultural areas and groundwater with hydrocarbons from the oil field mix to produce a zone of sulfate reduction that removes hydrocarbons and arsenic from groundwater but produces elevated sulfide (S 2- ) concentrations (maximum concentration 29 mg per liter, mg/L). In this study, multiple approaches were required to resolve the overlapping effects of land uses, oil-field infrastructure, and natural processes on the distributions of hydrocarbons and arsenic in groundwater. The combined use of geographic, historical, physical, chemical, isotopic, and other information to constrain processes could be a useful approach for studies in other hydrocarbon-production areas. This is particularly important where land uses affect aquifer hydrology to an extent that causes mixing of groundwaters with different chemical compositions.

California↗
USGS2025

Structural setting and geothermal potential of northeastern Reese River Valley, north-central Nevada: Highly prospective detailed study site for the INGENIOUS project

The northeastern part of the Reese River basin situated ~15 km southeast of Battle Mountain, Nevada, scored highly in the Nevada geothermal play fairway analysis (PFA) for hosting potential hidden geothermal systems. This site (also referred to as Argenta Rise) was therefore chosen for detailed study in the INGENIOUS project (INnovative Geothermal Exploration through Novel Investigations Of Undiscovered Systems). The high PFA scores resulted primarily from favorable structural settings (e.g., fault intersections and pull aparts) with relatively high slip rates on Quaternary faults. The INGENIOUS project is utilizing additional parameters and more rigorous analytical techniques to further advance exploration at this site. This includes integration of geological (e.g., Quaternary fault mapping) and new geophysical datasets (e.g., gravity, magnetics, MT data, and five reprocessed seismic reflection profiles) to build a structural model and to identify specific favorable sites for potential geothermal upwellings. Two-meter temperature surveys were also conducted in the area (139 measurements). This part of north-central Nevada is characterized by systems of intersecting northerly and ENE-striking faults within the broader Humboldt structural zone, a poorly understood belt of ENE-striking faults and relatively high heat flow extending across northern Nevada. Kinematic analysis of exposed fault surfaces shows that ENE-striking faults have accommodated sinistral-normal slip, and normal slip characterizes N- to NNE-striking faults. Northeastern Reese River Valley lies within a broad left step between major ENE-striking fault zones on the northern flanks of the Argenta Rim and Shoshone Range and thus corresponds to a broad pull-apart in the ENE-striking sinistral-normal fault system. Notably, the nearby Beowawe geothermal system in Whirlwind Valley (with abundant sinter, hot springs, and a geothermal power plant) occupies a fault intersection in a relatively small left step in a major ENE-striking sinistral-normal fault and may serve as an analogue for a potential hidden system in northeastern Reese River Valley. Existing geological maps, high-resolution lidar, and seismic reflection data demonstrate that northeastern Reese River Valley is structurally complex with multiple intersections between the ENE- and N- to NNE-striking fault systems. Some of these fault intersections correspond to low resistivity anomalies, magnetic lows, and/or very subtle 2-m temperature anomalies, which may indicate hidden geothermal upwellings. Three-dimensional modeling and temperature-gradient drilling are planned to further evaluate these sites for geothermal activity.

Nevada↗
USGS2025

Fossil diatoms in Arctic marine surface sediments

Diatoms are one of the main constituents of marine phytoplankton in the Arctic, and thanks to their siliceous skeletons, diatom fossils are relatively well preserved in sediments. Due to their species-specific sensitivity to different ocean conditions, their abundance and assemblages in sediments are routinely used by paleoceanographers to reconstruct the state of the surface ocean in the past (Koç 2007). By using statistical methods, we can obtain valuable knowledge about their ecological preferences (Oksman et al. 2019) and generate quantitative reconstructions of various parameters, such as sea-surface temperature and sea-ice concentration through time (Krawczyk et al. 2021; Sha et al. 2014). The Marine Arctic Diatoms (MARDI) working group (WG) ( pastglobalchanges.org/mardi ) aims to advance knowledge on marine-diatom ecology and diatom-based reconstructions by compiling and harmonizing data from surface-sediment samples across the Arctic.

Pages Magazine↗
USGS2026

Perspectives and insights gained from applying landscape ecological approaches to the Upper Mississippi riverscape

Advances in remote sensing, spatial data collection, hydraulic and geospatial models, and the maturation of landscape ecology during the early 2000s spurred the concept of rivers as landscapes or ‘riverscapes’. Since that time, the Upper Mississippi River Restoration (UMRR) Program - one of the world’s longest running large-river ecosystem restoration and monitoring programs - has applied landscape ecological principles to the Upper Mississippi River System (UMRS). I examine spatial patterns of aquatic nutrient concentrations, submersed aquatic vegetation, and fish and mussel communities, finding that within-river patterns are often patchy in nature, likely resulting from how nutrient processes and community traits interact with spatial variability in hydraulic exchange patterns. Local-scale patchy spatial patterns in the UMRS are embedded in system-scale longitudinal gradients and a hierarchical view of the river has provided a foundation for strategic river restoration planning. Key Words: Aquatic Habitat, hydraulic connectivity, patch, restoration, riverscape

BioScience↗
USGS2025

40Ar/39Ar geochronology supporting mineral resources research at USGS Denver

The 40 Ar/ 39 Ar geochronology method is used to date potassium-bearing rocks and minerals, based on the decay of 40 K to 40 Ar, which provides important temporal constraints for geological events. The USGS Denver Argon Geochronology Laboratory dates samples from a variety of projects, mainly in the USGS Mineral Resource Program and the National Cooperative Geologic Mapping Program, facilitating in-depth research into the timing of geological events and processes. Recent applications of this method include geochronological studies in the Yellow Pine district, revealing mineralization ages that range from 51 to 70 Ma, and investigations into the Bear Lodge alkaline complex, providing insights into the timing of REE carbonatite formation. Additionally, studies at Alunite Ridge highlight the method's potential for understanding concealed mineral deposits. Ongoing projects continue to leverage 40 Ar/ 39 Ar data to enhance understanding of geologic frameworks and mineral resources, underpinning its importance in modern geochronology and mineral resource assessment.

Conference Paper↗
USGS2025

Estimating the probability of export restrictions to inform mineral criticality

To assess risks associated with advanced technologies’ supply chain disruptions, governmental agencies and others have developed mineral “criticality” assessments, with criticality described using the economic impact and probability of supply chain disruptions. Previous work developed subjective supply risk indicators to approximate this probability, typically combining several factors such as supply diversity and trading partners’ political stability, where indicator weightings can substantially impact results. This work explicitly quantifies export barrier probability using an ensemble of machine learning classifiers, with probability estimates informed by exogenous variables, including prior barrier implementation and global export dominance. Major differences in high-probability countries and commodities are observed across models, but the ensemble method highlights Indonesia, China, Tanzania, and the United States as particularly high risk. The Supplementary Data File provides export barrier probability estimates for each analyzed country-commodity pair, enabling a direct, quantitative, objective contribution to assessing mineral criticality, enhancing risk identification and prioritization for policymakers.

Resources, Conservation, and Recycling↗
USGS2026

Predictable seismic cycles result from structural rupture barriers on oceanic transform faults

Earthquakes of magnitude ( M ) >5.5 on oceanic transform faults (OTFs) repeatedly rupture the same locked patches, sometimes quasiperiodically. These patches are separated by “barriers” that halt earthquake propagation and slip mostly aseismically. However, the physical processes governing this systematic behavior remain unclear. We analyzed two barriers along the Gofar transform fault that have arrested ~15 M 6 earthquakes over the past three decades. Ocean bottom seismometer data indicate that the barriers hosted intense microseismicity before the mainshocks and comprise multistrand faults and transtensional stepovers with 100- to 400-m lateral offset. These characteristics contradict earthquake rupture termination models invoking velocity-strengthening friction or large geometric steps and instead point to damage-enhanced porosity and dilatancy-strengthening mechanisms. By isolating rupture segments, the barriers regulate the quasiperiodic recurrence of OTF earthquakes.

Science↗
USGS2025

Flow heterogeneity controls dissolution dynamics in topologically complex rocks

Rock dissolution is a common subsurface geochemical reaction affecting pore space properties, crucial for reservoir stimulation, carbon storage, and geothermal energy. Predictive models for dissolution remain limited due to incomplete understanding of the mechanisms involved. We examine the influence of flow, transport, and reaction regimes on mineral dissolution using 29 time-resolved data from 3D rocks. We find that initial pore structure significantly influences the dissolution pattern, with reaction rates up to two orders of magnitude lower than batch conditions, given solute and fluid-solid boundary constraints. Flow unevenness determines the location and rate of dissolution. We propose two models describing expected dissolution patterns and effective reaction rates based on dimensionless metrics for flow, transport, and reaction. Finally, we analyze feedback between evolving flow and pore structure to understand conditions that regulate/reinforce dissolution hotspots. Our findings underscore the major impact of flow arrangement on reaction-front propagation and provide a foundation for controlling dissolution hotspots.

Geophysical Research Letters↗
USGS2024

Collision structures of the Prince William terrane and Chugach terrane docking along the Shumagin and Unimak convergent margins, Alaska, USA

Western Alaska’s convergent margins are composed of tectonostratigraphic terranes. On land, terrane assembly is recognized along boundaries or sutures between neighboring geologic elements with distinctly different origins. In marine areas where rock outcrops are covered by sediment, recognizing terrane sutures is problematic. A fault in seismic dip line 5 of the ALEUT project has been interpreted as a terrane suture. It is imaged intermittently down to the 30+-km-deep plate interface. Processing of ALEUT strike line 7 revealed the suture at ~18 km depths extending 300 km along the margin. Upper structures in line 5 are like the structures of adjacent seismic transects where imaging is only 8−10 km deep. They were previously not recognized as the upper reaches of terrane sutures and show structural details obscured at greater depths. The composite data are the basis for a simple tectonic model of terrane docking.

Alaska↗
USGS2025

Rupture process of the Mw7.0 December 5, 2024 Offshore Cape Mendocino earthquake

The Mw7.0 December 5, 2024 Offshore Cape Mendocino earthquake ruptured a km long portion of the east-west trending Mendocino fault zone (MFZ). In order to clarify the rupture process, we assemble three-component seismograms from regional seismic stations, horizontal coseismic displacement vectors derived from Global Navigation Satellite System (GNSS) time series, and a Sentinel-1 ascending interferogram. These data are interpreted with a model of slip distributed on two vertical fault planes representative of the eastern MFZ and spanning the ~70 km length of the aftershock zone. Assuming right-lateral strike slip, we find that the rupture initiates in the oceanic mantle at 20-30 km depth and proceeds unilaterally updip and toward the east. Early aftershocks locate adjacent to the peak slip areas, tracking the coseismic rupture propagation from oceanic mantle to shallower depth and implying a significant role of static stress transfer in driving aftershocks in an ocean plate environment.

California↗
USGS2025

Persistent deep long-period seismicity near the Lassen Volcanic Center

Deep long-period (DLP) earthquakes have been observed at many volcanic settings around the world and linked to the magmatic processes that drive volcanic unrest. At the Lassen Volcanic Center (LVC) of the Cascade arc, limited detection of DLP activity hinders classification of anomalous seismic behavior and its relationship to the LVC magmatic system. This study uses a template matching approach with seismic data from temporary nodal and permanent stations to detect and locate DLP earthquakes near the LVC between 2017 and 2024. Within the DLP catalog of 611 events, a transition occurs from scattered (2017–2020) to oscillatory (2021–2024) occurrence rates. During the oscillatory period, regional earthquakes with large amplitude velocity waveforms observed near the LVC are associated with abrupt changes in DLP occurrence rates. Continued monitoring of DLP activity has the potential to better define the processes that drive volcanic unrest at the LVC in the future.

California↗
USGS2026

Dynamic coupling between faulting, rifting and magmatism during 2021-2025 unrest on Reykjanes Peninsula, Iceland

Interactions among faulting, earthquakes, and eruptions are fundamental to plate tectonics and hazard forecasting yet rarely observed along mid-ocean ridges. On Iceland's Reykjanes Peninsula, seismotectonic–volcanic unrest resumed after nearly 800-year hiatus, providing an opportunity to observe these interactions during 2021–2025 activity. By integrating high-resolution seismicity, focal mechanisms, satellite geodesy, surface deformation, and eruption data, we document ∼4 m of total extension accommodated through 14 rifting episodes. The largest, in 2023, involved graben reactivation and diking, with seismic swarms and earthquake faulting that matched the surface ruptures, where strike-slip faulting preceded normal-faulting earthquakes and extension. The accrued extension was released by extension fractures triggered by magma accumulation. Long-term observations show no correlation between erupted magma volume, seismicity, and crustal extension. This highlights dynamic relation between rifting, faulting, and magmatism in transtensional settings and their implications for hazard assessment.

Reykjanes Peninsula↗
USGS2025

Reply to, “Comment on ‘The 1886 Charleston, South Carolina, earthquake: Relic railroad offset reveals rupture,’ by Roger Bilham and Susan E. Hough”

We welcome this opportunity to respond to Pratt et al. (2024) (hereinafter P24). Bilham and Hough (2023) proposed a “first-cut” elastic deformation model for the 1886 earthquake, a quantitative source model constrained by identified coseismic constraints. A key observation was the measurement of a lateral offset of a railroad line south of Summerville, leading to a model with predominately dextral slip and minor convergence, from which we concluded that active faulting had raised the Penholoway Marine Terrace >6 m since ∼770 ka. P24 questioned these constraints and proposed an alternative rupture model with predominantly reverse slip. This alternative model is neither consistent with coseismic constraints nor with other geophysical data. In a revised model presented here, we recognize that uplift of the Penholoway Terrace is confined to the eastern edge of the terrace, which we conclude results from active folding and tectonic transpression centered on the dextral fault that offset the railroad in 1886.

South Carolina↗
USGS2026

Generating geochemical and mineralogy distributions of soil in the conterminous United States using Bayesian hierarchical spatial models

Characterizing geochemical and mineralogical soil distributions across large spatial extents is essential for understanding mineral resources, ecosystem processes, and environmental risks. Rasters of soil geochemical distributions for the conterminous United States, however, are limited. We present a Bayesian modeling workflow and tool for generating predictive geochemical and mineralogy distribution maps for the conterminous United States using integrated nested Laplace approximation (INLA) with the stochastic partial differential equation approach. By modeling soil geostatistical data with environmental covariates (soil properties, topography, climate, and land cover), we generate predictive distributions of soil geochemistry that can be mapped or extracted for further analyses. As an example, we model the spatial distribution of trace elements in soil relevant to vertebrate health (cobalt, copper, iron, manganese, selenium, and zinc) and provide a workflow that can be used to generate and visualize predictive distributions of 39 other major and trace elements and 21 minerals of the soil survey, supporting a variety of ecological, environmental, and agricultural applications.

MethodsX↗
USGS2024

Summary of annual repeat magnetotelluric surveys of the Geysers geothermal field

As part of a multi-year project, funded by the California Energy Commission, annual magnetotelluric (MT) surveys have been collected at The Geysers geothermal field in northern California with the goal of measuring temporal changes within the steam field. The repeat surveys started in 2021 and repeated a 2017 survey (Peacock et al., 2020) with further extension to the southern part of the geothermal field. Temporal variations in the MT transfer functions are observed to be spatially coherent and compartmentalized. Mapping residual phase tensor ellipses demonstrates the direction of maximum change is often aligned with existing fracture orientations. Three dimensional inversion of the MT data, using the inversion results from previous years as the starting model, indicates that the steam reservoir has generally become more resistive over time (~10%), suggesting more steam in the field. A few pockets within the steam field have become more electrically conductive over time and are collocated with injection wells, suggesting either more fluid content in those zones, less steam, or more saline fluids.

California↗