USGS Science⌕ Search

SEARCH · USGS Science

Results for “Map”

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.

452 records · Page 19Linked to original sources
USGS2024

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

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

Idaho↗
USGS2025

Shallow geologic framework of the Mississippi Sound and the potential for sediment resources

The Mississippi Sound, an estuarine environment located between the mainland and barrier islands bordering the northern Gulf of America (formerly the Gulf of Mexico), serves as a vital ecosystem for the States of Mississippi and Alabama. Spanning approximately 100 kilometers from east to west and covering 1,400 square kilometers, the sound is home to marine industry and ports, and its shallow and brackish waters sustain a diverse array of marine life. Barrier islands along the southern edge of the sound separate the microtidal estuary from the Gulf of America. This protection from gulf wave action mediates current flow within the sound, resulting in predominantly fine-grained sediment deposition along the seafloor. This study, conducted by the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers, provides insight on fluvial and tidal processes spanning the past 5,000 years. The report synthesizes existing research to provide a comprehensive overview of the sound geology, from Pleistocene origins to present-day morphology, and utilizes high-resolution single channel seismic profiles and sediment data to identify and map sedimentary deposits and morphologic features at and below the seafloor. Despite its ecological significance, the Mississippi Sound faces environmental challenges, including water-quality issues, habitat degradation, storm-induced erosion, and the ongoing threats of sea-level rise and environmental changes. This study uses the present-day understanding of the sound's geology to inform coastal management decisions, hazard assessment, and potential mineral resources.

Louisiana, Mississippi↗
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↗
USGS2024

Preliminary observations of the April 5th, 2024, Mw4.8 New Jersey earthquake

On 5 April 2024, 10:23 a.m. local time, a moment magnitude 4.8 earthquake struck Tewksbury Township, New Jersey, about 65 km west of New York City. Millions of people from Virginia to Maine and beyond felt the ground shaking, resulting in the largest number (>180,000) of U.S. Geological Survey (USGS) “Did You Feel It?” reports of any earthquake. A team deployed by the Geotechnical Extreme Events Reconnaissance Association and the National Institute of Standards and Technology documented structural and nonstructural damage, including substantial damage to a historic masonry building in Lebanon, New Jersey. The USGS National Earthquake Information Center reported a focal depth of about 5 km, consistent with a lack of signal in Interferometric Synthetic Aperture Radar data. The focal mechanism solution is strike slip with a substantial thrust component. Neither mechanism’s nodal plane is parallel to the primary northeast trend of geologic discontinuities and mapped faults in the region, including the Ramapo fault. However, many of the relocated aftershocks, for which locations were augmented by temporary seismic deployments, form a cluster that parallels the general northeast trend of the faults. The aftershocks lie near the Tewksbury fault, north of the Ramapo fault.

New Jersey↗
USGS2026

Decadal-scale characteristics of natural and anthropogenic dune morphology along North Carolina Barrier Islands (SE Atlantic Coast)

Dunes offer protection to coastal communities from storms, yet quantitative research characterizing human influences on dune systems across decades and 10–100 km spatial scales constitutes a known data gap. This study analyzed a multidecadal archive of dune crest elevations and positions to map mesoscale spatiotemporal trends in foredune morphology across the barrier island coastline of North Carolina, SE Atlantic coast. Results are summarized within four categories of coastline representing a range of coastal management intensities, integrating beach width and dune lateral accommodation space to assess the comparative influence of differing intensities of anthropogenic intervention on dune morphology. These categories account for the presence of structures, roads, and alongshore connectivity to anthropogenic sediment sources (beach nourishment). These categories show distinct dune and beach system characteristics. Dunes are largest in the areas with the longest legacy of coastal management. As categories decrease in development intensity, dune and beach system characteristics generally depart from the large, static dunes and narrow beaches consistent with hold-the-line management. Where small or transgressive dunes exist on heavily developed coasts, they coincide with relatively strong hydrodynamic forcing or sediment-limiting geology. There is evidence of indirect coastal management influence on undeveloped barrier islands with alongshore connectivity to developed barriers. While this analysis was conducted at the mesoscale, the results provide researchers with the framework to make insights at smaller scales than those demonstrated in this study.

North Carolina↗
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↗
USGS2024

RegionGrow3D: A deterministic analysis for characterizing discrete three-dimensional landslide source areas on a regional scale

Regional-scale characterization of shallow landslide hazards is important for reducing their destructive impact on society. These hazards are commonly characterized by (a) their location and likelihood using susceptibility maps, (b) landslide size and frequency using geomorphic scaling laws, and (c) the magnitude of disturbance required to cause landslides using initiation thresholds. Typically, this is accomplished through the use of inventories documenting the locations and triggering conditions of previous landslides. In the absence of comprehensive landslide inventories, physics-based slope stability models can be used to estimate landslide initiation potential and provide plausible distributions of landslide characteristics for a range of environmental and forcing conditions. However, these models are sometimes limited in their ability to capture key mechanisms tied to discrete three-dimensional (3D) landslide mechanics while possessing the computational efficiency required for broad-scale application. In this study, the RegionGrow3D (RG3D) model is developed to broadly simulate the area, volume, and location of landslides on a regional scale (≥1,000 km 2 ) using 3D, limit-equilibrium (LE)-based slope stability modeling. Furthermore, RG3D is incorporated into a susceptibility framework that quantifies landsliding uncertainty using a distribution of soil shear strengths and their associated probabilities, back-calculated from inventoried landslides using 3D LE-based landslide forensics. This framework is used to evaluate the influence of uncertainty tied to shear strength, rainfall scenarios, and antecedent soil moisture on potential landsliding and rainfall thresholds over a large region of the Oregon Coast Range, USA.

Journal of Geophysical Research: Earth Surface↗
USGS2025

Characterizing the scale of regional landslide triggering from storm hydrometeorology

Rainfall strongly affects landslide triggering; however, understanding how storm characteristics relate to the severity of landslides at the regional scale has thus far remained unclear, despite the societal benefits that would result from defining this relationship. As mapped landslide inventories typically cover a small region relative to a storm system, here we develop a dimensionless index for landslide-inducing rainfall, A * , based on extremes of modeled soil water relative to its local climatology. We calibrate A * using four landslide inventories, comprising over 11 000 individual landslides over four unique storm events, and find that a common threshold can be applied to estimate regional shallow-landslide-triggering potential across diverse climatic regimes in California (USA). We then use the spatial distribution of A * , along with topography, to calculate the landslide potential area (LPA) for nine landslide-inducing storm events over the past 20 years, and we test whether atmospheric metrics describing the strength of landfalling storms, such as integrated water vapor transport, correlate with the magnitude of hazardous landslide-inducing rainfall. We find that although the events with the largest LPA do occur during exceptional atmospheric river (AR) storms, the strength of landfalling atmospheric rivers does not scale neatly with landslide potential area, and even exceptionally strong ARs may yield minimal landslide impacts. Other factors, such as antecedent soil moisture driven by storm frequency and mesoscale precipitation features within storms, are instead more likely to dictate the patterns of landslide-generating rainfall throughout the state.

California↗
USGS2026

Multi-site evaluation of a postfire debris-flow runout forecast method

Postfire debris flows pose a hazard to human life, property, and infrastructure when they travel from steep source areas to urbanized alluvial fans or other developed areas. Existing methods for rapid (<1 week) postfire debris-flow hazard assessment document the increase in the likelihood and size of debris flows as the magnitude of high-intensity rain necessary to initiate debris flows increases but do not indicate the extent of downstream debris-flow runout. Although many models for the simulation of debris-flow motion are available, there is no established approach for using these models to delineate locations susceptible to postfire debris-flow runout that (a) is feasible to use at the spatial scale of an entire fire; (b) is appropriate for runout onto unconfined areas; (c) reproduces observed relations between runout and rainfall intensity; and (d) characterizes inherent uncertainty in runout, even without spatiotemporally variable rainfall. We propose and evaluate a method for generating postfire debris-flow runout hazard maps that has all the above qualities. Selection of case studies prioritized events triggered by a range of rainfall intensities, locations within and outside of southern California, and observed runout onto unconfined topography. Qualitative and quantitative assessment of performance for four events indicate that simulation results broadly match observations albeit with some discrepancies at a scale larger than structure or land parcel level (approximately 20-m by 20-m). The method may be used to identify potentially hazardous areas immediately following a fire and to provide approximate runout forecasts when a storm is imminent.

Arizona, California↗
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↗
USGS2024

Substrate Enhancement Pilot Project—Monitoring summary and evaluation, Kootenai River, Idaho, 2012–22

To assess changes in substrate conditions and the efficacy of artificially placed substrates at select sites on the Kootenai River near Bonners Ferry, Idaho, the U.S. Geological Survey, in cooperation with the Kootenai Tribe of Idaho, completed repeat bathymetric, velocimetric, and underwater videography surveys. Collectively, three project sites throughout the Kootenai River make up the Substrate Enhancement Pilot Project (SEPP), an effort intended to improve spawning and egg incubation viability at locations identified to be aquatic habitat limited for the endangered Kootenai River white sturgeon (Acipenser transmontanus). Following the placement of coarse substrates at each site, bathymetric, velocimetric, and underwater videography data were collected from 2012 to 2022 to assess the role of deposition and erosion on maintaining suitable white sturgeon spawning and incubation substrates. Minimal erosion and deposition occurred at all Substrate Enhancement Pilot Project sites, according to interannual and intra-annual net volumetric changes between bathymetric surveys. Depending on the timing of bathymetric surveys relative to the annual peak streamflow conditions, isolated locations of deposition or erosion were observed at each site and the potential for deposition or erosion was supported by measured mean depth-averaged velocities. This study concluded that variability of deposition and scour were common at each site throughout the monitoring period and may be attributed to fluctuations in streamflow. Repeat bathymetric, underwater videography, and velocity mapping surveys were used to verify the interstitial spaces and surfaces of substrates at each SEPP site remained free of fine sediments for intervals longer than a year but were susceptible to deposition between high streamflow events.

Idaho↗
USGS2025

Beyond the mangroves: A global synthesis of tidal forested wetland types, drivers and future information opportunities

There is increasing awareness of the global diversity of tidal forested wetlands (TFWs) and their significance in the provision of ecosystem services. These ecosystems, including mangrove forests, tidal freshwater forested wetlands, supratidal forests and transitional forests together span multiple climatic zones, geomorphic settings, and inundation and salinity regimes. We utilise case studies across five continents to demonstrate the state of knowledge among TFWs. Intertidal mangroves are the best-defined of the TFWs thanks to decades of research on their geomorphology, hydrology and ecology across their broad distribution. Non-mangrove forest settings, however, demonstrate more diverse hydrological, biochemical and vegetation conditions. In many cases, non-mangrove forests are situated at upper intertidal or supratidal elevations, where surface waters and groundwater are subject to interactions between tides freshwater inputs. Salinity datasets show variations ranging from tidal freshwater forested wetlands and ‘low-salinity mangroves’ to mesohaline or marine salinities, often with high temporal variability. While the floristic composition of non-mangrove forests vary among biogeographic regions, locally dominant TFW species are commonly distributed beyond the tidal niche into non-tidal wetland and upland forests. This presents challenges for traditional remote sensing approaches to ecosystem mapping, which are mostly lacking for non-mangrove forests. Geomorphic approaches and developments in machine learning offer opportunities to address this.

Preprint↗
USGS2026

Software to support remote sensing of river discharge based on critical flow theory

Water resource management requires accurate observations of streamflow but standard field methods for measuring river discharge ( Q ) are costly and can be hazardous for equipment and personnel. Remote sensing has become a viable alternative, but many image-based techniques require field data for calibration and depth and velocity can seldom be mapped with a single sensor. A new approach based on critical flow theory, in contrast, allows both of these attributes to be inferred from readily available image data. This technique only pertains to sites with standing waves, called undular hydraulic jumps (UHJs), but a recent investigation demonstrated its potential to provide accurate discharge estimates. This paper introduces software designed to facilitate Inferring Q from UHJs Identified in River Images (InQUIRI). The package includes modules for retrieving data from image servers, making the measurements of wavelength and width required to calculate discharge, inferring a representative wavelength from a profile digitized along a wave train, combining multiple estimates to obtain an ensemble median discharge, and assessing accuracy via comparison to gage records from the U.S. Geological Survey. By making these steps easier to implement, InQUIRI enables users to apply the workflow to a variety of UHJ-containing images. Accumulating more case studies, some successful and others less so, would help constrain the range of applicability of the critical flow approach and foster development of refined guidelines for selecting and measuring waves. The software described herein could play an important role in promoting informed use of this new technique for non-contact streamflow measurement.

Arizona, Colorado, New Mexico, Utah↗
USGS2024

Constraining mean landslide occurrence rates for non-temporal landslide inventories using high-resolution elevation data

Constraining landslide occurrence rates can help to generate landslide hazard models that predict the spatial and temporal occurrence of landslides. However, most landslide inventories do not include any temporal data due to the difficulties of dating landslide deposits. Here we introduce a method for estimating the mean landslide occurrence rate of deep-seated rotational and translational slides derived solely from high-resolution (≤3 m) elevation data and globally available estimates of the diffusion coefficient for sediment flux. The method applies a linear diffusion model to the roughest landslide deposits until they reach a representative non-landslide roughness distribution. This estimates the time for a landslide deposit to be unrecognizable in high-resolution digital elevation data, which we term the mean lifetime of the landslide. Using the mean lifetime and number of landslides within an area of interest, we can estimate the mean occurrence rate of landslides over that domain. We validate this approach using a comprehensive temporal inventory of landslides in western Oregon created using age-roughness curves that are calibrated with high-resolution elevation data and radiocarbon data. We find good agreement between our diffusion method and the existing age-roughness-derived estimates, producing mean lifetimes of 4500 and 5200 years (4% difference), respectively. Hazard maps produced using the two methodologies generally agree, with the maximum differences in landslide probability reaching 0.1. Due to the relative abundance of high-resolution elevation data compared with age-dated landslides, our method could help constrain landslide occurrence rates in areas previously considered unfeasible.

Oregon↗
USGS2026

Hazard potential of compound flooding from rainfall, storm surge, and groundwater in coastal New York and Connecticut

Compound flood events, the co-occurrence of multiple flood drivers, can result in flood hazard potential exceeding that of any single driver alone. To evaluate compound flooding in a semi-urbanized coastal area, historical records dating back to 1970 are used to study the co-occurrences of high precipitation, storm surge, and shallow groundwater conditions along the coastlines of New York and Connecticut. Joint return periods for coincident precipitation-surge events were computed using statistical dependence models and compared to the assumption of independence as a ratio, referred to here as a return period adjustment. Results indicate distinct seasonality where compound events in the area disproportionately occur in the cold season between October and April. Return period adjustments range from a factor of 1 to almost 9, demonstrating the range in precipitation-storm surge dependence across the study area. Across all 24 station triad locations, groundwater levels were elevated during times of precipitation- surge co-occurrence, reflecting the tendency for coastal storms and shallow groundwater conditions to co-occur seasonally. The result is a pseudo-trivariate compound flood hazard score and corresponding hazard map that integrates dependence between daily precipitation-surge events and overall monthly groundwater levels (as a precondition) into a relative compound hazard score. The location with the highest compound flood hazard score is on the south shore of Long Island, as well as locations across coastal Connecticut where groundwater levels compound the co-occurrence of heavy precipitation and storm surge.

Connecticut, New York↗
USGS2026

Monitoring changes in Landsat thermal features in urban and non-urban interfaces from 1986 to 2023 in two international urban centers: Implications for climate and global issues

Rapid urbanization is reshaping thermal environments worldwide, with the strongest impacts occurring at the interface between urban and non-urban areas. Impervious surfaces, as key indicators of urban expansion, are critical for monitoring urban growth and assessing surface urban heat island (SUHI) effects. Land use and land cover change (LULCC) provides an essential link between urban dynamics and their environmental and societal consequences. Here, we integrated the U.S. Geological Survey (USGS) Climate Global Issues (CGI) Land Cover Product with Landsat thermal time-series to investigate SUHI evolution in two contrasting metropolitan regions: Wuhan, China, and Brasília, Brazil. Using data spanning 1986–2023, we analyzed the relationships between land cover, Landsat-based land surface temperature (LST), and SUHI intensity, and identified persistent thermal hotspots. Results demonstrate that the land cover data utilized increases the accuracy of impervious surface mapping along urban–rural gradients. Average SUHI intensities were 3.4 °C in Wuhan and 3.3 °C in Brasília, with statistically significant warming trends of 0.04 °C/year and 0.01 °C/year, respectively. Maximum temperature proved to be a robust indicator of SUHI intensification, capturing long-term upward trends. Our findings highlight the important role of urban land cover dynamics in shaping temporal SUHI variability and hotspot emergence. This prototype framework demonstrates the scientific and policy value of combining long-term land cover monitoring information with satellite thermal monitoring to quantify and track SUHI at city scale, supporting sustainable urban planning and climate adaptation strategies.

Remote Sensing↗
USGS2026

The geologic history of the Chehalis Forearc Basin, Washington State, USA

The Chehalis basin is located between the Cascade arc and the Coast Range in southwest Washington State. It consists of sedimentary and volcanic rocks deposited throughout the Cenozoic and is underlain by the Siletzia terrane, a thick, large igneous province accreted in the Eocene. Here, we constrain evolution of the Chehalis basin depocenter using isochore maps derived from isostatic gravity anomalies, borehole data, and stratigraphy for several time periods: the Eocene (47.6–36.8 Ma), the Oligocene (36.8–20 Ma), and the Neogene (20–0 Ma). Our results suggest that local subsidence is driven by deformation on faults that bound and intersect the basin. We see northward depocenter migration, interpreted as a shift in deformation from the northwest striking Cedar Creek fault to the west striking Doty fault. We interpret these data in terms of the long-term north-south shortening and clockwise rotation of the Cascadia forearc. During the Eocene, the Cedar Creek fault was preferentially aligned with north-south shortening, but became less active when clockwise rotation brought the Doty fault into east-west alignment with the stress field. Sediment accumulation rates decreased from 196 m/Myr in the Paleogene (∼40–20 Ma) to 27 m/Myr in the Neogene (∼20–0 Ma) as Miocene uplift of the Coast Range limited accommodation space and shifted the depositional environment from marine to fluvial. Our results are consistent with the geologic evolution of the Portland and Tualatin basins to the south and reveal a depositional and structural history uniquely shaped by clockwise rotation of the Cascadia forearc.

Washington↗
USGS2024

High resolution identification and quantification of diffuse deep groundwater discharge in mountain rivers using continuous boat-mounted helium measurements

Discharge of deeply sourced groundwater to streams is difficult to locate and quantify, particularly where both discrete and diffuse discharge points exist, but diffuse discharge is one of the primary controls on solute budgets in mountainous watersheds. The noble gas helium is a unique identifier of deep groundwater discharge because groundwater with long residence times is commonly enriched in helium. In this study, a portable mass spectrometer was used to measure longitudinal variation in dissolved helium concentrations in two mountainous rivers at high spatial resolution not feasible with traditional sampling techniques. Helium profiles were then simulated using a mass-balance model to quantify longitudinal variation in groundwater discharge to the receiving rivers. Results indicate helium concentrations were enriched by multiple orders of magnitude above atmospheric equilibrium in both rivers and that this persisted for up to 18 km below observed pulse inputs in the Colorado River. Helium mass-balance models match observed longitudinal patterns with the exception of sharp initial increases in helium observed in the rivers. Increased longitudinal groundwater discharge rates correspond to mapped geologic structures in both watersheds that likely transport deep geothermal water. Models show variable sensitivity to spatial assignment of input variables representing the groundwater source, illustrating the importance of collecting data from discrete groundwater discharges where possible. The methodology shows promise for field experiments designed to assess air–water exchange rates and to quantify total groundwater discharge from a combination of discrete and diffuse sources.

Colorado, Utah↗