SEARCH · USGS Science

Results for “Journal of Research of the U.S. Geological Survey”

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

SourceResearch connectionFrom yearThrough year

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

1,697 records · Page 62Linked to original sources

USGS2026

Insights into Mountain Pass carbonatite formation from in-situ sulfur isotopes and geochemistry of sulfate and sulfide minerals

The Mountain Pass carbonatite stock hosts a world-class rare earth element deposit and may be classified as a carbonate-sulfate igneous rock, as it contains on average > 50 volume percent carbonate minerals and 20 to 30 volume percent sulfate minerals. The sulfates range in composition from barite to celestine and locally occur with sparse sulfide minerals. We investigate the origin of sulfur enrichment and the occurrence of sulfur-bearing minerals in the Mountain Pass carbonatite with in-situ sulfur isotope and mineral chemistry. Barite cores with δ 34 S of 1 to 3‰ do not coexist with sulfides, whereas celestine rims with δ 34 S of > 3‰ are associated with sulfides with δ 34 S < -10‰. We propose a model in which sulfur-bearing sediments were subducted during episodes of plate convergence in the Mojave Province that preceded Mountain Pass magmatism. Metasomatism of the overlying mantle by melts derived from the subducted sediments generated an unusually carbon- and sulfur-rich source to yield carbonatite magmas. Sulfur from primary carbonatite magmas and ~ 1 to 7% sulfur from subducted sediment melts yielded a slightly enriched δ 34 S composition (relative to depleted mantle δ 34 S of -1‰) for early crystallizing barite. Celestine rims on magmatic barite cores formed at low, hydrothermal temperatures (< 350 °C) based on S isotope thermometry for equilibrium celestine-galena and celestine-pyrite pairs. The sparse sulfides in the carbonatite stock are not in equilibrium with the primary barite cores and therefore do not permit S isotope thermometry estimates of magmatic temperatures. The S/Se ratios of sulfide minerals (> 3,400) typically exceed primitive mantle values (S/Se of 3,340), also consistent with their derivation from hydrothermal fluids. Trace occurrences of sulfide and sulfate minerals in alkaline silicate stocks related to the carbonatite stock have similar δ 34 S compositions and yield similarly low formation temperatures, suggesting regionally extensive and chemically similar sulfur-bearing hydrothermal fluids that imparted lithologically diverse rocks with a consistent sulfur isotope fingerprint.

California, Nevada↗

USGS2026

Light-dependent activity of deepwater sculpin (Myoxocephalus thompsonii) across substrates with and without predation risk

Light availability strongly influences predator–prey interactions in deepwater ecosystems, where visual constraints shape both foraging success and prey behavior. The behavioral response of deepwater sculpin ( Myoxocephalus thompsonii ) to siscowet lake trout ( Salvelinus namaycush siscowet ) was studied under ecologically relevant light intensities spanning several orders of magnitude typical of daytime downwelling light in the 20–100 m depth range of Lake Superior. Trials were conducted over varying substrates (gravel, sand, and black fabric). Deepwater sculpin showed a significant preference for gravel over sand and black fabric. In the absence of siscowet, sculpin movement frequency increased as light intensity decreased. Sculpin reaction distance to siscowet was influenced by both light and substrate. Reaction distance was shortest at low light intensities, peaked at intermediate light intensities (3.05 × 10⁹ to < 6.0 × 10⁹ photons m⁻2 s⁻1), and declined again at the highest light intensities tested. In the presence of siscowet, sculpin activity was suppressed at the upper end of the light range (≥ 6.0 × 10⁹ photons m⁻2 s⁻1). The greatest increase in movement occurred between 6.0 × 10⁹ and 3.05 × 10⁹ photons m⁻2 s⁻1, which corresponded to the range where siscowet prey capture declined, suggesting sculpin exploit this low-light window to move with reduced risk. Reduced activity in the presence of predators is common among cryptic species, and our findings suggest that sculpin restrict movement at higher light levels to avoid detection by siscowet.

Wisconsin↗

USGS2024

Characterizing directivity in small (M 2.4-5) aftershocks of the Ridgecrest sequence

Directivity, or the focusing of energy along the direction of an earthquake rupture, is a common property of earthquakes of all sizes and can cause increased hazard due to azimuthally dependent ground‐motion amplification. For small earthquakes, the effects of directivity are generally less pronounced due to reduced rupture size, yet the directivity in small events can bias source property estimates and provide important insights into general regional faulting patterns. However, due to observational limitations, directivity is usually only measured and modeled for large events. As such, many studies of small earthquakes either ignore directivity altogether or assume a constant rupture direction for all events in a cluster. In our study, we apply a refined directivity fitting method constrained with two separate methods of source deconvolution to the dataset of aftershocks of the 2019 Ridgecrest earthquakes, which contain a large number of well‐recorded small‐to‐mid sized earthquakes occurring in close proximity to each other. The revealed directivity of 100+ small (M 2.4–5) earthquakes is highly heterogeneous and primarily oblique to and away from the main fault strike, suggesting a complex postseismic stress redistribution. In addition, the energy focusing effect of directivity appears to bias the selection of high‐quality data from stations in the direction of rupture, leading to average stress‐drop increases of 50% if directivity is not accounted for.

California↗

USGS2024

Probabilistic assessment of postfire debris-flow inundation in response to forecast rainfall

Communities downstream of burned steep lands face increases in debris-flow hazards due to fire effects on soil and vegetation. Rapid postfire hazard assessments have traditionally focused on quantifying spatial variations in debris-flow likelihood and volume in response to design rainstorms. However, a methodology that provides estimates of debris-flow inundation downstream of burned areas based on forecast rainfall would provide decision-makers with information that directly addresses the potential for downstream impacts. We introduce a framework that integrates a 24 h lead-time ensemble precipitation forecast with debris-flow likelihood, volume, and runout models to produce probabilistic maps of debris-flow inundation. We applied this framework to simulate debris-flow inundation associated with the 9 January 2018 debris-flow event in Montecito, California, USA. When the observed debris-flow volumes were used to drive the probabilistic forecast model, analysis of the simulated inundation probabilities demonstrates that the model is both reliable and sharp. In the fully predictive model, however, in which debris-flow likelihood and volume were computed from the atmospheric model ensemble's predictions of peak 15 min rainfall intensity, I 15 , the model generally under-forecasted the inundation area. The observed peak I 15 lies in the upper tail of the atmospheric model ensemble spread; thus a large fraction of ensemble members forecast lower I 15 than observed. Using these I 15 values as input to the inundation model resulted in lower-than-observed flow volumes which translated into under-forecasting of the inundation area. Even so, approximately 94 % of the observed inundated area was forecast to have an inundation probability greater than 1 %, demonstrating that the observed extent of inundation was generally captured within the range of outcomes predicted by the model. Sensitivity analyses indicate that debris-flow volume and two parameters associated with debris-flow mobility exert significant influence on inundation predictions, but reducing uncertainty in postfire debris-flow volume predictions will have the largest impact on reducing inundation outcome uncertainty. This study represents a first step toward a near-real-time hazard assessment product that includes probabilistic estimates of debris-flow inundation and provides guidance for future improvements to this and similar model frameworks by identifying key sources of uncertainty.

California↗

USGS2024

New 10Be-26Al isochron burial dating informs the Pliocene and Pleistocene evolution of the lower Colorado River, southwestern United States

Four new 10 Be- 26 Al isochron burial ages ranging from 4.4 to 2 Ma on ancestral Colorado River deposits in the lower Colorado River corridor (LCRC) help constrain the river’s evolution during the Pliocene and early Pleistocene. They help fill a gap between previous work that focused on older and younger deposits: Older dated deposits include the 5 Ma Bouse Formation, which records the integration of the Colorado River through a series of preexisting basins to the Gulf of California and the ca. 4.5 - 3.5 Ma Bullhead Alluvium, a 200 to 300 m thick aggregational package that immediately followed integration. The much younger, 100 - 70 ka, Chemehuevi Formation is another major aggradation package mapped throughout the LCRC. The new burial ages on the facies of Santa Fe Railway (4.37 ± 0.71 Ma), boulder conglomerate of Bat Cave Wash (2.12 ± 0.26 and 2.05 ± 0.31 Ma), and the Palo Verde alluvium (3.03 ± 0.26 Ma) partially fill in a 3.5 M.y. gap between the deposition of the Bullhead Alluvium and the Chemeheuvi Formation and document the timescales over which the Colorado River was able to remove the Bullhead aggradational package and initiate newer and smaller aggradational pulses.

Arizona, California, Nevada↗

USGS2025

Framework for assessing intranational risk of critical mineral sectors: Case study of India

Recent geopolitical actions regarding the supply of critical minerals have highlighted both the threats posed by countries dominating a critical mineral supply chain, and the potential mitigations served by alternative supplier countries. However, no systematic quantification exists for the intranational factors that determine such threats or mitigations. This study outlines a methodological framework for quantifying the intranational risk profile of critical mineral sectors and compares them with the sector's international potential for leverage. The four intranational factors of concentration, synchronization, private sector involvement, and regulatory obstacles were analyzed across multiple scales in order to determine an intranational risk score that quantifies the resilience or fragility of the critical mineral sector in question. Network analysis was also separately conducted to identify the most important entities operating across all critical mineral sectors. Finally, this study assessed the critical mineral sector's global share of reserves and production to determine an international risk score that quantifies the country's critical mineral sector's leverage potential. A case study of India showed that two critical mineral sectors with high leverage potential, chromite and barite, were respectively intranationally resilient and fragile. Of the remaining sectors with minimal leverage, five sectors such as aluminum were intranationally resilient and may serve as viable mitigation options for future critical mineral supply. Three sectors such as titanium were intranationally fragile, but recent regulatory developments may both decrease intranational risk and increase international risk, potentially making them major critical mineral sectors in the future.

Resources Policy↗

USGS2026

Orientation dependence of probabilistic seismic hazard estimates from CyberShake physics-based simulations

Earthquake ground‐motion intensities, such as pseudospectral accelerations (SAs), can vary significantly with horizontal orientation. However, conducting probabilistic seismic hazard analysis (PSHA) for each horizontal orientation is challenging because current ground‐motion models used in PSHA consider only a single horizontal intensity value, usually the median across all orientations, known as RotD50. To address this limitation, we employ physics‐based simulations for PSHA, which contain full waveforms from which ground‐motion intensities can be computed for all horizontal orientations to study directional seismic hazard. We apply our approach to the latest CyberShake study of the Greater Los Angeles metropolitan area, developed by the Statewide California Earthquake Center, finding that seismic hazard at a 2475‐yr return period, a common value used for earthquake‐resistant design, varies significantly with horizontal orientation. For instance, for SAs at 3 s, the maximum seismic hazard across all horizontal orientations is, on average, 15% higher than the median RotD50 hazard, with these differences becoming more pronounced at longer periods. These observed variations can generally be attributed to physical mechanisms that polarize seismic waves, such as the radiation pattern of the earthquake source and the influence of the subsurface structure. These results may have important implications for earthquake engineering applications, particularly for long‐period structures in areas with substantial horizontal variations in seismic hazard.

California↗

USGS2026

Towards a refined definition of “tidal influence” for the coastal wetland sciences

The astronomical tide relates to lunar (and to a lesser extent solar) influences on rising and falling water bodies, and is a key control of many coastal processes. However, tides manifest along coastal margins in different ways to muddle application beyond Sir Isaac Newton’s first descriptions. Here, we argue that tide, or specifically “tidal,” has not been adequately defined for the coastal wetland sciences to facilitate contemporary interdisciplinary distinction for rapidly changing coastlines and future climate adaptation planning. We provide a brief history of the tidal concept, and we describe tides in a diversity of uncommonly recognized marginal coastal environments. Common astronomical tide manifestations are described through hydrographic examples, then expanded to include those tidal environments exposed additionally to fluvial, meteorological (wind), or groundwater influences. Our expanded definition of “tidal” is intended to be a forward-looking construct as environments will inevitably change with sea-level rise, water extraction, and anthropogenic flow alteration. Our definition of wetland tidal influence, which we broaden, includes those “wetlands affected by astronomic water level fluctuations within coastally restricted environments causing frequent or infrequent surface inundation or groundwater tidal variation that affects the biogeochemistry of the soil (e.g., ion exchange capacity, oxygen state, mineralization processes)”. Astronomic water level fluctuations are often amplified or dampened by non-astronomical events, such as wind. Expanding the definition of what is considered tidal by the scientific community is critical for improved understanding of coastal wetland function and the suite of ecosystem services they provide people and society, including climate mitigation opportunity, management, and adaptation.

Wetlands↗

USGS2024

Mapping the resistivity structure of Walker Ridge 313 in the Gulf of Mexico using the marine CSEM method

A marine controlled source electromagnetic (CSEM) campaign was carried out in the Gulf of Mexico to further develop marine electromagnetic techniques in order to aid the detection and mapping of gas hydrate deposits. Marine CSEM methods are used to obtain an electrical resistivity structure of the subsurface which can indicate the type of substance filling the pore space, such as gas hydrates which are more resistive. Results from the Walker Ridge 313 study (WR 313) are presented in this paper and compared with the Gulf of Mexico Gas Hydrate Joint Industry Project II (JIP2) logging while drilling (LWD) results and available seismic data. The hydrate, known to exist within sheeted sand deposits, is mapped as a resistive region in the two dimensional (2D) CSEM inversion models. This is consistent with the JIP2 LWD resistivity results. CSEM inversions that use seismic horizons provide more realistic results compared to the unconstrained inversions by providing sharp boundaries and architectural control on the location of the resistive and conductive regions in the CSEM model. The seismic horizons include: 1) the base of the gas hydrate stability zone (BGHSZ), 2) the top of salt, and 3) the top and bottom of a fine grained marine mud interval with near vertical hydrate filled fractures, to constrain the CSEM inversion model. The top of salt provides improved location for brines, water saturated salt, and resistive salt. Inversions of the CSEM data map the occurrence of a ‘halo’ of conductive brines above salt. The use of the BGHSZ as a constraint on the inversion helps distinguish between free gas and gas hydrate as well as gas hydrate and water saturated sediments.

Louisiana↗

USGS2024

Application of non-stationary shear-wave velocity randomization approach to predict 1D seismic site response and its variability at two downhole array recordings

Accounting for uncertainties in seismic site response is crucial to improving the performance of one-dimensional (1D) ground response analyses (GRAs) at downhole array recording sites. In addition to site effects, uncertainties in 1D-GRAs can also be contributed from the seismic source and/or path. Though often representing not more than one percent of the distance (path) from the source, site conditions are known to have an enormous influence on ground shaking. In this study, we focus on the site shear-wave velocity ( V S ) structure, which is the main ingredient for estimating the variability of site response. As such, V S can manifest aleatory uncertainties related to the effects of small-scale spatial heterogeneities within the near surface, thus V S can substantially modify ground shaking during earthquakes. We apply a novel V S randomization approach to propagate the small-scale heterogeneities of V S to estimate seismic site response within a non-stationary probabilistic framework. The randomization approach generates samples of V S profiles that are used to perform several 1D-GRAs and obtain an averaged site response and related variability. The proposed method is implemented on data recorded at two downhole array sites with different subsurface soil conditions: a soft soil site on Treasure Island (California, United States of America) and a rock outcrop site in Cadarache (South-East France). We show that synthetic surface-to-borehole transfer functions from 1D-GRAs provide an acceptable fit to the empirical transfer functions from low-motion earthquake records and succeed in reproducing most of the site-specific seismic response variability. The remaining mismatch between transfer functions is likely due to insufficient precision on the seismic bedrock and the impedance contrast. The variability in site response is discussed with emphasis on the role of V S small-scale heterogeneities, attenuation, and input motion incidence angle in ground motion variability for the site and soil conditions at both locations.

California↗

USGS2025

Rare earth elements on the Moon

Rare earth elements (REEs) are a scarce but vital resource for our modern economies and lifestyles. Since the late 1990s, China has supplied the vast majority of the world’s refined REEs. Increasing global demand has broadened the search for REE deposits to unconventional places, including the Moon. Although most lunar rocks have very low REE concentrations, Apollo samples showed that one type of lunar rock containing potassium (K), REEs, and phosphorus (P)—known by the acronym KREEP—has high concentrations of REEs. Data from orbiting satellites have identified locations where substantial deposits of KREEP are likely. The viability of mining these deposits depends on the evolution of REE economics, the development of the Earth-Moon infrastructure, and the findings from future lunar mineral exploration missions.

Fact Sheet↗

USGS2024

Paleoseismology of the Sawtooth fault and implications for fault behavior in the epicentral region of the 2020 Mw 6.5 Stanley, Idaho, earthquake

The 2020 moment magnitude ( M w ⁠⁠ ) 6.5 Stanley, Idaho, earthquake raised questions about the history and extent of complex faulting in the northwestern Centennial Tectonic Belt (CTB) and its relation to the Sawtooth normal fault and Eocene Trans‐Challis fault system (TCFS). To explore faulting in this area, we excavated a paleoseismic trench across the Sawtooth fault along the western margin of the CTB, and compared an early Holocene (9.1 ± 2.1 ka, 1δ ⁠ ) rupture at the site with lacustrine paleoseismic data and fault mapping in the 2020 epicentral region. We find: (1) a history of partial to full rupture of the Sawtooth fault ( M w ⁠ 6.8–7.4), (2) that shorter ruptures ( M w ≤ ⁠⁠6.9 ) are likely along distributed and discontinuous faults in the epicentral region, (3) that this complex system that hosted the 2020 earthquake is not directly linked to the Sawtooth fault, (4) that the northeast‐trending TCFS likely plays a role in controlling fault length and rupture continuity for adjacent faults, and (5) that parts of the TCFS may facilitate displacement transfer between normal faults that accommodate crustal extension and rotation. Our results help unravel complex faulting in the CTB and imply that relict structures can help inform regional seismic hazard assessments.

Idaho↗

USGS2024

Predictions of groundwater PFAS occurrence at drinking water supply depths in the United States

Per- and polyfluoroalkyl substances (PFAS), known colloquially as “forever chemicals”, have been associated with adverse human health effects and have contaminated drinking water supplies across the United States owing to their long-term and widespread use. People in the United States may unknowingly be drinking water that contains PFAS because of a lack of systematic analysis, particularly in domestic water supplies. We present an extreme gradient boosting model for predicting the occurrence of PFAS in groundwater at the depths of drinking water supply for the conterminous United States. Our model results indicate that 71 to 95 million people in the conterminous United States potentially rely on groundwater with detectable concentrations of PFAS for their drinking-water supplies prior to any treatment.

conterminous United States↗

USGS2026

Diverging mineral chemistry of iron and nickel throughout Earth’s changing redox conditions reveals foundation for their evolution as protein cofactors

Iron (Fe) and nickel (Ni) were both foundational to early metabolism, yet their biological trajectories diverged as Earth’s surface redox state changed. Here, we integrate mineral chemistry network analysis, protein metal-site coordination-sphere analysis, and curated redox comparisons to test how geochemistry and metalloprotein architecture co-evolved. Mineral network analyses show broader electronegativity variation and network diversity for Fe-bearing minerals through time relative to Ni-bearing minerals. In structural analyses of protein metal centers in a combined Fe/Ni protein structure set, it is shown that Fe- and Ni-associated environments differ in amino-acid composition, hydropathy structure, and cysteine representation. The greater chemical diversity and electronegativity variation in Fe minerals mirror the higher redox and structural versatility of Fe-binding proteins. The presence of Fe in a broader range of mineral and protein environments demonstrates the chemical adaptability of the metal, from the anoxic Archean to oxidative Earth surface conditions following the Great Oxidation Event. Iron, with its broad redox potential range in Fe-oxidoreductases, has a central role in both anaerobic and aerobic metabolisms. Nickel, by contrast, is less widespread in biology. Today, Ni is predominantly employed in deeply branching anaerobic pathways and by proteins with narrower redox potential ranges. Our results show that evolutionary processes, constrained by metal chemistry, habitually utilize Fe as a redox generalist while retaining Ni in specialized roles. The divergent paths of Ni and Fe, from rocks to proteins, demonstrate the intimate relationship between planetary geochemistry and metabolic origins on Earth and suggest that Fe/Ni geochemistry may inform habitability assessments in extraterrestrial environments when interpreted within specific planetary environmental contexts.

Life↗

USGS2026

Temporal and spatial changes in seismic attenuation associated with inferred fluid migration in the 2016 central Apennines earthquake sequence

Prior work suggests that high‐frequency seismic attenuation acts as a highly sensitive proxy for crustal permeability and fluid mobility in fractured media. We test the hypothesis that the fault system responsible for the 2016–2017 Amatrice–Visso–Norcia–Capitignano sequence acted as an impermeable seal, compartmentalizing pressurized fluids until dynamic rupture triggered widespread fluid diffusion. By tracking across the sequence the spatiotemporal evolution of the S ‐wave anelastic attenuation parameter, we identify large, positive low‐frequency attenuation anomalies emerging within the hanging wall following the Amatrice mainshock and strictly preceding subsequent large ruptures. Conversely, we observe weaker, negative anomalies in the footwall, anticorrelated in time with those of the hanging wall, revealing a massive asymmetry in fluid redistribution and permeability evolution across the fault system. Furthermore, aftershock migration rates reveal distinct linear alignments in a distance‐reduced time space, allowing us to explicitly track and quantify episodes of lateral and upward fluid migration. These physically consistent patterns suggest that stress‐driven fluid diffusion directly weakens adjacent fault patches, dictating the spatiotemporal migration of seismicity. We conclude that near‐real‐time monitoring of seismic attenuation may help detect fluid redistribution in active fault systems and may provide useful information for time‐dependent seismic hazard assessment.

central Apennines↗

USGS2000

Natural attenuation of chlorinated volatile organic compounds in ground water at Area 6, Naval Air Station Whidbey Island, Washington

Natural attenuation is a viable alternative to pump and treat for meeting remediation objectives in the vicinity of the southern contaminant plume. The combination of historically low contaminant concentrations in ground water, a landfill cap that limits source area contributions, favorable conditions for degradation of VC, and a relatively long downgradient distance to potential receptors are all favorable for natural attenuation as a remediation alternative. Natural attenuation could effectively meet all but one remediation goal that extraction wells PW- 2, PW-4, PW-6, PW-7, PW-8 and PW-9 are currently being employed to meet. The goal of preventing migration of all VC across the site boundary could not be met by natural attenuation. Some VC would migrate south of the Navy boundary, but the potential for subsequent VC mineralization downgradient of the base and the existing institutional controls would result in minimal additional risk from using natural attenuation. In the western contaminant plume natural attenuation is not currently a viable alternative to pump and treat for meeting remediation objectives. There is a possibility that rates for reductive dechlorination of TCE and TCA could increase substantially if the plume was allowed to migrate beneath the Oak Harbor landfill, but there are not enough data to be certain of such an increase. TCA and TCE concentrations at the leading edge of the western plume need to be reduced to at least 25 and 40 µg/L, respectively, to be protective of potential downgradient receptors. Source area TCA and TCE concentrations have decreased substantially over the past ten years, and the extraction wells PW-3 and PW-5 in particular are removing a significant mass of contaminants from ground water, so natural attenuation may be a viable alternative for the western plume in the future. There would be some possible side benefits of using natural attenuation as an alternative to pump and treat in the southern contaminant plume. The first would be that the resulting decrease in the amount of treated water that would need to be recharged in the swale north of the landfill would result in less off-base migration of contamination across the western site boundary. The second benefit would be a substantial reduction in the amount of dissolved iron and manganese being extracted from the shallow aquifer and run through the treatment system. Removing that source of operation and maintenance problems would result in more effective containment and removal of contamination in the western contaminant plume. The most critical data gap identified in this evaluation is the paucity of contaminant chemistry information downgradient of the Navy boundary in the vicinity of the southern contaminant plume. Without such data, the behavior of the plume and the protectiveness of natural attenuation to downgradient receptors cannot be verified, and field attenuation rates for VC cannot be determined directly. The existing long-term monitoring plan would need to be reviewed and revised if natural attenuation is selected as a remedy for the contamination in the southern plume. In particular, additional performance monitoring wells may be required downgradient of the property.

Washington↗

USGS2024

Rainfall intensification amplifies exposure of American Southwest to conditions that trigger postfire debris flows

Short-duration, high-intensity rainfall can initiate deadly and destructive debris flows after wildfire. Methods to estimate the conditions that can trigger debris flows exist and guidance to determine how often those thresholds will be exceeded under the present climate are available. However, the limited spatiotemporal resolution of climate models has hampered efforts to characterize how rainfall intensification driven by global warming may affect debris-flow hazards. We use novel, dynamically downscaled (3.75-km), convection-permitting simulations of short-duration (15-min) rainfall to evaluate threshold exceedance for late 21st-century climate scenarios in the American Southwest. We observe significant increases in the frequency and magnitude of exceedances for regions dominated by cool- and warm-season rainfall. We also observe an increased frequency of exceedance in regions where postfire debris flows have not been documented, and communities are unaccustomed to the hazard. Our findings can inform planning efforts to increase resiliency to debris flows under a changing climate.

Arizona, California, Nevada, New Mexico, Utah↗

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↗