SEARCH · USGS Science

Results for “Earth System Science Data”

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

USGS2025

Map of topographic lineaments interpreted as recent surface ruptures along the Bennett Valley and Southern Maacama Fault Zones, Sonoma County, California

This study documents lidar-illuminated topographic lineaments interpreted as evidence of recent surface fault ruptures and surface ruptures related to distributed deformation along the Bennett Valley Fault Zone and the southernmost Maacama Fault Zone in the northern San Francisco Bay area (fig. 1, on map sheet). Together, these fault zones form a structural connection across a right stepover between the main Maacama Fault Zone and the Rodgers Creek Fault, overlapping principal strands of the San Andreas plate boundary system north of San Francisco Bay (figs. 1, 2, on map sheet) that accommodate about a quarter of the ~40 millimeters per year of regional dextral tectonic slip (Parsons and others, 2013). Although much of the Bennett Valley Fault Zone is seismically active (McLaughlin and others, 2012; Sowers and others, 2010), only the north end of the fault zone (the “Spring Valley strand,” indicated on southern map panel of map sheet and labeled “D” in figure 2, on map sheet) was previously known to have youthful geomorphic expression and to displace Holocene deposits (McLaughlin and others, 2008; Sowers and others, 2010; Sowers and others, 2016). Holocene activity along the length of the Bennett Valley Fault Zone and southernmost Maacama Fault Zone, with possible implications for rupture propagation, continuity, and slip transfer, had not been identified. However, fault splays that project eastward from the Rodgers Creek Fault have been mapped and hypothesized to accommodate slip transfer to the Bennett Valley Fault Zone (Hecker and Randolph Loar, 2018). The generally subtle and distributed nature of surface-rupture evidence along the Bennett Valley and southernmost Maacama Fault Zones, and extensive vegetation cover, had left recent faulting previously unmapped along most of the zone. The map presented here represents a new compilation of inferred surface-rupture features detected using high-resolution topographic lidar data from an airborne lidar survey of Sonoma County, California (OpenTopography, 2014). These data, which enable subtle topographic features to be discerned, indicate that recent (likely Holocene) surface ruptures extend throughout the Bennett Valley and southernmost Maacama Fault Zones. The rupture-lineament map was created using a GIS (geographic information system) framework and is included herein as an image map at a scale of 1:36,000 and as digital datasets (included as supplemental information to this report). The mapping is intended to lay the groundwork for future studies designed to better characterize how plate-boundary slip is accommodated on this important and complex system of faults.

California↗

USGS2024

Cross section N–N' through the Valley and Ridge province of the southern Appalachian basin, from Greene County, west-central Alabama, to Bibb County, central Alabama

Introduction Geologic cross section N–N′ is the sixth in a series of geologic cross sections constructed by the U.S. Geological Survey to document and improve understanding of the geologic framework and petroleum systems of the Appalachian basin. Cross section N–N′ provides a regional view of the structural and stratigraphic framework of the Appalachian basin in the Valley and Ridge province in western and central Alabama; it spans approximately 69 miles (mi) (111 kilometers [km]). This geologic cross section is a companion to geologic cross sections E–E′ , D–D′ , C–C′ , I–I′ , and A–A′ that are located approximately 350 to 550 mi (563 to 885 km) to the northeast. Cross section N–N' complements earlier geologic cross sections through the Alabama part of the Appalachian basin. Although some of the other cross sections show more structural and stratigraphic detail, they are of more limited extent geographically and stratigraphically. Cross section N–N′ contains information that is useful for evaluating energy resources in the Appalachian basin. Although the Appalachian basin petroleum systems are not shown on the cross section, many of their key elements (such as source rocks, reservoir rocks, seals, and traps) can be inferred from lithologic units, unconformities, and geologic structures shown on the cross section. Other aspects of petroleum systems (such as the timing of petroleum generation and petroleum migration pathways) may be evaluated by burial history, thermal history, and fluid flow models based on what is shown on the cross section. In addition, cross section N–N′ may be used as a reconnaissance tool to identify plausible geologic structures and strata for the subsurface storage of liquid waste or for the sequestration of carbon dioxide.

Alabama↗

USGS2025

Variable partitioning of lithium in rhyolitic melt during decompression and ascent

The partitioning behavior of Li in magmatic systems is increasingly being investigated due to the economic importance of Li in the transition to sustainable energy resources (e.g., batteries). However, at upper crustal pressures, it remains uncertain whether Li preferentially partitions into the vapor or liquid (brine) phase or remains in the silicate melt. This complicates our ability to determine where Li resides—silicate melt, minerals, or fluid phase—upon eruption, a crucial factor for understanding its postdepositional movement and concentration into a brine or volcano-sedimentary deposit. Here, we present a novel investigation into the behavior of Li within natural evolved melts during continuous magma decompression and ascent using melt embayments (open melt inclusions). Mineral-hosted melt embayments preserve records of the evolving composition of the exterior melt, including degassing pathways and ascent timescales, when paired with appropriate diffusion coefficients. Lithium concentration profiles were measured in quartz-hosted melt embayments from the rapidly quenched eruptive phases of five rhyolitic, caldera-forming eruptions to investigate the behavior of Li during magma decompression and ascent, where vapor partitioning and ascent dynamics were previously established by investigating H 2 O and CO 2 profiles. We find that in four systems, embayments contain lower interior Li concentrations than the coerupted melt inclusions; the fifth system contains the same Li concentrations in embayments and melt inclusions. However, many of these embayments contain gradients, with 84 % preserving Li enrichment near the melt-bubble interface, as compared to their interior concentration. We interpret these characteristics to represent two distinct stages of Li partitioning during magma decompression and ascent, in contrast to existing literature that proposes only one type of partitioning behavior. The first stage is interpreted as melt depletion of Li, likely driven by partitioning into an exsolved supercritical fluid phase, supported by the strong correlation between the extent of Li depletion and Cl concentration in the melt, as well as the decompression rate. This behavior then fundamentally shifts, where Li reenriches in the melt, postulated to be driven by the unmixing of the supercritical fluid phase at shallow pressures. For the one system that did not develop Li gradients through decompression, we attribute this to the lower values of Na and Cl in the melt, potentially inhibiting the partitioning of Li into a fluid phase. Importantly, the behavior of Li during decompression is not consistent within or between volcanic centers, highlighting the need for systematic experimental investigation in variable composition melts at pressures relevant to conduit dynamics. This knowledge would improve our ability to model Li profiles to understand magma decompression, and predict where Li resides (e.g., stored in volcanic glass, gas, or crystals) upon eruption prior to any later extraction.

Economic Geology↗

USGS2025

Simulating human behavior under earthquake early warning

Earthquakes are a rapid-onset hazard where advance planning and learning plays a key role in mitigating injuries and death to individuals. Recent advances in earthquake detection have resulted in the development of earthquake early warning (EEW) systems. These systems can provide advance warning to predetermined geographic regions that an earthquake is in progress, which may result in individuals receiving warning seconds before significant shaking is felt at their location. This additional time could allow individuals to take more effective protective actions during the immediate disaster. To demonstrate this, we created an agent-based simulation of a basic apartment that allowed us to randomly and repeatedly simulate an individual receiving and responding to an EEW message. The results of our preliminary simulation show that, in our study environment, earthquake early warning alerts have the potential to allow for sufficient time for individuals to take protective actions.

Washington↗

USGS2025

Assessment of coastal and fluvial morphodynamic changes using Structure-for-Motion: A case study of the Sfȃntu Gheorghe Mouth (Danube Delta, Romania)

The ability to accurately map erosion, flooding, and habitat loss in coastal environments is crucial for formulating national strategies aimed at preventing and mitigating the impacts of natural disasters. A fundamental component of this process is the implementation of coastal morphodynamics monitoring through Structure-from-Motion (SfM) techniques, utilizing high-resolution 2D/3D data obtained from aerial photogrammetry. To assess morphodynamic changes over a three-year period (2022 – 2024), several SfM-based photogrammetric studies were conducted, each year, in the Romanian sector of the Danube-Black Sea coastal zone, specifically at the mouth of one of the Danube River distributaries (Sf Gheorghe branch) into the Black Sea, and along the left bank, near Sf Gheorghe locality, located within the Danube Delta Biosphere Reserve (DDBR). The essential equipment for aerial photogrammetry comprises Unmanned Aerial Vehicles (UAVs) and Global Navigation Satellite Systems (GNSS). In this study, the UAV used was a DJI Mavic 3T (Enterprise/Thermal) drone, complemented by two Trimble R12i and R4 GNSS systems, as well as approximately 10 Ground Control Points (GCPs). Data acquisition and processing were carried out using specialized photogrammetric software (Agisoft Metashape) along with various GIS tools (e.g., Blue Marble Geographics Global Mapper and ESRI ArcMap). The photogrammetric products generated for the study, as detailed in this paper, include Digital Elevation Models (DEMs), Digital Terrain Models (DTMs), orthomosaics (orthophotos), and others. At Sfântu Gheorghe beach, a comparison between 2023 and 2024 photogrammetric surveys revealed that the left bank of the Sf. Gheorghe Arm, at the river mouth into the Black Sea, suffered from a twist (erosion) of up to 64 metres. Additionally, on the selected perimetre (total area of 31,910 square meters ) from the beach and dune zone of Sf. Gheorghe, an area of up to 16,202 square meters was eroded between 2023 and 2024. This contrasts with the period between 2022 and 2023, during which deposition predominated. Erosion at the Danube mouths and the adjacent Black Sea coastline is driven by a complex interaction of natural and anthropogenic factors. Natural processes, including subsidence, sea-level rise, and episodic extreme storm events, contribute significantly to coastal dynamics. Meanwhile, human-induced factors, such as upstream hydrotechnical works that limits sediment transport, cutting of navigation canals, as well as the exacerbating effects of climate change, further accelerate erosion. The recent Structure-from-Motion (SfM) surveys provide essential quantitative data, enabling a detailed analysis of both short-term and long-term morphodynamic changes influenced by seasonal variations and extreme hydrometeorological events in this highly dynamic coastal system.

Sfȃntu Gheorghe Mouth (Danube Delta)↗

USGS1994

Hydrogeology of the Chickasaw National Recreation Area, Murray County, Oklahoma

The Travertine District (Park) of the Chickasaw National Recreation Area, operated and maintained by the National Park Service, is near the City of Sulphur in south-central Oklahoma. The Park was established in 1902 because of its unique hydrologic setting, which includes Rock Creek, Travertine Creek, numerous mineralized and freshwater springs, and a dense cover of riparian vegetation. Since the turn of the century several flowing artesian wells have been drilled within and adjacent to the Park. Discharge from many of these springs and the numbers of flowing wells have declined substantially during the past 86 years. To determine the cause of these declines, a better understanding of the hydrologic system must be obtained. The U.S. Geological Survey, in cooperation with the National Park Service, has appraised hydrologic information obtained for the Park from several studies conducted during 1902-87. The principal geologic units referred to in this report are the Arbuckle Group and the overlying Simpson Group. These rocks are of Upper Cambrian to Middle Ordovician age and are composed of dolomitic limestone, with some sandstones and shales in the Simpson Group. Surface geologic maps give a general understanding of the regional subsurface geology, but information about the subsurface geology within the Park is poor. The Simpson and Arbuckle aquifers are the principal aquifers in the study area. The two aquifers are not differentiated readily in some parts of the study area because of the similarity of the Simpson and Arbuckle rocks; thus, both water-bearing units are referred to frequently as the Arbuckle-Simpson aquifer. The aquifers are confined under the Park, but are unconfined east and south of the Park. Precipitation on the outcrop area of the Arbuckle aquifer northeast and east of the Park recharges the freshwater springs (Antelope and Buffalo Springs) near the east boundary of the Park. The source of water from mineralized springs located in the central part of the Park, and flowing wells within and north of the Park, is believed to be a mix of waters from rocks of the Arbuckle and Simpson Groups. The source of water from two highly mineralized springs, Bromide and Medicine, that ceased to flow in the early 1970’s is believed to be from the Simpson Group. Water-quality characteristics reflect the sources of ground water in the study area. The highly mineralized springs near the western end of the Park are a sodium chloride type with dissolved solids greater than 4,500 mg/L. The freshwater springs near the eastern end of the Park are a calcium bicarbonate type with total dissolved solids of less than 400 mg/L. Flow from the artesian wells has declined substantially during the past 86 years and the wells are estimated to currently discharge only about 10 percent of the total flow reported in 1939. The depletion is believed to be caused by a gradual lowering of the hydraulic head within the aquifer. The influence on the hydrologic system of local municipal and industrial pumping from the Arbuckle-Simpson aquifer is difficult to discern because the system is much more sensitive to precipitation than to pumpage. Ground-water levels and spring flows in this region respond rapidly to precipitation. The effects of withdrawals from the City of Sulphur and Oklahoma Gas and Electric Company power-plant water-well fields are not discernible at wells and springs. The hydrologic system may be influenced by pumping, particularly during extended dry periods of several years, but the impact of pumping on the system cannot be determined without further investigation.

Oklahoma↗

USGS2025

Preliminary depth to basement modeling at Salton Sea, California

The San Andreas Fault – Imperial Fault (SAF-IF) transtensional step-over zone along the southern margin of the Salton Sea hosts substantial geothermal production and lithium brine resources. Recent volcanism at the Salton Buttes and active seismicity along the SAFIF fault system highlight active tectonic and magmatic processes that pose natural hazards and may impact energy and mineral production. Characterizing the subsurface architecture and extent of concealed alteration associated with this tectono-magmatic system enhances understanding of these active processes, associated hazards, and resources. We have compiled a gravity database, consisting of new and re-processed existing data, from which we have constructed a new isostatic residual gravity anomaly map of the Salton trough. We have used this new gravity dataset together with a compilation of publicly available borehole data to develop new depth to basement inversion models for the region. These depth to basement models help to constrain basin geometries, inform alteration mapping, and reveal variations in basement rocks. Due to the concealed nature of the complex tectonic framework at the Salton trough, it is necessary to utilize geophysical methods for subsurface characterization. These new depth to basement models are a first step toward constructing 2D and 3D geophysical and geologic models of the Imperial Valley and Salton Sea geothermal area. This analysis complements other geophysical initiatives, including magnetotelluric (MT) modeling (Tokmakoff et al., 2024), magnetic mapping (Glen and Earney, 2023, 2024) and potential field modeling, and seismic studies focused on hazard and resource investigations in the Imperial Valley.

California↗

USGS2024

Event classification, seismicity, and eruption forecasting at Great Sitkin Volcano, Alaska: 1999–2023

The frequency content of volcanogenic seismicity is often used to classify events and their spatial and temporal progression is then used to map subsurface volcanic processes. The progression of volcano-seismic events and associated source processes also plays a critical role in eruption forecasting. Here we develop and evaluate a computerized methodology for characterizing volcano-seismic event types using Frequency Index and Average Peak Frequency. We apply and test this technique at Great Sitkin Volcano, Alaska, classifying over 9000 hypocenters between 1999 and 2023. This 24-year time span covers periods of seismic quiescence, earthquake activity on nearby tectonic (bookshelf) faults, precursory unrest from 2016 to 2021, and the explosive onset in May 2021 of the ongoing effusive eruption. We use the spatial and temporal evolution of classified event types to map the active volcanic and tectonic processes, develop a conceptual model of the subsurface magmatic system, and perform a retrospective analysis of eruption forecasts at Great Sitkin Volcano between 2016 and the present. The classification and progression of hypocenters suggests the subsurface Great Sitkin Volcano magmatic system consists of a mid- to lower- crustal source zone between 10 and 40 km depth and an upper crustal magma storage area between −1 and 10 km depth (hypocenter depth is referenced to sea level and negative depths reflect height above sea level). The earliest precursors occurred in July 2016 and consisted of deep long-period and volcano-tectonic earthquakes at mid-crustal depths suggesting the subsequent unrest and eruption were triggered by a deeper intrusion of magma. This mid-crustal seismic activity was immediately followed by the onset upper-crustal long-period events and volcano-tectonic earthquakes VTs suggesting a strong linkage between the shallow and deeper portions of the magmatic system. The upper crustal area was likely capped by the 1974 lava dome until the magmatic explosion on May 26, 2021.

Alaska↗

USGS2026

Lunar analog study using portable gamma-ray neutron detector: Radiochemical mapping of silicic-to-basaltic volcanic terrains in the San Francisco Volcanic Field, Arizona

Studying lunar silicic volcanism provides key insights into the Moon’s crustal evolution, magmatic processes, and volcanic history. The Lunar-VISE (Lunar Vulkan Imaging and Spectroscopy Explorer) mission will investigate the Gruithuisen domes, a unique lunar region hypothesized to have formed through silicic volcanism. Using instruments on a Firefly Aerospace lander and a Honeybee Robotics rover, Lunar-VISE will analyze mineralogy, geochemistry, and surface properties to determine the origin and evolution of the domes, with a gamma-ray and neutron spectrometer (LV-GRNS) among its payload instruments. In preparation for this mission, we conducted preliminary fieldwork using a handheld gamma-ray neutron detector with NaI(Tl) and Cs LiYCl :Ce (CLYC) scintillators. This study focuses on various rhyolitic and basaltic volcanic centers in the San Francisco Volcanic Field (SFVF) near Flagstaff, Arizona—specifically Sugarloaf Peak, Bonito Lava Flow, and Robinson Mountain, a region containing several well-characterized lunar analog sites. The SFVF was selected for this study due to its broad compositional diversity spanning basaltic to rhyolitic compositions, its well-preserved volcanic morphologies, and its extensive use in previous NASA field campaigns and astronaut training exercises, making it an ideal terrestrial laboratory for testing planetary exploration techniques. We measured natural radioactivity, specifically from potassium (K), thorium (Th), and uranium (U) and other elements in their decay chains, which serve as diagnostic tracers of magmatic differentiation and crustal evolution processes, to assess geochemical variability across compositionally diverse terrains. The detector’s sensitivity was assessed across varying concentration levels. Regionally averaged concentrations of radioisotopes were determined by gamma-ray spectroscopy at selected sites. Our measurements reveal significant compositional variations between sites, with Sugarloaf Peak (rhyolitic, silicic dome) exhibiting the highest average radioisotope concentrations (K: 3.29 ± 0.24 wt%, U: 14.35 ± 1.81 ppm, Th: 27.14 ± 2.43 ppm), while Bonito Lava Flow (K: 1.13 ± 0.08 wt%, U: 3.86 ± 0.52 ppm, Th: 7.52 ± 1.02 ppm), and Robinson Mountain (K: 1.46 ± 0.16 wt%, U: 5.59 ± 1.06 ppm, Th: 11.75 ± 1.98 ppm) show lower concentrations aligned with basaltic compositions. Gamma-ray fluxes were elevated by a factor of approximately three to four at Sugarloaf Peak relative to nearby basaltic terrains, consistent with expected geochemical differentiation patterns. Furthermore, at meter scales, proximity to geological features significantly affects measurements. Th concentrations adjacent to a cliff face at Sugarloaf’s base were 28% higher than values measured 3-5 meters away from the same feature, demonstrating localized compositional heterogeneity. Analysis of station-by-station measurements reveals that K, U, and Th concentrations show an elevation trend at Robinson Mountain and generally higher concentrations near Sugarloaf Peak’s summit, suggesting progressive magmatic differentiation and/or the presence of more evolved lithologies at higher elevations. By mapping these elements on Earth using GRNS, we aim to optimize measurement techniques for Lunar-VISE and similar rover-borne missions by demonstrating the applicability of portable GRNS instruments for planetary surface exploration in an analog environment.

Arizona↗

USGS2024

Surficial geology and Quaternary fault map of the Las Vegas Valley, Clark County, Nevada

The Surficial Geology and Quaternary Fault Map of the Las Vegas Valley, Clark County, Nevada is a 1:50,000-scale compilation of published 1:24,000-scale geologic maps integrated with new field and desktop mapping. This geologic map compilation and GIS database are part of a broader study on the Quaternary faults in the Las Vegas Valley, the results of which are summarized in this report. The map compilation utilized a variety of digital base maps to evaluate the characteristics of surficial deposits and Quaternary fault scarps including lidar data in undeveloped parts of the valley and orthorectified historical aerial photos in the urbanized center of the valley. The map distinguishes twenty surficial map units including alluvial-fan and ephemeral-wash deposits, groundwater discharge deposits of the Las Vegas Formation, and alluvium deposits of Las Vegas Wash. Quaternary faults include the Quaternary Las Vegas Valley fault system in the center of the valley and the Frenchman Mountain fault system along the eastern margin. All Quaternary fault traces in the valley were evaluated, and the previously published fault mapping was modified as needed. The new mapping and uniform correlation of surficial deposits displaced by Quaternary faults yields a consistent valley-wide characterization of the recency of fault activity. The map and GIS database include a Quaternary geochronology compilation of 119 previously published surficial deposit ages of various types, 35 new luminescence ages, and 3 new radiocarbon ages. The new age data further bracket the ages of the Las Vegas basin Quaternary stratigraphy and provide new constraints on the timing of Quaternary fault activity.

Nevada↗

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↗

USGS2024

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

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

Idaho↗

USGS2023

Cursed? Why one does not simply add new data sets to supervised geothermal machine learning models

Recent advances in machine learning (ML) identifying areas favorable to hydrothermal systems indicate that the resolution of feature data remains a subject of necessary improvement before ML can reliably produce better models. Herein, we consider the value of adding new features or replacing other, low-value features with new input features in existing ML pipelines. Our previous work identified stress and seismicity as having less value than the other feature types (i.e., heat flow, distance to faults, and distance to magmatic activity) for the 2008 USGS hydrothermal energy assessment; hence, a fundamental question regards if the addition of new but partially correlated features will improve resulting models for hydrothermal favorability. Therefore, we add new maps for shear strain rate and dilation strain rate to fit logistic regression and XGBoost models, resulting in new 7-feature models that are compared to the old 5-feature models. Because these new features share a degree of correlation with the original relatively uninformative stress and seismicity features, we also consider replacement of the two lower-value features with the two new features, creating new 5-feature models. Adding the new features improves the predictive skill of the new 7-feature model over that of the old 5-feature model; albeit, that improvement is not statistically significant because the new features are correlated with the old features and, consequently, the new features do not present considerable new information. However, the new 5-feature XGBoost model has a statistically significant increase in predictive skill for known positives over the old 5-feature model at p = 0.06. This improved performance is due to the lower-dimensional feature space of the former than that of the latter. In higher-dimensional feature space, relationships between features and the presence or absence of hydrothermal systems are harder to discern (i.e., the 7-feature model likely suffers from the “curse of dimensionality”).

Geothermal Resources Council Transactions↗

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

Linking stream-reach nitrogen loads and groundwater “reachsheds” to inform wastewater-nitrogen management actions, Cape Cod, Massachusetts

Study Region Cape Cod, Massachusetts, U.S.A. Study Focus Anthropogenic nitrogen (N) is a key factor in degrading groundwater and surface-water quality, particularly in coastal New England where onsite wastewater systems are prevalent. This study evaluated whether direct N-load measurements in streams on Cape Cod, Massachusetts, coupled with flow-path information from groundwater-flow models, can effectively identify potential land areas where nitrogen mitigation could substantially reduce loads to receiving waters. Nitrogen fluxes were measured along 63 stream reaches during winter and summer and paired with simulated groundwater recharge areas to identify and rank potential areas for reduction of nonpoint-source N inputs. New Hydrologic Insights for the Region Reach-scale nitrate-N loads ranged from −39.1–1182 kg-N/yr per 100 m of stream, indicating spatially variable groundwater inputs across seasons. “Reachsheds” — areas contributing groundwater recharge to specific stream reaches — were delineated using a regional groundwater-flow model. Strong correlations were found between observed N loads and land-use characteristics, especially the number of septic systems and total N inputs from the sum of considered sources. Observed N loads were moderately correlated with recharge area size and wastewater flow estimates. Correlating reach-specific groundwater N loads with land use and parcel-scale nitrogen-yield data identified reachsheds with the highest potential for N load reduction. This approach enables targeted implementation of restoration efforts to optimize nutrient management and support regional load reduction.

Massachusetts↗

USGS2024

Volcanic gas monitoring

Introduction As magma rises through the crust, decreasing pressure conditions allow volatiles to exsolve from the magma. These volatiles then migrate upward through the crust, where they can be stored at shallower levels or escape to the atmosphere. Rising magma also heats rock masses beneath volcanic centers, causing water in shallow aquifers and hydrothermal systems to boil and release additional gases and steam (see chapter F , this volume; Ingebritsen and Hurwitz, 2024). The chemistry and quantity of gases that reach the surface during periods of quiescence or volcanic unrest can reveal that gas-rich magma is ascending, crystallizing, or alternatively stalling, with important implications for volcanic hazard (for example, Sutton and others, 1992; Aiuppa and others, 2007, 2021; Werner and others, 2009, 2011, 2012; Moretti and others, 2013; de Moor and others, 2016; Lewicki and others, 2019; Edmonds and others, 2022; Kern and others, 2022; Kunrat and others, 2022). Most volcanoes in Alaska and the western United States are characterized by weak degassing, with one or more low-temperature fumaroles (typically near the local boiling temperature of water) and connect to a deeper and sometimes extensive hydrothermal system (for example, McGee and others, 2001; Symonds and others, 2003a, b). Hydrothermal systems will affect the chemistry of rising gases exsolved from deeper magma (Symonds and others, 2001), including sulfur dioxide (SO 2 ), hydrogen chloride (HCl), and water vapor (for example, Doukas and Gerlach, 1995; Gerlach and others, 1998, 2008; Symonds and others, 2001; Werner and others, 2013). As an example, depending on factors such as temperature, pressure, and oxidation state, rising SO 2 will react with groundwater to form hydrogen sulfide (H 2 S) gas, dissolved sulfate (SO 4 2− ), or elemental sulfur (Christenson, 2000; Symonds and others, 2001; Werner and others, 2008). The reaction and dissolution of SO 2 into shallow groundwater is commonly referred to as scrubbing, and can reduce the likelihood that ascending, degassing magma can be detected. Carbon dioxide, however, in addition to exsolving from magma early in the ascent process, is not easily removed by hydrothermal fluids (Lowenstern, 2001). As scrubbing and other processes take place, the SO 2 /H 2 S, CO 2 /SO 2 , and CO 2 /H 2 S ratios may change. High rates of SO 2 emission indicate that magma has moved to relatively shallow levels in the volcano and that the system has heated up enough to establish dry pathways from depth to the surface. Monitoring multiple gas species and the total output of those species is thereby useful for volcano monitoring during both periods of quiescence, to establish background degassing conditions, and during unrest, when gas geochemistry and emission rates can provide information on changing conditions, such as magma ascent. To provide context for multidisciplinary volcano forecasts, we focus on the following two key required capabilities: (1) characterizing baseline geochemistry and gas discharge from volcanoes and volcanic regions and (2) monitoring changes in gas geochemistry and discharge to inform forecasts of volcanic eruptions and their effects. Sufficient baseline data must be collected to identify and interpret anomalous degassing associated with volcanic unrest (for example, Sorey and others, 1998; Rouwet and others, 2014). Differences in volcano type, baseline degassing rates, local hydrology, and geography (for example, high versus low latitude) will result in a different baseline for each volcano. Volcanoes of any threat level that exhibit one or more degassing phenomena would ideally be monitored by techniques needed to establish baseline degassing data, with the sampling frequency of baseline data dictated by the threat level (table E1). Additional monitoring techniques become necessary during periods of unrest. In general, three of the most important techniques for gas monitoring are (1) direct sampling of fumarole, spring, and soil gases for laboratory geochemical measurements, (2) measurements of the chemical composition of the volcanic plume and emission rates of major gas species (for example, H 2 O, CO 2 , SO 2 , and H 2 S) by satellite, airborne, or ground-based techniques, and (3) measurements of diffuse emissions of CO 2 and other gases through soils. Various methods and instruments may be useful both for baseline studies and during unrest.

Scientific Investigations 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↗

USGS2024

Near-real-time earthquake-induced fatality estimation using crowdsourced data and few-shot large-language models

When a damaging earthquake occurs, immediate information about casualties (e.g., fatalities and injuries) is critical for time-sensitive decision-making by emergency response and aid agencies in the first hours and days. Systems such as the Prompt Assessment of Global Earthquakes for Response (PAGER) by the U.S. Geological Survey (USGS) were developed to provide a forecast of such impacts within about 30 min of any significant earthquake globally. However, existing disaster-induced human loss estimation systems often rely on early casualty reports manually retrieved from global traditional media, which are labor-intensive, time-consuming, and have significant time latencies. Recent approaches use keyword matching and topic modeling to identify human casualty-relevant information from social media but tend to be error-prone when dealing with complex semantics in multi-lingual text data and parsing dynamically changing and conflicting human death and injury numbers shared by various unvetted sources in social media platforms. In this work, we introduce an end-to-end framework to significantly improve the timeliness and accuracy of global earthquake-induced human loss forecasting using multi-lingual, crowdsourced social media. Our framework integrates (i) a hierarchical casualty extraction model built upon large language models, prompt design, and few-shot learning to retrieve quantitative human loss claims from social media, (ii) a physical constraint-aware, dynamic-truth discovery model that discovers the truthful human loss from massive noisy and potentially conflicting human loss claims, and (iii) a Bayesian updating loss projection model that dynamically updates the final loss estimation using discovered truths. We test the framework in real-time on a series of global earthquake events in 2021 and 2022 and show that our framework effectively automates the retrieval of casualty information faster but with comparable accuracy to those now retrieved manually by the USGS. The code associated with this work is made available at: https://github.com/SusuXu-s-Lab/Hierarchical-Earthquake-Casualty-Information-Retrieval

International Journal of Disaster Risk Reduction↗