SEARCH · USGS Science

Results for “Environmental Modelling & Software”

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

SourceResearch connectionFrom yearThrough year

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

834 records · Page 30Linked to original sources

USGS2025

Sensitivity of benthic biota and toxicity of cadmium, cobalt, copper, nickel, lead, and zinc mixtures in Washington, United States, and British Columbia, Canada

Relative sensitivities and responses of juvenile white sturgeon ( Acipenser transmontanus ), Hyalella azteca , two families of mayfly (Ephemerellidae, Heptageniidae), one family of caddisfly (Brachycentridae), and a natural community of benthic macroinvertebrates (BMI) to multiple metals are predicted using previously collected laboratory and field samples and a metal mixture model. Biological responses in single metal exposures are used to parameterize toxicity functions, which include accumulations of hydrogen and selected metals on biological receptors, intrinsic potencies of hydrogen and metals, sensitivities of organisms, and times of exposure. The model then is used to predict responses in multiple metal laboratory exposures and field-collected porewater. The following sensitivity sequence in porewater was determined based on endpoints of survival or total abundance: juvenile white sturgeon greater than (>) Ephemerellidae family > Hyalella azteca > Heptageniidae family about equal to (~) benthic macroinvertebrate community > Brachycentridae family. The fraction of porewater samples that are predicted to have adverse impacts on benthic biota (20-percent or greater negative response) depends on organism sensitivities and metal toxicities, and ranges from 44 to 48 percent for juvenile white sturgeon, 23 to 26 percent for the Ephemerellidae family, 16 to 22 percent for Hyalella azteca , 5 to 8 percent for the Heptageniidae family and BMI community, and 0 percent for the caddisfly family. The most toxic porewater in the upper Columbia River Basin (UCR) is at the backwater bar site at Deadmans Eddy and China Bend. The model also indicates that the element responsible for the most toxic conditions in UCR porewater is copper for all organisms, except Hyalella azteca and the metal-insensitive Brachycentridae family. Copper and lead result in the most toxic conditions for Hyalella azteca . This approach and results can aid in assessing metal toxicity and its potential risk to aquatic biota in ecosystems impacted by historical mining activities.

British Columbia, Washington↗

USGS2024

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

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

California↗

USGS2026

Scenarios to assess the future water availability in the Mississippi River Valley Alluvial Aquifer for the Cache River and Grand Prairie Regions of Arkansas

The U.S. Geological Survey, as part of the Arkansas Groundwater Initiative, developed forecast scenarios using previously calibrated MODFLOW 6 groundwater models that focused on the Cache and Grand Prairie Critical Groundwater Areas to assess the impact of future climate and water management strategies on the Mississippi River Valley alluvial aquifer. A Soil Water Balance model was used to forecast recharge and irrigation water use. The forecast scenario period was from January 1, 2019, through December 31, 2055, with monthly stress periods. Twenty scenarios were simulated and included seven alternate climate forecasts, five 13 general groundwater pumping reduction scenarios (round 1), and groundwater pumping reduction scenarios by crop type and for the Bayou Meto Water Management Project and Grand Prairie Area Demonstration Project (round 2). Declines in saturated thickness within the Cache Critical Groundwater Area were larger for 18 of the 20 scenarios as compared to outside of the Critical Groundwater Area. The largest average increase in saturated thickness inside the Critical Groundwater Area was 6.4 m which occurred for the round 1, 50 percent reduction scenario. Automatic reductions in groundwater pumping by MODFLOW 6 in the Cache simulation ranged from 0.02 to 13.1 percent of total groundwater pumping. For the Grand Prairie model domain, the average change in saturated thickness of the Mississippi River Valley alluvial aquifer inside the Critical Groundwater Area for the forecast period ranged between -6.6 to 1.7 m. The average saturated thickness of the Mississippi River Valley alluvial aquifer inside the Grand Prairie Critical Groundwater Area declined for 16 of the 20 scenarios. The average reduction in requested groundwater pumping for all scenarios inside the Grand Prairie Critical Groundwater Area was 25.1 percent, and the largest reduction was 46.5 percent.

ESS Open Archive↗

USGS2025

Invited perspectives: Integrating hydrologic information into the next generation of landslide early warning systems

Although rainfall-triggered landslides are initiated by subsurface hydro-mechanical processes related to the loading, weakening, and eventual failure of slope materials, most landslide early warning systems (LEWS) have relied solely on rainfall event information. In previous decades, several studies demonstrated the value of integrating proxies for subsurface hydrologic information to improve rainfall-based forecasting of shallow landslides. More recently, broader access to commercial sensors and telemetry for real-time data transmission has invigorated new research into hydrometeorological thresholds for LEWS. Given the increasing number of studies across the globe using hydrologic monitoring, mathematical modeling, or both in combination, it is now possible to make some insights into the advantages versus limitations of this approach. The extensive progress demonstrates the value of in situ hydrologic information for reducing both failed and false alarms, through the ability to characterize infiltration during, as well as the drainage and drying processes between major storm events. There are also some areas for caution surrounding the long-term sustainability of subsurface monitoring in landslide-prone terrain, as well as unresolved questions in hillslope hydrologic modeling, which relies heavily on the assumptions of diffuse flow and vertical infiltration but often ignores preferential flow and lateral drainage. Here, we share a collective perspective based on our previous collaborative work across Europe, North America, Africa, and Asia to discuss these challenges and provide some guidelines for integrating knowledge of hydrology and climate into the next generation of LEWS. We propose that the greatest opportunity for improvement is through a measure-and-model approach to develop an understanding of landslide hydro-climatology that accounts for local controls on subsurface storage dynamics. Additionally, new efforts focused on the subsurface hydrology are complementary to existing rainfall-based methods, so leveraging these with near-term precipitation forecasts is a priority for increasing lead times.

Alaska↗

USGS2024

Simulating present and future groundwater/surface-water interactions and stream temperatures in Beaver Creek, Kenai Peninsula, Alaska

In many places, coldwater ecosystems are facing increasing pressure from anthropogenic warming. This study examined stream temperatures and the water balance in the Beaver Creek watershed on the Kenai Peninsula in south-central Alaska—an area that is experiencing rapid warming. Low-gradient streams near the Kenai coast provide important spawning and rearing habitat for salmon but may be especially vulnerable to rising temperatures, because of long residence times, inflows from abundant riparian wetlands, and reliance on groundwater discharge that may also warm, or decrease in volume with rising evapotranspiration. In recent decades, observed maximum 7-day temperatures have consistently exceeded statistical (regression-based) projections. Here we simulate total streamflows and temperatures with a physics-based model that links the Soil Water Balance, MODFLOW 6 and SNTEMP simulation codes on a 7-day timestep. The model is based on existing data and groundwater levels, instream flows, and stream temperatures collected during 2019–23. Future climate scenarios were developed for 2023–50 from downscaled climate projections. Results indicate that groundwater discharge is about 64 percent of the total streamflow during the months of May through September. Total streamflow and groundwater discharge are expected to remain similar to current conditions through 2050. Stream temperatures are expected to rise; by midcentury, near the Beaver Creek mouth the model predicts 34 to 63 additional days per year with average weekly temperatures above 13 degrees Celsius, 14 to 81 additional days with average weekly temperatures above 15 degrees Celsius, and routine exceedances of 20 degrees Celsius during the warmest periods. Projected stream temperatures vary spatially. Areas of high groundwater inflows in the lower main stem and some tributaries may be most resilient to warming air temperatures during dry conditions. During storm events, groundwater-dominated tributaries may have the coolest stream temperatures.

Alaska↗

USGS2026

Flood-inundation maps of the Blue River and major tributaries in and near Kansas City, Missouri, 2023–25

Digital flood-inundation maps for 35.5 miles of the Blue River, in and near Kansas City, Missouri, were created by the U.S. Geological Survey (USGS) in cooperation with the City of Kansas City. Maps were also developed for a combined 7.3 miles of Brush Creek (4.4 miles) and Indian Creek (2.9 miles), two primary tributaries to the Blue River in the study reach. The flood-inundation maps, available through the USGS Flood Inundation Mapping Program website at https://water.usgs.gov/osw/flood_inundation/ , depict estimates of the areal extent and depth of flooding, depth of water above first-floor elevations in selected buildings, water velocity, water-surface elevations, and time of travel corresponding to selected water levels or gage heights (hereafter referred to as “stages”) at nine reference USGS streamgages—five on the Blue River (USGS station numbers 06893590, 06893578, 06893500, 06893150, and 06893100), two on Brush Creek (USGS station numbers 06893562 and 06893557), one on Indian Creek (USGS station number 06893390), and one on the Missouri River (USGS station number 06893000). Near-real-time stages at these streamgages may be obtained from USGS Water Data for the Nation at https://doi.org/10.5066/F7P55KJN or the National Weather Service National Water Prediction Service at http://water.noaa.gov/ , which also forecasts flood hydrographs at seven of these sites (USGS station numbers 06893000, 06893590, 06893578, 06893500, 06893150, 06893557, and 06893390). Flood profiles were computed for eight map reaches by means of two-dimensional hydraulic models. The models were calibrated using the USGS-developed stage–streamflow relations at each reference streamgage. Two map reaches on the lower Blue River also include the effects of backwater from the Missouri River on flood inundation extent and hydraulic characteristics. The hydraulic models were used to compute water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datums. The profile stages ranged from the National Weather Service “Action Stage” or near bankfull to a stage exceeding the highest recorded water level at each streamgage. The simulated water-surface profiles were then combined with a digital elevation model (derived from light detection and ranging data having a nonvegetated vertical accuracy of a maximum 10-centimeter root-mean-square error) to delineate the area flooded at each water level and the associated hydraulic characteristics. The availability of these maps, along with information regarding current stage from the USGS streamgage and forecasted high-flow stages from the National Weather Service, will provide emergency management personnel, resource managers, and residents with information that could be critical for flood-response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Kansas, Missouri↗

USGS2025

Favorability mapping for hydrothermal power resource assessments of the Great Basin, USA

The U.S. Geological Survey (USGS) is updating the 2008 assessment of conventional hydrothermal resources for the Great Basin in the western United States. As part of this work, the workflow for hydrothermal resource favorability maps is being modified to integrate modern data-driven machine learning (ML) methods. Improvements include: [1] using new and refined evidence layers (features); [2] using an order of magnitude more training sites (labeled examples); [3] utilizing simple but non-linear supervised ML algorithms; [4] representing positive training sites (wells with measured heat flow) with their ordinal value proportional to the magnitude of convective upflow (i.e., low, high, or very high convective signals instead of past strategies using positive-negative labels); [5] supplementing training sites with additional sites with low convective signals to represent diverse under-sampled areas where hydrothermal systems are unlikely to exist; [6] comparing with competing approaches; and [7] utilizing Monte Carlo cross-validation to estimate and evaluate prediction uncertainty. For the new favorability map, over half of the power-producing systems (i.e., 15 of 28) are predicted in the 99th percentile of most favorable locations (i.e., the highest 1 % of favorability, corresponding to 1 % of the map area), exceeding the performance of past models that have explicitly used power plants as training sites. Previous favorability maps predicted approximately half of the power-producing hydrothermal systems above the 80th percentile (i.e., 20 % of the map area). For the new favorability map, 93 % of power-producing systems (i.e., 26 of 28) are above the 80th percentile. The power-producing systems for which the new model does not perform well are either comparatively small, low-temperature systems or systems also not predicted well by prior modeling approaches, suggesting that these few systems are unusual when compared with most power-producing systems. Focusing research on these known, seemingly different systems may yield new insights and subsequent discovery of new prospects.

California, Idaho, Nevada, Oregon, Utah↗

USGS2024

Conceptualization and simulation of groundwater flow and groundwater availability in the Boone and Roubidoux aquifers in northeastern Oklahoma, 1980–2017

Oklahoma Groundwater Law (Oklahoma Statute § 82-1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations to determine the maximum annual yield for the State’s groundwater basins. The Boone and Roubidoux aquifers (also known as the Springfield Plateau aquifer and Ozark aquifer, respectively) are bedrock aquifers that extend from northeastern Oklahoma into Kansas, Arkansas, and Missouri. At present (2024), the Oklahoma Water Resources Board has yet to legally issue orders for the final determination of maximum annual yields for the Boone and Roubidoux aquifers. To support determination of a maximum annual yield, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, developed a hydrogeologic framework, a conceptual groundwater-flow model, and a calibrated numerical groundwater-flow model for the Boone and Roubidoux aquifers. Three types of groundwater-availability scenarios were simulated by using the calibrated numerical model. These scenarios were used to (1) estimate equal-proportionate-share groundwater withdrawal rates (groundwater withdrawal applied equally over the aquifer), (2) quantify the potential effects of projected groundwater withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical 10-year drought. For the Boone aquifer, equal-proportionate-share groundwater withdrawal rates were 1.10, 0.98, and 0.96 acre-feet per acre per year for the 20-, 40-, and 50-year scenarios, respectively. For the Roubidoux aquifer, equal-proportionate-share groundwater withdrawal rates were 1.76, 1.34, and 1.25 acre-feet per acre per year for the 20-, 40-, and 50-year simulations, respectively. For the 50-year scenarios, stream seepage was minimally affected. Over the 10-year drought scenario, groundwater storage in the Boone and Roubidoux aquifers decreased by 660,451 acre-feet (6.7 percent) and 508,472 acre-feet (1.0 percent), respectively.

Oklahoma↗

USGS2025

Magnitude, depth and methodological variations of spectral stress drop within the SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest Earthquake Sequence

We present the first ensemble analysis of the 56 different sets of results submitted to the ongoing Community Stress Drop Validation Study using the 2019 Ridgecrest, California, earthquake sequence. Different assumptions and methods result in different estimation of the source contribution to recorded seismograms, and hence to the source parameters (principally corner frequency, f c ⁠ , spectral stress drop, Δσ, and seismic moment, M 0 ⁠ ) obtained from modeling calculated source spectra. For earthquakes smaller than magnitude (M) 2.5 there is negligible correlation between the f c values obtained by different studies, implying that no present method is reliable using available data. For larger magnitude events, correlation between f c measurements of different studies, within even a small M range is always higher than spectral ⁠Δσ , because the f c measurements simply reflect the underlying physical decrease in f c with increasing M. We model the observed trends of submitted f c with both magnitude and depth. Most methods report an increase in spectral Δσ with M, although a magnitude‐invariant spectral Δσ is within the confidence limits. The depth dependence is smaller and depends on whether a study allows attenuation to vary with source depth; a combination of depth‐dependent attenuation correction, and depth‐dependent shear‐wave velocity can compensate for reported depth trends. We model the submitted values to remove differing M and depth variation to investigate the relative interevent variability. We find consistent relative variation between individual events, and also lower relative spectral Δσ in the northwest of the aftershock sequence, and higher on the cross fault and in the region of main fault intersection. This large‐scale comparison implies that absolute spectral Δσ estimates are dependent on the methods used; studies of different regions or using different methods should not be directly compared and improved constraints on path and site corrections are needed to resolve these absolute spectral Δσ differences.

California↗

USGS2024

The impact of source time function complexity on stress drop estimates

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

Bulletin of the Seismological Society of America↗

USGS2025

Slow rupture, long rise times, and multi-fault geometry: The 2020 M6.4 southwestern Puerto Rico mainshock

The M 6.4 mainshock of the southwestern Puerto Rico seismic sequence on 7 January 2020, was one of the most impactful modern earthquakes in the northeastern Caribbean. Due to its offshore location and complex aftershock distribution, its source kinematics remain poorly constrained. This active sequence illuminated a complex set of previously unrecognized structures that indicate multiple causative faults may have slipped during its rupture. Here, we utilize seismic and geodetic observations to enhance model resolution, estimate the finite slip of the mainshock, and test a multi-segment, geologically realistic fault geometry. Our refined model finds a lower rupture velocity and longer rise times than typical for an event of this magnitude. This indicates a slow-evolving rupture process that resembles characteristics of a tsunami earthquake. Although this normal/strike-slip faulting event was not tsunamigenic, these qualities, if pervasive for this region, have important implications for future seismic monitoring and hazards in southwestern Puerto Rico.

Puerto Rico↗

USGS2025

Rebuilding a volcano one lava flow at a time—Visualizing the lava dome-building eruption in the crater of Mount St. Helens, 1982–1986

Between 1980 and 1986, the U.S. Geological Survey made a series of 1:2,000-scale topographic contour maps from aerial photographic surveys to monitor the eruption. These maps were made for operational purposes and were not intended for publication. Since then, advances in technology made it possible to digitize the original, highly detailed hardcopy maps and derive new digital data elevation models of the surface of the lava dome. These digital elevation models allow for the visualization of the progression of the eruption and reveal the rubbly, chaotic surface of the lava flows and dome. Additionally, these new data help fill gaps in the long-term record of topographic changes that have occurred at the volcano since the May 18, 1980, eruption.

Washington↗

USGS2025

Avian navigation: Comparing the olfactory navigational “map” and the infrasound direction-finding hypotheses to aeronautics

Animal navigation has long been a fascinating but bewildering subject. Humans and animals might well share similar navigational strategies because they developed within the same physical environments. A “map-and-compass” model has been proposed to explain the two-step avian navigational process, but the “map” step has remained elusive. Although scalar values from bicoordinate geomagnetic or atmospheric olfactory gradients have been considered foundational to the avian map, neither has proved convincing engendering decades of controversy. The olfactory map, and an alternative infrasound direction-finding (IDF) hypothesis, are discussed in this review. The olfactory map hypothesis currently requires extensive stable gradients of trace-odor ratios, but such gradients are highly unlikely within a turbulent and rapidly mixed lower atmosphere. The IDF hypothesis, on the other hand, postulates a two-step navigational model analogous to the maritime and aeronautical radio direction-finding technique. This review was also written to encourage further investigation, and direct testing, of the acoustic navigational process. The IDF hypothesis, at present, appears the better explanation of observed avian navigational behavior and accuracy within the atmosphere’s physical environment.

Journal of Comparative Physiology A↗

USGS2024

Assessment and characterization of ephemeral stream channel stability and mechanisms affecting erosion in Grand Valley, western Colorado, 2018–21

The Grand Valley in western Colorado is in the semiarid Southwest United States. The north side of the Grand Valley has many ungaged ephemeral streams, which are of particular interest because (1) the underlying bedrock geology, Late Cretaceous Mancos Shale, is a sedimentary rock deposit identified as a major salinity contributor to the Colorado River and (2) despite infrequent streamflows of short duration, monsoon-derived floods in these ephemeral streams can carry substantial amounts of sediment downstream, affecting upstream and downstream banks and channel cross sections. The study area is of interest, because salinity, or the total dissolved solids concentration, in the Colorado River causes an estimated $300 million to $400 million per year in economic damages in the United States, and it is estimated 62 percent of the Upper Colorado River Basin’s total dissolved solid loads originate from geologic sources. In an effort to minimize salt contributions to the Colorado River from public lands administered by the Bureau of Land Management, a comprehensive salinity control approach is typically used to reduce nonpoint sources of salinity through land management techniques and practices. In 2018, the U.S. Geological Survey, in cooperation with the Bureau of Land Management, began an assessment of ephemeral streams located on the north side of the Grand Valley, western Colorado, to characterize stream channel stability and identify mechanisms affecting erosion. The U.S. Geological Survey developed a method for automatically extracting channel cross-section geometry from existing remotely sensed terrain models. Based on estimated flood stage and surrogate streamflows, hydraulic characteristics were calculated. Furthermore, the channel geometries and hydraulic characteristics were used to estimate channel stability using a statistical model. Cross-section stabilities were determined from a stream channel stability assessment for a subset of 1,406 visited (field observed) locations out of 13,415 cross sections, which were delineated from remotely sensed terrain models. The application of Manning’s resistance equation in combination with multiple logistic regression models demonstrated channel stability can be estimated with a 0.845 goodness of fit for a validation dataset when using a combination of drainage area, width-to-depth ratio, sinuosity, and shear stress as the explanatory variables. Stream channel stability was extrapolated for 13,415 unvisited (not field observed) cross sections using the multiple logistic regression model and defined explanatory variables. Mapping of the ephemeral streams and their associated stabilities may be used by the Bureau of Land Management to prioritize areas for remediation or changes in management strategies to reduce sediment and salinity loading to the Colorado River. The study found channel stability within the ephemeral streams to be spatially variable, longitudinally discontinuous, and dictated by changes in channel bed slope. The stable ephemeral streams were relatively wide and shallow and often had smaller drainage areas with less potential for producing shear stresses capable of overcoming channel adhesion. A change in channel bed slope can provide the means necessary to generate shear stresses appropriate to initiate erosion and a subsequent stability transition to incising channels. Channel widening happens when either or both banks of an incising channel reach a critical height for mass wasting, or when channel curvature causes higher sidewall stress. Regardless, widening channels can promote increases in sinuosity and subsequently reduce steep channel bed slopes. Consequently, stable and widening channels can have comparable bed slopes, making channel bed slope a poor explanatory variable to predict channel stability overall, despite its function to initiate channel instability. The results were based on a surrogate 0.10 annual exceedance probability (AEP; return period equal to the 10-year flood) interval streamflow, although it was recognized fluctuations in streamflow would also affect channel stability. Past and current changes within the study area affect streamflow; therefore, mechanisms affecting erosion include land use disturbances, soil compaction, loss of vegetation cover, drought, less frequent and more extreme precipitation, and fires—which all intensify the potential runoff and erosion within the study area.

Colorado↗

USGS2026

Emergent headwater types for watershed monitoring, research, and planning

Headwaters—that is headwater streams and the watersheds they drain—strongly influence water quality, ecosystem services, and hydrologic connectivity, yet remain poorly mapped, monitored, and understood at large spatial scales. Using available spatial data to classify headwater “types” could therefore enhance the utility of sparse monitoring data to advance research and watershed planning—and help update literature-based perceptual models of headwaters as steep, forested, and minimally disturbed watersheds. Here, we analyze headwaters across the contiguous United States to quantify how their attributes (climate, topography, land cover, soils) vary regionally and differ from their downstream watersheds. Our analysis (a) identified seven dominant headwater types with distinct combinations of biophysical and climatic attributes and (b) revealed systematic contrasts with downstream watersheds, with strongest differences in attributes representing potential flowpath depth. Our findings offer a transferable framework for classifying, modeling, and managing these critical yet understudied components of the hydrologic network.

conterminous United States↗

USGS2026

A case study of maximum depositional ages from terrestrial sandstones near the Cretaceous—Paleogene transition, western Williston Basin, USA

We present a new Bayesian method for deriving maximum depositional ages from detrital K-feldspar using total-fusion 40 Ar/ 39 Ar data. Individual analyses weighted by K/Ca ratio, age uncertainty, and percent radiogenic argon prioritize a result to come from accurately measured volcanic sources. Results from this method applied to sandstones from the Hell Creek region of northeastern Montana show that detrital K-feldspar maximum depositional ages align closely with detrital U-Pb zircon ages as well as tephra based chronostratigraphic constraints from both 40 Ar/ 39 Ar and U-Pb data. An age model informed by maximum depositional ages and available tephra data provides an estimate of 1.74 ± 1.04 Ma for the duration of deposition for the Hell Creek Formation. Combined age model and maximum depositional age data suggest ≤ 1.84 Ma of missing time is represented by the Hell Creek and Fox Hills formational disconformity, coinciding with the final regression of the Western Interior seaway.

Montana↗

USGS2025

Using probability difference to compare streamflow information of alternatives for efficient operation of monitoring networks

Efficient operation of streamflow monitoring networks requires investments in technology and labor that provide the greatest benefits from available resources. Economic analyses comparing the costs and benefits from different types of alternatives for monitoring have not been practical to implement. Streamflow information provides a generic measure of benefits that can be incorporated into operational decisions as an objective for monitoring networks. A methodology for comparing how accuracy, monitoring period, and monitoring instead of modeling affects streamflow information is developed from information-theoretic approaches for network design but contributes three novel features: (1) a probability-difference model for conditional probability of monotonically paired variables, (2) explicit discounting of unverified information that may exceed the accuracy of streamflow records, and (3) run analysis to account for non-stationarity in streamflow probabilities. Application of the methodology to the U.S. Geological Survey streamflow monitoring network indicates the value of monitoring period to reduce the uncertainty of streamflow probabilities and, thus, increase streamflow information. The methodology has important limitations, particularly for sites with non-perennial streamflow, but demonstrates that probability difference could be used to evaluate operational alternatives to increase the efficiency of monitoring networks.

PLOS Water↗

USGS2025

Comparison of hydrologic data and water budgets between 2003–08 and 2018–23 for the eastern part of the Arbuckle-Simpson aquifer, south-central Oklahoma

The Arbuckle-Simpson aquifer is divided spatially into three parts (eastern, central, and western). The largest groundwater withdrawals are from the eastern part of the Arbuckle-Simpson aquifer, which provides water to approximately 39,000 people in Ada and Sulphur, Oklahoma, and surrounding areas. The Arbuckle-Simpson aquifer, including the eastern part, is designated a sole source aquifer for its service area. Based primarily on data collected between 2003 and 2008, a series of comprehensive hydrologic studies of the Arbuckle-Simpson aquifer was published to provide the information necessary to perform groundwater-flow model simulations so that the Oklahoma Water Resources Board could determine how much water could be withdrawn from the aquifer while maintaining flow to springs and streams. As part of the Phase 1 studies, an aquifer water budget was developed from a numerical model for the period 2003–08. For this report, Phase 1 refers to the 2003–08 data collection period, although for some of the analyses, data collected prior to 2003 were used to inform model development work. Allocation of water from this aquifer was then established by the Oklahoma Water Resources Board in 2013. Additional well-spacing rules were also established by the Oklahoma Water Resources Board for sensitive sole source groundwater basins. To determine how the water budget for the eastern part of the Arbuckle-Simpson aquifer has changed over time, recently collected hydrologic data (2018–23) were compared to data collected during 2003–08. The analysis of changes in the aquifer water budget from 2003–08 to 2018–23 could help resource managers better understand changes in the overall balance of water in storage and the potential effects on streamflow, changes in groundwater levels, and the effects of different water uses in the aquifer area on available water in the eastern part of the Arbuckle-Simpson aquifer and streams overlying the eastern part of the Arbuckle-Simpson aquifer.

Oklahoma↗