USGS Science⌕ Search

SEARCH · USGS Science

Results for “Environmental Science and Technology - Water”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,665 records · Page 68Linked to original sources
USGS2026

High‐resolution surface deformation and slip distribution observations for the 2023 Kahramanmaraş, Türkiye, earthquake sequence help constrain the rupture process

Splay, or branch, faults are a common geometric feature of earthquake surface ruptures and may provide constraints on the rupture behavior of an earthquake. The 2023 M w 7.8 Pazarcık and M w 7.5 Elbistan, Türkiye, earthquakes are examples of ruptures with multiple small splays, and the Pazarcık earthquake nucleated on a splay fault, the Narlı fault, before rupturing bilaterally on the East Anatolian fault (EAF). Here, we present 3‐m‐resolution surface displacement from subpixel correlation of Planet Dove optical images for the entirety of both ruptures with corresponding surface slip distributions. For a 30‐km‐long study region spanning the Narlı‐EAF intersection, we compare surface slip derived from five data sets with different resolutions (on‐the‐ground, WorldView, Planet Dove, Sentinel‐2, and Sentinel‐1) to elucidate complementary information. In addition, we integrate information from the surface expression of faulting with published dynamic rupture simulations and rupture process studies to constrain a rupture evolution for the Pazarcık earthquake that is consistent across data sets. This work highlights the complementary nature of disparate surface slip data sets and the role that high‐resolution surface displacement information, including from fault splays, can play in constraining nonunique rupture models and refining understanding of the earthquake rupture process.

Bulletin of the Seismological Society of America↗
USGS2026

A diatom-based quantitative sea-ice proxy for the Bering and Chukchi seas

Sea ice affects Earth's climate system on both regional and global scales. Its incorporation into climate can be used to achieve more accurate predictions of future climate. However, instrumental records of sea-ice concentration do not extend earlier than 1978. In an effort to extend this record, we constructed a proxy using the generalized additive model based on relative abundances of five easy-to-identify diatom species found in sediment samples across the Bering and Chukchi seas. Here we present the first quantitative diatom-based sea-ice proxy developed for Beringia. The developed proxy has been applied to two sediment cores in the Bering Sea ranging from 0 to 25.7 ka (HLY0204 51JPC) and 369 to 430 ka (IODP Exp 323 Site U1345) and one in the Chukchi Sea ranging from 2.7 to 10 ka (HLY0204 24JPC). The obtained reconstructions of sea-ice concentrations are similar, but not identical to previously published qualitative and nearby records based on other proxies. Because our results are quantitative, they can be incorporated into regional climate models. The proxy is publicly available as an R Shiny application (app) and can be applied to any diatom count from marine sediments in the region.

Bering Sea, Chukchi Sea↗
USGS2025

Geologic framework and Holocene sand thickness offshore of Seven Mile Island, New Jersey

The U.S. Geological Survey assessed the Quaternary evolution of Seven Mile Island, New Jersey, to quantify coastal sediment availability, which is crucial for establishing sediment budgets, understanding sediment dispersal, and managing coastlines. This report presents preliminary interpretations of seismic profiles, maps of Holocene sand thickness from the shoreline to 2 kilometers offshore, and tables quantifying the volume of available sediment along the coastal margin based on data collected during 2021 and 2022. The results reveal spatial variability in the thickness and cross-shore extent of Holocene sand. The study area was separated into northern, central, and southern zones by using underlying stratigraphy and geomorphic features. The characteristics and spatial extent of the Holocene sand deposit indicate that hydrodynamic processes contribute to its spatial variability. Northern Seven Mile Island contains the thickest deposits of Holocene sand that were formed by sediment bypass around the Townsends Inlet ebb-tidal delta. Specifically, swash bars have welded to the updrift end of Seven Mile Island and have formed thick deposits of Holocene sand that thicken landward and taper seaward. Despite their thickness, these deposits have the smallest cross-shore extent; therefore, northern Seven Mile Island has the smallest volume of Holocene sand of the three geomorphic zones. Central Seven Mile Island has the thinnest Holocene sand deposits because this section of the barrier island is outside the influence of ebb-tidal deltas. Southern Seven Mile Island has the greatest volumes of Holocene sand because of increased accommodation and deposition adjacent to the Hereford Inlet ebb-tidal delta. Even though tidal inlets exert variable influence on the three geomorphic zones, sediment is distributed fairly uniformly within each geomorphic zone; each of the three zones contains 31.05–36.48 percent of the volume of available Holocene sand.

New Jersey↗
USGS2024

Application of non-stationary shear-wave velocity randomization approach to predict 1D seismic site response and its variability at two downhole array recordings

Accounting for uncertainties in seismic site response is crucial to improving the performance of one-dimensional (1D) ground response analyses (GRAs) at downhole array recording sites. In addition to site effects, uncertainties in 1D-GRAs can also be contributed from the seismic source and/or path. Though often representing not more than one percent of the distance (path) from the source, site conditions are known to have an enormous influence on ground shaking. In this study, we focus on the site shear-wave velocity ( V S ) structure, which is the main ingredient for estimating the variability of site response. As such, V S can manifest aleatory uncertainties related to the effects of small-scale spatial heterogeneities within the near surface, thus V S can substantially modify ground shaking during earthquakes. We apply a novel V S randomization approach to propagate the small-scale heterogeneities of V S to estimate seismic site response within a non-stationary probabilistic framework. The randomization approach generates samples of V S profiles that are used to perform several 1D-GRAs and obtain an averaged site response and related variability. The proposed method is implemented on data recorded at two downhole array sites with different subsurface soil conditions: a soft soil site on Treasure Island (California, United States of America) and a rock outcrop site in Cadarache (South-East France). We show that synthetic surface-to-borehole transfer functions from 1D-GRAs provide an acceptable fit to the empirical transfer functions from low-motion earthquake records and succeed in reproducing most of the site-specific seismic response variability. The remaining mismatch between transfer functions is likely due to insufficient precision on the seismic bedrock and the impedance contrast. The variability in site response is discussed with emphasis on the role of V S small-scale heterogeneities, attenuation, and input motion incidence angle in ground motion variability for the site and soil conditions at both locations.

California↗
USGS2025

Quantitative evaluations of earthquake early warning performance using “Did You Feel It?” and post-alert surveys

We examine responses to the U.S. Geological Survey’s “Did You Feel It?” (DYFI) survey and its companion earthquake early warning (EEW) questionnaire to assess the performance of the U.S. ShakeAlert EEW system directly from the alert recipients’ perspectives. ShakeAlert rapidly detects earthquakes and develops alert information, but as official alert delivery partners issue these alerts, it is thus difficult to determine how many people were alerted and when. We investigate DYFI reports for six California earthquakes that had EEW alerts and substantial responses to the DYFI EEW questionnaire. Comparisons of ShakeAlert predictions to reported intensities demonstrate that magnitude estimation accuracy is not necessarily indicative of ground-motion prediction accuracy. Perceived warning time distributions indicate that estimating maximum-expected warning times using the S-wave arrival is a reasonable assumption when discussing public EEW performance. However, we also find many reports of shorter warning times, late alerts, and missed alerts than expected based on ShakeAlert publication times, suggesting alert delivery latencies are substantial and highly variable. The novelty of our analysis is that we demonstrate that the DYFI EEW survey provides useful EEW efficacy information—independent of the specific alerting pathway—that can be used to inform our choices for conveying EEW performance.

California↗
USGS2024

Testing rate‐and‐state predictions of aftershock decay with distance

We analyze aftershocks of the 2019 M 7.1 Ridgecrest mainshock and isolated M 5–6 mainshocks in southern California to test predictions made by the rate‐and‐state friction model of Dieterich (1994) . Rate‐and‐state friction predicts that the seismicity rate after a stress step follows Omori decay, where the Omori c ‐value, which is the saturation in aftershock rate observed at small times, is larger for smaller stress steps. Put in the context of an aftershock sequence, this predicts that the Omori c ‐value will be systematically larger at greater distances from the mainshock. To our knowledge, this predicted effect has not been observed. In part this may be because the Omori c ‐value is difficult to measure because it often reflects short‐term catalog incompleteness rather than a true saturation in aftershock rate. We explore the dependence of the Omori c ‐value on the distance to the mainshock by applying the “ a ‐positive” method ( van der Elst and Page, 2023 ). This method is insensitive to short‐term aftershock incompleteness and allows resolution of the true aftershock rate deep into the mainshock coda. For aftershocks of the Ridgecrest mainshock and stacked M 5–6 mainshocks, we observe systematic differences in early aftershock rates, relative to mainshock distance, consistent with the predictions of rate‐and‐state friction. Furthermore, for the larger Ridgecrest dataset, we observe that aftershocks nearer to the mainshock start earlier, and we resolve a flattening of the Omori curve consistent with a larger Omori c ‐value for the farthest aftershocks, as predicted by Dieterich (1994) .

Seismological Research Letters↗
USGS2026

Surface-wave relocation and characterization of the 2024–2025 dyking episode along the Fentale–Dofen segment of the Ethiopian rift

From late-December 2024 to mid-March 2025, a 50-km-long dyke intrusion triggered over 300 earthquakes (magnitude 4 to 5.9) between Fentale and Dofen volcanoes along the Northern Main Ethiopian Rift. Dyke intrusions periodically occur along the Fentale–Dofen magmatic segment and are an expression of ongoing rift extension. Preliminary analyses using interferometric synthetic aperture radar revealed extensive ground deformation (up to 60 cm), which closely matched the temporal and spatial evolution of surface manifestations and earthquake locations from global catalogues. While global catalogues are critical for real-time monitoring, the precision of locations in remote and or sparsely instrumented regions can be low. In this investigation, we present surface-wave relocation results of the dyking episode that began near Fentale volcano in December 2024. We estimate relative locations using differential traveltimes measured from regional-to-teleseismic distance surface-wave observations of earthquakes reported by the U.S. Geological Survey. Relative relocations reduce the initial region of diffuse seismicity to a 50-km-long narrow band bounding the strike of surface manifestations and the zone of maximum surface deformation. We demonstrate the precision of surface-wave relocations over incremental time periods, capturing the progression of dyking from seismic onset through seismic migration and caldera subsidence. Results showcase the utility of surface-wave relocations in the characterization of dyking episodes and provide complementary insights into the current understanding of the Fentale–Dofen volcanic plumbing system.

Geophysical Journal International↗
USGS2025

A global database of soil microbial phospholipid fatty acids and enzyme activities

Soil microbes drive ecosystem function and play a critical role in how ecosystems respond to global change. Research surrounding soil microbial communities has rapidly increased in recent decades, and substantial data relating to phospholipid fatty acids (PLFAs) and potential enzyme activity have been collected and analysed. However, studies have mostly been restricted to local and regional scales, and their accuracy and usefulness are limited by the extent of accessible data. Here we aim to improve data availability by collating a global database of soil PLFA and potential enzyme activity measurements from 12,258 georeferenced samples located across all continents, 5.1% of which have not previously been published. The database contains data relating to 113 PLFAs and 26 enzyme activities, and includes metadata such as sampling date, sample depth, and soil pH, total carbon, and total nitrogen. This database will help researchers in conducting both global- and local-scale studies to better understand soil microbial biomass and function.

Scientific Data↗
USGS2026

Regional characterization of coal resources in the Gulf Coast Coal Region—In support of coal as a potential source for critical minerals and rare earth elements

There is increasing interest in extracting critical minerals (CM), including rare earth elements (REE), from coals in the United States to address the overreliance on imported REE. The Gulf Coast and the Fort Union regions are the two major lignite-bearing regions within the United States. Recent REE and CM studies of the lignite in these regions have indicated that the coals may be a viable source material for REE and CM extraction. To evaluate in-place coal as a potential source of REE and CM, the coal resources need to be quantified. This study presents the results of a regional analysis of Gulf Coast Coal Region lignite and bituminous coal resources that might be available as potential sources of REE and CM. The resource analysis used kriging methods to develop isopleth maps of cumulative coal thickness throughout the region using data from drill holes and measured mine and outcrop sections for a total of 31,181 data points. The estimated total remaining coal resource in the Gulf Coast Coal Region is about 83 billion metric tons in the upper 90 meters (approximately 300 feet) of the subsurface. The current estimate of remaining coal resources includes approximately 1.88 billion metric tons of coal that have been extracted in the region since 1983. Texas accounted for 39 percent (32 billion metric tons) of the remaining total resource, followed by Mississippi (24 percent, 20 billion metric tons ), Louisiana (14 percent, 12 billion metric tons), Tennessee (10 percent, 8.5 billion metric tons), and Arkansas (6 percent, 5.1 billion metric tons). Alabama, Georgia, Kentucky, and Missouri each accounted for less than 5 percent of the total resource; Georgia had the smallest resource estimated at 7 million metric tons. Here we report the first known statewide lignite resource estimates for Georgia, Kentucky (820 million metric tons), and Missouri (1,800 million metric tons). A comparison of the results of this study with those of previous Gulf Coast Coal Region and Fort Union Coal Region resource studies is difficult because each study used different data sources, assessment methodologies, overburden depths, and qualifying coal thicknesses.

Alabama, Arkansas, Florida, Georgia, Illinois, Lou↗
USGS2025

The influence of pre-existing structures on geothermal springs: Inferences from potential field mapping in Surprise Valley, CA and other sites In the northwestern Great Basin

Surprise Valley, located in the northwestern Great Basin, is an asymmetric extensional basin that marks a major tectonic transition between the relatively un-extended volcanic Modoc Plateau to the west, and the Basin and Range to the east that has undergone 10-15% extension. In addition, it sits just north of the Walker Lane which accommodates up to 20% of dextral slip associated with Pacific-North American plate interactions. Thermal springs issue from eight areas within Surprise Valley. Most of these occur within the basin and are not situated on the main basin forming range-front faults. As a result, efforts to resolve the structural setting of the valley’s hydrothermal system have relied on geophysics to characterize basin structure and geology. Extensive efforts to map the basin with ground and airborne magnetics have revealed a >35 km-long linear, intra-basin magnetic high, interpreted as a buried dike swarm. Geothermal springs on the eastern side of the valley, including Seifert hot springs, Leonards hot springs, and Surprise Valley hot springs (SVHS), are all situated along the magnetic high and occur at local breaks and bends in the anomaly, suggesting that fracture permeability is enhanced along the feature and particularly at these discontinuities. Recent studies, including drilling over the anomaly near SVHS that likely intersected dike material, as well as subsequent mapping and sampling of dikes outcropping along the anomaly on the playa surface south of SVHS, confirm (as previously inferred) that mafic intrusives are the principal source of the anomaly. Similar interpretations made in two other valleys (in southern Oregon and northwestern Nevada), where inferred intra-basin dikes appear to be spatially correlated with hot springs or prospective geothermal resource areas, suggest that the impact of pre-existing basement structure on hydrothermal activity may pertain more generally to other hydrothermal settings throughout the Great Basin. If so, efforts to map basement may enhance understanding structural controls on some geothermal systems. Furthermore, similarities across these disparate sites suggest that magmatism may play a much larger role in accommodating extension and influencing basin evolution across the western Great Basin than previously recognized.

Conference Paper↗
USGS2024

Post-glacial stratigraphy and late Holocene record of great Cascadia earthquakes in Ozette Lake, Washington, USA

Ozette Lake is an ~100-m-deep coastal lake located along the outer coast of the Olympic Peninsula (Washington, USA); it is situated above the locked portion of the northern Cascadia megathrust but also relatively isolated from active crustal faults and intraslab earthquakes. Here we present a suite of geophysical and geological evidence for earthquake-triggered mass transport deposits (MTDs) and related turbidite deposition in Ozette Lake since ca. 14 ka. Comprehensive high-resolution bathymetry data, seismic reflection profiles, and sediment cores are used to characterize the post-glacial stratigraphic framework and examine paleoseismic evidence in the lacustrine sediments. Stacked sequences of MTDs along the steep eastern flanks of the lake appear to grade basin-ward from thick, chaotic, blocky masses to thin, parallel-bedded turbidite beds. The discrete turbidite event layers are separated by fine-grained (silt and clay) lake sedimentation. The event layers are observed throughout the lake, but the physical characteristics of the deposits vary considerably depending on proximity to primary depocenters, steep slopes, and subaqueous deltas. A total of 30–34 event deposits are observed in the post-glacial record. Radiometric dating was used to reconstruct a detailed sedimentation history over the last ~5.5 k.y., develop an age model, and estimate the recurrence (365–405 yr) for the most recent 12 event layers. Based on sedimentological characteristics, temporal overlap with other regional paleoseismic chronologies, and recurrence estimates, at least 10 of the dated event layers appear to be sourced from slope failures triggered by intense shaking during megathrust ruptures; the recurrence interval for these 10 events is 440–560 yr. Thus, Ozette Lake contains one of the longest and most robust geological records of repeated shaking along the northern Cascadia subduction zone.

Washington↗
USGS2026

Small earthquake moment magnitude and implications for frequency–magnitude scaling of injection induced earthquakes of the Raton Basin

Accurate estimation of earthquake source parameters—such as moment magnitudes, corner frequencies, and stress drops—is essential for improving seismic hazard assessments and understanding earthquake physics. In this study, moment magnitudes ( M W ) are calculated for 31,581 earthquakes associated with wastewater injection in the Raton Basin (located along the border between northern New Mexico and southern Colorado) between 2016 and 2024 using radiative transfer theory to fit coda decay envelopes. Our results show that it is feasible to estimate moment magnitudes down to M W ~1 with coda envelopes from a small local monitoring network. Significant differences were found between M W and local magnitudes ( M L ) for small earthquakes ( M < 3.0). A linear relationship was optimized to convert M L to M W : M W = 0.7 M L + 0.96 and M W = 0.73 M L + 0.99 (for the events reported by the U.S. Geological Survey), which can be applied in future studies of Raton Basin seismicity. We find that b -values calculated employing different methods and using M L are approximately 1.0, while those using M W range from 1.2 to 1.4. A larger estimate of the b -value could influence interpretations of the statistical behavior of earthquakes associated with injection and consequently seismic hazard assessments based on a magnitude–frequency distribution. The potential differences between local versus moment magnitude-based earthquake statistics should be considered in other seismically active regions.

Colorado, New Mexico↗
USGS2021

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

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

Pasco↗
USGS2025

Origin and evolution of mafic volcanism associated with 3 m.y. of andesite production at the Goat Rocks volcanic cluster, southern Washington Cascade Range

More than 3 m.y. of mafic volcanism near the Goat Rocks volcanic cluster in the southern Washington Cascade Range, USA, lends insight into the evolution of basalts and the subarc mantle at a long-lived, major arc volcanic locus. We contribute field observations, 40 Ar/ 39 Ar dates, paleomagnetic directions, and bulk rock and mineral compositions to characterize nine mafic units that erupted in association with the Goat Rocks volcanic cluster. The time frame of mafic volcanism, ca. 3.6 Ma to 60 ka, encompasses the lifespan of the central volcanic cluster (3.1 Ma to 115 ka), with a lull from ca. 2.7 Ma to 1.4 Ma. A climactic period of voluminous mafic activity and far-traveled lava flows, including construction of the Hogback Mountain shield volcano, coincided with voluminous andesite eruptions from the central volcanic cluster. The basaltic rocks in the Goat Rocks area are calc-alkaline to barely tholeiitic and have high field strength element depletion relative to large-ion lithophile elements characteristic of calc-alkaline basalts (CAB) of the Cascade volcanic arc. Unlike at neighboring andesitic volcanic centers (Mounts Adams, St. Helens, and Rainier), no other mafic end members such as high-aluminum olivine tholeiite (HAOT) or intraplate-type basalt (IPB) are present at or near the Goat Rocks volcanic cluster, although some of the calc-alkaline basalts in this study have IPB-like affinities. The Goat Rocks mafic units exhibit two main temporal trends in composition: (1) the most primitive basalts erupted earlier, compared to less primitive and more evolved compositions later, and (2) high field strength element concentrations are higher in the younger basalt units relative to the oldest two. In contrast to these temporal trends, the mafic units define two compositional groups that recur through time, a low-Sr and a high-Sr group, each with distinct trace element and Sr and Nd isotope ratios. Although radiogenic isotope ratios are generally aligned with High Cascades CAB and HAOT, some extend toward IPB of Mount Adams and Simcoe Mountains volcanic field. Olivine-dominated crystal fractionation at shallow pressure from a small range of parent magma compositions accounts for much of the variation among the basalts and basaltic andesites. A high-pressure fractionation model is plausible for only one of the youngest basalt units (basalt of Walupt Lake volcano). Mafic recharge and crustal assimilation accounts for the incompatible-element enriched composition of basaltic andesites erupted during construction of the largest andesitic centers, further supporting sustained basalt mass flux and thermal energy driving andesite genesis. We model the most primitive members of the Goat Rocks mafic units as partial melts of successively less depleted mantle in time. Variable degrees of fluxing with fluids and melts from subduction explain the distinction between high-Sr and low-Sr groups. We propose that mantle metasomatism by ancestral subduction and fluid-flux melting is heterogeneously distributed through the local subarc mantle and played a greater role in the genesis of the high-Sr basalt group. The limited range of primitive basalt types around the Goat Rocks volcanic cluster contrasts with the much greater diversity of basalts throughout the southern Washington to northern Oregon Cascade arc. On the other hand, the central volcanic cluster encompasses nearly the entire diversity observed at neighboring composite volcanoes. In the case of the Goat Rocks area at least, and perhaps attributable to the entire region, this means that the genesis of diverse intermediate magmas is independent from and does not require vastly different parental basalt compositions.

California, Oregon, Washington↗
USGS2025

Constraining large magnitude event source and path effects using ground motion simulations

The purpose of this study is to use ground motion simulations to investigate ways in which source and path effects for large magnitude events can be represented in non-ergodic GMMs. While we initially developed computation techniques using CyberShake simulations, the range of magnitudes and source-site combinations is not adequate to replicate what is observed empirically. We therefore designed a new ground motion simulation study, which includes earthquakes with a large range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of simulations (M4-M7), we then develop a non-ergodic GMM with the simulation data. We find that the within-site residuals are dominated by the radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within-site residuals among multiple source realizations. Finally, after removing the source effects from the within-site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site with the large events, while the second approach considers all small events on the fault plane. The results indicate that the path effects of large events cannot be satisfactorily approximated with that of small events using either approach.

Conference Paper↗
USGS2026

Neutron scattering reveals fractionation of natural gas mixtures in unconventional petroleum reservoir pores: Perspectives on energy resource recovery and storage

In unconventional petroleum reservoirs hydrocarbon fluids are hosted by both mineral and organic matter pores. These pores can have diameters that range from microns to less than a single nanometer and, for unconventional reservoirs, there is evidence that small pores ( <20 nm diameter) may constitute a large proportion of the available space. Understanding subsurface volumes and how fluids behave in them can be helpful for predicting hydrocarbon production and storage in the subsurface. One area with knowledge gaps regarding hydrocarbon behavior in small pores is the possibility for mixtures to fractionate (i.e., unmix) based on pore size or pore type. Mixture fractionation as a function of pore size could impact recovery of hydrocarbons, drive compositional shifts during production, and limit fluid storage within candidate reservoirs. To investigate natural gas fractionation in small geologic pores, we applied total neutron scattering to probe methane-ethane mixtures at reservoir pressures (up to ≈30 MPa) and temperature (60°C) within a sample from the Upper Cretaceous Niobrara Formation. Neutron scattering data reveal only minor fractionation occurs between methane and ethane in 20-nm diameter sample mesopores. Increased fractionation is observed for sample micropores, with up to 72% (±1% at 1-sigma) methane found in 2 nm diameter pores following injection of a 50%-50% methane-ethane mixture. These data provide rarely available direct experimental observations of hydrocarbon mixture behavior under nanoconfinement in a sample from an important unconventional petroleum reservoir. Our results are discussed in the context of evaluating hydrocarbon resources in unconventional reservoir meso- and micropores, reconciling observed gas composition changes during production, and more broadly, understanding subsurface pore volumes within an energy storage framework.

Fuel↗
USGS2025

Induced earthquakes are generally not tidally triggered in Oklahoma and Kansas

Human-induced earthquakes occur along critically stressed faults as injected wastewater simultaneously heightens fluid pressure and pushes faults to failure. We investigate the possibility that small stresses imposed by Earth tides could trigger earthquakes in the induced seismicity region of Oklahoma and Kansas from 2011 to 2018. We decluster a catalog consisting of ∼110,000 earthquakes using three methods (Reasenberg, nearest-neighbor distance, and phase-bin). We find no significant tidal earthquake triggering using Schuster's p -value test for the declustered catalogs as a whole. We search for localized triggering using discretized space-time cells and find ∼0–6% of cells have significant tidal triggering which is close to what is randomly expected (5%) and indicates there is an insignificant amount of tidal triggering for the full study region. One area that has significant p -values across multiple time windows, ∼2014–2016 is ∼15 km from a region of large wastewater injection volume. It is possible that localized tidal triggering occurs for this time and area because faults remain critically stressed and are particularly susceptible to slip under the small stress load from the semidiurnal tide. Possible explanations for the lack of tidal triggering in our broader study are that the pre-seismic stressing rate in the earthquake nucleation area is faster than the tidal stressing rate (∼3 kPa/day), faults are not close enough to critically stressed to be affected by tidal forcing, and that nucleation occurs over longer periods than the tides considered in this study (∼1, ∼14 days). Fluid injection could be the source of a higher pre-seismic stress rate.

Kansas, Oklahoma↗
USGS1976

Transport and dispersion of fluorescent tracer particles for the dune-bed condition, Atrisco Feeder Canal near Bernalillo, New Mexico

A tracer technique in which mineral particles were coated with fluorescent dyes was used to study the rates of transport and dispersion of sediment particles`of various diameters and specific gravities for a dune-bed condition in an alluvial channel. The experiment was conducted in the Atrisco Feeder Canal near Bernalillo, N. Mex., between May 1 and July 14, 1967. A continuous point source of tracers, approximated by injections at 5- or 10-minute intervals, was maintained for 7 days so that the steady-dilution procedure could be used to calculate the transport rate of bed material. After termination of the injection process, the spatial-integration procedure was used to follow the movement of the tracers downstream and to calculate the transport rate. Samples of the bed material in transport and the accompanying tracers moving along the surface of the dune bed were obtained periodically throughout the study with the "dustpan" sampler especially designed for fluorescent tracer studies. In addition, the spatial distributions of the tracers in the dune bed were determined three times during the study by core sampling.

New Mexico↗