SEARCH · USGS Science

Results for “Transactions of the American Fisheries Society”

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,688 records · Page 47Linked to original sources

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↗

USGS2026

Suitability of ground-motion models for seismic hazard assessment in Puerto Rico and the U.S. Virgin Islands

We perform linear mixed-effects analyses with a ground-motion dataset to evaluate how well ground-motion models (GMMs) fit active crustal, subduction interface, and subduction intraslab earthquakes in Puerto Rico and the US Virgin Islands. Most of the GMMs reproduce the spatial variation in peak ground motions with earthquake magnitude and rupture distance but predict ground motions 0.3–1.0 natural log units (35%–270%) greater than observations. Two GMMs developed for Puerto Rico that are based on ground-motion records from mostly magnitude 4–5 earthquakes do not perform as well. We attribute the period-dependent overprediction in observed ground motions to differences between the observed site response and the linear site response in the GMMs. Consequently, we developed region-specific GMMs by adjusting the period-dependent linear site response coefficients and period-dependent constant coefficients to remove most of the bias between predicted ground motions and observations. For the analysis, we compile ground-motion records and process waveforms to build a dataset with 10,127 records at 72 stations from 849 magnitude 4.0–6.4 earthquakes between 1 January 2006 and 31 March 2024. The earthquakes include active crustal, subduction interface, and subduction intraslab events. We evaluate the GMMs using the time-averaged shear wave speed in the top 30 m ( ), which we compile from site surveys and proxy values computed from horizontal to vertical spectral ratios. Site terms exhibit strong consistency across GMMs and crustal and subduction earthquakes, indicating that the linear mixed-effects analysis successfully isolates the effects of local site response. The event terms show little spatial correlation and more substantial variability than in other regions, which we attribute to uncertainties in the earthquake magnitudes. This analysis guides the selection of the GMMs for the 2025 update of the National Seismic Hazard Model for Puerto Rico and the US Virgin Islands.

Puerto Rico, U.S. Virgin Islands↗

USGS2025

Avian navigation: Comparing the olfactory navigational “map” and the infrasound direction-finding hypotheses to aeronautics

Animal navigation has long been a fascinating but bewildering subject. Humans and animals might well share similar navigational strategies because they developed within the same physical environments. A “map-and-compass” model has been proposed to explain the two-step avian navigational process, but the “map” step has remained elusive. Although scalar values from bicoordinate geomagnetic or atmospheric olfactory gradients have been considered foundational to the avian map, neither has proved convincing engendering decades of controversy. The olfactory map, and an alternative infrasound direction-finding (IDF) hypothesis, are discussed in this review. The olfactory map hypothesis currently requires extensive stable gradients of trace-odor ratios, but such gradients are highly unlikely within a turbulent and rapidly mixed lower atmosphere. The IDF hypothesis, on the other hand, postulates a two-step navigational model analogous to the maritime and aeronautical radio direction-finding technique. This review was also written to encourage further investigation, and direct testing, of the acoustic navigational process. The IDF hypothesis, at present, appears the better explanation of observed avian navigational behavior and accuracy within the atmosphere’s physical environment.

Journal of Comparative Physiology A↗

USGS2024

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

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

Alaska↗

USGS2024

Investigation of land cover within wetland complexes at Dixie Meadows, Churchill County, Nevada, from October 2015 to January 2022

The U.S. Geological Survey investigated land cover at subannual time steps within six wetland areas in Dixie Valley, Churchill County, Nevada, from October 2015 to January 2022. As requested by the U.S. Fish and Wildlife Service, we used aerial photography and satellite remote sensing data to map surface water and other land cover types within the wetland complexes. We identified five land cover classes using the green normalized difference vegetation index (gNDVI) and its inverse relationship to the normalized difference water index (NDWI) within three U.S. Department of Agriculture National Agriculture Imagery Program aerial images (acquired in 2015, 2017, and 2019) and 110 European Space Agency Sentinel-2 satellite images (acquired 2015–2022). The relative wetness of soil conditions within each land cover class is estimated by comparison to previously published observations of relative conductivity measured by 79 field-based sensors within the wetlands from 2019 to 2021. We mapped the areal coverage of the five land cover classes for approximately 385 acres (1,559,000 square meters [m²]) comprising six individual wetland complexes as well as a larger 1,298- acre (5,254,000-m 2 ) area of interest inclusive of the wetland complexes and adjacent landscape. Land cover of open water (Class 5) primarily within ponds at one of the wetland complexes comprised 8,333 m 2 , on average, of the wetland complexes. Land cover of mixed shallow surface water, saturated soil, and vegetation (Class 4) comprised 111,723 m 2 on average of the wetland complexes. Land cover of dense green vegetation canopy cover (Class 3) that often (46 percent of observations) had underlying surface water or saturated soil conditions comprised 592,522 m 2 on average of the wetland complexes. The remaining areas of the wetland complexes not mapped as these three land cover types (Classes 2 and 1) had sparse vegetation or bare soil cover and commonly (greater than or equal to 67 percent of observations) had dry soil conditions. The investigation of land cover detailed in this report could inform future efforts to map land cover more precisely via higher resolution remote sensing or ground-based surveying or could be incorporated with other environmental monitoring data to characterize habitat and hydrology of the wetland complexes at Dixie Meadows.

Nevada↗

USGS2025

Performance analysis of oil recovery and CO2 retention in a greenfield residual oil zone: CO2-EOR in Tall Cotton Field (Permian Basin, West Texas, USA)

Residual oil zones (ROZs) can offer significant oil resources via enhanced oil recovery (EOR) as well as subsurface carbon dioxide (CO 2 ) retention during injection. If injected CO 2 is anthropogenic, the ROZs can offer a substantial geologic storage potential. The ROZs below the oil/water contact (OWC) of main pay zones (MPZ) in conventional reservoirs or brownfields, are more commonly developed for CO 2 injection and oil production and reported in the literature. However, CO 2 -EOR in greenfield ROZs, reservoirs without a MPZ present, have rarely been developed for CO 2 -EOR operation. The Tall Cotton Field of West Texas, Permian Basin, which started production in 2015 (Phase 1) and expanded in 2017 (Phase 2) from the San Andres Limestone, is one of the first examples of greenfield ROZs developed for EOR by injecting CO 2 . This paper analyses EOR and CO 2 retention performance of Tall Cotton Field using allocated injection and production data from inverted 5-spot well patterns of Phase-1 and -2 developments. Production and injection data allocated to each of the 28 identified patterns (nine 20-acre patterns for Phase-1, three 20-acre and sixteen 10-acre patterns for Phase-2) were analyzed for historical and forecasted oil recovery using ratio-trend decline analysis, and for CO 2 retention performance of the patterns. The allocated data were further used to calculate injected reservoir pore volume and void replacement ratios (VRR) for the analysis period. Quantitative results indicated that oil recovery factors of the 5-spot patterns varied between 4–10 %, and 5–30 % between the end of injection and the forecast periods, respectively. Storage of CO 2 , on the other hand, increased to a mean value of ∼7130 MMscf per pattern in Phase-1 and to a mean storage of 3700 MMscf per pattern in Phase-2 until the end of injection, followed by a decline after the end of injection and into the forecast period. Resulting CO 2 utilization factors ∼6–50 Mscf/bbl were estimated at the end of injection. Overall, presented results suggested that developing greenfield ROZs for CO 2 -EOR can be as promising as brownfield ROZs and mature MPZs for EOR and underground storage of injected CO 2 . For Tall Cotton Field, results suggest that Phase-2 patterns generally outperformed Phase-1 for oil recovery factors, while Phase-1 performed better in CO 2 retention performance metrics. This is the first study in the literature that reports a detailed CO 2 -EOR performance analysis of a greenfield ROZ in the Permian Basin, which can potentially allow for comparison with MPZs and brownfield ROZs.

Texas↗

USGS2026

Decadal trends in the quality of groundwater used for public drinking-water supply in California, 2004–2023, California groundwater ambient monitoring and assessment program, priority basin project

This study provides a comprehensive assessment of decadal changes in the quality of groundwater used for public drinking-water supply at 444 monitoring sites across California during 2004–2023. We assessed decadal step trends in groundwater quality for 145 water-quality constituents and geochemical indicators statewide and across geographic and land-use based network groups. We evaluated the statistical significance of directional changes (predominant increase or decrease of constituent concentrations) and the magnitude of those changes across all network groups. Uranium showed the most widespread directional and high-magnitude increases of all constituents with regulatory benchmarks statewide, particularly in the agriculture-dominated Central Valley as well as urban- and desert-dominated regions of Southern California. Fluoride and perchlorate showed the most widespread directional and high-magnitude decreases of all constituents with regulatory benchmarks statewide, which were also most pronounced in Southern California. Although arsenic and nitrate did not often register significant directional changes across network groups, they showed widespread, high-magnitude changes in both directions (increase and decrease) at levels often exceeding 10 percent of respective regulatory benchmarks statewide. Triazine herbicides (atrazine and simazine) and the gasoline oxygenate methyl tert-butyl ether (MTBE) showed significant directional decreases statewide, but not at levels considered to be of high magnitude compared to respective regulatory benchmarks. We observed significant directional and high-magnitude increases of total dissolved solids (TDS) statewide, which were most pronounced in agricultural areas. Analysis of explanatory geochemical indicators indicated that prevalent statewide increases of alkalinity and calcium were the predominant components of the observed statewide increases in TDS by mass. Widespread increases in groundwater alkalinity and calcium across agricultural and urban areas may be related, in part, to warm-season irrigation and other anthropogenic factors that have shifted soil weathering dynamics over the long term. Increasing alkalinity concentrations were related to increasing uranium concentrations, particularly in areas with aquifer materials derived from granitic rocks. Conversely, increasing calcium concentrations were related to decreasing fluoride concentrations, particularly in areas where fluoride occurred naturally at elevated concentrations. Decrease of perchlorate, triazine herbicides, and MTBE are likely related to decreased anthropogenic source inputs over time and natural attenuation in aquifers.

EarthArXiv↗

USGS2024

Assessing locations susceptible to shallow landslide initiation during prolonged intense rainfall in the Lares, Utuado, and Naranjito municipalities of Puerto Rico

Hurricane Maria induced about 70 000 landslides throughout Puerto Rico, USA, including thousands each in three municipalities situated in Puerto Rico's rugged Cordillera Central range. By combining a nonlinear soil-depth model, presumed wettest-case pore pressures, and quasi-three-dimensional (3D) slope-stability analysis, we developed a landslide susceptibility map that has very good performance and continuous susceptibility zones having smooth, buffered boundaries. Our landslide susceptibility map enables assessment of potential ground-failure locations and their use as landslide sources in a companion assessment of inundation and debris-flow runout. The quasi-3D factor of safety, F 3 , showed strong inverse correlation to landslide density (high density at low F 3 ). Area under the curve (AUC) of true positive rate (TPR) versus false positive rate (FPR) indicated success of F 3 in identifying head-scarp points (AUC = 0.84) and source-area polygons (0.85 ≤ AUC ≤ 0.88). The susceptibility zones enclose specific percentages of observed landslides. Thus, zone boundaries use successive F 3 levels for increasing TPR of landslide head-scarp points, with zones bounded by F 3 at TPR = 0.75, very high; F 3 at TPR = 0.90, high; and the remainder moderate to low. The very high susceptibility zone, with 118 landslides km −2 , covered 23 % of the three municipalities. The high zone (51 landslides km −2 ) covered another 10 %.

Puerto Rice↗

USGS2021

Geology, mineralogy, and cassiterite geochronology of the Ayawilca Zn-Pb-Ag-In-Sn-Cu deposit, Pasco, Peru

The Ayawilca deposit in Pasco, Peru, represents the most significant recent base-metal discovery in the central Andes and one of the largest undeveloped In resources globally. As of 2018, it hosts an 11.7 Mt indicated resource grading 6.9% Zn, 0.16% Pb, 15 g/t Ag, and 84 g/t In, an additional 45.0 Mt inferred resource grading 5.6% Zn, 0.23% Pb, 17 g/t Ag, and 67 g/t In, and a separate Sn-Cu-Ag inferred resource of 14.5 Mt grading 0.63% Sn, 0.21% Cu, and 18 g/t Ag. Newly obtained U–Pb dates for cassiterite by LA-ICP-MS (22.77 ± 0.41 and 23.05 ± 2.06 Ma) assign the Ayawilca deposit to the Miocene polymetallic belt of central Peru. The polymetallic mineralization occurs as up to 70-m-thick mantos hosted by carbonate rocks of the Late Triassic to Early Jurassic Pucará Group, and subordinately, as steeply dipping veins hosted by rocks of the Pucará Group and overlying Cretaceous sandstones-siltstones of the Goyllarisquizga Group. Relicts of a distal retrograde magnesian skarn and cassiterite (stage pre-A) were identified in the deepest mantos. The volumetrically most important mineralization at Ayawilca comprises a low-sulfidation assemblage (stage A) with quartz, pyrrhotite, arsenopyrite, chalcopyrite, Fe-rich sphalerite, and traces of stannite and herzenbergite. Stage A sphalerite records progressive Fe depletion, from 33 to 10 mol% FeS, which is compatible with the observed transition from low- to a subsequent intermediate-sulfidation stage (B) marked by the crystallization of abundant pyrite and marcasite. Finally, during a later intermediate-sulfidation stage (C) sphalerite (up to 11 mol% FeS), galena, native bismuth, Cu-Pb-Ag sulfosalts, siderite, Mn-Fe carbonates, kaolinite, dickite, and sericite were deposited. This paragenetic evolution shows striking similarities with that at the Cerro de Pasco Cordilleran-type polymetallic deposit, even if at Ayawilca stage C did not reach high-sulfidation conditions. The occurrence of an early retrograde skarn assemblage suggests that the manto bodies at Ayawilca formed at the transition between distal skarn and skarn-free (Cordilleran-type) carbonate-replacement mineralization. Mineral assemblages define a T- f S 2 evolutionary path close to the pyrrhotite-pyrite boundary. Buffering of hydrothermal fluids by underlying Devonian carbonaceous phyllites of the Excelsior Group imposed highly reduced conditions during stage A mineralization (log f O 2 < − 30 atm). The low f O 2 favored efficient Sn mobility during stages pre-A and A, in contrast to other known ore deposits in the polymetallic belt of central Peru, in which the occurrence of Sn minerals is minor. Subsequent cooling, progressive sealing of vein walls, and decreasing buffering potential of the host rocks promoted the shift from low- (stage A) to intermediate-sulfidation (stages B and C) states. LA-ICP-MS analyses reveal significant In contents in Fe-rich sphalerite (up to 1.7 wt%), stannite (up to 1908 ppm), and chalcopyrite (up to 1185 ppm). The highest In content was found in stage A sphalerite that precipitated along with chalcopyrite and stannite, thus pointing to the early, low-sulfidation assemblage as prospective for this high-tech metal in similar mineral systems. Indium was likely incorporated into the sphalerite crystal lattice via Cu + + In 3+ ↔ 2 Zn 2+ and (Sn, Ge) 4+ + (Ga, In) 3+ + (Cu + Ag) + ↔ 4 Zn 2+ coupled substitutions. Indium incorporation mechanisms into the stannite and chalcopyrite crystal lattices remain unclear.

Pasco↗

USGS2026

Forecasting volcanic activity in Germany—A multi-criteria approach

Igneous activity, including shallow intrusions and volcanism, has the potential to disrupt underground critical infrastructure. Notably, future underground infrastructure projects like high-level radioactive waste repositories must be sited in areas of extremely low disruption probability by igneous activity. In Germany, according to the Repository Site Selection Act of 2017 (Standortauswahlgesetz, or StandAG), areas in which Quaternary volcanism is either present or future volcanic activity is expected within the next 1 million years (m.y.) must be excluded from the site selection process. Although the locations of regions with Quaternary volcanism are reasonably well known in Germany, forecasting potential igneous activity at intraplate volcanic fields is challenging, as many processes and their interactions control the spatial distribution of volcanic centers. Here, a semi-quantitative, multi-criteria approach is proposed for a regional evaluation of the relative potential of future igneous activity in Germany. A variety of geoscientific indicators are used, including seismic anomalies in Earth’s mantle, gravity data, tectonic activity, sutures, ground motion, earthquakes, mantle degassing centers, and geochronological data of volcanic rocks. The indicators describe the sequence of processes from potential melt generation in Earth’s mantle, through ascent and accumulation of melt within the lithosphere, to eruption at Earth’s surface. In total, 15 out of 30 proposed geoscientific indicators are selected and quantified using 20 total assigned parameters. Defined threshold values are used to spatially delimit relevant parameter properties to focus on areas with higher potential of future magmatic activity. To consider uncertainties of parameters and their underlying processes, which are usually more spatially extensive below ground, buffer zones are defined in which values of relevance decrease with increasing distance from the initial lateral shape of a parameter. Normalized parameters are combined into an index, whose spatial value distribution is used to differentiate the relative potential of future igneous activity (within the next 1 m.y.). The sensitivity of the results is shown by varying the weighting factors for the relevant parameters in country-wide index maps. Thereby, profiles illustrate the distribution of the resulting index values and respective index fractions of various parameters. Different index maps for the relative potential of future igneous activity are presented and can be used for hazard assessments.

Professional Paper↗

USGS2024

Shallow storage of the explosive Earthquake Flat Pyroclastics magma body, Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand: Evidence from phase-equilibria experiments

Rhyolitic tuffs range widely in their crystal contents from nearly aphyric to crystal-rich, and their crystal cargoes inform concepts of upper crustal magma reservoirs. The Earthquake Flat pyroclastics (Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand) are 10 km 3 of rhyolitic tuffs with abundant (~ 40 vol.%) plagioclase and quartz, minor biotite, hornblende, and orthopyroxene, and accessory Fe-Ti oxides, apatite, and zircon, set in high-silica rhyolitic glass. Major minerals form large, euhedral phenocrysts and abundant glomerocrysts with few disequilibrium textures excepting some faintly resorbed quartz. Plagioclase phenocrysts have thick rims of nearly constant composition near An 30 , and hornblende is weakly zoned or unzoned. The abundant and texturally complex mineral assemblage contrasts with the nearby (~ 25 km), nearly synchronous, but more voluminous and crystal-moderate rhyolite tuffs from Rotoiti caldera. New H 2 O-saturated phase-equilibria results on the erupted Earthquake Flat melt (glass) determine its co-saturation with the partial phenocryst assemblage of plagioclase, quartz, biotite, and Fe-Ti oxides at: 140 MPa, 755 ºC. These closely approximate the conditions of the pre-eruptive magma body assuming it was saturated with nearly pure H 2 O and at an f O 2 of ~ Ni–NiO. Absence of hornblende and orthopyroxene from the synthesized assemblages may result from those minerals being in a peritectic reaction relation with melt to produce biotite, so they would not grow from the liquid used as starting material. Experimental results on Rotoiti rhyolite (Nicholls et al. 1992) show that the two bodies resided at similar pressures, temperatures, and f O 2 s. Lower crystal abundance of the Rotoiti tuffs may result from slight compositional differences. We interpret that the Earthquake Flat pyroclastics were sourced from the crystal-rich periphery of a mushy reservoir system with the Rotoiti occupying a more melt-rich central location. Uncertain is whether this was a single intrusion zoned continuously in crystallinity, or discrete adjacent intrusions, but our results illustrate and quantify complexities of magma storage across relatively short distances.

North Island, Taupo Volcanic Zone↗

USGS1990

Floods of February 1989 in Tennessee

Rainfall amounts of over 5 inches the night of February 13 and the morning of February 14, 1989, caused flooding in areas of Middle and West Tennessee. The towns of Lebanon in Middle Tennessee and Obion in West Tennessee were most severely affected. Most of the business district in Lebanon and many residential areas in Obion were flooded. Recurrence intervals for 24-hour rainfall totals were as high as 25 years at some sites but most peak discharges had recurrence intervals of less than 10 years. Rainfall amounts for the period February 13-20, 1989, peak stages and discharges for this flood, the peak of record, and a list of discharge measurements made during the flood are documented in this report.

Tennessee↗

USGS2025

Tracing fluid evolution during batholith and deposit formation: In-situ and bulk rock Mo isotopes from the Questa porphyry Mo deposit, New Mexico, USA

The Oligocene Questa porphyry Mo deposit, located in New Mexico, USA, is a well-characterized system that formed over 400 kyr, punctuating the 8 Myr formation of the Questa-Latir batholith. Detailed sampling of the intrusions and minerals associated with magmatichydrothermal mineralization, as well as bookending barren intrusions, allows for examination of the Mo isotope record of the deposit and its host batholith. Through this established chronologic framework, such data can yield insights into the metallogenesis of deposit formation and long-term magma genesis evolution for upper crustal magmatism. The Mo isotopic compositions of barren intrusions yield a younging trend towards lighter δ 98 Mo. However, intrusions associated with Mo deposit formation yield heavier δ 98 Mo relative to barren magmatism, by up to 2.0‰. This indicates that they were affected by an isotopically heavy and Mo-rich fluid. The in-situ data from molybdenite covers the range observed in mineralized bulk rocks. However, individual samples consistently yield Mo isotopic variations of 0.4‰ to 0.5‰. The consistency of the δ 98 Mo variation in in-situ analyses of molybdenite within hand samples and individual mineral grains is likely the result of MoS 2 crystallization in a closed system and reflects an almost complete precipitation of Mo.

New Mexico↗

USGS2025

Estimating the probability of export restrictions to inform mineral criticality

As demand for advanced technologies rises, mineral commodities will increase in geopolitical importance. To assess risks associated with mineral commodity supply chain disruptions, governmental agencies and others have developed "criticality" assessments, with criticality described using the economic impact and probability of supply chain disruptions. In previous work, subjective supply risk indicators were developed to approximate this probability, typically combining several factors such as supply diversity and political stability of trading partners, where indicator weightings can substantially impact results. This work explicitly quantifies trade barrier probability using an ensemble of several machine learning classifiers, with probability estimates informed by exogenous variables such as prior trade barrier implementation and global export dominance. Major differences in the 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. This approach enables a direct, quantitative, objective approach to assessing trade barrier probability, enhancing risk identification and prioritization for policymakers.

SSRN↗

USGS2025

Petrogenesis and mineralization potential of spinifex komatiitic basalts in the Bradley Peak greenstone terrane, Wyoming Province

Komatiitic volcanic rocks are important hosts of Ni sulfide mineralization and record early Earth evolution; however, those in the well-studied Archean Wyoming Province have received little attention. Here, we elucidate the timing and petrogenesis of the Bradley Peak komatiitic volcanic rocks using field and textural observations, geochronology, and geochemistry. Detrital and igneous zircon U-Pb ages for two samples from previously undated units support published age determinations, placing the eruption age at 2.72 Ga. Stratigraphy of the volcanic flows was mapped and 36 samples including cumulates, greenschists, and spinifex-textured rocks were collected. Whole-rock geochemistry was used to classify the spinifex-textured samples as Al-undepleted komatiitic basalts (11–17 wt% MgO). Platinum-group element concentrations (n = 25) are like those in global Al-undepleted komatiitic basalts, and PGE/Ti ratios do not indicate the volcanic flows likely host sulfide mineralization. Initial ε Nd values of −0.5 to +4.7 (n = 16), indicate that these lavas were derived from a depleted mantle source and have negligible evolved crust contamination. The primary magma to the komatiitic basalt flows is estimated to have had 19 wt% MgO and be derived from ∼15 to 25 % mantle partial melting at 3–4 GPa. Trace element chemistry and thermodynamic modeling suggest the primary melt assimilated local banded iron formation. Although the Bradley Peak komatiitic basalts do not contain positive evidence of magmatic sulfide deposits, depleted Au in the flows suggests they could be source rocks for nearby orogenic gold deposits.

Wyoming↗

USGS2025

Airborne geophysics for geologic mapping of critical mineral systems in the United States southern midcontinent

The increased demand for clean energy technology and a significant reliance on foreign supply chains have given impetus to understanding critical mineral systems and locating potential resources within the United States. At least thirteen critical mineral-bearing systems have been identified throughout the U.S. southern Midcontinent (Hofstra and Kreiner, 2020) but much of the region’s geologic framework is concealed by vegetation and sedimentary cover that hinder traditional geologic mapping efforts. Airborne geophysical data provide an effective way to overcome these obstacles and to provide additional insight into the deeper structures that underlie shallow mineralization. However, legacy airborne magnetic and radiometric data were collected using now-outdated instruments and methods, inconsistent survey parameters, and large flight-line spacings resulting in low-resolution data that present challenges to regional-scale study and interpretation. Over the last decade, the U.S. Geological Survey Earth Mapping Resources Initiative (EMRI) and National Cooperative Geologic Mapping Program have conducted a series of high-resolution airborne magnetic and radiometric surveys across the southern Midcontinent (Fig. 1) as part of an effort to improve understanding of the geophysical framework and natural resource potential in the region. These surveys are designed using modern survey methods and instruments with consistent parameters for flight-line spacing and flight height relative to magnetic sources. The EMRI airborne surveys are planned in collaboration with State geological surveys based on focus areas (Dicken et al., 2022) according to the presence of or potential for critical mineral deposits. High-resolution airborne magnetic and radiometric data cover focus areas such as the southeast Missouri iron metallogenic province and South-Central iron-oxide-apatite (IOA) – iron-oxide-copper-gold (IOCG) province, the Magnet Cove alkaline-carbonatite complex, the Midwest Permian ultramafic dike district, the Illinois-Kentucky fluorspar district, and several Mississippi Valley-type lead-zinc deposits and districts (Fig. 1). These focus areas represent known deposits or prospective host systems of critical minerals including rare earth elements (REEs), platinum-group elements (PGEs), cobalt, lithium, fluorspar, niobium, titanium, vanadium, lead, zinc, gallium, germanium, and many more. Other significant geologic and geophysical features covered include the Reelfoot rift, the New Madrid seismic zone, the Illinois basin, the Arkoma basin, the South-Central magnetic lineament, and the Kentucky-Tennessee magnetic anomaly (Fig. 1). This presentation focuses on new airborne magnetic and radiometric data with continuous coverage across parts of six states, preliminary interpretations, examples of geologic mapping applications, and discussion of newly discovered magnetic anomalies and follow-up investigations.

Conference Paper↗

USGS2025

Systematic approach to prioritize wells for effective groundwater monitoring and management in the Arkansas Headwaters Basin, Colorado, USA

Study region The Arkansas Headwaters Basin, an intermountain basin in the Southern Rocky Mountains of North America. Study focus Our specific focus is choosing a set of wells to support a possible future regional groundwater-surface water model that would support water management. We present a three-step process using multiple criteria to score, predict, and choose prioritized wells that capture the full distribution of data including extremes. The three-step process provides accessible visualizations, fiscally efficient well prioritization, and screening useful for subsequent groundwater modeling. The novelty of the proposed methodology is the systematic approach integrating a scoring and a predictive approach to support a selection path. The systematic approach may be broadly adapted for other basins. New hydrological insights for the region Understanding regional hydrology hinges on efficient collection of hydrologic data that captures the relevant dynamics including extremes. The present study, a case study for a particular basin in the Southern Rocky Mountains, is the first use of a scripted (R software) strategy to select an economical and representative set of monitoring wells. Our findings suggest caution when using proximity as a proxy for correlation, because proximal wells in the same geologic formation and similar depths are not always correlated. In the Arkansas Headwaters Basin, subsurface geology may be less influential on groundwater elevations than broader hydrologic influences, such as regional drought.

Colorado↗

USGS2025

Uncertainty reduction for subaerial landslide-tsunami hazards

Subaerial rock slopes may generate a tsunami by rapidly moving into the water. Large uncertainty in landslide characteristics propagates into large uncertainty in tsunami hazard, making hazard assessment more difficult for land and emergency managers. Once a potentially tsunamigenic landslide is identified, it may not be clear which landslide characteristics contribute most significantly to uncertainty in the tsunami hazard. Our aim is to document the relative worth of different landslide characteristics (e.g., size, material properties) for reducing uncertainty in landslide-tsunami hazard assessments. Isolating the relative importance of specific landslide characteristics may inform prioritization of data collection and improve efficiency in understanding hazard. To accomplish this, we generated a set of 288 landslide-tsunami simulations in which we systematically varied the size and material properties of possible failure extents at the Barry Arm landslide complex in northwestern Prince William Sound, Alaska, USA. We find that for landslides smaller than the receiving waterbody, the landslide volume has the strongest effect on resulting wave characteristics and thus the highest leverage on reducing uncertainty in tsunami hazard. In contrast, for landslides substantially larger than the waterbody, the duration of rapid movement of the landslide has the highest leverage. Based on our results, we propose a classification scheme for subaerial landslides based on the relative size of the landslide and waterbody. Additionally, our results support the generation of a tsunami height transfer function between existing tide gages and a nearby coastal city. These results have direct implications for the practice of operational early warning.

Alaska↗