Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Mining, Metallurgy and Exploration”

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

Regional characterization of coal resources in the U.S. Gulf Coast

There is increasing interest in extracting critical minerals (CM), including rare earth elements (REE), from coals in the United States to address the overreliance on imported REE. The U.S. Gulf Coast and the Williston basins are the two major lignite-bearing basins within the country. Recent REE and CM studies of the lignite in these basins have indicated that the coals may be a viable source material for REE and CM extraction. To evaluate in-place coal as a potential source of REE and CM, the coal resources need to be quantified. This study presents the results of a regional analysis of the U.S. Gulf Coast lignite and bituminous coal resources that might be available as potential sources of REE and CM. The resource analysis used kriging methods to develop isopleth maps of cumulative coal thickness throughout the region using data from 31,181 drill holes and other data points. The estimated total coal resource in the Gulf Coast is about 83 billion metric tons in the upper 90 m (~ 300 ft) of the subsurface. Texas accounted for 40 percent (32 billion metric tons) of the total resource, followed by Mississippi (24 %, 20 billion metric tons), Louisiana (14 %, 12 billion metric tons), Tennessee (10 %, 8.5 billion metric tons), and Arkansas (6 %, 5.1 billion metric tons). The remaining states each accounted for less than 5 percent of the total resource. Georgia had the smallest resource estimated at 7 million metric tons. Here we report the first known state-wide lignite resource estimates for Georgia, Kentucky (820 million metric tons), and Missouri (1,800 million metric tons). A comparison of the results of this study with those of previous Gulf Coast and Williston Basin resource studies is difficult because each study used different data sources, assessment methodologies, overburden depths, and qualifying coal thicknesses. Coal-power electric generation has sharply decreased in past decades and mining of these coals for CM and REE could provide additional co-products such as activated carbon and other uses such as fertilizer (soil enhancer).

EarthArXiv↗
USGS2024

Long-term communication of aftershock forecasts: The Canterbury earthquake sequence in New Zealand

On 14 February 2016, a magnitude (M)5.7 earthquake struck in Christchurch New Zealand (Aotearoa in the Maori language). The shaking caused damage to historic facades, power outages, cliff collapses, rock falls, and liquefaction but no reported injuries or fatalities. This Valentine’s Day earthquake was an aftershock in the Canterbury earthquake sequence (CES), which began on 4 September 2010 with the M7.1 Darfield Earthquake and included the destructive and fatal M6.2 Christchurch aftershock on 22 February 2011. This study, eight months after the Valentine’s Day earthquake and six years after the initiation of the CES, is the first to explore long-term aftershock forecast information and communication needs. The exploratory study also aimed to gather feedback on aftershock scenarios, an alternative form for communicating the forecast. The qualitative study involved workshops with emergency managers, public health officials, and members of the public in Christchurch. Key findings for long-term communication throughout an earthquake sequence include: 1. divergent earthquake experiences affect aftershock communication response and information needs; 2. understanding aftershock sequence behavior is foundational to sense-making when large aftershocks occur; 3. strategic earthquake sequence updates from the trusted science agency and local agencies could serve as important reminders for earthquake preparedness; 4. communication of aftershock forecast uncertainty could aid with both the credibility of the information and living with uncertainty, and 5. inclusion of impact information and preparedness advice into aftershock forecast scenarios could provide links to actionable information. The paper derives implications for research and practice of long-term communications during an aftershock sequence.

Christchurch↗
USGS2025

Structural analysis of brittle-plastic shear zones in the Sangre de Cristo Range, southern Colorado USA: Superposition of Rio Grande rift extension on Laramide contraction

The Sangre de Cristo Range in southern Colorado exposes some of the deepest Cenozoic structural levels in the Rocky Mountain region, including mylonitic shear zones associated with both the Laramide orogeny and Rio Grande rift. We investigated the relation between Laramide contraction and Rio Grande rift extension with detailed geologic mapping, kinematic analysis, and geochronometry in a 50 km 2 area centered on the Independence Mine shear zone (IMSZ). The 15−100-m-thick IMSZ is one of several shallowly to moderately (∼45° ± 20°) W-SW−dipping brittle-plastic shear zones along the western flank of the range. These shear zones display microstructural evidence of initiation as top-NE contractional mylonite zones, consistent with regional Laramide kinematics, which have been pervasively overprinted by shear fabrics indicating top-SW extensional reactivation. Both top-NE and top-SW shear fabrics involve cataclasis and quartz dislocation creep, although top-SW shear is more commonly localized along phyllosilicate-lined shear bands. Shear zones are hosted predominately within Proterozoic gneiss, and contain abundant chlorite and white mica derived from alteration of hornblende and feldspar, which indicates that weakening driven by fluid reactions played an important role in localizing strain. Extensional overprinting appears to be most pervasive along more steeply dipping portions of shear zones and where secondary phyllosilicates form an interconnected weak phase, which suggests that reactivation was primarily controlled by geometry and rheological contrasts inherited from contraction. One top-SW shear zone adjacent to the IMSZ cuts a late Oligocene gabbro stock, and monazite grains synkinematic with top-SW shear in the IMSZ yielded late Oligocene to Early Miocene U-Th-Pb dates that correspond with initiation of the Rio Grande rift. Reactivation of weak reverse faults may represent an important structural control during initial extension in the middle crust, prior to slip along the high-angle Sangre de Cristo normal fault system.

Colorado↗
USGS2024

A decade of shaking in the Garden City: The dynamics of preparedness, perceptions, and beliefs in Canterbury, New Zealand, and implications for earthquake information

This study explored earthquake preparedness over time - before, during, and 10 years after the Canterbury Earthquake Sequence (CES) in Aotearoa New Zealand (NZ). Surveys of Canterbury residents were conducted in 2009, 2013 and again in 2021, using variables derived from Community Engagement Theory (CET). The surveys measured earthquake perceptions and beliefs, participation and engagement, and preparedness actions. Results were compared across the three samples. Findings indicate that perceptions and beliefs (e.g. risk perception, outcome expectancy beliefs), and types of preparedness actions taken (e.g. collection of survival items, structural preparedness, community and agency relationships), differed over time, depending on people’s experiences before, during, and after the CES. For example, during and after the CES people were more likely to believe that preparing provided a benefit to daily life, but less likely to think it could reduce property damage, perhaps due to people’s experiences of disruption and damage during the earthquakes. An understanding of such dynamics can assist with the provision and timing of risk and preparedness information. This study highlights the importance of providing applicable and actionable preparedness information, that is relevant to people’s experiences, throughout an earthquake sequence. Such information might evolve and change in focus over time depending on risks and needs. Focus could also be given to information that builds peoples beliefs and capacities to undertake preparedness in evolving situations. Understanding preparedness in the context of different experiences and timeframes is useful in helping update models such as the CET, where the dynamics of time might be better incorporated.

Canterbury↗
USGS2025

High frequency and region-scale simulations of large (Mw7+) earthquakes on the southern Whidbey Island fault, Washington, USA

We simulate ground shaking in western Washington State from hypothetical M w 7.0–7.5 earthquakes on the southern Whidbey Island fault (SWIF). Ground motions are modeled considering kinematic source distributions on a complex fault plane, a 3D seismic velocity model, and region‐specific soil velocity models. We run simulations with varying model resolutions, including regional‐scale simulations with a maximum‐modeled frequency of ∼1 Hz and local‐scale simulations with a maximum‐modeled frequency of ∼2.5 Hz. Additional local‐scale simulations are run considering high‐resolution surface topography. We explore how source parameters (i.e., magnitude, hypocenter location, and dip direction) and 3D velocity structure impact peak shaking intensity and its variability. In particular, we find that earthquakes on the SWIF would likely produce strong shaking throughout the populated Puget Lowland, including in the cities of Everett, Seattle, Bellevue, and Tacoma, Washington. Simulated short‐period ( T ≤ 2 s) spectral accelerations are strong throughout the Puget Lowland, and long‐period shaking ( T ≥ 5 s) is strong in the deep regional sedimentary basins, especially the Everett and Seattle basins. Source parameters strongly influence intra‐ and interevent variability in response, primarily through changes in source and site geometry, as well as rupture directivity. We also note a potential coupling between rupture directivity and basin effects, wherein directivity pulses are seemingly guided through the region’s deep, interconnected sedimentary basins. Overall, this work highlights the impacts of 3D source, path, and site effects on seismic hazard in the U.S. Pacific Northwest and substantially expands the catalog of simulated ground motions for Puget Sound area crustal faults.

Washington↗
USGS2026

Lithium brine potential of North Dakota: Results from the 2025 sampling program

Lithium (Li), classified as a critical mineral and key in energy storage applications because of its high energy density, faces unprecedented demand. This has driven interest in unconventional sources of Li, especially brines from sedimentary basins. Accordingly, the spatial distribution of basinal brinehosted Li resources across North America is becoming better defined, including the Williston Basin underlying Saskatchewan, Canada. This has led to growing investment and exploration by industry throughout the 2020s; however, despite this interest, the Li potential of basinal brines of the Williston Basin in North Dakota remain underexplored. This work presents the results from the 2025 joint brine sampling campaign between the North Dakota Geological Survey, the United States Geological Survey (USGS), and The University of Texas at Austin. Twenty-seven fluid samples were collected from producing oil and water source wells from eight stratigraphic intervals across western North Dakota and analyzed for major cation and anion concentrations. Similar to southeastern Saskatchewan, the highest Li concentrations were identified in the Devonian Duperow Formation, with concentrations up to 179 mg/L being observed in a well commingled between the Duperow and Red River Formations. Some samples from the overlying Birdbear Formation exceeded 75 mg/L, whereas Li concentrations were lower in the Madison Group (33-44 mg/L) and the Three Forks (50 mg/L), Dawson Bay (43 mg/L), Winnipegosis (23 mg/L), and Red River (40.0 mg/L) Formations. This work indicates Li resource potential exists across the Williston Basin and advances our understanding of the Li distribution in basinal brines which could provide a new domestic source of this critical mineral.

North Dakota↗
USGS2026

Uncertainty quantification of geologic energy storage in depleted gas reservoirs using material balance equations embedded in a hierarchical errors-in-variables model

The storage potential for gas in geologic settings, such as depleted hydrocarbon reservoirs and solution-mined salt caverns, is becoming salient to future energy infrastructure planning. Technologies such as carbon capture, utilization, and storage, carbon dioxide-enhanced oil recovery, and natural gas and hydrogen storage help meet growing energy demands, reduce carbon emissions to meet climate goals, and provide energy security amid geopolitical uncertainties. Therefore, estimates of underground gas storage capacity could be useful for efficiently navigating the energy transitions. Material balance is a fundamental method in reservoir engineering for estimating original gas in place and potential storage capacity at the scale necessary for national assessments of subsurface pore space resources. However, the deterministic method cannot accommodate multiple data sources or quantify uncertainty in predictions. In this study, a method that embeds material balance equations within a hierarchical errors-in-variables model is proposed which allows the estimation of the distributions of reservoir properties needed for assessments. Uncertainties associated with these reservoir properties have traditionally been expert-elicited, whereas the uncertainty estimates from the proposed models are data-driven. Capacity and uncertainty estimates can be used in a probabilistic resource assessment, supplementing information traditionally used by assessors or even replacing this expert elicitation step when data are unavailable. Various regression models are compared in a case study of the Michigan Basin, a large contributor to the United States’ current natural gas storage capacity. In particular, errors-in-variables models help ameliorate regression dilution and can quantify uncertainty in predictions of pressure in addition to storage capacity. Overfitting is addressed by quantifying generalization error and model averaging in simple and stratified cross-validation against reported working gas capacity, representing the varying quality and quantity of available data. Incorporating a statistical framework into existing numerical methods in reservoir engineering can improve the quality of estimation, and in particular, this method brings rigor to uncertainty quantification as part of a larger effort by the U.S. Geological Survey to assess domestic energy gas storage resources in depleted hydrocarbon reservoirs.

Mathematical Geosciences↗
USGS2026

The Great Salt Lake (Utah, USA) acoustic and sedimentary archive of Wasatch Front earthquakes

As the largest low-relief basin in the western United States, the Great Salt Lake offers an opportunity to explore acoustic and sedimentary archives of earthquake ruptures in a shallow terminal basin. Here, we use subbottom compressed high intensity radiated pulse (Chirp) profiles and sedimentary cores to address questions related to how Great Salt Lake sediments respond to surface rupture and earthquake shaking and whether the lake is sensitive to both Great Salt Lake fault (GSLF) and Wasatch fault zone (WFZ) earthquakes. Chirp profiles reveal laterally continuous seismic stratigraphy with vertically displaced strata, onlap surfaces, and growth stratigraphy produced by GSLF surface ruptures. We identify and date two ruptures on the Fremont Island section and three on the Antelope Island section; temporal overlap of these ruptures supports synchronous rupture of both sections. Combined with previous data, our results yield a record of four Holocene surface-rupturing (M W 6.9−7.3) earthquakes on the GSLF with a mean recurrence of ∼3.2 k.y. Sediment cores include ≤41-cm-thick, massive, homogeneous to normally graded disturbance beds that interrupt laminated bioclastic sediment. Disturbance bed correlations yield four lacustrine earthquake deposits since ca. 1.7 ka; the youngest possibly records the most recent GSLF surface rupture whereas older deposits were likely produced by WFZ earthquakes. Modeled historical earthquakes yield a Great Salt Lake sediment disturbance shaking threshold of MMI ∼VI, controlled in part by the lake-bottom expression of the GSLF scarp. Our integration of primary (fault-related) and secondary (shaking-related) earthquake evidence improves our understanding of sediment disturbance processes and the archival potential of low-relief lake basins.

Utah↗
USGS2024

Bioconcentration of per- and polyfluoroalkyl substances and precursors in fathead minnow tissues environmentally exposed to aqueous film-forming foam-contaminated waters

Exposure to per- and polyfluoroalkyl substances (PFAS) has been associated with toxicity in wildlife and negative health effects in humans. Decades of fire training activity at Joint Base Cape Cod (MA, USA) incorporated the use of aqueous film-forming foam (AFFF), which resulted in long-term PFAS contamination of sediments, groundwater, and hydrologically connected surface waters. To explore the bioconcentration potential of PFAS in complex environmental mixtures, a mobile laboratory was established to evaluate the bioconcentration of PFAS from AFFF-impacted groundwater by flow-through design. Fathead minnows ( n = 24) were exposed to PFAS in groundwater over a 21-day period and tissue-specific PFAS burdens in liver, kidney, and gonad were derived at three different time points. The ∑PFAS concentrations in groundwater increased from approximately 10,000 ng/L at day 1 to 36,000 ng/L at day 21. The relative abundance of PFAS in liver, kidney, and gonad shifted temporally from majority perfluoroalkyl sulfonamides (FASAs) to perfluoroalkyl sulfonates (PFSAs). By day 21, mean ∑PFAS concentrations in tissues displayed a predominance in the order of liver > kidney > gonad. Generally, bioconcentration factors (BCFs) for FASAs, perfluoroalkyl carboxylates (PFCAs), and fluorotelomer sulfonates (FTS) increased with degree of fluorinated carbon chain length, but this was not evident for PFSAs. Perfluorooctane sulfonamide (FOSA) displayed the highest mean BCF (8700 L/kg) in day 21 kidney. Suspect screening results revealed the presence of several perfluoroalkyl sulfinate and FASA compounds present in groundwater and in liver for which pseudo-bioconcentration factors are also reported. The bioconcentration observed for precursor compounds and PFSA derivatives detected suggests alternative pathways for terminal PFAS exposure in aquatic wildlife and humans.

Environmental Toxicology and Chemistry↗
USGS2025

Geochemistry and radiogenic isotopes constrain the mantle source region of the Mountain Pass Intrusive Suite, California

The Mountain Pass carbonatite stock is the largest rare earth element (REE) deposit and only active REE mine in the United States. The carbonatite intrusion and spatially associated alkaline silicate intrusions constitute the Mountain Pass Intrusive Suite, which is located within the Mojave Province in California. Both the carbonatite and the alkaline silicate rocks are enriched in large ion lithophile elements and light REEs, and less enriched to depleted in high field strength elements, indicating the mantle source region was metasomatically enriched in incompatible trace elements. The cause of this metasomatic mantle enrichment and the genetic relationship between the carbonatite and the alkaline silicate stocks are poorly understood. In this study, major and trace element geochemical data and isotopic (Rb-Sr, Sm-Nd, and Lu-Hf) data are presented to constrain genesis of the Mountain Pass Intrusive Suite, from mantle source region to the intrusion of the stocks. Our geochemical data are consistent with derivation of the alkaline silicate and carbonatite melts through partial melting from a shared mantle source region rather than through liquid immiscibility or fractional crystallization and separation of a carbothermal fluid. Although the Rb-Sr isotopic system in the Mountain Pass Intrusive Suite is disturbed at the whole-rock scale, the isotopic systems for whole-rock Sm-Nd (εNd i = ‐2.2 ± 0.8) and zircon Lu-Hf (εHf i = 0.1 ± 1.1) are robust and support mantle derivation of the magmas. Geochemical modeling using experimentally derived partition coefficients was used to identify possible causes of enrichment in incompatible elements through metasomatism in the mantle source region. Modeling of metasomatism by melts derived by partial melting of deeply subducted carbonated sediments approximates observed Mountain Pass Intrusive Suite trace element chemistry. Scattered εHf i in inherited zircon (2.8 ± 2.6) is consistent with derivation from an arc-related environment with substantial crustal contamination. Paleotectonic studies in the Mojave Province indicate that regional subduction preceded emplacement of the Mountain Pass Intrusive Suite by ∼300 Ma. Melting of the Mountain Pass source region may have been caused by post-collisional thermal relaxation and extension.

California↗
USGS2025

Dynamic baseflow storage estimates and the role of topography, geology and evapotranspiration on streamflow recession characteristics in the Neversink Reservoir Watershed, New York

Estimates of dynamic groundwater volumes supplying baseflow to streams are important for water availability projections during extended periods of drought. The primary goals of this study were to provide dynamic storage volume estimates, inferred from streamflow recession analysis, for baseflow regimes within seven gaged catchments within the Neversink Reservoir Watershed (NRW), a critical municipal water source for New York City. Additionally, geomorphological properties, surficial geology and hydro-meteorological processes were quantified and described in relation to time and spatially variable recession behaviour and storage estimates across the NRW. To explore these relationships, we (1) evaluated seasonal trends in streamflow recession behaviour in relation to modelled potential evapotranspiration (PET) and catchment runoff rates, (2) derived empirical streamflow models for cool-season runoff using both linear and nonlinear reservoir assumptions for baseflow and (3) calculated metrics related to the geology and geomorphology of each catchment and compared these metrics to area normalised baseflow dynamic storage estimates. Results show that baseflow recession behaves as a nonlinear reservoir, and applying linear groundwater reservoir assumptions may underestimate the total dynamic storage volumes compared to what would be predicted for a nonlinear reservoir. Increases in PET caused decreases in storage conditions that resulted in increased recession rates and nonlinearity in streamflow recession during the growing season. Additionally, we found that while no single physical catchment characteristic solely predicted catchment storage dynamics, sediment volume and stream gradients were stronger predictors of normalised storage volumes than catchment surface area or surface topography alone. Within the NRW, catchments with the highest sediment volume exhibited the lowest recession rates and higher dynamic storage volumes, while the smallest catchment, mostly devoid of sediment, had the fastest recession rate and lowest dynamic storage volume.

New York↗
USGS2026

Tracking baseflow supply dynamics using SWOT data from small groundwater-dominated lakes

In situ surface-water monitoring strategies are biased towards larger perennial streams and lakes and are generally not designed to track mechanisms of baseflow supply contributed by the dynamic storage of aquifers. Additionally, small (< 1 km 2 ) groundwater-influenced lakes and wetlands globally have little in situ monitoring infrastructure. We explored the utility of remotely sensed Surface Water Ocean Topography Satellite (SWOT) data, collected from 2023 onward, to characterise the seasonal and multi-year water-level trends of groundwater flow-through kettle lakes distributed across the permeable sediments of eastern Massachusetts, USA. This analysis indicated that water levels for kettle lakes with areas down to approximately 0.05 km 2 are resolvable in the study area. Our examination of 17 kettle lakes found that SWOT water-surface elevation data closely tracked groundwater levels in adjacent monitoring wells where available, including the timing of seasonal patterns (highest levels generally in late spring), although there was some variation between years and there was a substantial lag in the timing of high water levels for a lake located downgradient from a 30-m-thick vadose zone. Furthermore, SWOT-observed water-level increases in kettle lakes tracked with baseflow increases in two adjacent groundwater-dominated streams, as would be expected from increased hydraulic gradients. Unlike spectral remote sensing, SWOT data are generally not affected by cloud cover, resulting in a potential for groundwater-dominated lakes to be sentinels of dynamic storage patterns, including identification of baseflow drought lags, which are currently ill-defined hydrological processes. SWOT monitoring of groundwater-influenced surface waters shows potential for augmenting existing monitoring wells and streamgages as continuous monitors of groundwater levels and baseflow supply in permeable terrain.

Massachusetts↗
USGS2025

Fiber-optic sensing for earthquake hazards research, monitoring and early warning

The use of fiber‐optic sensing systems in seismology has exploded in the past decade. Despite an ever‐growing library of ground‐breaking studies, questions remain about the potential of fiber‐optic sensing technologies as tools for advancing if not revolutionizing earthquake‐hazards‐related research, monitoring, and early warning systems. A working group convened to explore these topics; we comprehensively examined the application of fiber optics in various aspects of earthquake hazards, encompassing earthquake source processes, crustal imaging, data archiving, and technological challenges. There is great potential for fiber‐optic systems to advance earthquake monitoring and understanding, but to fully unlock their capabilities requires continued progress in key areas of research and development, including instrument testing and validation, increased dynamic range for applications focused on larger earthquakes, and continued improvement in subsurface and source imaging methods. A key current stumbling block results from the lack of clear data archiving requirements, and we propose an initial strategy that balances data volume requirements with preserving key data for a broad range of future studies. In addition, we demonstrate the potential for fiber‐optic sensing to impact monitoring efforts by documenting the data completeness in a number of long‐term experiments. Finally, we outline the features of a instrument testing facility that would enable progress toward reliable and standardized distributed acoustic sensing data. Overcoming these current obstacles would facilitate progress in fiber‐optic sensing and unlock its potential application to a broad range of earthquake hazard problems.

California, Oregon, Washington↗
USGS2024

Wildfire, extreme precipitation and debris flows, oh my! Channel response to compounding disturbances in a mountain stream in the Upper Colorado Basin, USA

Compounding changes in climate and human activities stand to increase sediment input to rivers in many landscapes, including via discrete perturbations such as post-fire debris flows. Because sediment supply is a dominant control on river morphology, understanding mountain river responses to sediment regime perturbations is critical to predicting and addressing downstream effects to infrastructure, water security and aquatic habitat. A growing body of literature explores the causes, likelihood, size and composition of post-fire debris flows, but the channel response to these disturbances remains poorly studied. This study used repeat field surveys, time-lapse photographs and pre- and post-disturbance remote sensing datasets to document and analyse space- and time-varying channel response to post-fire debris flows along a steep mountain stream in the Upper Colorado River Basin, USA. Specifically, we evaluated channel morphology and bed composition changes, correlations between channel changes and valley and channel attributes, and the relative importance of spring snowmelt versus summer monsoon events. Several cross-sectional channel change types were observed from lidar a month after post-fire debris-flow events, including channelized and braided incision into deposits, incision into the pre-fire channel bed, bank erosion and no change. Channel changes were most correlated with pre-fire channel width, valley width and unit stream power, and these relationships could be tested in other burned locations to evaluate their transferability. Repeat channel surveys before and after snowmelt indicate rapid recovery and channel narrowing following major sediment disturbances, although sediment deposits remained in the channel margins. Together, these results highlight the importance of field and remote sensing-based channel surveys to improve understanding of, and potential to predict, mountain channel response to compounding climate disturbances.

Earth Surface Processes and Landforms↗
USGS2026

Mid-Cretaceous porphyritic magmatism in the Beartooth Mountains of the northern Laramide foreland and its connection to Frontier Formation sediment dispersal in Bighorn basin, Montana-Wyoming (USA), and Cordilleran tectonism

Porphyritic intrusions in the southeast Beartooth Mountains and porphyritic cobbles from the Torchlight conglomerate of the Cenomanian Frontier Formation in Bighorn basin (Montana-Wyoming, USA) bear striking resemblance. We utilize geologic mapping, petrography, geochemistry, and geochronology to characterize rocks from both locations and test whether the Beartooth Mountains area was the provenance for Bighorn basin cobbles. Mapping on Line Creek Plateau in the Beartooth Mountains shows three units with porphyritic phenocrysts consisting dominantly of plagioclase, plagioclase + potassium feldspar, and plagioclase + quartz. Petrography shows embayed quartz, opacitic rims, and zoned plagioclase. Geochemistry characterizes the intrusions as trachydacite to trachyte and trace elements exhibit enrichment in light rare earth elements and slight depletion in heavy rare earth elements relative to primitive mantle. Zircon U-Pb dating shows substantial Archean inheritance and ca. 104−92 Ma dates; some samples display multiple mid-Cretaceous date populations. Bighorn basin porphyritic cobbles exhibit similar modal mineralogy, near-identical geochemistry, and Archean inheritance and mid-Cretaceous geochronology. We suggest the porphyritic cobbles were sourced from basement-hosted porphyritic intrusions now exposed in the Beartooth Mountains, given that geochemistry does not agree with coeval volcanic sources along the Frontier Formation depositional fairway (Idaho batholith) or sedimentary-hosted Beartooth Mountains porphyritic sills intruded into Cambrian shale. Stratigraphic evidence does not support kilometer-scale mid-Cretaceous Beartooth Mountains exhumation to expose basement-hosted porphyritic intrusions, and we therefore present a volcanic hypothesis to bring rocks to the surface. This mid-Cretaceous multimillion-year magmatism is unexplained by tectonic models for the north-central Laramide foreland, and thus we further explore the emplacement processes, spatiotemporal significance, and tectonic implications of these rocks.

Montana, Wyoming↗
USGS2026

Site response models based on geometric parameters for southern California sedimentary basins

Site response in sedimentary basins is influenced by complex three-dimensional (3D) features, including trapping of seismic waves, focusing of seismic energy and basin resonance. Current ground motion models (GMMs) incorporate basin effects using one-dimensional parameters like V S30 and shear wave velocity isosurface depths, which are limited in capturing lateral and 3D effects. To address these limitations, we develop seismic site response models based on novel parameters that represent multi-dimensional properties of the Los Angeles Basin (LAB) geometry and shear wave velocity. We define a basin shape for the LAB using depth to subsurface geologic interfaces associated with the oldest sedimentary deposits (depth to a particular shear wave velocity horizon, i.e., 1.5 km/s - z 1.5 ) and the depth to the crystalline basement ( z cb ) which are determined using geologic cross sections and community seismic velocity model profiles. We explore a suite of geometric descriptors computed for the LAB and southern California, from which three parameters with the greatest predictive potential are selected and evaluated using empirical ground motion residual analyses in combination with the Boore et al. GMM. The results demonstrate that the zonal heterogeneity index ( ), standard deviation of the absolute difference between z 1.5 and z cb ( ) and standard deviation of z cb ( ) each provide a reduction in site-to-site variability ( ϕ S2S ) of empirical GMMs. The reduction in ϕ S2S is period-dependent, with average decreases of 3%, 26% and 6% for , , and , respectively. Although these reductions are modest from an engineering application perspective, they are statistically significant, underscoring the inherent difficulty in fully characterising complex basin effects. Collectively, these findings indicate that the inclusion of basin-specific geometric parameters yields measurable, albeit incremental, improvements in site response prediction and establishes a framework for the progressive refinement of seismic hazard characterisation within sedimentary basins.

California↗
USGS2025

Evolution of permeability and strength recovery of shear fracture under hydrothermal conditions

Geothermal energy is a clean and renewable resource that depends on the ability to move water through hot rock. In many locations, the ability to move water through rock requires the presence of extensive natural or human-made systems of fractures. However, these fracture systems are influenced by a variety of complex processes that occur at the temperature and pressure conditions found in geothermal reservoirs. To make geothermal systems more efficient, it is important to understand how these fractures evolve over time in response to these processes. This study explored how fractures in rock change under high-temperature and pressure conditions, like those found in geothermal reservoirs, through laboratory experiments and numerical simulations. In general fractures, tend to close due to pressure from the surrounding rock preventing the flow of water. One set of laboratory experiments revealed that fractures in granitic rock can weaken over time at elevated temperatures. These weak fractures can then slip slowly potentially keeping the fracture open for a longer period. This study also showed that at the highest temperature examined (250 °C) large differences between the chemical composition of the rock and the fluid moving through it could improve water flow temporarily likely due to the water dissolving parts of the rock. Flow rate predicted by models shows overall good agreement with the flow rate measured in our experiments but with some long time-scale variations that may be due to mechanical wear of the fracture surface. A second set of experiments found that, over the temperature range examined, the closure of fractures was caused by brittle failure driven by high stresses at points where fracture surfaces were in contact. This result suggests that fracture closure in geothermal reservoirs could be counteracted by maintaining high pore pressures. Additionally, we also tested a new method using electrical resistance to measure fracture closure in controlled systems with high precision.

Final Technical Report↗
USGS2025

Noble and base metal distribution and processes affecting ore tenors in the disrupted lower stratigraphy of the Stillwater Complex, USA

Exploration continues for contact-style Ni-Cu sulfide and chromitite-associated PGE mineralization in ultramafic rocks of the Stillwater Complex. At the Iron and Chrome Mountain areas, massive sulfides occur along the complex’s footwall contact and anomalous concentrations of PGE+Au are associated with the three lowermost chromitite seams. Southeast of Chrome Mountain, magmatic layering is highly disrupted by the presence of faults, magmatic breccias, serpentinized discordant dunites, pyroxenite pegmatoids, and disaggregated chromitite seams. The bulk rock chemistry, sulfide chemistry, and noble metal mineralogy of samples from this area were examined to determine the deportment of PGE and processes that led to enrichments in PGE, Au, Cu, Co and Ni. Results show that a sulfide liquid was the principal collector of PGE. If sulfide liquid was initially deposited with chromite, it was disaggregated or redistributed by subsequent melt or fluid infiltration, which may have resulted in the offset of peak PGE(+Cu, Ni) from peak Cr 2 O 3 concentrations, and upgraded PGE tenors. Upon cooling, PGE exsolved from sulfides to form discrete bismuth tellurides, arsenides, arsenic sulfides, antimonides, and alloys, commonly along the margins of sulfide globules. Calculated metal tenors are highest in the disseminated sulfides southwest of Chrome Mountain, whereas massive and net-textured sulfides near the Iron Mountain-Camp zone represent monosulfide solid solution cumulates. At progressively shallower levels, higher metal tenors combined with lower S/Se ratios are consistent with increasing R-factors from 100 to 100,000. Serpentinization and talc-tremolite alteration resulted in S loss through partial replacement of sulfides by secondary silicate+carbonate+magnetite+sulfide assemblages, further upgrading Ni-Cu-PGE tenors. The present work shows that processes responsible for the disruption of magmatic layering and post-magmatic fluid alteration along the intrusion’s lower contact led to noble and base metal enrichments.

Montana↗