SEARCH · USGS Science

Results for “Rapid Communications in Mass Spectrometry”

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

USGS2025

Chemical and isotopic fractionation during melt inclusion formation

Melt inclusions are used to study the origin and evolution of magmas. The extent to which they represent equilibrium melt compositions, however, critically hinges on the ratio of crystal growth rate to diffusion rate in melts. If the rate of crystal growth is limited by the supply of nutrients and the dissipation of unwanted components, the trapped melt will be depleted in compatible elements and enriched in incompatible elements. Despite widespread recognition of the potential for melt inclusions to trap disequilibrium compositions, their identification remains a challenge. Here, we demonstrate how stable isotopes provide a solution to this problem. Melt inclusions in basaltic tephra from Kīlauea volcano, Hawai‘i, display Mg isotope fractionation up to 0.49‰ ( 26 Mg/ 24 Mg) relative to bulk rock. These observations indicate that kinetic processes, such as chemical diffusion, are at play, as equilibrium mineral-melt isotope fractionation at high temperatures is an order of magnitude smaller. We propose that the heavy Mg isotopic compositions of the melt inclusions were generated by the faster diffusion and thus preferential incorporation of 24 Mg by the growing olivine, leaving the slower diffusant ( 26 Mg) behind. The incompatible elements (e.g., Ca, Ti, and K) are predicted to display similarly large fractionations when boundary layer effects are significant. These findings show that stable isotopes can be a useful tool in identifying disequilibrium compositions in melt inclusions.

Geology↗

USGS2026

Magnetic storms and geoelectric hazards

Magnetic storms induce geoelectric fields at Earth's surface that can interfere with grounded long-line systems. The September 1859 storm disrupted global telegraph operations, the March 1989 storm caused a blackout in Canada and interfered with electric-power-transmission systems in the United States, and other storms have had related impacts. The geographic and temporal dependence of geoelectric fields are functions of both geomagnetic variation and local surface impedance, which differ considerably across different geological regions. These dependencies can be mapped across the contiguous United States by combining magnetotelluric impedance tensors with ground magnetometer time series. This review illustrates such mapping for the 1989 storm and shows that power-system interference was experienced where surface impedance is high, and when and where geoelectric fields were intense. Statistical analyses indicate that storms comparable to that of March 1989 occur roughly once every four solar cycles. Ongoing developments in numerical modeling and real-time monitoring are anticipated to enable prediction of geoelectric hazards. ▪ Magnetic storms can induced electric fields in the solid Earth that interfere with electric-power-transmission systems. ▪ Geoelectric hazards depend on the storm-time geomagnetic disturbance and the electrical conductivity structure of Earth. ▪ Historically, impacts on telecommunication and power-transmission systems in the United States have been concentrated in the East and Midwest. ▪ The future occurrence of a magnetic superstorm could cause widespread disruption of electric-power-transmission systems.

Annual Review of Earth and Planetary Sciences↗

USGS2025

Hydrogeologic framework and conceptual model of the Red River alluvial aquifer east of Lake Texoma, southeastern Oklahoma, 1980–2022

The 1973 Oklahoma Groundwater Law (Oklahoma Statutes §82-1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin. At present (2025), the Oklahoma Water Resources Board has not established a maximum annual yield for the Red River alluvial aquifer east of Lake Texoma. To support the evaluation and determination of a maximum annual yield, a hydrogeologic framework and conceptual groundwater-flow model were developed to assess groundwater availability in the Red River alluvial aquifer east of Lake Texoma. The scope of this hydrologic investigation is the alluvium and terrace containing the Red River alluvial aquifer in Oklahoma between Lake Texoma, the Texas State line, and the Arkansas State line, an extent referred to in this report as “the eastern part of the Red River alluvial aquifer.” Parts of the alluvium and terrace extent in Arkansas and Texas are included in some analyses to address hydrologic influences from outside the aquifer’s boundaries in Oklahoma. The eastern part of the Red River alluvial aquifer in southeastern Oklahoma consists of approximately 401,280 acres of Quaternary alluvium and terrace deposits associated with the Red River and its major tributaries. Mean annual recharge to the aquifer for the 1980–2022 study period was estimated to be 8.62 inches per year, or 17.98 percent of the mean annual precipitation over the same period (47.94 inches). This mean annual recharge rate is equivalent to an inflow of approximately 288,250 acre-feet per year for the eastern part of the Red River alluvial aquifer. Recharge estimated using the Soil-Water-Balance code accounts for 98.7 percent of the conceptual-model inflows to the eastern part of the Red River alluvial aquifer. Saturated-zone evapotranspiration accounts for 11.9 percent and net streambed seepage accounts for 87.4 percent of the outflows in the conceptual model.

Arkansas, Oklahoma, Texas↗

USGS2024

The impact of source time function complexity on stress drop estimates

Earthquake stress drop—a key parameter for describing the energetics of earthquake rupture—can be estimated in several different, but theoretically equivalent, ways. However, independent estimates for the same earthquakes sometimes differ significantly. We find that earthquake source complexity plays a significant role in why theoretically (for simple rupture models) equivalent methods produce different estimates. We apply time‐ and frequency‐domain methods to estimate stress drops for real earthquakes in the SCARDEC (Seismic source ChAracteristics Retrieved from DEConvolving teleseismic body waves, Vallée and Douet, 2016 ) source time function (STF) database and analyze how rupture complexity drives stress‐drop estimate discrepancies. Specifically, we identify two complexity metrics—Brune relative energy (BRE) and spectral decay—that parameterize an earthquake’s complexity relative to the standard Brune model and strongly correlate with the estimate discrepancies. We find that the observed systematic magnitude–stress‐drop trends may reflect underlying changes in STF complexity, not necessarily trends in actual stress drop. Both the decay and BRE parameters vary systematically with magnitude, but whether this magnitude–complexity relationship is real remains unresolved.

Bulletin of the Seismological Society of America↗

USGS2026

Earthquake catalog for the Fairbanks region of central Alaska, 2014–2024, based on waveform cross-correlation

The Fairbanks region of central Alaska is part of a broad zone of intraplate crustal deformation, situated north of the Denali fault and north of the ongoing collision and flat‐slab subduction of the Yakutat oceanic plateau. Seismicity in the Fairbanks region occurs both in diffuse areas as well as in well‐defined lineaments, such as the left‐lateral Salcha fault, which hosted the 1937 M 8 7.3 earthquake. Starting with the regional seismicity catalog, we perform waveform cross‐correlation, network‐matched filtering, and relative relocation to obtain an enhanced seismicity catalog over the time period 2014–2024. Based on the relocated catalog, we interpret a set of 15 fault segments, including two conjugate faults and two new faults east of the previously documented fault system. Considering the combined seismicity in the Minto and Fairbanks regions, the median depth of seismicity decreases from east (6 km) to west (20 km). Our interpreted faults provide guidance for future tectonic modeling and assessment of seismic hazards in this region.

Alaska↗

USGS2025

Using public participatory geographic information systems (PPGIS) to explore uses and values for Mojave Trails National Monument, California

Many people ascribe a variety of values to public lands and waters, but some values are more difficult to assess and quantify than others. Public participatory geographic information systems (PPGIS) are tools that have been used to help quantify and map the public’s diverse values for a landscape. This work describes the first known Office of Management and Budget–approved use of PPGIS by a Department of the Interior bureau. The U.S. Geological Survey developed an internet-based application to aid in gathering PPGIS data, called Values Mapping for Planning in Regional Ecosystems (VaMPIRE). Further, this work describes the first pilot of the VaMPIRE application in coordination with the Bureau of Land Management to collect spatial data and other survey data regarding the public’s uses of and values for locations within Mojave Trails National Monument. We emailed the link to the VaMPIRE application to an interested party email list in 2024 with 207 valid emails and received 74 responses; we also received 47 responses from members of an off-roading social media group. Of the list of 16 value options, recreation was the most popular value for the monument, followed by wilderness and inspirational. Over 1,000 points were placed throughout the monument, indicating locations people use or value, with the locations spread throughout the entire monument. Additionally, most survey respondents stated their ability to receive benefits in locations they mapped would not change in response to a hypothetical scenario related to recreational facility development. This report describes exploratory results from the first use of the VaMPIRE tool in Mojave Trails National Monument and includes reflections on how the process went and considerations for future use of VaMPIRE.

California↗

USGS2025

Forecast, monitor, adapt: A multi-agency strategy to protect people from postfire debris flows

In 2020, a wildfire burned across Glenwood Canyon in Colorado, USA. A history of postfire debris flows in the region and a hazard assessment for the burn area indicated that potentially life-threatening debris flows could be triggered by rainfall within months of a wildfire. As a result, four government agencies evaluated strategies to help mitigate hazards, including the loss of human life, that may be associated with debris-flow events. After the fire, 26 large debris flows occurred in the summer of 2021 and three sediment-laden flows occurred in the summer of 2023, but there were no major injuries or fatalities reported. We found that integrating hazard assessment/ forecasting, monitoring, and adaptation scenarios was a successful strategy for reducing postfire debris-flow risks to human life (including injuries and fatalities). Weather forecasts and estimates of debris-flow triggering rainfall thresholds, likelihood, and volume were used to anticipate the timing, location, and magnitude of debris-flow events. Rainfall monitoring and detailed recordkeeping of storms that triggered debris flows were used to validate and update debris-flow warning thresholds that varied with time following the wildfire. Although the governmental agencies working in this burn area had distinct and differing agency mandates, they were able to integrate information to reduce the risk of debris-flow events to human life.

Colorado↗

USGS2025

Using gridded seismicity to forecast the long-term spatial distribution of earthquakes for the 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

Gridded (or background) seismicity models are a critical component of probabilistic seismic hazard assessments, accounting for off‐fault and smaller‐magnitude earthquakes. They are typically developed by declustering and spatially smoothing an earthquake catalog to estimate a long‐term seismicity rate that can be used to forecast future earthquakes. Here, we present new gridded seismicity models for use in the 2025 National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The previous PRVI NSHM was released in 2003, and our new models incorporate updates to both data and methodology. We utilize an updated earthquake catalog based on improved Puerto Rico Seismic Network data with newly characterized completeness epochs. The catalog is divided into crustal, subduction interface, and intraslab seismicity using new methods and Slab2 subduction zone geometries. To forecast the long‐term spatial distribution of earthquakes, we use an updated methodology developed for the 2023 U.S. 50‐state NSHM, considering three declustering methods and two spatial smoothing methods based on 2D Gaussian kernels. To adapt it for the complex seismotectonics of the region, we also adopt probabilistic methods to account for events with unknown depths and uncertainties in tectonic classification, and develop a new method for spatial scaling to counteract the effects of spatial variability in network coverage while maintaining the use of smaller events. Finally, we test the performance of these spatial models in forecasting the location of M w ≥ 5earthquakes in the region. Our updated methodology improves the representation of epistemic uncertainty relative to the 2003 model, and our results demonstrate the effectiveness of the new measures we have introduced to address heterogeneities in network detection and systematically evaluate forecast performance.

Puerto Rico, U.S. Virgin Islands↗

USGS2025

Geologic map of the Greater Antilles and Virgin Islands

Introduction The geologic map of the Greater Antilles and Virgin Islands is a compilation of information from the literature, integrated to provide a seamless geologic map of the region. This map was prepared to serve as a base map for a mineral resource assessment of the region. Several small-scale regional geologic maps of the region have been prepared in the past. This report supersedes an earlier version of the geologic map of the Greater Antilles and the Virgin Islands that was released as U.S. Geological Survey (USGS) Open-File Report 2019–1036. For this report, the regional geologic overview shown on sheet 1 covers Cuba; the island of Hispaniola, which includes Haiti and Dominican Republic; Jamaica; the Cayman Islands; Puerto Rico; and the U.S. and British Virgin Islands. Sheet 2 shows the geology of Cuba and Cayman Islands, sheet 3 the geology of Jamaica, sheet 4 the geology of Hispaniola, and sheet 5 the geology of Puerto Rico and Virgin Islands. Accompanying the maps are three tables: table 1 lists the map units (in alphabetical order by map-unit label), their ages, and the countries or territories in which they can be found. These are listed in alphabetical order for ease in locating unit names from the symbols on the maps. Table 2 is a list of sources for the geologic map compilation of Puerto Rico. Table 3 shows formations in the Devil’s Race Course Group of Jamaica.

Greater Antilles, Virgin Islands↗

USGS2026

GST-1: A high-resolution global sediment thickness model

Global Sediment Thickness 1 (GST-1) is a high-resolution sedimentary thickness model calculated on a 0.125° x 0.125° grid. It modifies the sediment thickness of the 1° x 1° Earth Crustal Model 1 (ECM1) by means of 3D inversions of free air gravity anomalies. GST-1 is calculated by performing structural inversions on high-density contrasts across two crustal boundaries: the sediment – basement interface and the crystalline crust – upper mantle interface. The inversions are calculated in each of ten overlapping 3D models that span the globe. These ten models are merged to obtain the GST-1 global model, providing an eight-fold increase in lateral spatial resolution in comparison with ECM1 and CRUST 1.0. Our sediment thickness model exploits the nearly continuous sampling of gravity data when compared to the irregular, sparse sampling of seismic refraction data. Sediment thickness values in GST-1 are in excellent agreement with independently derived cross sections from well-studied sedimentary basins, and within expected resolution limits of seismic refraction data. GST-1 offers a robust, high resolution global model of sedimentary thickness to support studies of sedimentary basins.

Tectonophysics↗

USGS2025

Geologic map of MTM −15032 and −20032 quadrangles, western Ladon basin, Mars

Our mapping efforts focused on the geomorphic terrains and geologic units contained within Mars Transverse Mercator (MTM) −15032 and −20032 quadrangles. These two quadrangles are located along the west side of Ladon basin and span lat −12.5° N. to −22.5° N. and long 325° E. to 330° E. The western part of Ladon basin and its bounding basin ring structures to the west preserved features that help to understand the long history of drainage across the Margaritifer Terra region of Mars. Our geologic map provides new insight into the extent of aqueous activity across western Ladon basin and the adjoining highlands, which includes the northern part of Ladon Valles, Arda Valles, numerous small valleys in the western highlands, and light-toned layered deposits associated with these fluvial features. The Mars Odyssey Thermal Emission Imaging System (THEMIS) infrared (IR) daytime mosaic (100 meters per pixel) was used as the primary base map. We constructed the geologic map of western Ladon basin at 1:1,000,000 scale. We identified 20 geologic units in the map area, which we divided into the following groups: crater units; volcanic units; chaotic units; basin fill units; crater, valley, and channel units; and plateau and highlands units. We observed the following sequence of events: (1) formation of Ladon and Holden impact basins in the middle Noachian, producing a mountainous unit; (2) landscape degradation and infilling during the Late Noachian to Late Hesperian forming terra units; (3) formation of Ladon Valles in the Late Noachian to Early Hesperian by catastrophic flooding, producing channel units; (4) accumulation of sediments in the Late Noachian to Early Amazonian, forming basin fill units; (5) formation of a smooth crater fill unit during the Early to Late Hesperian; (6) deposition during the Late Hesperian to Early Amazonian, creating light-toned layered units within Ladon Valles, Ladon basin, and other smaller valley networks along the western uplands, as well as formation of a light-toned unit along some crater floors; (7) formation a of chaotic unit in the Early to Middle Amazonian, formation of alluvial fans in the Early Amazonian, and eruption of a volcanic unit in the Middle Amazonian; and (8) formation of craters throughout the geologic history of the map region.

Scientific Investigations Map↗

USGS2024

Are the horizontal-to-vertical spectral ratios of earthquakes and microtremors the same?

We consider the similarities and differences between earthquake and microtremor horizontal‐to‐vertical spectral ratios (eHVSR and mHVSR, respectively) using a dataset of 161 sites in southern California. Quantitative comparisons are made in terms of the eHVSR and mHVSR lognormal median curves, as well as the frequencies and amplitudes associated with the fundamental‐ and higher‐mode resonances where present. The results show only 58% of the eHVSR–mHVSR pairs agree in terms of their median curve and only 25% of the eHVSR–mHVSR pairs agree in terms of shared resonances, which increases to 68% if flat HVSRs are considered equivalent. Furthermore, while the shared resonances match very well in terms of frequency (root mean square error, RMSE, <0.11 Hz), the amplitudes of those resonances do not agree (RMSE >1.6). These findings demonstrate that while eHVSR and mHVSR agree at some sites, they are not equivalent at all sites. To investigate if the agreement between eHVSR and mHVSR could be related to features of the microtremor data, earthquake recordings, and/or the site conditions, three machine learning (ML) models at varying levels of interpretability are presented. The ML models—which include multivariate logistic regression, gradient‐boosted trees, and support vector machines—show only partial success at using site‐specific data to predict whether eHVSR and mHVSR will likely agree in terms of their median curve (accuracy of 78%) and number of resonances (accuracy of 84%). Therefore, we conclude that while eHVSR and mHVSR can be quite similar in terms of resonant frequencies at some sites, they are not identical at all sites. Furthermore, preliminary evidence shows that the agreement of eHVSR and mHVSR can be predicted a priori given features of the microtremor measurements, earthquake recordings, and site conditions, although a larger dataset will be necessary for developing a robust predictive model.

California↗

USGS2024

Slope Unit Maker (SUMak): An efficient and parameter-free algorithm for delineating slope units to improve landslide modeling

Slope units are terrain partitions bounded by drainage and divide lines. In landslide modeling, including susceptibility modeling and event-specific modeling of landslide occurrence, slope units provide several advantages over gridded units, such as better capturing terrain geometry, improved incorporation of geospatial landslide-occurrence data in different formats (e.g., point and polygon), and better accommodating the varying data accuracy and precision in landslide inventories. However, the use of slope units in regional ( > 100 km 2 ) landslide studies remains limited due, in part, to the large computational costs and/or poor reproducibility with current delineation methods. We introduce a computationally efficient algorithm for the parameter-free delineation of slope units that leverages tools from within TauDEM and GRASS, using an R interface. The algorithm uses geomorphic laws to define the appropriate scaling of the slope units representative of hillslope processes, avoiding the often ambiguous determination of slope unit size. We then demonstrate how slope units enable more robust regional-scale landslide susceptibility and event-specific landslide occurrence maps.

Natural Hazards and Earth Systems Sciences (NHESS)↗

USGS2025

Preface to the focus section on earthquake-induced ground deformation

Ground failure due to earthquake-induced liquefaction, surface fault rupture, and landsliding has caused substantial damage in the past and continues to pose risks to critical infrastructure in the future. The consideration of ground failure is central to engineering analysis and design, ranging from site-specific engineering assessments of ground failure potential to broader regional-scale ground failure assessments of near-real-time earthquake impacts or longer-term portfolio risk. Uncertainties in all aspects of the problem complicate the prediction of earthquake-induced ground failure, including the identification of hazard-susceptible geomaterials, characterizing their spatial extent and physical properties, spatio-temporal variability in groundwater conditions, characterization of earthquake ground motions, assessing ground failure severity, and linking ground failure severity to infrastructure damage and loss.

Seismological Research Letters↗

USGS2025

Pathways for potential exposure to onshore oil and gas wastewater: What we need to know to protect human health

Produced water is a chemically complex waste stream generated during oil and gas development. Roughly four trillion liters were generated onshore in the United States in 2021 (ALL Consulting, 2022, https://www.gwpc.org/wp-content/uploads/2021/09/2021_Produced_Water_Volumes.pdf ). Efforts are underway to expand historic uses of produced water to offset freshwater needs in water-stressed regions, avoid induced seismic activity associated with its disposal, and extract commodities. Understanding the potential exposures from current and proposed produced water uses and management practices can help to inform health-protective practices. This review summarizes what is known about potential human exposure to produced water from onshore oil and gas development in the United States. We synthesize 236 publications to create a conceptual model of potential human exposure that illustrates the current state of scientific inquiry and knowledge. Exposure to produced water can occur following its release to the environment through spills or leaks during its handling and management. Exposure can also arise from authorized releases, including permitted discharges to surface water, crop irrigation, and road treatment. Knowledge gaps include understanding the variable composition and toxicity of produced water released to the environment, the performance of treatment methods, migration pathways through the environment that can result in human exposure, and the significance of the exposures for human and ecosystem health. Reducing these uncertainties may help in realizing the benefits of produced water use while simultaneously protecting human health.

Environmental Health Perspectives↗

USGS2025

An exploration of the relative influence of physical models for Omori’s law

Omori’s law states that the rate of aftershocks decays as a function of inverse time. There are multiple physical explanations that we reduce into a nonlinear mixed effects relation of three terms: (1) a Rate/State expression that can account for static/dynamic and viscoelastic triggering caused directly by the mainshock, (2) a fluid diffusion triggering term, and (3) a randomized secondary triggering (cascade) term. We fit free physical-model parameters to an observed aftershock sequence through two nonlinear regression methods to find the relative contributions of physics-based models in an observed aftershock sequence. Results from both methods show that Rate/State models overpredict aftershock rates by ∼0–30%. Secondary aftershocks cause a net negative contribution (seismicity rate reduction that corrects overprediction by other terms) ranging between ∼0 and 30%. All regression solutions yield negative secondary triggering contributions without being guided to do so. A physical explanation for this is that aftershock occurrence relieves stress from the crust, ultimately causing the sequence to extinguish itself. Fluid diffusion triggering contributions range from ∼0 to 20%. Diffusion processes are observed to be shorter in time than the full duration of an aftershock sequence and they are also spatially limited, diminishing their influence. Our results apply to an aftershock decay curve from the 2016 Central Apennines earthquake sequence, meaning that our specific results may not be general. Our primary conclusion is that any one physical model cannot alone fit the observed sequence as well as the combination of three we investigated.

Frontiers in Earth Science↗

USGS2026

Stepovers and beyond: Structural control of The Geysers geothermal system and the broader Clear Lake region

Fault geometry exerts a first-order control on geothermal systems by governing stress localization, fracture development, and permeability, yet in complex fault networks or broader shear zones, the relative influence of individual geometric features is often difficult to resolve. In the northern California Coast Ranges, The Geysers geothermal field is commonly interpreted to occur within a releasing stepover, although no single, clearly defined stepover is identified in published studies. To investigate the structural controls on The Geysers and the broader Clear Lake region, a two-dimensional elastic boundary element model is developed to evaluate spatial patterns of dilational strain associated with progressively more complete fault geometries. Model results show that dilation in the region is not controlled by a single structure but instead reflects the combined effects of multiple interacting fault elements. Three primary controls are identified: (1) opposing bends in the regional strike-slip fault system, including a releasing bend along the Maacama fault; (2) the southern fault tip of the Collayomi fault, which generates a prominent dilational lobe beneath the southern Geysers; and (3) a releasing stepover between the Collayomi fault and the Geyser Peak–Mercuryville–Big Sulphur Creek fault system, inferred to collectively behave as a right-lateral shear zone bounding the western margin of The Geysers. Predicted dilational strain magnitudes are sufficient to localize permeability between faults. These results highlight that incorporating complete fault networks and bedrock geological mapping can enhance geothermal assessments and provide a transferable framework for evaluating structurally controlled permeability in tectonically active regions.

California↗

USGS2025

Evidence for fluid pressurization of fault zones and persistent sensitivity to injection rate beneath the Raton Basin

Subsurface wastewater injection has increased the seismicity rate within the Raton Basin over more than two decades, with the basin-wide injection rate peaked between 2009-2015. To understand the evolution of injection-induced earthquakes, we systematically analyzed 2016-2024 broadband recordings with a machine-learning-based phase picker and constructed a catalog with 95,993 earthquakes (-1≤ M L ≤4.3). We then inverted for full centroid moment tensors (CMT) for 90 M L ≥ 2 events, with a special interest in constraining the non-double-couple components via probabilistic metrics. Both relocations and CMT solutions support basement-rooted normal faults, including graben and half-graben structures. Furthermore, we observe the non-double-couple components that imply elevated pore pressure in the fault zones. An earthquake cluster emerged in the north-central basin in 2023, preceded by ~1-yr of increased injection volume from wells within 15km. Despite a basin-wide decrease in the injection volume, we highlights the persistence of seismicity that remains to sensitive to injection rates within the Raton Basin.

Colorado, New Mexico↗