USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Open Source Education”

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,683 records · Page 27Linked to original sources
USGS2024

Characterizing directivity in small (M 2.4-5) aftershocks of the Ridgecrest sequence

Directivity, or the focusing of energy along the direction of an earthquake rupture, is a common property of earthquakes of all sizes and can cause increased hazard due to azimuthally dependent ground‐motion amplification. For small earthquakes, the effects of directivity are generally less pronounced due to reduced rupture size, yet the directivity in small events can bias source property estimates and provide important insights into general regional faulting patterns. However, due to observational limitations, directivity is usually only measured and modeled for large events. As such, many studies of small earthquakes either ignore directivity altogether or assume a constant rupture direction for all events in a cluster. In our study, we apply a refined directivity fitting method constrained with two separate methods of source deconvolution to the dataset of aftershocks of the 2019 Ridgecrest earthquakes, which contain a large number of well‐recorded small‐to‐mid sized earthquakes occurring in close proximity to each other. The revealed directivity of 100+ small (M 2.4–5) earthquakes is highly heterogeneous and primarily oblique to and away from the main fault strike, suggesting a complex postseismic stress redistribution. In addition, the energy focusing effect of directivity appears to bias the selection of high‐quality data from stations in the direction of rupture, leading to average stress‐drop increases of 50% if directivity is not accounted for.

California↗
USGS2024

Recent applications of the USGS National Crustal Model for Seismic Hazard Studies

The U.S. Geological Survey is developing the National Crustal Model (NCM) for seismic hazard studies to facilitate modeling site, path, and source components of seismic hazard across the conterminous United States. The NCM is composed of a 1km grid of geophysical profiles, extending from the Earth’s surface into the upper mantle. It is constructed from a threedimensional (3D) geologic framework and geophysical rules that use (1) a petrologic and mineral physics database; (2) a 3D temperature model; and (3) a calibrated rock type- and age-dependent porosity model. Parameters needed to estimate site response for existing ground motion models (GMMs), including the time-averaged velocity in the upper 30 meters (VS30), the depths to 1.0 and 2.5 km/s shear-wave velocity (Z1.0 and Z2.5), and sediment thickness, can be computed from the NCM. As GMMs continue to improve in the future, other metrics could also be extracted or derived from the NCM, such as fundamental period, site attenuation (ko), a fully frequency-dependent site response function, or 3D geophysical volumes for wavefield simulations. Application of the NCM may also benefit other aspects of seismic hazard analysis, including better accounting for path-dependent attenuation and geometric spreading, more accurate estimation of earthquake source properties such as hypocentral location and stress drop, and calculation of crustal strength profiles that inform estimates of the base of seismicity.

conterminous United States↗
USGS2025

Towards mobile wind measurements using joust configured ultrasonic anemometer for applications in gas flux quantification

Small uncrewed aerial systems (sUASs) can be used to quantify emissions of greenhouse and other gases, providing flexibility in quantifying these emissions from a multitude of sources, including oil and gas infrastructure, volcano plumes, wildfire emissions, and natural sources. However, sUAS-based emission estimates are sensitive to the accuracy of wind speed and direction measurements. In this study, we examined how filtering and correcting sUAS-based wind measurements affects data accuracy by comparing data from a miniature ultrasonic anemometer mounted on a sUAS in a joust configuration to highly accurate wind data taken from a nearby eddy covariance flux tower (aka the Tower). These corrections had a small effect on wind speed error, but reduced wind direction errors from 50° to >120° to 20–30°. A concurrent experiment examining the amount of error due to the sUAS and the Tower not being co-located showed that the impact of this separation was 0.16–0.21 ms − 1 "> ms − 1 , a small influence on wind speed errors. Lower wind speed errors were correlated with lower turbulence intensity and higher relative wind speeds. There were also some loose trends in diminished wind direction errors at higher relative wind speeds. Therefore, to improve the quality of sUAS-based wind measurements, our study suggested that flight planning consider optimizing conditions that can lower turbulence intensity and maximize relative wind speeds as well as include post-flight corrections.

Alaska↗
USGS2024

Short-term sediment dispersal on a large retreating coastal river delta via 234Th and 7Be sediment geochronology: The Mississippi River Delta Front

Many Mississippi River Delta studies have shown recent declines in fluvial sediment load from the river and associated land loss. In contrast, recent sedimentary processes on the subaqueous delta are less documented. To help address this knowledge gap, multicores were collected offshore from the three main river outlets at water depths of 25–280 m in June 2017 just after the peak river discharge period, with locations selected based on 2017 U.S. Geological Survey seabed mapping. The coring locations included the undisturbed upper foreset, mudflow lobes, gullies, and the undisturbed prodelta. Nine multicores were analyzed for Beryllium-7 activity, and four cores were analyzed for excess Thorium-234 activity via gamma spectrometry, granulometry and X-radiography. Our results indicate a general trend of declining 7 Be and 234 Th activities and inventories with increasing distance from sources and in deeper water. The core X-radiographs are graded from the predominantly physically stratified nearshore to the more bioturbated offshore, consistent with the sedimentation patterns. Sediment focusing assessed via the 7 Be and 234 Th sediment inventories shows preferential sedimentation in gully and lobe environments, whereas the upper foreset and prodelta focusing factors are relatively depleted. Overall, short-term sediment deposition from the main fluvial source remains active offshore from all three major river outlets, despite the overall declining river load.

Louisiana↗
USGS2025

Base-flow sampling to enhance understanding of the groundwater flow component of nitrogen loading in small watersheds draining into Long Island Sound

Excessive nitrogen discharge is a major concern for the Long Island Sound. Programs have been implemented to reduce point sources of nitrogen to the sound, but little is known about the nonpoint sources. This study aims to better understand the current groundwater contributions of nitrogen from nonpoint sources in the Long Island Sound watershed. During the spring and summer of 2022, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, collected water-quality samples to analyze nutrients (nitrogen and phosphorus), chloride, and bromide at 45 stations in the Long Island Sound watershed in Connecticut, New York, and Rhode Island. The stations were in small drainage watersheds (5 to 30 square kilometers) in the southern part of the Long Island Sound watershed. During two separate synoptic sampling events, water-quality samples and instantaneous streamflow measurements were collected under base-flow conditions (where the streamflow is dominated by groundwater inputs rather than overland flow or runoff flow). One sampling event was in the nongrowing season (April 24–25, 2022), and the other was in the growing season (June 30–July 1, 2022). To calculate instantaneous nitrogen loads and yields, streamflow was measured at the time of sample collection. Nitrogen concentrations, loads, and yields varied among sampling stations and by season. Total filtered nitrogen concentrations were generally lower in the nongrowing season (from less than 0.14 to 1.9 milligrams per liter) than in the growing season (from less than 0.23 to 3.0 milligrams per liter). Nitrate plus nitrite concentrations showed little variation between the nongrowing and growing seasons. Unfiltered ammonia plus organic nitrogen concentrations were generally lower in the nongrowing season (from less than 0.07 to 0.83 milligram per liter) than in the growing season (from 0.11 to 0.98 milligram per liter). In contrast, total filtered and unfiltered nitrogen loads and yields were higher in the nongrowing season than during the growing season, likely because streamflows were higher during the nongrowing season. Total unfiltered nitrogen yields during the nongrowing season ranged from less than 0.15 to 5.0 kilograms per square kilometer per day. Total unfiltered nitrogen yields during the growing season ranged from less than 0.12 to 2.5 kilograms per square kilometer per day. Total filtered nitrogen yields during the nongrowing season ranged from less than 0.13 to 5.2 kilograms per square kilometer per day. Total filtered nitrogen yields during the growing season ranged from less than 0.06 to 2.5 kilograms per square kilometer per day.

Connecticut, New York, Rhode Island↗
USGS2024

Assessing locations susceptible to shallow landslide initiation during prolonged intense rainfall in the Lares, Utuado, and Naranjito municipalities of Puerto Rico

Hurricane Maria induced about 70 000 landslides throughout Puerto Rico, USA, including thousands each in three municipalities situated in Puerto Rico's rugged Cordillera Central range. By combining a nonlinear soil-depth model, presumed wettest-case pore pressures, and quasi-three-dimensional (3D) slope-stability analysis, we developed a landslide susceptibility map that has very good performance and continuous susceptibility zones having smooth, buffered boundaries. Our landslide susceptibility map enables assessment of potential ground-failure locations and their use as landslide sources in a companion assessment of inundation and debris-flow runout. The quasi-3D factor of safety, F 3 , showed strong inverse correlation to landslide density (high density at low F 3 ). Area under the curve (AUC) of true positive rate (TPR) versus false positive rate (FPR) indicated success of F 3 in identifying head-scarp points (AUC = 0.84) and source-area polygons (0.85 ≤ AUC ≤ 0.88). The susceptibility zones enclose specific percentages of observed landslides. Thus, zone boundaries use successive F 3 levels for increasing TPR of landslide head-scarp points, with zones bounded by F 3 at TPR = 0.75, very high; F 3 at TPR = 0.90, high; and the remainder moderate to low. The very high susceptibility zone, with 118 landslides km −2 , covered 23 % of the three municipalities. The high zone (51 landslides km −2 ) covered another 10 %.

Puerto Rice↗
USGS2025

Petrogenesis and mineralization potential of spinifex komatiitic basalts in the Bradley Peak greenstone terrane, Wyoming Province

Komatiitic volcanic rocks are important hosts of Ni sulfide mineralization and record early Earth evolution; however, those in the well-studied Archean Wyoming Province have received little attention. Here, we elucidate the timing and petrogenesis of the Bradley Peak komatiitic volcanic rocks using field and textural observations, geochronology, and geochemistry. Detrital and igneous zircon U-Pb ages for two samples from previously undated units support published age determinations, placing the eruption age at 2.72 Ga. Stratigraphy of the volcanic flows was mapped and 36 samples including cumulates, greenschists, and spinifex-textured rocks were collected. Whole-rock geochemistry was used to classify the spinifex-textured samples as Al-undepleted komatiitic basalts (11–17 wt% MgO). Platinum-group element concentrations (n = 25) are like those in global Al-undepleted komatiitic basalts, and PGE/Ti ratios do not indicate the volcanic flows likely host sulfide mineralization. Initial ε Nd values of −0.5 to +4.7 (n = 16), indicate that these lavas were derived from a depleted mantle source and have negligible evolved crust contamination. The primary magma to the komatiitic basalt flows is estimated to have had 19 wt% MgO and be derived from ∼15 to 25 % mantle partial melting at 3–4 GPa. Trace element chemistry and thermodynamic modeling suggest the primary melt assimilated local banded iron formation. Although the Bradley Peak komatiitic basalts do not contain positive evidence of magmatic sulfide deposits, depleted Au in the flows suggests they could be source rocks for nearby orogenic gold deposits.

Wyoming↗
USGS1992

Public water-supply systems and water use in Tennessee, 1988

This report summarizes the results of a study conducted by the U.S. Geological Survey, in cooperation with the Tennessee Department of Environment and Conservation (TDEC), Division of Water Supply in 1988. Data gathered during an inventory by the TDEC were collated to determine water use, supply sources, population served, and design and storage capacities of the systems. The inventory was limited to systems that were active on June 30, 1988. Results of a survey of the systems conducted by the Tennessee Department of Health and Environment during 1988 were a primary source of data for this report. Data from computer and manual files maintained by the Tennessee Department of Health and Environment and the U.S. Geological Survey also were used. The Division of Water Supply, TDEC, surveyed 541 public water-supply systems. These systems served 81 percent of the population of the State, or 3.95 million people. The gross per capita use statewide for public-supplied water was 179 gallons per day. Total water withdrawals for public supply increased about 39 percent from 510 million gallons per day (Mgal/d) in 1980, to 708 Mgalld in 1988. During the same period, the population increased about 7 percent. Surface-water withdrawals accounted for 63 percent (446 Mgal/d) of the total water withdrawn in the State. All of these withdrawals occurred in the Tennessee (56 percent or 249 Mgal/d) and the Ohio (44 percent or 197 Mgalld) hydrologic regions. Ground water supplied 262 Mgal/d or 37 percent of the total water withdrawn by public-supply systems statewide. Of that amount, 79 percent, or 208 Mgalld, was used in western Tennessee.

Tennessee↗
USGS2024

New methodology for assessing underground natural gas storage resources – Example from Michigan Basin, United States

Energy consumption in the United States (U.S.) and across the world is shifting away from traditional fossil fuels like coal and oil, and towards natural gas and renewable sources, including hydrogen. Because gas demand is typically greatest during cold seasons and renewable sources sometimes produce variable supplies, it is important to store energy for use when demand exceeds supply. Whereas batteries and tanks typically store energy above ground, geologic (underground) storage may be able to retain much greater quantities of energy over much longer durations (e.g., Matos et al. 2019; Buursink et al. 2023). Consequently, the U.S. Geological Survey (USGS) is developing a methodology to locate new geologic energy storage sites with an initial focus on depleted hydrocarbon reservoirs followed by estimates of associated pore space or gas storage capacity. Depleted hydrocarbon fields in Michigan Basin in the U.S. already host large quantities of injected natural gas (about 681 billion cubic feet of seasonal gas capacity in 42 gas storage facilities; U.S. Energy Information Administration 2024). Therefore, we used this basin as an example to assess potential new storage resources alongside existing storage capacity and carbon dioxide sequestration resources (e.g., Katz and Coats 1968; Haagsma et al. 2020). In Michigan Basin gas storage is taking place in reservoirs within several key hydrocarbon producing formations with different lithologies, including the informal Stray sandstone of the Michigan Formation and the pinnacle reefs of the Niagara Group (Gautier et al. 1995; Swezey et al. 2015). The assessment methodology introduced here consists of two steps relying on both well information and reservoir data.

Michigan↗
USGS2026

Earthquake ground-motion model adjustments for the San Francisco Bay area

We develop adjustments to ergodic ground‐motion models (GMMs) to improve their performance in the San Francisco Bay Area (SFBA). GMMs are widely used in hazard assessments to estimate characteristics of ground shaking based on known properties of the source, path, and site. Such models are often developed using datasets containing records from various regions, resulting in models that represent median ground‐motion behavior, which may not adequately represent ground motions within subregions. This is true for the SFBA, where ground motions attenuate more rapidly with distance than in many other parts of California that dominate GMM databases. To support improved seismic hazard estimates in the SFBA, we calculate regional constants and anelastic attenuation coefficient adjustments relative to two commonly used ergodic GMMs: BSSA14 ( Boore et al. , 2014 ) and ASK14 ( Abrahamson et al. , 2014 ). These adjustments are obtained for a suite of ground‐motion intensity measures (peak ground acceleration, peak ground velocity, and 5%‐damped pseudospectral acceleration at oscillator periods ranging from 0.075 to 10 s) using mixed‐effects regression. Use of the regionally adjusted models reduces the overall bias by up to 0.5 natural log units for BSSA14 and up to 0.6 natural log units for ASK14. We demonstrate one application of our attenuation adjustments and their implications in an earthquake early warning case study of the 2014 M 6.0 South Napa earthquake. The predicted extent of shaking using the adjusted models better matches observed shaking at large source‐to‐site distances, especially for lower shaking intensities, thus potentially reducing overalerting. We encourage the use of our model adjustments when ergodic models are considered for seismic hazard studies in the SFBA.

California↗
USGS2025

To heal or not to heal?: 2. The moment-recurrence time behavior of repeating earthquakes in the 2011 Prague, Oklahoma aftershock sequence is consistent with laboratory healing rates

The timing and failure conditions of an earthquake are governed by the interplay between fault reloading and restrengthening. The moment-recurrence time behavior of repeating earthquakes can give observational estimates of fault healing rates; however, it is difficult to link these observed healing rates to laboratory studies of frictional healing in part because of uncertainty in lithology. Here, we study the 2011 Prague earthquake sequence, which includes repeating earthquakes in the Arbuckle group and the granitic basement, and compare them to laboratory experiments on samples of the Arbuckle and Troy granite (representative of the basement rock) (Okamoto et al., 2025, https://doi.org/10.1029/2024JB030573 ). We find three spatially distinct groups of repeating earthquakes with different moment-recurrence behavior: (a) constant moment-recurrence time in the Arbuckle group, (b) scattered moment-recurrence time at the intersection of the foreshock-mainshock fault in the granitic basement, and (c) moment-predictable behavior outside of the foreshock-mainshock fault intersection also in the granitic basement. Our observation of stagnant healing for repeating sequences in the Arbuckle group is consistent with laboratory observations of low healing rates for moderately high pore fluid pressures in Arbuckle samples. For the moment-predictable group, the source radius that is required in order to match healing rates is consistent with source radius estimations when taking into account reasonable attenuation of the P- pulse width. Overall, we observe diverse healing behaviors in the seismic families that are consistent with laboratory healing rates, providing seismic evidence that contact-scale frictional mechanisms are relevant to large-scale earthquake dynamics.

Oklahoma↗
USGS2021

Critical review of mercury methylation and methylmercury demethylation rate constants in aquatic sediments for biogeochemical modeling

Mercury is a toxin that causes neurological impairments in adults, is particularly harmful for fetuses and children, and is deadly in severe cases, making it a worldwide health concern. Methylmercury (MeHg) is the environmentally relevant form of mercury (Hg) because it biomagnifies along the food chain. Methylmercury is mainly produced in aquatic sediments via methylation of inorganic Hg (Hg(II)) and transformed back via demethylation. Because transformation rates determine MeHg concentrations, quantification of methylation and demethylation rates is needed to inform management of MeHg. Published rate constants for Hg(II) methylation ( 𝑘 𝑚 ) and MeHg demethylation ( 𝑘 𝑑 ) vary greatly, stemming partly from differences in experimental methods. We conducted a comprehensive review of rate laws, evaluated published rate constants, and performed biogeochemical simulations to assess variability in reported 𝑘 𝑚 and 𝑘 𝑑 . Based on selected studies employing the same pseudo-first-order rate law and similar experimental methods, we found that 𝑘 𝑚 = 0.04 ± 0.03 d −1 is a reasonable range for wetland sediments. Over a number of environments, maximum 𝑘 𝑑 was smaller at sites without Hg source ( 𝑘 𝑑 = 0.5 d −1 ) than at sites with identified Hg source ( 𝑘 𝑑 = 1.8 d −1 ). Larger variability and higher uncertainty in 𝑘 𝑑 compared to 𝑘 𝑚 highlight the need for more research on MeHg demethylation rates. This critical review: (a) aids the design of future experimental studies of 𝑘 𝑚 and 𝑘 𝑑 ; (b) provides guidance for comparing rate constants from different studies; (c) presents a biogeochemical reaction model to assess rate constants; and (d) informs selection of 𝑘 𝑚 and 𝑘 𝑑 values from the literature for use in model simulations.

Critical Reviews in Environmental Science and Tech↗
USGS2025

Evidence of mineral alteration in a salt marsh subterranean estuary: Implications for carbon and trace element cycling

Subterranean estuaries (STE) in salt marshes are biogeochemically active zones where interactions between terrestrial groundwater and seawater drive complex cycling of carbon and trace elements, influenced by mineral dissolution. These systems, characterized by fine-grained organic-rich peat overlying permeable coastal aquifers, play a crucial role as a blue carbon sink, yet their geochemical dynamics remain poorly understood. We investigated dissolved trace elements, carbon, silica, and radium isotopes in a salt marsh STE (Sage Lot Pond, Waquoit Bay, MA) over seasonal and annual cycles. Our results reveal that groundwater and estuarine water circulation through marsh peat and aquifer sediments leads to enrichments of dissolved organic and inorganic carbon (DOC and DIC), Si, Ba, and Mn, with variable source/sink behavior of Fe and net removal of U. Submarine groundwater discharge dominated Ba fluxes, whereas pore water drainage from marsh peat acted as the main sink for U and source of Si. Fe cycling was variable, with terrestrial Fe largely removed as groundwater passed through the STE, consistent with Fe-sulfide and amorphous phase formation. Radium isotope ratios identified two distinct subsurface flow pathways, influenced by metal-oxide cycling and organic matter breakdown. Si production was decoupled from DIC, suggesting Si originates from mineral alteration, whereas DIC results from both mineral weathering and microbial respiration. Silicate mineral alteration, coupled with marsh pore water drainage, accounts for up to 16% of annual DIC exports (66 g C m −2 y −1 ), highlighting the importance of STEs in coastal carbon and trace element cycling, especially as marshes face environmental change.

Massachusetts↗
USGS2026

Selected water-quality data from the Cedar River and Cedar Rapids well fields, Cedar Rapids, Iowa, 2017–22

The Cedar River alluvial aquifer is the source of drinking water in Cedar Rapids, Iowa. Production wells are completed in the alluvial aquifer approximately 40 to 80 feet below land surface. The City of Cedar Rapids and the U.S. Geological Survey have studied the groundwater-flow system and water quality of the aquifer in the vicinity of Cedar Rapids since 1992. Results of these studies documented hydrologic conditions, water quality, and geochemistry of the alluvial aquifer and interactions with the Cedar River. Water-quality samples were collected for studies involving well field monitoring, trends, source-water protection, groundwater geochemistry, surface-water–groundwater interaction, and pesticides in groundwater and surface water. Water quality was analyzed for dissolved major ions (boron, bromide, calcium, chloride, fluoride, iron, magnesium, manganese, potassium, silica, sodium, sulfate, and total dissolved solids), dissolved nutrients (ammonia as nitrogen, ammonia plus organic nitrogen as nitrogen, nitrite plus nitrate as nitrogen, nitrite as nitrogen, orthophosphate as phosphorus, and phosphorus), dissolved organic carbon, and selected pesticides. Physical characteristics (alkalinity, dissolved oxygen, pH, specific conductance, and water temperature) were measured on site and recorded for each water sample collected. This report presents the results of routine water-quality data-collection activities from October 2017 through September 2022. Methods of data collection, quality assurance, water-quality analyses, and statistical procedures are presented. Data include the results of water-quality analyses from quarterly sampling from monitoring wells, production wells, two water treatment plants, and the Cedar River at Blairs Ferry Road at Palo, Iowa, streamgage (U.S. Geological Survey station number 05464420), as well as monthly nutrient sampling from the Cedar River and Morgan Creek near Covington, Iowa, streamgage (U.S. Geological Survey station number 05464475).

Iowa↗
USGS2025

Fiber-optic sensing for earthquake hazards research, monitoring and early warning

The use of fiber‐optic sensing systems in seismology has exploded in the past decade. Despite an ever‐growing library of ground‐breaking studies, questions remain about the potential of fiber‐optic sensing technologies as tools for advancing if not revolutionizing earthquake‐hazards‐related research, monitoring, and early warning systems. A working group convened to explore these topics; we comprehensively examined the application of fiber optics in various aspects of earthquake hazards, encompassing earthquake source processes, crustal imaging, data archiving, and technological challenges. There is great potential for fiber‐optic systems to advance earthquake monitoring and understanding, but to fully unlock their capabilities requires continued progress in key areas of research and development, including instrument testing and validation, increased dynamic range for applications focused on larger earthquakes, and continued improvement in subsurface and source imaging methods. A key current stumbling block results from the lack of clear data archiving requirements, and we propose an initial strategy that balances data volume requirements with preserving key data for a broad range of future studies. In addition, we demonstrate the potential for fiber‐optic sensing to impact monitoring efforts by documenting the data completeness in a number of long‐term experiments. Finally, we outline the features of a instrument testing facility that would enable progress toward reliable and standardized distributed acoustic sensing data. Overcoming these current obstacles would facilitate progress in fiber‐optic sensing and unlock its potential application to a broad range of earthquake hazard problems.

California, Oregon, Washington↗
USGS2024

Radiogenic strontium- and uranium-isotope tracers of water-rock interactions and hydrothermal flow in the Upper Geyser Basin, Yellowstone Plateau Volcanic Field, USA

Natural radiogenic isotopes (primarily 87 Sr/ 86 Sr) from hot springs in the Upper Geyser Basin of the Yellowstone Plateau volcanic field and associated rocks were used to evaluate groundwater flow patterns, water-rock reactions, and the extent of mixing between various groundwater sources. Thermal waters have very low uranium concentrations and 234 U/ 238 U activity ratios near 1.0, which limit their utility as tracers in this reducing setting. Thermal waters have higher Sr concentrations (<22 ng/g) and a wide range of 87 Sr/ 86 Sr values that vary both temporally at individual discharge sites and between adjacent springs, indicating that conduits tap different subsurface reservoirs to varying degrees. Sr from local rhyolites have 87 Sr/ 86 Sr compositions that bound the range of values observed in groundwater throughout the basin. Non-boiling springs on the west flank of the basin discharge water with low 87 Sr/ 86 Sr consistent with flow through young volcanic rocks exposed at the surface. Boiling springs in the central basin have higher 87 Sr/ 86 Sr values reflecting interactions with older, more radiogenic volcanic rocks. Variability in upwelling thermal waters requires mixing with a low 87 Sr/ 86 Sr component derived from young lava or glacial sediments, or more likely, from deeper sources of hot groundwater circulating through buried Lava Creek Tuff having intermediate 87 Sr/ 86 Sr. Isotope data constrain basin-wide output of thermal water to 110–140 kg·s −1 . Results underscore the utility of radiogenic Sr isotopes as valuable tracers of hydrothermal flow patterns and improve the understanding of temperature-dependent water-rock reactions in one of the largest continental hydrothermal systems on Earth.

Wyoming↗
USGS2026

Mid-Cretaceous porphyritic magmatism in the Beartooth Mountains of the northern Laramide foreland and its connection to Frontier Formation sediment dispersal in Bighorn basin, Montana-Wyoming (USA), and Cordilleran tectonism

Porphyritic intrusions in the southeast Beartooth Mountains and porphyritic cobbles from the Torchlight conglomerate of the Cenomanian Frontier Formation in Bighorn basin (Montana-Wyoming, USA) bear striking resemblance. We utilize geologic mapping, petrography, geochemistry, and geochronology to characterize rocks from both locations and test whether the Beartooth Mountains area was the provenance for Bighorn basin cobbles. Mapping on Line Creek Plateau in the Beartooth Mountains shows three units with porphyritic phenocrysts consisting dominantly of plagioclase, plagioclase + potassium feldspar, and plagioclase + quartz. Petrography shows embayed quartz, opacitic rims, and zoned plagioclase. Geochemistry characterizes the intrusions as trachydacite to trachyte and trace elements exhibit enrichment in light rare earth elements and slight depletion in heavy rare earth elements relative to primitive mantle. Zircon U-Pb dating shows substantial Archean inheritance and ca. 104−92 Ma dates; some samples display multiple mid-Cretaceous date populations. Bighorn basin porphyritic cobbles exhibit similar modal mineralogy, near-identical geochemistry, and Archean inheritance and mid-Cretaceous geochronology. We suggest the porphyritic cobbles were sourced from basement-hosted porphyritic intrusions now exposed in the Beartooth Mountains, given that geochemistry does not agree with coeval volcanic sources along the Frontier Formation depositional fairway (Idaho batholith) or sedimentary-hosted Beartooth Mountains porphyritic sills intruded into Cambrian shale. Stratigraphic evidence does not support kilometer-scale mid-Cretaceous Beartooth Mountains exhumation to expose basement-hosted porphyritic intrusions, and we therefore present a volcanic hypothesis to bring rocks to the surface. This mid-Cretaceous multimillion-year magmatism is unexplained by tectonic models for the north-central Laramide foreland, and thus we further explore the emplacement processes, spatiotemporal significance, and tectonic implications of these rocks.

Montana, Wyoming↗
USGS2026

Long-term intermittent connection between the western Snake River Plain and Columbia basin: A two-phased incision history of Hells Canyon

For more than a century, researchers have debated the route of the Snake River across the northern Cordillera and U.S. Pacific Northwest, including the associated incision history of Hells Canyon. Here, we use detrital zircon U-Pb provenance analysis of Miocene strata upstream and downstream of Hells Canyon to constrain the evolution in fluvial pathways over time. Downstream of Hells Canyon, we show that the ca. 10−6 Ma Clarkston Heights gravel was dominantly sourced from the nearby Salmon and Clearwater rivers but also with substantial input from drainages that now flow into the western Snake River Plain. These same western Snake River Plain drainages were contributing to time-equivalent Lake Idaho strata, upstream of Hells Canyon, indicating a Hells Canyon fluvial connection between the western Snake River Plain and Columbia basin before 6 Ma. Results from Lake Idaho strata suggest a two-phased lacustrine history in the western Snake River Plain: Before ca. 6.0 Ma, lake strata were derived from local tributaries, suggesting that the western Snake River Plain was isolated from the modern upper Snake River system. Between ca. 4.3 Ma and 2.2 Ma, the source area expanded to include eastern tributaries in conjunction with migration of high-standing topography of the Yellowstone hotspot. Our results challenge the long-held hypothesis of “capture” of the Snake River through Hells Canyon. Instead, we infer that Hells Canyon was a long-established route for outflow of lakes occupying the western Snake River Plain, possibly intermittently, during a ca. 10−2 Ma phase of slow regional incision, followed by rapid incision in Hells Canyon starting ca. 2 Ma in conjunction with erosion of the topographic barrier impounding Lake Idaho.

Idaho, Oregon↗