USGS Science⌕ Search

SEARCH · USGS Science

Results for “The Planetary Science Journal”

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

Reduction of red bed sedimentary rocks in connection with energy metal ore formation: A case study from the Sinbad seep, Mesa County, Colorado

The Paradox Basin’s Sinbad seep is a modern analog for ancient bleaching of red bed sediments by introduced alkaline, reducing brines. This bleaching, involving reductive alteration of former red beds, is essential ground preparation that enables the altered rocks to trap Cu, U, and V from later oxidized fluids, forming ore deposits. Study of Sinbad thus offers insights into these metallic mineralization processes in the Paradox and other sedimentary basins. The Sinbad seep occurs where shallow groundwater interacts with organic-rich petroleum source rocks, flows up a fault, and discharges into Salt Creek. Ratios of Na/Cl, Na/Br, and Cl/Br, plus high sulfate concentrations, indicate that the salinity of the seep water originated from dissolution of halite and gypsum during topographically driven flow of meteoric water across a diapir of the Pennsylvanian Paradox Formation, rather than from deep basinal brines. In contact with the organic-rich shales of the Paradox, these SO 4 -rich waters become reduced through bacterial reduction, producing H 2 S. The waters react with red beds of the Permian Cutler Formation, causing pervasive bleaching. The bleaching at Sinbad is characterized by iron-conservative reduction of ferric iron in diagenetic hematite and detrital ilmenite and magnetite to Fe sulfides. The ferrous iron is retained in these sulfides likely as sorbed Fe. Associated alteration includes precipitation of quartz and feldspar overgrowths, partial dissolution of detrital quartz and feldspar caused by pressure solution, formation of clay and authigenic rutile, and precipitation of carbonate and gypsum. The Fe sulfides rapidly degenerate to a mixture of jarosite and iron oxides on weathering at the surface. Among the major-oxide elements and most trace elements, there is no statistical difference between unbleached and bleached samples. The exceptions are uranium and sulfur, which are somewhat greater in bleached samples. Fission track radiography illuminates that uranium is preferentially concentrated in iron oxide cements, iron sorbed into detrital clasts, and finely disseminated iron oxide in illite cement in bleached samples. A leaching experiment suggests that uranium may be more easily available for mobilization from rocks that have undergone bleaching alteration, making them potential U sources. In addition, bleached rocks form effective traps for U, V, and some Cu mineralization.

Colorado↗
USGS2024

Introduction to recommended capabilities and instrumentation for volcano monitoring in the United States

Introduction The National Volcano Early Warning System (NVEWS) was authorized and partially funded by the U.S. Government in 2019. In response, the U.S. Geological Survey (USGS) Volcano Hazards Program asked its scientists to reflect on and summarize their views of best practices for volcano monitoring. The goal was to review and update the recommendations of a previous report (Moran and others, 2008) and to provide a more detailed analysis of capabilities and instrumentation for monitoring networks for U.S. volcanoes. This Scientific Investigations Report and its chapters reflect those USGS scientists’ views and summaries and will serve as a guide for future network upgrades funded through NVEWS. Given the well-documented hazards posed by volcanoes to population centers and aviation (for example, Blong, 1984; Scott, 1989; Neal and others, 1997, 2019; Guffanti and others, 2010; Shroder and Papale, 2014; Prata and Rose, 2015; Palmer, 2020), volcano monitoring is critical for ensuring public safety and for mitigating the impacts of volcanic activity. Accurate and timely forecasts are facilitated by well-designed monitoring networks that are in place long enough to allow for background behavior to be recognized and understood. Because precursory signals may be limited and unrest may progress rapidly to an eruption, our goal is to deploy monitoring systems that enable detection of the reactivation of dormant volcanoes as early as possible, allowing for public safety and risk mitigation. NVEWS planning is also informed by the results of Ewert and others (2005, 2018), whereby 161 U.S. volcanoes are currently categorized and ranked commensurate with their relative threat. In each chapter, author(s) considered the need for some redundancy of instrumentation and telemetry, given the likelihood of occasional equipment failure, particularly in extreme and remote environments. Establishing digital telemetry networks requires advanced planning, sighting, radio-shot testing, and, inevitably, troubleshooting in the field. This is harder to achieve rapidly during a crisis; thus, an important goal for monitoring U.S. volcanoes is to establish digital telemetry backbones with redundancy and extra capacity to absorb additional instruments should a volcano begin to exhibit signs of unrest (fig. A1). The National Telecommunications and Information Administration (NTIA) imposed new regulations in the United States, eliminating the use of older analog radios for many purposes, which had been one previous means for redundant data delivery. However, the resulting conversion from analog to digital systems usefully enables stations to accommodate new and multivariate real-time data streams (for example, Global Navigation Satellite System [GNSS] receivers, infrasound arrays, gas spectrometers, visible and infrared cameras, and broadband seismometers). We note that other USGS and broader national and international hazard programs can leverage NVEWS instrumentation plans. Examples of this include the following: Improved seismic coverage of volcanoes will increase the capability of the USGS Earthquake Hazards Program to detect and locate earthquakes, estimate ground shaking, and provide timely early warnings through the ShakeAlert Earthquake Early Warning System (Given and others, 2018). The National Oceanic and Atmospheric Administration’s Tsunami Program will benefit from additional seismic stations, particularly within the sparsely instrumented Aleutian Islands, Northern Mariana Islands, and American Samoa. Infrasound stations can detect signals from landslides, debris flows and lahars, floods, and weather events, providing benefits to the National Weather Service and the USGS Landslide Hazards Program.

Scientific Investigations Report↗
USGS2024

Are field observations of surface rupture useful? An example from the 2023 Mw 7.8 Pazarcık, Turkey (Türkiye), earthquake

Field investigations have long been an important component of the scientific response to surface‐faulting earthquakes. However, in light of advances in remote data and models, the question arises whether field‐based observations of surface rupture remain useful for understanding rupture processes and seismic hazards. We approach this question using a field‐based study of the central 2023 M w 7.8 Pazarcık, Turkey (Türkiye), earthquake rupture, at the intersection of the east Anatolian fault (EAF) and Narlı fault. Our field observations include the surface rupture expression and extent of the central EAF and northernmost Narlı fault in generally forested and steep terrain and 68 measurements of left‐lateral surface displacement. These data improve our understanding of the Pazarcık rupture complexity, resolve the surface geometry of the Narlı fault–EAF connection, and exhibit a clear (>2 m) change in surface displacement across this intersection zone that confirms remote‐based coseismic slip models. Our study shows that focusing field efforts in areas of obscured or low‐resolution remote data can yield essential data for refining rupture extent, documenting perishable on‐fault displacement, and improving postearthquake situational awareness. A comparison of similarly large‐magnitude continental surface‐rupturing earthquakes indicates that displacement uncertainties relate to a complex set of factors, including measurement methods, rupture complexity, and displacement magnitude. Our study validates the need for postearthquake field observations, which, when driven by clear motivating questions and knowledge of methodological strengths and limitations, provide high‐resolution rupture data that complement remote‐based models.

Pazarcık↗
USGS2024

Ungulate migrations of the Western United States, volume 4

Broadly distributed across the Western United States, ungulates (hooved mammals) play an important role in ecosystem function by affecting vegetation communities and forming the prey base for large carnivores. Additionally, ungulates provide economic benefits to regional communities through tourism and hunting and hold cultural significance for many Tribal communities. Many ungulates migrate seasonally between distinct summer and winter ranges to take advantage of spatially and temporally variable food sources and avoid threats such as predators and deep snow. Increasingly, these migrations are threatened by the growing human footprint and associated subdivisions, energy development, and increased traffic volume. Efforts to study ungulate populations and conserve their migrations received support in recent years from the U.S. Department of the Interior Secretarial Order No. 3362, which provided Federal support for enhancing habitat quality for ungulates across the Western States. In response to Secretarial Order No. 3362, the U.S. Geological Survey (USGS) established the Corridor Mapping Team, a collaboration among USGS and participating State and Federal wildlife management agencies and numerous Tribal Nations. Together, the Corridor Mapping Team maps ungulate migrations throughout the Western United States in the USGS “Ungulate Migrations of the Western United States” report series. This report (volume 4) details migrations and seasonal ranges from 31 new herds throughout nine Western States. Additionally, this report includes updates to two herds published in previous reports. Including this report, the report series has provided the mapped migrations and seasonal ranges of 182 unique herds and has provided a map-based inventory of the documented ungulate migrations across the Western United States for biologists, managers, policy makers, and conservation practitioners. This report also discusses how the mapping efforts associated with the Corridor Mapping Team can be used to guide management and policy regarding renewable energy development and ungulate disease, specifically chronic wasting disease, in the Western United States.

Arizona, California, Colorado, New Mexico, Nevada,↗
USGS2026

Microclimate mediates the strength and direction of avian biotic interactions

Theory predicts that that the strength and direction of species interactions can shift from being competitive in benign environments toward being facilitative in stressful environments. However, the environmental context dependency of species interactions has rarely been tested in animal communities. We capitalized on a 15-year, landscape-scale dataset, collected annually in a relatively stable old-growth forest environment to test the long-held hypothesis that the strength and direction of species interactions might be mediated by climatic conditions. It is generally accepted that competitive and facilitative interactions drive the distributions of many species. Using multi-species dynamic occupancy models applied to long-term data, we tested whether annual settlement by bird species could affect either the persistence or settlement by other phylogenetically related species, and whether these interactions are mediated by microclimate. We found that species interactions were influenced by microclimate for some, but not all avian species pairs. Related species pairs more often showed settlement dynamics that were indicative of attraction rather than repulsion. In some cases, competitive interactions at warmer microclimates flipped to become facilitative in colder ones. However, the reverse was also true: facilitative interactions amplified for some species under warm conditions. Furthermore, species pairs that were closely related were more likely to exhibit competitive relationships along at least part of the microclimatic gradient. Our results highlight the importance of using long-term data to incorporate competitive and facilitative interactions into species distribution models and support the notion that the strength and direction of species interactions can be dependent on microclimatic environmental conditions.

Oregon↗
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↗
USGS2023

The Saudi Geological Survey-U.S. Geological Survey northern Harrat Rahat project—Styles, rates, causes, and hazards of volcanism near Al Madīnah al Munawwarah, Kingdom of Saudi Arabia

Active volcanic systems pose serious hazards to people and property including inundation and incineration by lava, blanketing by tephra (volcanic ash), exposure to noxious volcanic gases, and damage from shallow earthquakes triggered by ascending molten material (magma). To improve understanding of volcanism and associated seismicity on the western Arabia Plate, the Saudi Geological Survey and the U.S. Geological Survey conducted a multi-year investigation of the northern Harrat Rahat volcanic field adjacent to the city of Al Madīnah al Munawwarah, Kingdom of Saudi Arabia. Project components included creation of a high-resolution digital topographic base; interpretation of eruptive history supported by detailed geologic mapping, paleomagnetism, and abundant high-precision geochronology of volcanic deposits; assessments of eruptive styles and volcanic hazards by physical volcanology; investigation of the origins of magmas in the mantle and of their differentiation in the crust revealed by chemical and isotopic petrology; gravity and magnetotelluric surveys to reveal crustal structures and to search for magma reservoirs; and regional and local seismic tomography and analyses of seismic hazards. Project results are presented in this Professional Paper as chapters written for technical scientific audiences. This initial chapter introduces the project and briefly summarizes results in plain language for readers who have more general backgrounds.

Al Madīnah al Munawwarah↗
USGS2025

Assessing flood water infiltration and storage in a restored floodplain

In urban areas, floodplain restoration is gaining prominence as a strategy for restoring the natural functions of floodplain ecosystems and reducing flood risk. This has spurred research into potential interactions between floodwaters, the hyporheic zone, and the floodplain aquifer. An urban restored stream in Wisconsin, USA, was used as a case study to examine four methods to estimate floodplain infiltration and storage during overbank floods. We characterised flood-related infiltration over a 4-year period from 2018 through 2021 by simultaneously and continuously measuring groundwater levels and vertical temperature profiles with stream water levels linked to high-resolution flood inundation maps. High-resolution topographic data helped to quantify surface floodplain storage and the unsaturated soil volume relative to flood stage. Infiltration estimates from the simple methods align well with those from the more complex methods; however, the complex methods provide additional insights about the factors influencing infiltration. Results from all methods indicate that the volume of water that vertically infiltrates during floods is likely small relative to the total volume of the flood, with 0.08%–0.52% of flood water infiltrating into the floodplain, on average. Spatially variable vertical hydraulic gradients, driven by flood depth, groundwater level, and permeability, imply heterogeneous patterns of infiltration across the floodplain. Gradients favourable for infiltration typically occurred during the onset of flooding but, over the study period, were mostly (98% of the time) favourable for groundwater discharge to the channel (non-flood periods). These findings highlight the importance of considering surface-groundwater dynamics, floodplain soils, and unsaturated floodplain volume in defining the benefits of floodplain infiltration for flood attenuation.

Wisconsin↗
USGS2026

Magnitude conversion relations create substantial differences in seismic hazard models

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Seismological Research Letters↗
USGS2025

A soil velocity model for improved ground motion simulations in the U. S. Pacific Northwest

Near-surface seismic velocity structure may significantly impact the intensity, duration, and frequency content of ground shaking during an earthquake. In this study, we compile 649 shear wave velocity (Vs) profiles throughout the U.S. Pacific Northwest and southern British Columbia (PNW) and use these measured profiles to develop a representative soil velocity model for four major Holocene soil provinces: Puget Lowlands, Willamette Valley, fill and alluvium, and `other' soils. The resulting soil velocity model shows good agreement to measured data for a wide range of site conditions, with variability between different geologic domains reflecting fundamental differences in depositional environments. We then show that using this regional soil velocity model in simulations of the 2001 M6.8 Nisqually, Washington earthquake improves the fit to observed high-frequency (≥ 0.5 Hz) ground motions in the Puget Sound region compared to simulations that do not incorporate shallow (≤ 200 m) seismic velocity structure. Overall, this work shows that incorporating localized soil velocity profiles into seismic velocity models is important for accurately estimating high-frequency ground motion and regional seismic hazard in earthquake simulations. Future earthquake simulations and hazard studies in the PNW could incorporate these soil velocity profiles to capture the region's distinct site response characteristics.

Washington↗
USGS2024

Temperature

The “National Field Manual for the Collection of Water-Quality Data” (NFM) provides guidelines and procedures for U.S. Geological Survey (USGS) personnel who collect data used to assess the quality of the Nation’s surface-water and groundwater resources. This chapter, NFM A6.1, provides guidance and protocols for the measurement of temperature of air, of a surface-water body or in groundwater, which include the scientific basis of the measurement, selection and maintenance of equipment, calibration verification, troubleshooting, and procedures for measurement and reporting. It updates and supersedes USGS Techniques of Water-Resources Investigations, book 9, chapter A6.1, version 2.0, by Franceska D. Wilde. Temperature of air and water is routinely measured when water samples are collected, is often measured continually at USGS streamgages, and is a parameter regularly measured during laboratory and field experiments. The field method for measuring temperature described in this chapter is applicable to air and most natural waters. Before 2017, the NFM chapters were released in the USGS Techniques of Water-Resources Investigations series. Effective in 2018, new and revised NFM chapters are being released in the USGS Techniques and Methods series; this series change does not affect the content and format of the NFM. More information is in the general introduction to the NFM (USGS Techniques and Methods, book 9, chapter A0) at https://doi.org/10.3133/tm9A0 . The authoritative current versions of NFM chapters are available in the USGS Publications Warehouse at https://pubs.usgs.gov/ . Comments, questions, and suggestions related to the NFM can be addressed to nfm@usgs.gov .

Techniques and Methods↗
USGS1997

Nutrients and pesticides in ground water of the Ozark Plateaus in Arkansas, Kansas, Missouri, and Oklahoma

A total of 229 ground-water samples were collected from 215 sites as part of the Ozark Plateaus study unit of the National Water-Quality Assessment Program. These samples were collected from 1993 through 1995 using a network of springs and wells with three scale-dependent components. The first component, the study-unit survey, consisted of 99 randomly selected springs and domestic wells in the Springfield Plateau and Ozark aquifers. The second component, two land-use studies, consisted of 42 springs and domestic wells in a poultry-dominated agricultural area and 40 springs and domestic wells in a cattle-dominated agricultural area overlying the Springfield Plateau aquifer. The third component, the small-watershed study, consisted of 4 springs, 18 domestic wells, and 11 monitoring wells in a small basin within the poultry land-use study area. Samples were analyzed for major ions, nutrients, dissolved organic carbon, methylene blue active substances, tritium, and 88 pesticides and metabolites.The water-quality data from these samples were analyzed with descriptive and statistical methods. Nitrite plus nitrate, which was detected more often and in greater concentrations than any of the other nutrients, ranged from less than 0.05 to 25 milligrams per liter as nitrogen. Nitrite plus nitrate concentrations positively correlated to percent agricultural land use around each site. Median nitrite plus nitrate concentrations generally were greater in samples from springs than in samples from wells. Concentrations of nitrite, ammonia, and ammonia plus organic nitrogen were also affected by land use and also by concentrations of dissolved oxygen in the ground water. Concentrations of phosphorus and orthophosphate probably were affected by land use and also by phosphorus solubility. Pesticides were detected in 80 of 229 samples from 73 of 215 sites. A total of 20 pesticides were detected with a maximum of 5 pesticides detected in any 1 sample. The most commonly detected pesticides were tebuthiuron, atrazine, prometon, desethylatrazine, and simazine. Maximum concentrations ranged from 0.003 to 1.0 microgram per liter. The occurrence and distribution of pesticides were related to land use. Percent agricultural land use was greater for samples with pesticides detected than for samples with no pesticides detected. Pesticides were detected more often in samples from springs than in samples from wells. The occurrence of pesticides also was related to seasonality and chemical characteristics, such as solubility and persistence, of the compounds.

Arkansas, Kansas, Missouri, Oklahoma↗
USGS2026

Bathymetric survey and storage capacity of Upper Lake Mary near Flagstaff, Arizona in 2024

The U.S. Geological Survey (USGS), in cooperation with the city of Flagstaff, collected bathymetric, light detection and ranging (lidar), and land-survey data of Upper Lake Mary in Arizona during the months of April and October 2024. The city of Flagstaff uses a combination of groundwater from well fields throughout the Flagstaff area and surface water, mainly from Upper Lake Mary, for its potable water supply. The purpose of the survey is to update previous surveys using new technology and compare the results to previous surveys to determine if there was a decrease in storage capacity that could affect the city’s water supply. The lakebed was mapped in April 2024 using a vessel equipped with a multibeam echosounder (MBES) and mobile lidar scanner with positioning captured using a real-time kinematic global navigation satellite system (RTK GNSS) base and receivers. In October 2024, areas of the reservoir that were too shallow for the boat and shoreline that were not captured by the vessel-based lidar were surveyed on foot using hand-held RTK GNSS receivers. At full pool (spillway elevation of 6,831.82 feet above NAVD 88 [2,082.34 meters (m)], Upper Lake Mary has a storage capacity of 16,449.80 acre-feet (20,290,611.73 cubic meters) and a surface area of 953.57 acres (3,860,926.075 square meters). The reservoir is 5.7 miles (9.7 kilometers) long and varies in width from 326 feet (99.36 m) near the central, narrow portion of the reservoir to 2,613 feet (796.44 m) in the upper portion. Comparisons between this survey and the previous two surveys from the 1950s and 2006 indicate no apparent decrease in reservoir area or storage capacity. Results of the 2024 survey indicate that Upper Lake Mary’s storage capacity increased by 0.9 percent from the 2006 survey and a 1.6 percent increase in surface area from the 2006 survey.

Arizona↗
USGS2025

Assimilation of reduced carbon triggers platinum alloy saturation in mafic and ultramafic magmas

It is generally observed that magmatic sulfide ores have higher ratios of Pd/Pt than the mantle-like values of their parental magmas. This discrepancy has defied simple explanation because the partitioning behavior of both elements between sulfide and silicate liquids is very similar. Assimilation of sulfur- and carbon-rich country rocks by mafic and ultramafic magmas is considered a critical, if not essential, step in the formation of magmatic base metal sulfide deposits. Although there is general consensus that the assimilation of external sulfur and carbon promotes sulfide saturation, the effect of carbon assimilation on the solubilities of platinum-group elements in natural S-bearing silicate melt has been overlooked. In this study, we investigate the variations of platinum and palladium solubilities during assimilation of graphite and methane through thermodynamic modeling, in comparison with data from an array of highly distinctive magmatic sulfide ore systems representing ages from Archean to Paleozoic, melt compositions from komatiite to basalt, and magmatic settings including lavas, hypabyssal intrusions, plutonic continental arc roots, and plutonic layered intrusions, namely: Raglan, Norilsk-Talnakh, Lac des Iles, and the J-M Reef of the Stillwater Complex. We model assimilation-fractional crystallization processes to estimate the reduction of oxygen fugacity ( ⁠fO2⁠ ) of the melt due to incorporation of graphite and methane. The simulations show that although Pd remains highly soluble during the progressive assimilation of reduced carbon, Pt solubility decreases significantly as the silicate melt becomes increasingly reduced. With less than 8 % of sediment assimilation, Pt alloy may saturate and then deviate from sulfide-undersaturated silicate melts, concomitantly increasing the Pd/Pt value of the remaining melts of the Raglan and Norilsk-Talnakh systems. For the Lac des Iles and Stillwater systems, a higher extent of assimilation is needed to reach Pt saturation because of the relatively carbon-poor nature of the lower crustal rocks. The assimilation of methane volatiles is shown to be more effective than graphite assimilation, and it provides a pathway to Pt alloy fractionation in the absence of detectable amounts of bulk host-rock assimilation. High Pd/Pt values have been documented in many world-class magmatic sulfide deposits whose parental magmas have demonstrably experienced crustal contamination. Our model suggests that although anomalous Pd/Pt values may be explained by other mechanisms such as incongruent melting of preexisting sulfide or differences in the diffusivities of the metals within achieving equilibration, the assimilation of graphite or methane may play an important role in the global occurrence of magmatic sulfide ores with elevated Pd/Pt values.

Economic Geology↗
USGS2025

Groundwater salinity: Applying the specific conductance and water type proxy

Groundwater is increasingly needed for water supplies but may have limited utility in some locations because of its salinity. Salinity, often expressed as total dissolved solid (TDS), is frequently estimated using specific conductance (SC) measurements. However, the commonly used proxy (0.65 multiplied by SC to indicate TDS, common in many handheld meters) can result in inaccurate TDS estimates. First, the TDS–SC relationship is not linear over the entire concentration range of groundwater. Furthermore, the TDS (and salinity)–SC relationships vary substantially depending on the major-ion composition. Here we develop a proxy method utilizing SC and major-ion water type to estimate TDS and salinity specifically for groundwaters. Compared to most surface waters, groundwater tends to have a wider range of salinity (fresh to highly saline) and higher concentrations of bedrock-derived solutes such as carbonate ions, silica, and many other ions. The dataset used to develop the proxies includes water chemistry data from 149,059 discrete groundwater samples. The groundwater proxies, which employ nonlinear log–log relations, utilize five water types (HCO 3 , Cl, Ca-Mg-SO 4 , Na-K-SO 4 , and mixed waters), are accurate (median percent difference between TDS and salinity determined using the proxy compared to discrete measurements was <±0.8%) over a wide range of SC (up to 200 mS/cm), rapid, cost-effective, and can be measured on-site.

Groundwater↗
USGS2025

Mine waste as a potential source of critical minerals and other commodities: Examples from the Four Corners states, USA

The growing demand for critical minerals and other mineral resources has raised concerns about possible supplies of these essential commodities. Mine waste is a potential source of these essential commodities. We compiled a geospatial database of publicly available data of the largest non-fuel mine waste features (>200,000 m 2 in areal extent) in the Four Corners states of the United States, where most of those features are from porphyry Cu deposits. The combined volume of those large porphyry Cu mine waste features is approximately 17 cubic kilometers, 60% of which is ore-related material such as tailings. Using publicly available data on density, grade, and previous recovery values, we estimate the contained endowments of Ag, As, Au, Bi, Cu, Mo, Re, S, Sb, Se, Te, and W. These estimates indicate endowments within ore-related mine wastes are collectively comparable to endowments of giant to supergiant deposits. If fully recovered, these commodities could meet current global demand from less than 1 year (Sb) to more than 200 years (Re), underscoring the enormous untapped resource potential of mine waste.

Arizona, Colorado, New Mexico, Utah↗
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↗
USGS2025

Deterministic physics-based earthquake sequence simulators match empirical ground-motion models and enable extrapolation to data poor regimes: Application to multifault multimechanism ruptures

We use the deterministic earthquake simulator RSQSim to generate complex sequences of ruptures on fault systems used for hazard assessment. We show that the source motions combined with a wave propagation code create surface ground motions that fall within the range of epistemic uncertainties for the Next Generation Attenuation‐West2 set of empirical models. We show the model is well calibrated where there are good data constraints, and has good correspondence in regions with fewer data constraints. We show magnitude, distance, and mechanism dependence all arising naturally from the same underlying friction. The deterministic physics‐based approach provides an opportunity for better understanding the physical origins of ground motions. For example, we find that reduced stress drops in shallow layers relative to constant stress drop with depth lead to peak ground velocities in the near field that better match empirical models. The simulators may also provide better extrapolations into regimes that are poorly empirically constrained by data because physics, rather than surface shaking data parameterizations, is underlying the extrapolations. Having shown the model is credible, we apply it to a problem where observations are lacking. We examine the case of crustal faults above a shallow subduction interface seen to break coseismically in simulations of the New Zealand fault system. These types of events were left out of consideration in the most recent New Zealand national seismic hazard model due to the modeling complexity and lack of observational data to constrain ground‐motion models (GMMs). Here, we show that in the model, by breaking up the coseismic crustal and interface rupturing fault motions into two separate subevents, and then recombining the resulting ground‐motion measures in a square‐root‐of‐sum‐of‐squares incoherent manner, we reproduce well the ground‐motion measures from the full event rupture. This provides a new method for extrapolating GMMs to more complex multifault ruptures.

Seismological Research Letters↗