USGS Science⌕ Search

SEARCH · USGS Science

Results for “Science of the Total Envionrment”

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

Peak streamflow trends in Montana and northern Wyoming and their relation to changes in climate, water years 1921–2020

Frequency analysis on annual peak streamflow (hereinafter, peak flow) is essential to water-resources management applications, including critical structure design (for example, bridges and culverts) and floodplain mapping. Nonstationarity is a statistical property of a peak-flow series such that the distributional properties (the mean, variance, or skew) change either gradually (monotonic trend) or abruptly (shift, step change or change point) through time. Not incorporating or accounting for observed nonstationarity into peak-flow frequency analysis might result in a poor representation of the true probability of large floods and thus misrepresent the actual flood risks to life and property. This report summarizes how hydroclimatic variability might affect the temporal and spatial distributions of peak-flow data in the State of Montana (and northern Wyoming) and is part of a larger study to document peak-flow nonstationarity and hydroclimatic changes across a nine-State region consisting of Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flows. This study was completed in cooperation with the Montana Department of Natural Resources and Conservation. The purpose of this report is to characterize temporal and spatial patterns of nonstationarity in peak flows and hydroclimatology in Montana and northern Wyoming. In this evaluation, peak-flow, daily streamflow, and model-simulated gridded climatic data were examined for monotonic trends, change points, and other statistical properties that might indicate changing climatic and environmental conditions. This report includes background information on the study area, the history of U.S. Geological Survey peak-flow data collection and frequency analysis in Montana, and the review of research relating to hydroclimatic variability and change in Montana. This study might help provide a framework for addressing potential nonstationarity issues in peak-flow frequency updates that commonly are completed by the U.S. Geological Survey in cooperation with other agencies throughout the Nation. The analytical structure of this study includes analyses of monotonic trends and change points in numerous hydroclimatic variables in assigned 30-, 50-, 75-, and 100-year analysis periods. For Montana and part of Wyoming, the 30-, 50-, 75, and 100-year analyses included 157, 70, 48, and 12 streamgages, respectively. For those streamgages, nonstationarities were analyzed in the following variables: (1) climatic variables, including annual and seasonal (winter, spring, summer, and fall) temperature and precipitation; (2) daily streamflow variables, including the annual center of volume duration, annual center of volume median, and peaks over threshold with a mean of four events per year; and (3) annual peak-flow variables, including peak-flow timing and magnitude. A likelihood approach was used to express statistical confidence and assign the nonstationarity results as likely upward or downward (highest statistical confidence), somewhat likely upward or downward (less statistical confidence), or about as likely as not (little statistical confidence; hereinafter, neutral). For the nonstationarity analyses of the climatic, daily streamflow, and peak-flow variables, the results are presented in detail and discussed with respect to statewide patterns and geographic variability. For each of the 30-, 50-, and 75-year analyses, peak-flow change-point and monotonic trend analyses were compiled for streamgages classified with likely downward or likely upward trends. For those streamgages, the associated basin characteristics and nonstationarity results for peak-flow timing, daily streamflow, and climatic variables were investigated and statistically compared to discern associations among other variables that might contribute to the peak-flow nonstationarity results. The 50- and 75-year peak-flow nonstationarities identified in this study are mostly downward, in association with mostly upward temperature and potential evapotranspiration:precipitation monotonic trends. For the 50-, 75-, and 100-year analyses, the peak-flow change points are predominantly downward and are concentrated in the 1970s and 1980s, which indicates general consistency among the longer trend periods. These findings are in association with substantial research documenting globally rising temperature and atmospheric greenhouse gas concentrations that might be largely attributed to anthropogenic activities. Anthropogenic effects might represent long-term (on the order of several decades to more than a century) climate changes that might happen within highly variable natural climate fluctuations. Several paleo studies in the north-central United States have indicated that hydroclimatic extremes (that is, low- and high-streamflow conditions) before European settlement have been outside of extremes since the 1900s. Depending on the interactions of anthropogenic effects and natural climate variability, extreme high-streamflow conditions might occur in the future, even in the presence of long-term downward peak-flow trends.

Montana, Wyoming↗
USGS2026

From collision to subduction: Thermal-kinematic inversions constrain plate boundary structure and dip-slip activity in southwestern New Zealand

Constraints on three-dimensional fault system geometry and kinematics are required to understand both long-term tectonics and contemporary seismic hazard. In many cases, these parameters must be inferred from a combination of different data or analyses. Here, we invert low-temperature thermochronological data to quantitatively parameterize both the crustal-scale geometry and late Cenozoic dip-slip activity of the plate boundary system in southwest New Zealand. We find that the modern transpressive plate boundary system has accommodated the bulk of Australian-Pacific Plate convergence since establishment in the middle Miocene-early Pliocene. Shortening has occurred via (a) dextral-reverse slip on the moderate angle, listric central Alpine Fault décollement, and (b) reverse-slip on either the sharply curved, anti-listric northern Puysegur subduction interface or low-angle offshore thrust fault system that splays from the high-angle southern Alpine Fault at depth. We resolve late Cenozoic shortening rates on these structures that are indistinguishable, and interpret these to reflect the transfer of convergent deformation outboard of the predominantly strike-slip southern Alpine Fault. Our work helps elucidate the transition from strain localized oblique continental collision to strain partitioned subduction. We also provide an estimate of the across-strike width of the central-southern Alpine Fault section boundary at seismogenic depths. This is the first quantitative constraint on the geometry of a structural feature inferred to have arrested several major prehistoric earthquakes. This study furthers understanding of fault structure and dip-slip kinematics along the greatest sources of seismic hazard in and offshore of southwest New Zealand.

South Island↗
USGS2026

USGS Geochron—A database of geochronological and thermochronological dates and data—Technical documentation

Geochronological and thermochronological data are essential for constraining the timing and rates of geological processes, supporting geologic mapping, natural hazard assessment, and resource exploration. The U.S. Geological Survey (USGS) Geochron Database is a centralized, relational database that integrates USGS and State geological survey data in accordance with the National Geologic Mapping Act of 1992 and its reauthorizations. This report documents the structure and development of the database, including its standardized schema, controlled vocabularies, and protocols for compiling legacy and newly published data. The database is designed to adhere to FAIR (Findable, Accessible, Interoperable, and Reusable) data principles and currently (2026) includes data from a wide range of geochronological and thermochronological methods. Public access is provided through the USGS Geochron Database Explorer, a browser-based geographic information system (GIS) interface, and through versioned data releases available on ScienceBase, which provide sample-level summary data and detailed analytical information. The USGS Geochron Database is structured to be extensible, supporting the inclusion of additional methods and data types as they are compiled. This format ensures long-term utility and aligns with the Mapping Act’s directive to create a national archive of geochronological information that adds interpretive value to geologic map data. Its design reflects a commitment to data transparency, scientific reproducibility, and national geoscience priorities.

Data Report↗
USGS2025

ShakeAlert® version 3: Expected performance in large earthquakes

The ShakeAlert earthquake early warning (EEW) system partners along with U.S. Geological Survey (USGS) licensed operators deliver EEW alerts to the public and trigger automated systems when a significant earthquake is expected to impact California, Oregon, or Washington. ShakeAlert’s primary goal is to provide usable warning times before the arrival of damaging shaking. EEW is most likely to achieve this goal in large‐magnitude earthquakes. In recent years, ShakeAlert has gone through a series of upgrades to its underlying scientific algorithms aimed at improved performance during large earthquakes. Version 3 of this software recently went live in the production system and includes improvements to all algorithms. The main seismic algorithms that detect an earthquake and characterize its location, magnitude, and fault rupture orientation are faster than older versions. Other key changes include: using real‐time geodetic data to characterize the magnitude growth in large earthquakes; the introduction of an alert pause procedure to compromise between speed near the epicenter and improved accuracy at larger distances; and the inclusion of a nonergodic site‐response model in the ground‐motion predictions. ShakeAlert has achieved its primary goal of usable warning times before strong shaking at some locations in real‐time operations in recent M 6 earthquakes. Using offline tests, we demonstrate usable warning times are possible for many sites with peak shaking values of modified Mercalli intensity (MMI) 7–8 in M 7+ earthquakes and also for many MMI 8–9 sites in M 8+ earthquakes. ShakeAlert partners use a variety of MMI and magnitude thresholds in deciding when to alert their users within bounds set by the USGS. Our study shows that there is room to raise the magnitude thresholds up to about M 5.5 without adversely affecting performance in large earthquakes. The ground‐motion criteria are more complex owing to a significant drop‐off in warning times between the MMI 4 and 5 levels of predicted shaking. However, widely used ShakeAlert products, such as the MMI 3 and 4 contour products, can provide sufficiently long warning times before strong shaking in moderate‐to‐great earthquakes to enable a range of protective actions.

Bulletin of the Seismological Society of America↗
USGS2017

Field-trip guide to Mount Hood, Oregon, highlighting eruptive history and hazards

This guidebook describes stops of interest for a geological field trip around Mount Hood volcano. It was developed for the 2017 International Association of Volcanology and Chemistry of the Earth’s Interior (IAVCEI) Scientific Assembly in Portland, Oregon. The intent of this guidebook and accompanying contributions is to provide an overview of Mount Hood, including its chief geologic processes, magmatic system, eruptive history, local tectonics, and hazards, by visiting a variety of readily accessible localities. We also describe coeval, largely monogenetic, volcanoes in the region. Accompanying the field-trip guidebook are separately authored contributions that discuss in detail the Mount Hood magmatic system and its products and behavior (Kent and Koleszar, this volume); Mount Hood earthquakes and their relation to regional tectonics and the volcanic system (Thelen and Moran, this volume); and young surface faults cutting the broader Mount Hood area whose extent has come to light after acquisition of regional light detection and ranging coverage (Madin and others, this volume). The trip makes an approximately 175-mile (280-kilometer) clockwise loop around Mount Hood, starting and ending in Portland. The route heads east on Interstate 84 through the Columbia River Gorge National Scenic Area. The guidebook points out only a few conspicuous features of note in the gorge, but many other guides to the gorge are available. The route continues south on the Mount Hood National Scenic Byway on Oregon Route 35 following Hood River, and returns to Portland on U.S. Highway 26 following Sandy River. The route traverses rocks as old as the early Miocene Eagle Creek Formation and overlying Columbia River Basalt Group of middle Miocene age, but chiefly lava flows and clastic products of arc volcanism of late Miocene to Holocene age.

Oregon↗
USGS1997

Hydrologic and water-quality conditions in the Horse Creek basin, west-central Florida, October 1992-February 1995

A baseline study of the 241-square-mile Horse Creek basin was undertaken from October 1992 to February 1995 to assess the hydrologic and water-quality conditions of one of the last remaining undeveloped basins in west-central Florida. During the period of the study, much of the basin remained in a natural state, except for limited areas of cattle and citrus production and phosphate mining. Rainfall in 1993 and 1994 in the Horse Creek basin was 8 and 31 percent, respectively, above the 30-year long-term average. The lowest and highest maximum instantaneous peak discharge of the six daily discharge stations occurred at the Buzzard Roost Branch and the Horse Creek near Arcadia stations with 185 to 4,180 cubic feet per second, respectively. The Horse Creek near Arcadia station had the lowest number of no-flow days with zero days and the Brushy Creek station had the highest number with 113 days. During the study, the West Fork Horse Creek subbasin had the highest daily mean discharge per square mile with 30.6 cubic feet per second per square mile, and the largest runoff coefficient of 43.7 percent. The Buzzard Roost Branch subbasin had the lowest daily mean discharge per square mile with 5.05 cubic feet per second per square mile, and Brushy Creek and Brandy Branch shared the lowest runoff coefficient of 0.6 percent. Brandy Branch had the highest monthly mean runoff in both 1993 and 1994 with 11.48 and 19.28 inches, respectively. During the high-baseflow seepage run, seepage gains were 8.87 cubic feet per second along the 43-mile Horse Creek channel. However, during the low-baseflow seepage run, seepage losses were 0.88 cubic foot per second. Three methods were used to estimate average annual ground-water recharge in the Horse Creek basin: (1) well hydrograph, (2) chloride mass balance, and (3) streamflow hydrograph. Estimated average annual recharge using these three methods ranged from 3.6 to 8.7 inches. The high percentage of carbonate plus bicarbonate analyzed at the Carlton surficial aquifer well could indicate an upward ground-water flow from the underlying intermediate aquifer system. Based on constituent concentrations in water samples from the six daily discharge stations, concentrations generally are lower in the upper three subbasins, West Fork Horse Creek, Upper Horse Creek, and Brushy Creek than in the lower three subbasins. Typically, concentrations were highest for major ions at Buzzard Roost Branch and nutrients at Brushy Creek.

Florida↗
USGS2025

Observations of tear-drinking by lepidopterans on moose ( Alces alces americana ) in northeastern North America

Lepidoptera have long been known to feed on the tears of vertebrates as a presumed source of minerals or nutrients. While this unusual behavior has been observed in a variety of species, only a single previous record has been documented outside of the tropics. Here, we present the first documentation of moths visiting the eyes of a bull moose ( Alces americanus americanus ), captured via trail camera in Green Mountain National Forest, Vermont, United States. We discuss the biogeography of this behavior, how it may differ between tropical and temperate climates, and its potential impact on moose health.

Vermont↗
USGS2025

AAPG Energy and Minerals Division Tight Oil and Gas Committee Activities and Commodity Report for 2021-2022: Tuscaloosa Marine Shale, Gulf Coast basin, Louisiana and Mississippi

The Upper Cretaceous Tuscaloosa marine shale (TMS) potential production area encompasses 20.4 million acres across central Louisiana (LA), southern Mississippi (MS), southwestern Alabama (AL), and a small southwestern section of the Florida panhandle (Hackley et al., 2018). It remains a minor and largely undeveloped unconventional shale oil play with production from the TMS confined along the east-west LA-MS State boundary (Fig. 1). The mean undiscovered, technically recoverable resources in the TMS are estimated at 1.5 billion barrels of oil and 4.6 trillion cubic feet of gas (Hackley et al., 2018). Geochemical analyses of source rock solvent extracts and oil samples indicate that, in the play area, the TMS is the primary source of shale oil produced from the TMS as well as of accumulated oils in the underlying conventional reservoirs of the lower Tuscaloosa (Hackley et al., 2020). The API gravity of TMS oils ranges from approximately 34 to 46 degrees (Hackley et al., 2020; Croke et al., 2020).

AAPG Wiki↗
USGS2025

Leveraging wildfire to augment forest management and amplify forest resilience

Successive catastrophic wildfire seasons in western North America have escalated the urgency around reducing fire risk to communities and ecosystems. In historically frequent-fire forests, fuel buildup as a result of fire exclusion is contributing to increased fire severity. The probability of high-severity fire can be reduced by active forest management that reduces fuels, prompting federal and state agencies to commit significant resources to increase the pace and scale of fuel reduction treatments. However, lower severity areas of wildfires also have the potential to act as “treatments,” and even catastrophic fires with large areas of high severity can still have substantial areas of lower severity fire that may be improving forest conditions locally. We quantified active management and wildfire severity across yellow pine and mixed conifer (YPMC) forests in the Sierra Nevada of California over a 22-year period (2001–2022). We did not detect increases in the area treated through time, but the area of beneficial wildfire (low to moderate severity) increased substantially, exceeding active treatment area in 8 of 22 years. Overall, beneficial wildfire treated ~17% more area than all treatments combined, and roughly four times more area than fire-related treatments alone. We then used disturbance history to evaluate resistance to high-severity wildfire and forest loss across the YPMC range. Of the 2.3 million ha YPMC of forests in 2001, 20% lost mature forests due to high-severity fire by 2022, which is nearly half of all YPMC area burned. Most of the landscape (47%) remains at risk of high-severity fire because it had no restorative disturbances, but 33% of the study area has some level of resistance to high-severity wildfire. In these areas, resistance will need to be enhanced and maintained over time via active management or managed wildfire, but these treatment needs will likely outpace capacity even under optimistic implementation scenarios. Given limited resources for implementing active management and the likelihood of a more fiery future, incorporating beneficial wildfire into landscape-level treatment planning has the potential to amplify the impact of active management treatments.

California, Nevada↗
USGS2025

REDPy: A Python tool for automated repeating earthquake detection and visualization

Detecting and cataloging seismic events are among the most fundamental tasks in seismology. Many standardized tools for these tasks exist, including the open‐source package repeating earthquake detector in Python (REDPy). REDPy generates an organized catalog of seismic events from continuous waveform data, in which events are automatically separated into groups (“families”) by their waveform similarity through cross‐correlation. REDPy also automatically generates various outputs that allow a user to visualize important trends in the catalog, which may be used in real time or in retrospective analyses to allow rapid identification of interesting features. The code was designed for near‐real‐time volcano monitoring but is applicable across a broad range of use cases in seismology and seismoacoustics. In this article, the utility and performance of REDPy are demonstrated on two highly seismogenic volcanic eruption sequences: the onset of the dome‐building eruption of Mount St. Helens, Washington, from 2004 to 2005, and the entirety of the summit caldera collapse sequence of Kīlauea, Hawai‘i, in 2018. This article is meant to be a companion to the documentation of the code; in addition to detailing the basic required inputs, script functionality, and resulting outputs, the reasonings behind several important design decisions are also discussed.

Seismological Research Letters↗
USGS2024

Earthquake relocations delineate discrete a fault network and deformation corridor throughout Southeast Alaska and Southwest Yukon

Deformation in southeastern Alaska and southwest Yukon is governed by the subduction and translation of the Pacific-Yakutat plates relative to the North American plate in the St. Elias region. Despite notable historical seismicity and major regional faults, studies of the region between the Fairweather and Denali faults are complicated by glacial coverage and the remote setting. In the last decade, significant improvements have been made to the density of regional broadband seismometer networks. We relocate more than 5,000 earthquakes between 2010 and 2021 in the region of southeastern Alaska and southwestern Yukon utilizing these improved seismic networks. With reductions in catalog uncertainty, particularly in depth, we quantify the thickness of the seismogenic layer in the crust throughout the region and locate seismicity on a shallow network of upper-crustal faults. Relocated earthquakes, combined with an updated focal-mechanism catalog, permit estimating and classifying motion of active faults. This includes mapping the Totschunda-Fairweather “Connector” fault, which plays an important role in explaining regional deformation, and identifying new faults like the Kathleen Lake fault. We draw similarities between our seismic observations and simplified conceptual models of regional tectonics, which describe a dominant transpressional regime and localized slip partitioning. Our results support a hypothesis where current deformation is taking place on a well-defined and evolved network of shallow faults in the corridor between the Totschunda-Fairweather “Connector” and Denali faults.

Alaska, Yukon↗
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↗
USGS2025

Young explosive eruptions from the Clear Lake volcanic field

The Clear Lake volcanic field is the northernmost and youngest field in a chain of volcanic fields in and near the California Coast Range mountains. For 2 million years, numerous eruptions have happened around (and through) Clear Lake. The most recent period of activity in the Clear Lake volcanic field probably started around 40,000 years ago and was mainly explosive eruptions concentrated on and near faults in and around the lake. The combination of hot magma and groundwater created explosive releases of hot steam. The craters, called maars, threw out volcanic ash, pumice, and lava fragments, which fell back to the ground and draped the surrounding landscape. The deposits are visible in many places around the southeast end of the lake.

California↗
USGS2026

Distributed volcanic fields—Patterns, products, processes, and plumbing

Distributed volcanic fields—also referred to as “monogenetic fields” or simply “volcanic fields”—are found on every continent, on the seafloors, and in every tectonic setting on Earth. Yet, our understanding of what constitutes a distributed volcanic field, and our knowledge of the common products, spatiotemporal patterns in activity, and plumbing systems of these fields, are commonly poor and (or) debated. Here we provide an overview of the current understanding of distributed volcanic fields (DVFs) through review and compilation of published data, with an emphasis on understanding the characteristics of fields that are essential to prepare for future hazardous volcanic activity. Distributed fields can erupt restricted to variable magma compositions. Activity in these fields produces nearly all types of volcanoes, including scoria and spatter cones, maars, low shields, lava fields, stratovolcanoes, and calderas, but monogenetic volcanoes are the defining characteristic of DVFs. The number of vents and lifespans of DVFs vary greatly; vent numbers range from 2 to greater than 1,000, and fields can remain active for 150,000 to 10 million years. Some characteristics of DVFs appear field dependent. Although the volumes of single monogenetic eruptions within DVFs globally are variable (from less than 0.1 to a few cubic kilometers), each field may have a characteristic range in eruption volumes. Eruptions within distributed volcanic fields can have a random distribution and timing, but spatial and (or) temporal vent clustering and alignment are common, and eruptions may be time or volume predictable. Average eruption recurrence intervals in DVFs span over an order of magnitude, from less than 3 to greater than 100 thousand years. Exposures of the shallow portions of older systems reveal that many magma batches end their ascent before reaching the surface, with some dikes stopping within only meters of the surface. Petrologic insights into the deeper plumbing systems of DVF volcanoes indicate variable depths of magma crystallization and (or) stalling, and that magma transport and storage systems may be unique to a particular DVF and its underlying crustal structure. Both seismic and petrologic data for mafic DVF eruptions corroborate intrusions of magma into storage systems in the months to years before eruption, as well as variable ascent rates through the crust (days to weeks). Few distributed fields are closely monitored, as they erupt infrequently, but these observations imply that the amount of time available to install monitoring instrumentation after the first felt seismicity may be short, as precursory activity is frequently on the order of days to months for modern DVF eruptions. Satellite-based geodesy may provide useful warning in such cases. This chapter concludes with outstanding questions and areas for future research that will aid in our understanding and preparedness for future DVF eruptions.

Professional Paper↗
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

Tectonic controls on volcanism and associated hydrothermal activity in a sediment-dominated mid-ocean ridge; Escanaba Trough

Mid-ocean ridges, the Earth's most extensive volcanic system, exhibit unique characteristics in sediment-dominant environments. Thick sediment cover insulates the crust and channels fluid along pathways that can lead to the formation of distinct crustal alteration patterns, exceptionally large mineral deposits, and specialized chemosynthetic ecosystems. This study presents an interdisciplinary investigation into the tectonics of the Escanaba Trough, a heavily sedimented axial valley at the southern Gorda Ridge in the Northeast Pacific Ocean. A primary challenge in such environments is overcoming the masking effect of thick sediments on basement structures that control magmatic and hydrothermal activity. We address this by employing three-dimensional (3D) magnetic modeling of high-resolution near-seafloor magnetic data collected by an autonomous underwater vehicle (AUV). The 2022 surveys with AUV Sentry provided data for 3D magnetic susceptibility models, refining our understanding of the geometry of sub-sediment laccoliths/saucer-shaped sills and hydrothermal alteration. In conjunction with a new 1:100,000 scale lithostratigraphic map, we outline the tectonic controls on the emplacement of Escanaba Trough's three main volcanic centers, characterize the geometry of its spreading segments, and provide volumetric data on the distribution of sub-sediment volcanism in the southern Gorda Ridge.

Escanaba Trough, Pacific Ocean↗
USGS2025

Applying U-Th disequilbrium for dating siliceous sinters

Continental hydrothermal systems are critical avenues for the crustal transport of heat and mass captured for geothermal energy and mineral exploration. Thus, understanding their temporal evolution and longevity is important for resource characterization. Deposits of microlaminated siliceous sinter, common surface expressions of high temperature reservoirs (> 170 °C), have the potential to trace hydrothermal histories. Geothermal reservoirs are often located within uranium-bearing silicic volcanic rock where subsurface fluid-rock interactions extract U into hydrothermal fluids. U incorporated in the surface sinter deposit has the potential to provide a U— Th disequilibrium dating option. We focus on samples from El Tatio geyser field in the Altiplano of northern Chile, the largest geothermal system in the Andes. Our resulting 230 Th/U ages, along with the water and deposit elemental compositions, suggest concentrations of U and Th vary predictably along the sinter apron. While distal facies containing the highest U concentrations (> 50 μg/g) are least affected by detrital Th corrections, they can display suspected open-system behavior. In contrast, more medial facies, where bacterial mats and other porous textures are co mmonly concentrated, have only trace amounts of U (< 0.1 μg/g), which leads to unreliable or geologically improbable dates. Proximal facies tend to date most consistently. By comparing existing 14 C ages with 230 Th/U results, 230 Th/U ages tend to be younger than the 14 C ages, supporting the presence of a 14 C-dead carbon influence. New data confirm that the onset of geothermal activity at El Tatio goes back to the late Pleistocene.

Altiplano, El Tatio geyser field↗
USGS2026

Detecting snow avalanche activity using infrasound: Hooker Valley, New Zealand

Snow avalanches pose considerable hazards to people and infrastructure in alpine environments. Traditional avalanche monitoring relies on meteorological data and visual observations, which can be limited in scope and timeliness. Infrasound offers a promising complementary monitoring tool by detecting the low-frequency sound waves generated by avalanches. Here, we present infrasound and camera observations during a 50-day field campaign in the Hooker Valley of Aoraki/Mount Cook National Park, New Zealand. Our study detected seven avalanches with the cameras, whereas the infrasound system identified only one of these events, which was the largest and occurred under conditions that likely favoured infrasound propagation. The infrasound system recorded numerous other events not captured by the cameras, indicating the benefit of further investigation to determine their sources. These findings highlight the potential of infrasound technology for detecting avalanches and providing broad spatial coverage, capturing events in areas not monitored by cameras, while also showcasing limitations in infrasound capabilities. The limited detection of smaller avalanches underscores the opportunity for further research to enhance detection capabilities and understand environmental influences such as snow cover and wind noise. Overall, this study emphasises the utility of multidisciplinary monitoring techniques to improve avalanche detection in alpine environments.

Hooker Valley↗