Geology topics

Browse USGS reports

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

At least 757 records · Page 42Linked to original sources

USGS2026

Flowering and flower maturation in valley oak Quercus lobata

We investigated factors driving flowering and flower maturation and their effects on acorn production at both the annual and individual trees levels in a population of valley oak ( Quercus lobata ) in central coastal California, USA. Among years, female flower production was negatively correlated with rainfall, while warm spring temperatures and a large acorn crop the previous year decreased the proportion of flowers maturing into acorns. Among trees, individuals with greater access to ground water growing in warmer microclimates produced more catkins and more female flowers. In addition, male and female flower abundance were positively correlated even after controlling for tree size and resources, countering the hypothesis of a life-history tradeoff between investment in male and female flowering. Factors limiting pollen availability include weather-driven effects (both the pollination Moran effect and environmental vetos), phenological synchrony, and pollen coupling, factors that were difficult to distinguish. All three flower indices explained a significant proportion of variance in acorn production at both the annual and individual tree levels. Most important, however, was variation in female flower maturation, while both female and male flower abundance explained between one-quarter and one-third of variance in acorn production. Quercus lobata is not only a flower maturation masting species, but also a female flower and male flower masting species. The underlying drivers of flower initiation, fertilization, and flower maturation yield important insights into the factors influencing masting behavior in perennial plants and represent an area that can benefit from more empirical study.

California↗

USGS2026

Assessment of extreme subsurface hydrologic conditions captured during atmospheric river storms in the San Francisco Bay area (California, USA) with applications to shallow landslide early warning

An increase in soil pore water pressure is the typical trigger for the majority of landslides caused by rainfall. Atmospheric river storms, common to the west coast of North America during the winter season, can deliver landslide triggering rainfall resulting in severe impacts to coastal communities. Using a network of hydrologic monitoring stations situated within landslide-prone terrain in the San Francisco Bay area of California (USA), we assess the meteorologic conditions and resulting hydrologic and landslide response resulting from eight consecutive storm events that caused thousands of shallow landslides during the winter of 2022–2023. We find disparate hydrological responses and resultant degrees of observed landsliding ranging from < 1 landslide/km2 to 18 landslides/km2 at the monitoring sites that reflect the interplay and differences between rainfall delivery, subsurface hydrological characteristics, and geotechnical properties at each site. Antecedent soil moisture from both early season rainfall and the first storm in the sequence played a critical role in setting up some hillslopes for failure. Subsequent storms then generated elevated pore water pressures for several hours with associated landsliding. However, we find that the occurrence of widespread landsliding required not only sufficient pore water pressure magnitude in susceptible hillslopes, but also full and prolonged effective soil saturation throughout hillslope profiles. Landslides may still occur at lower values and durations of effective saturation but are likely to be less extensive regionally. We present these findings within the context of research directions and improvements to landslide early warning systems first suggested by researchers 40 years ago.

California↗

USGS2026

Geothermal district energy systems coupled with seasonal underground thermal energy storage: A U.S. techno-economic screening by climate and geology

In the United States, cooling-dominated commercial building loads can cause geothermal heat pump-based district energy systems to accumulate a long-term subsurface thermal imbalance, motivating the incorporation of seasonal underground thermal energy storage. We developed a transferable workflow to evaluate geothermal district systems that pair ground heat exchangers with seasonal underground thermal energy storage. Using standardized hourly loads for seven commercial buildings and a uniform cost framework, we simulated ten U.S. cities with a physics-based ground heat exchanger model, subsurface storage simulations, and economic assessment to isolate the roles of climate and hydrogeology. In cooling-dominated cities, underground thermal energy storage supplied the majority of annual cooling, cutting electricity use and summer peaks substantially while achieving levelized costs comparable to or below conventional chiller-boiler plants. In cooler climates, the storage share shrunk, required borefield size and costs rose, and levelized cost of energy increased nearly linearly with declining underground thermal energy storage fraction, indicating storage fraction as the primary economic lever. Sensitivity analysis showed capital risk dominated by borefield drilling and surface heating, ventilation, and air-conditioning and piping, with underground thermal energy storage costs secondary. This workflow provides a transparent foundation for site-specific design and screening of next-generation geothermal district energy systems.

Arizona, Illinois, Maryland, Michigan, Mississippi↗

USGS2026

Society’s struggle with unwanted species: What centuries of history reveal about species bounty programs

Species bounty programs, much like bounty hunters charged to bring fugitives to justice, enlist the public to locate and remove unwanted species through financial incentives. With the goal of reducing population sizes, these programs address perceived ecological and economic damage caused by target species. In this study, we provide the first global assessment of species bounty programs, drawing on evidence from both historical and contemporary efforts across diverse regions and cultural contexts over the past eight centuries. We uncovered a long history of bounty programs involving at least 283 species—mammals, birds, fish, plants, reptiles, mollusks, insects, amphibians, and crustaceans—across 449 programs in 60 countries. Using this collective knowledge, we offer five perspectives on species bounty programs. First, bounty programs are launched for a variety of reasons, including economic (livestock, crops, fisheries, infrastructure), ecological (species, ecosystems), and social (human health) considerations related to unwanted species. Second, bounty programs vary in their design and implementation, ranging from well-planned operations with clear management and conservation objectives to ad hoc operations with limited articulation and investigation of project outcomes. Third, evidence points to unintended consequences, in which bounty programs result in the incidental removal of non-target species or in effects that may inadvertently benefit target species. Fourth, while not always the case, fraudulent activities have been reported, compromising the management outcomes of some programs. Fifth, public perception of bounty programs is highly dynamic and ensuring program engagement remains a persistent challenge. By reviewing the scattered narratives of past and present bounty programs globally, this review seeks to inform the evolving role of this management strategy.

npj Biodiversity↗

USGS2026

Mangrove microbiomes as drivers of ecosystem recovery and restoration success

The microbes found in the rhizosphere, roots, leaves and stem surfaces and within the internal tissues of mangrove vegetation and their environment constitute the microbiome of the ecosystem. The organisms in the microbiome include bacteria, protozoa, fungi, algae, amoebas, and slime molds, which assist in maintaining and restoring mangrove ecosystems. This review explores the role of microbiomes in the maintenance of healthy mangrove ecosystems and in the successful restoration of degraded mangrove ecosystems. Microbes have important roles in several geomicrobiological cycles shaping mangrove ecosystems, including transforming nitrogen, phosphorus, carbon, sulfur and iron in biogeochemical cycles. Mangrove microbiomes contribute to the adaptation of vegetation to the harsh abiotic conditions in coastal areas, enhance nutrient uptake, produce plant-growth-promoting substances, and degrade the mangrove litter and the pollutants that can hinder restoration. Soil microbes function as biofertilizers, biopesticides, and bioremediation agents. The microbial diversity, composition, and functional capacity are important in the restoration of mangroves through their influence on voluntary recruitment following hydrologic restoration, on the establishment success of planted seeds and propagules, and on the survival of transplanted saplings and nursery-raised seedlings. The knowledge of the beneficial attributes of the microbiome can enhance the overall success of mangrove restoration. Identifying future needs, such as microbial inoculant validation, field-scale trials, and integration with hydrological restoration, are essential.

Microorganisms↗

USGS2026

Evaluating the U.S. Geological Survey’s earthquake shaking hazard forecasts and their implications on seismic risk

We analyze the last six update cycles of the long-term probabilistic earthquake shaking hazard forecast from the U.S. Geological Survey (USGS) and discuss the changes in hazard estimates from the 1996 to the latest 2023 update of the National Seismic Hazard Model (NSHM) for the conterminous U.S. We summarize how our understanding of earthquake hazards has evolved over the last two decades and quantify the implications of changing hazards on people, buildings, and lifeline infrastructure. The net effect of changes in hazard estimates, with reference to the mean hazard from the 2023 NSHM update that has 2% probability of exceedance in 50 years, suggests an overall increase in total geographic areas facing very strong shaking levels (modified Mercalli Intensity (MMI) of VII or more) when compared with any of the previous cycles. This increase puts ~159.2 million people (an increase from 111.8 million when using the 2018 NSHM), ~47.3 million residential buildings (an increase from 32.5 million with the 2018 NSHM), and ~$25.6 trillion of economic exposure of all buildings (an increase from $17.4 trillion when using the 2018 NSHM) at risk from earthquakes. We demonstrate that small changes in hazard estimates do not necessarily imply small changes in risk estimates. Therefore, changes in risk estimates can be used to highlight key USGS NSHM updates for improving risk mitigation efforts.

Earthquake Spectra Journal↗

USGS2026

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

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

Fuel↗

USGS2026

Co-occurrence of pesticides and pharmaceuticals and personal care products (PPCPs) across Zostera marina (common eelgrass) communities

Anthropogenic pressures are driving changes in eelgrass communities, which are altering baseline conditions in estuarine environments. Field detections have validated the transport of land-sourced pollutants to aquatic systems; however, studies rarely sample concurrently for pesticides, and pharmaceuticals and personal care products (PPCPs) across environmental compartments. Moreover, studies on contaminant uptake by eelgrass and associated species are even more limited. In collaboration with the Confederated Tribes of the Coos, Lower Umpqua and Siuslaw Indians (CTCLUSI), this study collected samples of water, eelgrass, clams, and sediment at sites of Tribal significance in Southern Oregon to test for organic contaminants (i.e., herbicides and pharmaceuticals). Paired sampling was conducted for analysis by the CTCLUSI in tandem with the United States Geological Survey (USGS) in order for the Tribe to develop analytical standards for future sampling efforts. Ten pesticides and eight pharmaceuticals were detected across the four sites, with the highest number of overall detections (27) at the Florence Marina site. The insecticide bifenthrin was most frequently detected across all media (0.012–1.565 μg/g organic carbon in sediment, 2.7–30 ng/g in organismal tissue) and the anti-diabetic agent metformin was the most detected PPCP in clam tissues (1.33–3.78 ng/g). Pesticides and PPCPs were observed to co-occur in eelgrass habitats, with numerous pesticide detections across media types. These findings demonstrate numerous routes of exposure for estuarine organisms which could be addressed with pharmaceutical disposal strategies or pesticide use restrictions near these habitats.

Oregon↗

USGS2026

Localization of spatiotemporally heterogeneous subsurface flows using autoencoder-based deep learning framework for time-lapse self-potential tomography

Self-potential (SP) monitoring has emerged as a valuable method for characterizing subsurface hydrogeological features and processes due to its sensitivity to fluid-induced electrokinetic effects. Despite advancements in SP inversion, challenges remain in imaging groundwater dynamics from SP activities due to complex hydrological settings and transient noise. In this study, a deep learning autoencoder (AE)-based framework is proposed for the spatiotemporal localization of subsurface fluid movement from time-lapse SP tomography. Temporal segments of time-lapse numerical inversions were first derived from long-term SP monitoring conducted from a floodplain site in Oak Ridge, Tennessee, known for active hyporheic exchange. Subsequently, AE models based on vision transformer (ViT), convolutional long short-term memory (ConvLSTM), convolutional neural network, and temporal convolutional network were individually trained and compared on the SP tomography segments for reconstruction performance. Finally, the reconstruction error over time serves as an anomaly score to identify moments of active SP variation, whereas spatial distributions of errors within these moments are analyzed to image and localize regions associated with anomalous subsurface fluid movement. The results demonstrate that ConvLSTM- and ViT-AE are most capable for the localization task with contrasting error distributions and consistent delineation of anomalies. Applying the method to both SP arrays parallel and perpendicular to the stream produced consistent anomaly zones near a fault or karst feature, validating the robustness and generalization of the approach. These results demonstrate the potential of the proposed framework as a scalable and interpretable tool for spatiotemporal analysis of subsurface flow dynamics in complex hydrogeological systems.

Tennessee↗

USGS2026

Future water constraints on United States lithium mining under climate change

Lithium is necessary for low-carbon technologies that combat climate change, but lithium extraction is water-intensive. Changes in temperature and precipitation arising from climate change are altering water distribution, which could further strain supplies for new mines and industry, farms, and households. Here we explored how climate change, water use, and mining siting could impact lithium mining in the United States. We analyzed whether there would be sufficient water available to support the single existing and 22 proposed U.S. lithium mines at mid-century under four socioeconomic-climate scenarios and five climate models. Though dependent on socioeconomic-climate scenario, climate model, and lithium deposit type, available water supply in most subbasins would likely be unable to support new mines’ water demands, or even non-mining water demands from other sectors. Water scarcity could hinder the ability of the United States to produce enough lithium to meet domestic demand thereby necessitating higher imports.

conterminous United States↗

USGS2026

Growth patterns, sexual dimorphism, and seasonal body mass variation in Plains Zebra of the Serengeti region, Tanzania

Knowledge of growth patterns and physical properties of animals helps researchers assess population characteristics and health in the field. Information on age-specific body measurements and growth patterns is rare in ungulates because examining a large portion of a population is time-, budget-, and personnel-intensive in the wild. We collected age- and sex-specific body measurements from 442 Plains Zebra ( Equus quagga boehmi ) culled in the Serengeti region, Tanzania, between 1969 and 1971. The dataset contained 307 males and 133 females. We found that males were on average 2.8% larger than females in body measurements of heart girth, height at shoulder, neck-body length, hind foot length, and the body mass, indicating minimal sexual dimorphism in adult Plains Zebra. We used polynomial regression to explore the relationship between heart girth and body mass and applied third-degree polynomial regression models to predict body mass for 34 females and 115 males with recorded girth but missing body mass records. We applied 5 non-linear models—Gompertz, von Bertalanffy, Bertalanffy-Pütter, logistic, and Richards—to characterize sex-specific growth trajectories for each trait. Our results showed that different traits followed different growth patterns, with size asymptotes reached by ∼2 to 5 yr depending on trait and sex. Seasonal trends in adult body mass were also evaluated using generalized linear models incorporating sex and cyclical effects of day of year. While female body mass showed little seasonal variation, adult males exhibited a 16.9% increase in mass during the wet season, followed by a sharp decline, suggesting greater seasonal plasticity and energetic investment. Our work is one of the few studies on body size and growth in wild equids, which provides new information on age- and sex-specific growth trajectories for comparison across zebra subspecies over time and space.

Serengeti region↗

USGS2026

Effects of repeat prescribed burning in dry coniferous forests in national parks of California

Background Prescribed fire is a common approach to reduce fuels and mitigate fire hazards. The accumulation of live and dead fuels following initial treatment means that repeated application of prescribed fire could be used to maintain this benefit. However, the effect of repeated prescribed fires is not well documented in many dry coniferous forests in the western United States. Here, we present observations of changes in live trees and surface fuels following two prescribed fires in dry coniferous forests in national parks of California. Results Changes in forest structure and accumulation of surface fuels were similar over time following initial-entry and second-entry fires. An exception was that repeated fires were associated with substantial reductions in stem density. There were smaller changes in live tree basal area and stem biomass. Conclusions Our results indicate that following initial-entry fires, subsequent burning maintained reductions in surface fuel loads without major inadvertent losses of live tree basal area and stem biomass, implying the survival of large trees.

California↗

USGS2026

Range-wide relative abundance of the Appalachian grizzled skipper ( Pyrgus centaureae wyandot ) in the Eastern United States

The U.S. Fish and Wildlife Service has designated the Pyrgus centaureae wyandot (Appalachian Grizzled Skipper [AGS]) to be at-risk, based on its declining populations and the lack of information on its status. The objective of this study was to complete range-wide surveys to locate extant AGS colonies and to quantify the number of AGS observed at each location. From 2021–24, 284 surveys were done in 25 unique (that is, distinct) counties in 8 States in the Eastern United States — Maryland, Michigan, New York, North Carolina, Ohio, Pennsylvania, Virginia, and West Virginia. We found AGS in only two counties: Alleghany County, Virginia, and Greenbrier County, West Virginia. AGS were observed 180 times in these two counties. Our results can inform U.S. Fish and Wildlife decisions about where and how future AGS conservation efforts can be implemented.

Maryland, Michigan, New York, North Carolina, Ohio↗

USGS2026

Assessing environmental influences on gene flow for a migratory ungulate in a fire-prone landscape

Context Migratory species traverse long distances across complex landscape mosaics. Shifts in environmental attributes across space and time could affect movement paths and genetic connectivity of migratory species. Objectives We identified boundaries to genetic connectivity based on population structure and evaluated environmental influences on gene flow in a migratory ungulate. We projected how genetic connectivity might be affected by actual wildfire events. Methods We collected over 168,400 GPS locations for 97 bighorn sheep ( Ovis canadensis ) and estimated kinship with 5,140 single nucleotide polymorphisms for 95 individuals. We identified subpopulations and employed machine-learning optimization to evaluate how environmental attributes affect gene flow between them. We projected future genetic connectivity based on recent landscape changes. Results We observed that water bodies and > 70% canopy cover may limit effective dispersal, suggesting that shifts in these landscape attributes, influenced by climate and vegetation management, may affect future gene flow. We also predicted recent wildfire events will change gene flow between three subpopulations, predicting more flow where fire reduced canopy cover at a pinchpoint near a linear lake and less flow where fire was internal to the subpopulation. Conclusions Modifications in canopy cover could be applied to facilitate or impede connectivity and gene flow, depending on management goals. Large, semi-linear lakes could form movement barriers in other bighorn populations. Recent wildfires will likely increase gene flow between semi-isolated subpopulations, with implications for disease management. As wildfire and tree encroachment and densification increase, gene flow between other spatially structured populations in similar landscapes may shift.

British Columbia, Montana↗

USGS2026

Geochemical, mineralogical, and isotopic evidence for multi-stage genesis of the Hicks Dome REE + Y-HFSE-fluorite deposit, Illinois, USA

Hicks Dome hosts breccias enriched in rare earth elements (REE), Y, Th, F, Ba, Ti, Nb, and Be, alongside spatially associated lamprophyre dikes (ca. 271 Ma). Hicks Dome is located within the Illinois–Kentucky Fluorspar District, which hosts fluorite, Pb–Zn, and barite resources. This study investigates the genetic relationships between Hicks Dome mineralization in breccias, alkaline magmatism, and Illinois–Kentucky Fluorspar District mineralization. Lamprophyre dikes are light REE–enriched with chondrite-normalized abundances decreasing from La to Lu. The Host Breccia exhibits middle and heavy REE–enriched patterns that mirror those of the principal REE–Th host minerals, including fluorapatite, xenotime, and thorite. Textural evidence suggests recrystallization of phosphates, sulfates, and Ti–Nb oxides in the Host Breccia. U–Pb geochronology constrains multiple mineralizing events, with ages of 277 ± 18 Ma from low-Th apatite interpreted as main-stage mineralization, and 121.6 ± 9.7 Ma from high-Th apatite indicating later overprinting. O–H–C stable isotope data provide evidence for multiple stages of fluid-rock interaction and fluid mixing: (1) early magmatic fluids dissolved limestone country rock, (2) mixing between magmatic fluids and basinal brines led to main-stage mineralization in the Host Breccia, and (3) late-stage mineralization occurred following mixing of meteoric water and basinal brine. These results indicate that heavy REEs, high field strength elements, and fluorine precipitated proximal to its alkaline magmatic source because of fluid–rock interactions and fluid mixing. Subsequent fluid mixing drove late-stage recrystallization and additional fluorite formation, a process that may be similar to mineralization in the Illinois-Kentucky Fluorspar District.

Illinois↗

USGS2026

Moment magnitude for small earthquakes in the Delaware basin of west Texas and southeast New Mexico, USA

The Delaware Basin region of west Texas and southeast New Mexico has become one of the most prolific regions of seismic activity in the continental United States due to widespread hydraulic fracturing and wastewater disposal injection. In response to the increased number of earthquakes in this region, rapid and accurate characterization of earthquake sources is necessary to understand the evolution of seismic activity and level of seismic hazard associated with these earthquakes. This study re-evaluates earthquake magnitudes, estimating moment magnitude (MW) for small earthquakes in the Delaware Basin using 1) moment-rate spectra derived from S-wave coda envelopes, and 2) a relative magnitude method that relies exclusively on the ratio of waveform amplitudes between highly correlated waveform pairs. The coda-envelope method produces accurate M W estimates for small earthquakes ( M 1.5 – 3) that are consistent with independent, waveform modeled moment magnitudes for events with M W > 3 . Using the relative amplitudes method to extend these M W magnitudes to many other events, we successfully provide relative moment magnitude ( M W,rel ) values for 81% of the Texas Seismological Network catalog in the Delaware Basin region, and 45% of the USGS Induced Seismicity Project’s catalog of events in southeast New Mexico. The adoption and integration of the calibrated M W,rel method with current magnitude estimation methods offers valuable insights into the relationships between local and moment magnitude and will contribute to improved characterization of widespread induced seismicity.

New Mexico, Texas↗

USGS2026

Demographic responses to experimental manipulation of an enclosed population of invasive brown treesnakes (Boiga irregularis)

Understanding how invasive species populations respond to perturbations can give insight into effective management. However, opportunities to test theoretical population responses are rare, as experimental perturbations often run counter to management goals (e.g., suppression). On the island of Guam, a 5-ha outdoor enclosure (with no immigration or emigration) constructed in 2003 serves as a unique population-scale study area for invasive brown treesnakes ( Boiga irregularis ). Capture-mark-recapture studies (CMR) and several experiments have taken place since 2004. Here, we analyze a subset of the CMR data collected between 2008 and 2012, a time period encompassing two experimental perturbations: (1) supplemental feeding to determine if the population is food limited and (2) targeted removal of adult and large-subadult cohorts to quantify the timing and magnitude of the population's demographic response to a simulated control effort. We used the CMR data to fit a Pradel model in a Bayesian framework to compare the population growth rate, survival probability, and recruitment rate prior to, during, and post perturbations. We found that increasing food resources increased population growth and recruitment rates, but that survival remained relatively constant. The removal perturbation resulted in a decrease in population growth rate and survival immediately after the removal, but population growth rate and recruitment began increasing 1 year after the removal perturbation. Our results suggest that brown treesnakes in the enclosed population are prey limited, likely constraining population growth, and that the suppression of population growth after removal of adults and large sub-adults is short-lived (i.e., less than a year).

Guam↗

USGS2026

Streamflow and surface-water presence data availability across the conterminous United States: A review for headwater systems

Water is essential for life on Earth, supporting ecosystems, human health, and economic activities. Hydrology relies on observational data, and this paper discusses regional and national datasets for the conterminous United States (CONUS) publicly available as of 2023, focusing on headwaters, defined as first- and second-order streams at 1:24000 scale. It identifies 72 primary and secondary datasets and 11 repositories and argues how better integration and accessibility of hydrological data can improve research. The paper distinguishes between datasets where streamflow was the primary data collection objective and those where it was secondary. This distinction highlights opportunities to consider data from efforts peripheral to hydrology but is still useful for understanding hydrologic conditions. The analysis reveals that out of about 118 000 active and inactive stream observation sites, about 6.6% and 25% are located on first- and second-order streams, respectively. This indicates a substantial data gap for headwater systems, which account for over 77% of stream length in CONUS. Federal agencies manage 72% of hydrologic monitoring sites across all stream orders, but only 34% of these are in headwater systems. Academic institutions operate about 2% of sites, with almost half (48%) in headwater systems, focusing on ecosystem research. State agencies also operate about 2% of sites, primarily on larger systems, with 19% on headwaters. Additionally, 23% of sites are managed by multiple agencies. Spatial patterns further reveal pronounced disparities among physiographic regions. Eastern and coastal provinces show relatively dense monitoring, while central and western regions show sparse coverage. These gaps reflect historical priorities, logistical constraints, funding limitations, and the high cost of continuous instrumentation. To address biases in monitoring networks, data collection could be enhanced with low-cost monitoring, community science, and remote sensing technologies. This study also notes the benefits of long-term monitoring and prioritizing retention of streamgages with longer records.

conterminous United States↗