SEARCH · USGS Science

Results for “Developments in Water Science”

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.

At least 1,225 records · Page 68Linked to original sources

USGS2022

Long-term impacts of impervious surface cover change and roadway deicing agent application on chloride concentrations in exurban and suburban watersheds

Roadway deicing agents, including rock salt and brine containing NaCl, have had a profound impact on the water quality and aquatic health of rivers and streams in urbanized areas with temperate climates. Yet, few studies evaluate impacts to watersheds characterized by relatively low impervious surface cover (ISC; < 15 %). Here, we use long-term (1997-2019), monthly streamwater quality data combined with daily streamflow for six exurban and suburban watersheds in southeastern Pennsylvania to examine the relations among chloride (Cl − ) concentrations and ISC. Both flow-normalized Cl − concentrations and ISC increased over time in each of the six watersheds, consistent with changes in watershed management (e.g., ISC, road salt application, etc.). The watersheds that experienced the greatest changes in percent ISC (e.g., agriculture replaced by residential and commercial development) experienced the greatest changes in flow-normalized Cl − concentrations. We also utilized a comprehensive mass-balance model (2011–2018) that indicated Cl − inputs exceeded the outputs for the study watersheds. Road salt applied to state roads, non-state roads, and other impervious surfaces accounted for the majority of Cl − inputs to the six watersheds. Furthermore, increasing Cl − concentrations during baseflow conditions confirm impacts to shallow groundwater. Although flow-normalized Cl − concentrations are below the U.S. Environmental Protection Agency's chronic threshold value for impacts to aquatic organisms, year-round exceedances may result before the end of this century based on current trends. Though reduced Cl − loading to streams may be achieved by limiting the expansion of impervious surfaces in exurban and suburban watersheds, changes in baseflow concentrations are likely to be gradual because of the accumulated Cl − in groundwater.

Pennsylvania↗

USGS2012

Variance partitioning of stream diatom, fish, and invertebrate indicators of biological condition

Stream indicators used to make assessments of biological condition are influenced by many possible sources of variability. To examine this issue, we used multiple-year and multiple-reach diatom, fish, and invertebrate data collected from 20 least-disturbed and 46 developed stream segments between 1993 and 2004 as part of the US Geological Survey National Water Quality Assessment Program. We used a variance-component model to summarize the relative and absolute magnitude of 4 variance components (among-site, among-year, site × year interaction, and residual) in indicator values (observed/expected ratio [O/E] and regional multimetric indices [MMI]) among assemblages and between basin types (least-disturbed and developed). We used multiple-reach samples to evaluate discordance in site assessments of biological condition caused by sampling variability. Overall, patterns in variance partitioning were similar among assemblages and basin types with one exception. Among-site variance dominated the relative contribution to the total variance (64–80% of total variance), residual variance (sampling variance) accounted for more variability (8–26%) than interaction variance (5–12%), and among-year variance was always negligible (0–0.2%). The exception to this general pattern was for invertebrates at least-disturbed sites where variability in O/E indicators was partitioned between among-site and residual (sampling) variance (among-site = 36%, residual = 64%). This pattern was not observed for fish and diatom indicators (O/E and regional MMI). We suspect that unexplained sampling variability is what largely remained after the invertebrate indicators (O/E predictive models) had accounted for environmental differences among least-disturbed sites. The influence of sampling variability on discordance of within-site assessments was assemblage or basin-type specific. Discordance among assessments was nearly 2× greater in developed basins (29–31%) than in least-disturbed sites (15–16%) for invertebrates and diatoms, whereas discordance among assessments based on fish did not differ between basin types (least-disturbed = 16%, developed = 17%). Assessments made using invertebrate and diatom indicators from a single reach disagreed with other samples collected within the same stream segment nearly ⅓ of the time in developed basins, compared to ⅙ for all other cases.

Freshwater Science↗

USGS2007

High nutrient pulses, tidal mixing and biological response in a small California estuary: Variability in nutrient concentrations from decadal to hourly time scales

Elkhorn Slough is a small estuary in Central California, where nutrient inputs are dominated by runoff from agricultural row crops, a golf course, and residential development. We examined the variability in nutrient concentrations from decadal to hourly time scales in Elkhorn Slough to compare forcing by physical and biological factors. Hourly data were collected using in situ nitrate analyzers and water quality data sondes, and two decades of monthly monitoring data were analyzed. Nutrient concentrations increased from the mid 1970s to 1990s as pastures and woodlands were converted to row crops and population increased in the watershed. Climatic variability was also a significant factor controlling interannual nutrient variability, with higher nutrient concentrations during wet than drought years. Elkhorn Slough has a Mediterranean climate with dry and rainy seasons. Dissolved inorganic nitrogen (DIN) concentrations were relatively low (10-70 ??mol L-1) during the dry season and high (20-160 ??mol L-1) during the rainy season. Dissolved inorganic phosphorus (DIP) concentrations showed the inverse pattern, with higher concentrations during the dry season. Pulsed runoff events were a consistent feature controlling nitrate concentrations during the rainy season. Peak nitrate concentrations lagged runoff events by 1 to 6 days. Tidal exchange with Monterey Bay was also an important process controlling nutrient concentrations, particularly near the mouth of the Slough. Biological processes had the greatest effect on nitrate concentrations during the dry season and were less important during the rainy season. While primary production was enhanced by nutrient pulses, chlorophyll a concentrations were not. We believe that the generally weak biological response compared to the strong physical forcing in Elkhorn Slough occurred because the short residence time and tidal mixing rapidly diluted nutrient pulses. ?? 2006 Elsevier Ltd. All rights reserved.

Estuarine, Coastal and Shelf Science↗

USGS2016

Climate change and the Delta

Anthropogenic climate change amounts to a rapidly approaching, “new” stressor in the Sacramento–San Joaquin Delta system. In response to California’s extreme natural hydroclimatic variability, complex water-management systems have been developed, even as the Delta’s natural ecosystems have been largely devastated. Climate change is projected to challenge these management and ecological systems in different ways that are characterized by different levels of uncertainty. For example, there is high certainty that climate will warm by about 2°C more (than late-20th-century averages) by mid-century and about 4°C by end of century, if greenhouse-gas emissions continue their current rates of acceleration. Future precipitation changes are much less certain, with as many climate models projecting wetter conditions as drier. However, the same projections agree that precipitation will be more intense when storms do arrive, even as more dry days will separate storms. Warmer temperatures will likely enhance evaporative demands and raise water temperatures. Consequently, climate change is projected to yield both more extreme flood risks and greater drought risks. Sea level rise (SLR) during the 20th century was about 22cm, and is projected to increase by at least 3-fold this century. SLR together with land subsidence threatens the Delta with greater vulnerabilities to inundation and salinity intrusion. Effects on the Delta ecosystem that are traceable to warming include SLR, reduced snowpack, earlier snowmelt and larger storm-driven streamflows, warmer and longer summers, warmer summer water temperatures, and water-quality changes. These changes and their uncertainties will challenge the operations of water projects and uses throughout the Delta’s watershed and delivery areas. Although the effects of climate change on Delta ecosystems may be profound, the end results are difficult to predict, except that native species will fare worse than invaders. Successful preparation for the coming changes will require greater integration of monitoring, modeling, and decision making across time, variables, and space than has been historically normal.

San Francisco Estuary and Watershed Science↗

USGS2024

Rapid simulation of wave runup on morphologically diverse, reef-lined coasts with the BEWARE-2 (Broad-range Estimator of Wave Attack in Reef Environments) meta-process model

Low-lying, tropical, coral-reef-lined coastlines are becoming increasingly vulnerable to wave-driven flooding due to population growth, coral reef degradation, and sea-level rise. Early-warning systems (EWSs) are needed to enable coastal authorities to issue timely alerts and coordinate preparedness and evacuation measures for their coastal communities. At longer timescales, risk management and adaptation planning require robust assessments of future flooding hazard considering uncertainties. However, due to diversity in reef morphologies and complex reef hydrodynamics compared to sandy shorelines, there have been no robust analytical solutions for wave runup to allow for the development of large-scale coastal wave-driven flooding EWSs and risk assessment frameworks for reef-lined coasts. To address the need for fast, robust predictions of runup that account for the natural variability in coral reef morphologies, we constructed the BEWARE-2 (Broad-range Estimator of Wave Attack in Reef Environments) meta-process modeling system. We developed this meta-process model using a training dataset of hydrodynamics and wave runup computed by the XBeach Non-Hydrostatic process-based hydrodynamic model for 440 combinations of water level, wave height, and wave period with 195 representative reef profiles that encompass the natural diversity in real-world fringing coral reef systems. Through this innovation, BEWARE-2 can be applied in a larger range of coastal settings than meta-models that rely on a parametric description of the coral reef geometry. In the validation stage, the BEWARE-2 modeling system produced runup results that had a relative root mean square error of 13% and relative bias of 5% relative to runup simulated by XBeach Non-Hydrostatic for a large range of oceanographic forcing conditions and for diverse reef morphologies (root mean square error and bias 0.63 and 0.26 m, respectively, relative to mean simulated wave runup of 4.85 m). Incorporating parametric modifications in the modeling system to account for variations in reef roughness and beach slope allows for systematic errors (relative bias) in BEWARE-2 predictions to be reduced by a factor of 1.5–6.5 for relatively coarse or smooth reefs and mild or steep beach slopes. This prediction provided by the BEWARE-2 modeling system is faster by 4–5 orders of magnitude than the full, process-based hydrodynamic model and could therefore be integrated into large-scale EWSs for tropical, reef-lined coasts and used for large-scale flood risk assessments.

Natural Hazards and Earth System Sciences↗

USGS2013

Volcano crisis response at Yellowstone volcanic complex - after-action report for exercise held at Salt Lake City, Utah, November 15, 2011

A functional tabletop exercise was run on November 14-15, 2011 in Salt Lake City, Utah, to test crisis response capabilities, communication protocols, and decision-making by the staff of the multi-agency Yellowstone Volcano Observatory (YVO) as they reacted to a hypothetical exercise scenario of accelerating volcanic unrest at the Yellowstone caldera. The exercise simulated a rapid build-up of seismic activity, ground deformation, and hot-spring water-chemistry and temperature anomalies that culminated in a small- to moderate-size phreatomagmatic eruption within Yellowstone National Park. The YVO scientific team's responses to the unfolding events in the scenario and to simulated requests for information by stakeholders and the media were assessed by (a) the exercise organizers; (b) several non-YVO scientists, who observed and queried participants, and took notes throughout the exercise; and (c) the participants themselves, who kept logs of their actions during the exercise and later participated in a group debriefing session and filled out detailed questionnaires. These evaluations were tabulated, interpreted, and summarized for this report, and on the basis of this information, recommendations have been made. Overall, the YVO teams performed their jobs very well. The exercise revealed that YVO scientists were able to successfully provide critical hazards information, issue information statements, and appropriately raise alert levels during a fast-moving crisis. Based on the exercise, it is recommended that several measures be taken to increase YVO effectiveness during a crisis: 1. Improve role clarification within and between YVO science teams. 2. Improve communications tools and protocols for data-sharing and consensus-building among YVO scientists, who are geographically and administratively dispersed among various institutions across the United States. 3. Familiarize YVO staff with Incident Command System (ICS) procedures and protocols, and provide more in-depth training to appropriate staff members, as needed. 4. Train all science team members in the use of all analytical and computational tools available to them, in order to maximize effectiveness of teams in tracking and interpreting possible accelerating unrest at Yellowstone. Desirable pre-crisis preparations include: (a) updating a catalog of existing map and information products (and identifying additional products) that would be helpful during a crisis; (b) creating "to do" lists of early-crisis tasks for each scientific team; (c) coordinating radio frequencies among partner agencies; and (d) brief training on and promotion of the internal YVO Web log as a repository for scientific observations, data, photographs, and other material to be shared among YVO scientific teams during a crisis. This exercise was designed as an opportunity to practice response to a fast-developing volcano crisis and to test for organizational and procedural weaknesses that could emerge during a real crisis. This report is based upon the observations of the exercise organizers during the one-day exercise and upon written evaluations by the participants. It does not attempt to evaluate any other aspect of YVO or the scientific expertise of any of the highly competent YVO staff. Participants unanimously found the exercise to be helpful for improving their response capabilities, and it is our hope that the report will be a starting point for internal discussions that will make YVO even better-prepared for some future volcano crisis.

Utah↗

USGS2003

Escherichia coli at Ohio bathing beaches—Distribution, sources, wastewater indicators, and predictive modeling

Results of studies during the recreational seasons of 2000 and 2001 strengthen the science that supports monitoring of our Nation's beaches. Water and sediment samples were collected and analyzed for concentrations of Escherichia coli ( E. coli ). Ancillary water-quality and environmental data were collected or compiled to determine their relation to E. coli concentrations. Data were collected at three Lake Erie urban beaches (Edgewater, Villa Angela, and Huntington), two Lake Erie beaches in a less populated area (Mentor Headlands and Fairport Harbor), and one inland-lake beach (Mosquito Lake). The distribution of E. coli in water and sediments within the bathing area, outside the bathing area, and near the swash zone was investigated at the three Lake Erie urban beaches and at Mosquito Lake. (The swash zone is the zone that is alternately covered and exposed by waves.) Lake-bottom sediments from outside the bathing area were not significant deposition areas for E. coli . In contrast, interstitial water and subsurface sediments from near the swash zone were enriched with E. coli . For example, E. coli concentrations were as high as 100,000 colonies per 100 milliliters in some interstitial waters. Although there are no standards for E. coli in swash-zone materials, the high concentrations found at some locations warrant concern for public health. Studies were done at Mosquito Lake to identify sources of fecal contamination to the lake and bathing beach. Escherichia coli concentrations decreased with distance from a suspected source of fecal contamination that is north of the beach but increased at the bathing beach. This evidence indicated that elevated E. coli concentrations at the bathing beach are of local origin rather than from transport of bacteria from sites to the north. Samples collected from the three Lake Erie urban beaches and Mosquito Lake were analyzed to determine whether wastewater indicators could be used as surrogates for E. coli at bathing beaches. None of the concentrations of wastewater indicators of fecal contamination, including 3b-coprostanol and cholesterol, were significantly correlated (a=0.05) to concentrations of E. coli . Concentrations of the two compounds that were significantly correlated to E. coli were components of coal tar and asphalt, which are not necessarily indicative of fecal contamination. Data were collected to build on an earlier 1997 study to develop and test multiple-linear-regression models to predict E. coli concentrations using water-quality and environmental variables as explanatory variables. The probability of exceeding the single-sample bathing-water standard for E. coli (235 colonies per 100 milliliters) was used as the model output variable. Threshold probabilities for each model were established. Computed probabilities that are less than a threshold probability indicate that bacterial water quality is most likely acceptable. Computed probabilities equal to or above the threshold probability indicate that the water quality is most likely not acceptable and that a water-quality advisory may be needed. Models were developed at each beach, whenever possible, using combinations of 1997, 2000, and (or) 2001 data. The models developed and tested in this study were shown to be beach specific; that is, different explanatory variables were used to predict the probability of exceeding the standard at each beach. At Mentor Headlands and Fairport Harbor, models were not developed because water quality was generally good. At the three Lake Erie urban beaches, models were developed with variable lists that included the number of birds on the beach at the time of sampling, lake-current direction, wave height, turbidity, streamflow of a nearby river, and rainfall. The models for Huntington explained a larger percentage of the variability in E. coli concentrations than the models for Edgewater and Villa Angela. At Mosquito Lake, a model based on 2000 and 2001 data contained the explanatory variables rainfall, number of dry days preceding a rainfall, date, wind direction, wind speed, and turbidity. Additional research could include testing the threshold probabilities assigned for these models in subsequent years and comparing the models’ ability to predict recreational water quality to results from the current method—using antecedent E. coli concentrations. Each year the model is tested, new data can be added and model variables can be recalculated to determine whether the predictive ability improves with additional data.

Ohio↗

USGS2025

Mapping global coral vulnerability to stony coral tissue loss disease: Implications for biosecurity and conservation

Stony Coral Tissue Loss Disease (SCTLD) has devastated Caribbean coral reefs since 2014, but its potential for global impact remains uncertain. We developed predictive models to assess the worldwide vulnerability of coral reefs to SCTLD under different origin and spread hypotheses. Using random forest regression models incorporating coral taxonomy and zooxanthellae clade associations from 52 taxa, we projected SCTLD susceptibility and mortality patterns globally using six indices: Mean susceptibility per genus per location, Summed susceptibilities across genera per location, Summed susceptibilities across genera per realm, Mean mortality per genus per location, Summed mortalities across genera per location, and Summed mortalities across genera per realm. Models demonstrated strong predictive performance (R² = 0.57 for susceptibility; R² = 0.73 for mortality) and revealed that about 7% of coral genera per location are potentially susceptible to SCTLD. While mean susceptibility and mortality per genus were highest in the Tropical Atlantic, the summed susceptibility and mortality across genera were much higher in the biodiverse Central Indo-Pacific. Natural barriers could limit SCTLD’s spread, including the mid-Atlantic gap and the low diversity of the Tropical Eastern Pacific, supporting the contained disease hypothesis. However, the widespread distribution of susceptible genera across coral reef realms indicates significant vulnerability should SCTLD circumvent these barriers through human-mediated transport, particularly via ballast water or the aquarium trade. If SCTLD is an invasive pathogen originating in the Pacific, as shipping patterns for the aquarium trade suggest, mortality in its native range would likely be lower than our projections. These findings point to targeted intervention strategies, including enhanced monitoring at key locations, assessment of biosecurity needs in high-risk areas, and prioritized conservation efforts in vulnerable high-diversity regions to prevent SCTLD from spreading globally.

Frontiers in Marine Science↗

USGS2003

Physical stream habitat dynamics in Lower Bear Creek, northern Arkansas

We evaluated the roles of geomorphic and hydrologic dynamics in determining physical stream habitat in Bear Creek, a stream with a 239 km2 drainage basin in the Ozark Plateaus (Ozarks) in northern Arkansas. During a relatively wet 12-month monitoring period, the geomorphology of Bear Creek was altered by a series of floods, including at least four floods with peak discharges exceeding a 1-year recurrence interval and another flood with an estimated 2- to 4-year recurrence interval. These floods resulted in a net erosion of sediment from the study reach at Crane Bottom at rates far in excess of other sites previously studied in the Ozarks. The riffle-pool framework of the study reach at Crane Bottom was not substantially altered by these floods, but volumes of habitat in riffles and pools changed. The 2- to 4-year flood scoured gravel from pools and deposited it in riffles, increasing the diversity of available stream habitat. In contract, the smaller floods eroded gravel from the riffles and deposited it in pools, possibly flushing fine sediment from the substrate but also decreasing habitat diversity. Channel geometry measured at the beginning of the study was use to develop a two-dimensional, finite-element hydraulic model at assess how habitat varies with hydrologic dynamics. Distributions of depth and velocity simulated over the range of discharges observed during the study (0.1 to 556 cubic meters per second, cms) were classified into habitat units based on limiting depths and Froude number criteria. The results indicate that the areas of habitats are especially sensitive to change to low to medium flows. Races (areas of swift, relatively deep water downstream from riffles) disappear completely at the lowest flows, and riffles (areas of swift, relatively shallow water) contract substantially in area. Pools also contract in area during low flow, but deep scours associated with bedrock outcrops sustain some pool area even at the lowest modeled flows. Modeled boundary shear stresses were used to evaluate which flows are responsible for the most mobilization of the bed, and therefore, habitat maintenance. Evaluation of the magnitude and frequency of bed-sediment entrainment shows that most of the habitat maintenance results from flows that occur on average about 4 to 7 days a year. Our analysis documents the geomorphic and hydrologic dynamics that form and maintain habitats in a warmwater stream in the Ozarks. The range of flows that occurs on this stream can be partitioned into those that sustain habitat by providing the combinations of depth and velocity that stream organisms live with most of the time, and those flows that surpass sediment entrainment thresholds, alter stream geomorphology, and therefore maintain habitat. The quantitative relations show sensitivity of habitats to flow variation, but do not address how flow may vary in the future, or the extent to which stream geomorphology may be affected by variations in sediment supply.

Arkansas↗

USGS2024

Seasonally varying contributions of contemporaneous and lagged sources of instream total nitrogen and phosphorus load across the Illinois River basin

Quantifying nutrient sources in streams, their temporal and spatial variability, and drivers of that variability can support effective water resources management. Yet a lack of data and modeling capabilities has previously prevented comprehensive quantification across both space and time. Here a dynamic SPARROW (Spatially Referenced Regressions on Watershed attributes) model that accounts for a lagged delivery of nutrients to streams was developed and applied to simulate seasonal and source-specific total nitrogen (TN) and total phosphorus (TP) loads in streams across the Illinois River basin (IRB). Dynamic load predictions from 2000 through 2020 revealed that a third of the TN and a quarter of the TP instream load originated from non-point sources that were lagged in their delivery from land-application to streams by more than a season. This lagged mass was the largest overall TN source—which was estimated as a lagged expression of previous seasonal non-point sources including fertilizer, manure, atmospheric deposition and fixation, and urban land use. Treated wastewater effluent was the largest TP source exported from the basin, contributing 39 % of the TP load and 15 % of the TN load, and dominated the load in the upper Illinois River near Chicago. Loads in the lower river during this period, conversely, were attributed primarily to a mix of agricultural sources and their lagged fractions from headwater tributaries. Instream processes removed 10 % of the TN load while only 4 % of the TP load was removed during instream transport. With appropriate datasets, the models could be extended to other basins or time periods and used to forecast future seasonal nutrient loads.

Illinois↗

USGS2024

Streams, springs, and volcanic lakes for volcano monitoring

Introduction Volcanic unrest can trigger appreciable change to surface waters such as streams, springs, and volcanic lakes. Magma degassing produces gases and soluble salts that are absorbed into groundwater that feeds streams and lakes. As magma ascends, the amount of heat and degassing will increase, and so will any related geochemical and thermal signal. Subsurface magma movement can cause pressurization that alters hydrostatic head and may induce groundwater discharge. Fluid-pressure changes have been linked to distal volcano-tectonic earthquakes (White and McCausland, 2016; Coulon and others, 2017) and phreatic eruptions (for example, Yamaoka and others, 2016). Clearly, changes in groundwater and surface waters are both indicators of unrest and clues to how and where magma is rising toward the surface. Where possible, it is prudent to incorporate real-time hydrologic data into multiparameter monitoring of restless volcanoes. Hydrologic dynamics can also be tracked by changes in groundwater levels that are commonly measured in shallow boreholes (see of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume, on boreholes; Hurwitz and Lowenstern, 2024). Although inferred to be common, relatively few volcano-hydrology anomalies are well documented, and many are essentially anecdotal (Newhall and others, 2001), reflecting the fact that high-resolution time series remain rare. Extreme examples include the 2008 eruption of Nevado del Huila, Colombia, where relatively minor phreatomagmatic eruptions were accompanied by expulsion of as much as 300 million cubic meters of groundwater from fissures high on the volcano (Worni and others, 2011), generating large lahars. Substantial decreases in flow rate from springs about 8 kilometers from the summit of Mayon Volcano, Philippines, have been noted before most eruptions in the 20th century (Newhall and others, 2001). Stream monitoring at Redoubt Volcano in 2009 allowed Werner and others (2012) to recognize that groundwater was unable to absorb (or scrub) the high flux of volcanic gas and that a high CO 2 /SO 2 precursor signal had been evident for 5 months prior to the eruption. A key to better interpreting hydrologic anomalies—or even identifying them—is therefore obtaining adequate baseline data. Most hydrologic monitoring at U.S. volcanoes has been accomplished by intermittent sampling surveys with annual or less frequent sampling (for example, https://hotspringchem.wr.usgs.gov/index.php ). More frequent sampling, however, generally is needed to establish reliable baselines. A recent hydrologic and hydrothermal monitoring experiment at 25 sites and 10 of the 12 level 4 (very high threat) volcanoes in the U.S. portion of the Cascade Range demonstrated that there is sufficient temporal variability in hydrothermal fluxes, even during quiescent periods, that one-time measurements will commonly have limited interpretive value (Crankshaw and others, 2018). Thus, surveys are best augmented with data from streamgages (for example, Evans and others, 2004; Bergfeld and others, 2008). Streamflow (water discharge) data allow measured temperature and specific conductance to be converted to heat and solute mass fluxes, which could be insightful parameters for detecting anomalous activity (McCleskey and others, 2012). At the Yellowstone Caldera, long-term monitoring of river solutes has allowed calculation of the chloride flux, a proxy for heat discharge (Hurwitz and others, 2007; McCleskey and others, 2016) from the subsurface magma. This is readily accomplished because data from streamgages are continuously recorded and archived by the U.S. Geological Survey (USGS) National Water Information System (NWIS) (USGS, 2024). Similar studies on stratovolcanoes or shield volcanoes would be scientifically useful, and yet are logistically challenging, requiring streamgages on numerous radial drainages complemented by either frequent manual sampling or numerous deployments of equipment to measure water temperature and specific conductance as a proxy for water chemistry. Another challenge is that some volcanic areas, especially shield volcanoes, are characterized by near-surface porous rocks and soils, such that surface streams are rare and replaced by distant, dilute large-volume springs with only a trace of any original volcanically sourced water (Manga, 2001; Hurwitz and others, 2021). Volcanic lakes are worthy of special attention for monitoring efforts, as their temperature and composition can provide evidence of increased flux of volatile-rich fluids from below. Quantifying changes in volatile and heat release from magma can be simpler in lakes than for volcanoes with radial drainages and no major lakes. Moreover, volcanic lakes pose a range of hazards themselves, including phreatomagmatic eruptions, debris flows, flank collapse, tsunamis, and toxic gas release (Mastin and Witter, 2000; Delmelle and others, 2015; Manville, 2015; Rouwet and others, 2015)—hazards that have historically been responsible for substantial loss of life at many volcanoes worldwide (Manville, 2015). Catastrophic CO 2 release at Lake Nyos, Cameroon, in 1986 suffocated about 1,750 people and about 3,500 livestock and was probably triggered by a large landslide into the gas-saturated lake (Kling and others, 1987; Evans and others, 1993). Gas-charged springs in Soda Bay within Clear Lake (California) have caused almost a dozen deaths to bathers in the past hundred years (ABC News, 2000). A 2005 example of lake overturn and abundant gas release was documented at Mount Chiginagak in Alaska (Schaefer and others, 2008) but did not result in any human casualties. Although thermally stratified lakes, which promote trapping of exsolved magmatic gas, tend to develop in tropical regions, the phenomenon can also arise where salinity creates meromixis (a condition in which a lake does not mix completely), as occurs in Mono Lake, California (Jellison and Melack, 1993; Jellison and others, 1998). If magma erupts or flows into a lake, the interaction between hot magma and cold water can be explosive (Mastin and others, 2004; Zimanowski and others, 2015) and substantially expand the area affected by the eruption. Another hazard is the breaching of crater rims by landslides triggered by volcanic and (or) seismic activity. Under some circumstances, substantial volumes of water can be displaced, leading to large floods and lahars. Late Holocene lake flooding from Aniakchak Crater in the Alaska Peninsula (Waythomas, 2022) and from Paulina Lake in Newberry Crater, Oregon (Chitwood and Jensen, 2000), caused by the failure of outlet sills, testify to the substantial hazards at lake-filled calderas. Several volcanic systems in the United States host lakes known to receive heat and gas from underlying magma. These lakes vary widely in area, depth, and chemical composition. Lakes are present at level 4 volcanoes, including Crater Lake and Newberry Volcano in Oregon; Yellowstone Caldera in Wyoming; Long Valley Caldera, Clear Lake volcanic field, Medicine Lake, and Salton Buttes in California; and Aniakchak Crater, Mount Katmai, Fisher Caldera, Mount Okmok, and Kaguyak Crater, among others, in Alaska. A water lake was present in Halemaʻumaʻu, the crater of Kīlauea, Hawai‘i (fig. F1), from October 2019 to December 2020. Level 3 volcanoes with lakes include Mono Lake volcanic field (Calif.), Mount Bachelor (Ore.), Ukinrek Maars and Mount Chiginagak (Alaska), and Soda Lake (Nevada). In addition, there are lakes at many levels 1 and 2 volcanoes. In the United States, there are no strongly acidic lakes that receive abundant input of magmatic gas, such as those found at Mount Ruapehu (New Zealand), Ijen and Kelud (Indonesia), and Poás (Costa Rica). Nevertheless, many contain fluids that provide clues to magmatic processes below. Since publication of a previous report on recommended instrumentation for volcano monitoring (Moran and others, 2008), continuous hydrologic monitoring has become increasingly feasible. However, changes in water pressure, temperature, and chemistry remain, in general, poorly studied phenomena at volcanoes (Sparks, 2003; National Academies of Sciences, Engineering, and Medicine, 2017). Recent efforts by the USGS have included the temporary study of Cascade Range volcanoes, which included frequent (15 minute to hourly) temporal sampling of temperature, depth, and conductivity (Crankshaw and others, 2018; Ingebritsen and Evans, 2019). At Yellowstone Caldera, many streamgages have now added thermistors and specific conductance sensors, allowing estimation of time-dependent chloride flux as a proxy for variations in subsurface heat flux (McCleskey and others, 2012, 2016). Efforts to better understand lakes have also accelerated, with bathymetric mapping and sampling carried out at several locations in the United States. Especially thorough work was done at Yellowstone Lake thanks to the Hydrothermal Dynamics of Yellowstone Lake (HD-YLAKE, https://hdylake.org ) project, funded primarily by the National Science Foundation. In addition to geophysical surveys and recovery of cores and other samples, HD-YLAKE investigations included remotely operated vehicle (ROV) investigations of hydrothermal vents on the lake floor (fig. F2). Data collected by the ROV provided a better understanding of the thermal and chemical influx from lake-bottom hydrothermal systems (Sohn and others, 2017). In this chapter, we focus on detecting changes in the chemistry, temperature, discharge, or water levels of streams, springs, and lakes that can be caused by seismicity, volumetric strains, or increases in gas flux associated with ascending magma. There is unavoidable overlap with other chapters of this report. Samples of water and gas can also be obtained in boreholes ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume; Hurwitz and Lowenstern, 2024), both shallow and deep. Gas monitoring ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–E, 11 p., https://doi.org/10.3133/sir20245062e.">chapter E , this volume; Lewicki and others, 2024) relies in part on samples from springs and wells, particularly where measurable gas plumes are absent. Water acts as a trigger and lubricant for landslides and sediment-rich floods, and so hydrology has obvious relevance for lahar monitoring, as discussed in of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–H, 6 p., https://doi.org/10.3133/sir20245062h. ">chapter H (this volume; Thelen and others, 2024). Shared situational awareness among scientists engaged in geophysical, gas, and hydrologic monitoring will improve overall understanding of the volcanic hazard.

Scientific Investigations Report↗

USGS2026

The surface elevation table and marker horizon technique: A protocol for measuring wetland elevation dynamics, narrative (Version 2.0)

The National Park Service (NPS), in response to the growing evidence and awareness of the effects of climate change on federal lands, determined that monitoring wetland elevation change is a top priority in North Atlantic Coastal parks (Stevens et al. 2010). As a result, the NPS Northeast Coastal and Barrier Network (NCBN) in collaboration with colleagues from the U.S. Geological Survey (USGS) and the National Oceanic and Atmospheric Administration (NOAA) have developed a protocol for monitoring wetland elevation change and other processes important for determining the viability of wetland communities. Although focused on North Atlantic Coastal parks, this document is applicable to all coastal and inland wetland regions. Wetlands exist within a narrow range of elevation that is influenced by local hydrologic conditions. For tidally influenced coastal wetlands, local hydrologic conditions may be changing as sea levels continue to rise. As sea level rises, coastal wetland systems may either a) build elevation to maintain favorable hydrologic conditions for their survival or b) become increasingly inundated beyond their physiological tolerance resulting in an eventual conversion to open water. A better understanding of these processes will help to determine the present and future viability of coastal wetlands managed by the NPS and can help address measures to ensure these communities exist into the future. This protocol provides the reader with instructions and guidelines on designing a monitoring plan or study to: Quantify elevation change in wetlands with the Surface Elevation Table (SET); understand the processes that influence elevation change, including vertical accretion (SET and Marker Horizon methods); survey the wetland surface and SET mark to a common reference datum to allow for comparing sample stations to each other and to local tidal datums; survey the SET mark to monitor its relative stability; and establish water level equipment to characterize the hydrology and tidal datums for a particular wetland site or sites. The protocol is divided into two documents: The narrative (this document) presents an overview of all aspects of monitoring wetland elevation dynamics, and a document of Standard Operating Procedures (SOPs) with detailed instructions on the design, installation, data collection, and data management in support of monitoring changes in wetland elevation.

Science Report↗

USGS2023

Benchmarking high-resolution hydrologic model performance of long-term retrospective streamflow simulations in the contiguous United States

Because use of high-resolution hydrologic models is becoming more widespread and estimates are made over large domains, there is a pressing need for systematic evaluation of their performance. Most evaluation efforts to date have focused on smaller basins that have been relatively undisturbed by human activity, but there is also a need to benchmark model performance more comprehensively, including basins impacted by human activities. This study benchmarks the long-term performance of two process-oriented, high-resolution, continental-scale hydrologic models that have been developed to assess water availability and risks in the United States (US): the National Water Model v2.1 application of WRF-Hydro (NWMv2.1) and the National Hydrologic Model v1.0 application of the Precipitation–Runoff Modeling System (NHMv1.0). The evaluation is performed on 5390 streamflow gages from 1983 to 2016 ( ∼ 33 years) at a daily time step, including both natural and human-impacted catchments, representing one of the most comprehensive evaluations over the contiguous US. Using the Kling–Gupta efficiency as the main evaluation metric, the models are compared against a climatological benchmark that accounts for seasonality. Overall, the model applications show similar performance, with better performance in minimally disturbed basins than in those impacted by human activities. Relative regional differences are also similar: the best performance is found in the Northeast, followed by the Southeast, and generally worse performance is found in the Central and West areas. For both models, about 80 % of the sites exceed the seasonal climatological benchmark. Basins that do not exceed the climatological benchmark are further scrutinized to provide model diagnostics for each application. Using the underperforming subset, both models tend to overestimate streamflow volumes in the West, which could be attributed to not accounting for human activities, such as active management. Both models underestimate flow variability, especially the highest flows; this was more pronounced for NHMv1.0. Low flows tended to be overestimated by NWMv2.1, whereas there were both over and underestimations for NHMv1.0, but they were less severe. Although this study focused on model diagnostics for underperforming sites based on the seasonal climatological benchmark, metrics for all sites for both model applications are openly available online.

Hydrology and Earth System Sciences↗

USGS2022

Water availability drives instream conditions and life-history of an imperiled desert fish: A case study to inform water management

In arid ecosystems, available water is a critical, yet limited resource for human consumption, agricultural use, and ecosystem processes—highlighting the importance of developing management strategies to meet the needs of multiple users. Here, we evaluated how water availability influences stream thermal regimes and life-history expressions of Lahontan cutthroat trout ( Oncorhynchus clarkii henshawi ) in the arid Truckee River basin in the western United States. We integrated air temperature and stream discharge data to quantify how water availability drives stream temperature during annual spawning and rearing of Lahontan cutthroat trout. We then determined how in situ stream discharge and temperature affected adult spawning migrations, juvenile growth opportunities, and duration of suitable thermal conditions. Air temperatures had significant, large effects (+) on stream temperature across months; the effects of discharge varied across months, with significant effects (−) during May through August, suggesting increased discharge can help mitigate temperatures during seasonally warm months. Two models explained adult Lahontan cutthroat trout migration, and both models indicated that adult Lahontan cutthroat trout avoid migration when temperatures are warmer (~ > 12 °C) and discharge is higher (~ > 50 m 3 *s −1 ). Juvenile size was best explained by a quadratic relationship with cumulative degree days (CDD; days>4 °C) as size increased with increasing CDDs but decreased at higher CDDs. We also found an interaction between CDDs and discharge explaining juvenile size: when CDDs were low, higher discharge was associated with larger size, but when CDDs were high, higher discharge was associated with smaller size. Stream temperatures also determined the duration of juvenile rearing, as all juvenile emigration ceased at temperatures >24.4 °C. Together, our results illustrated how stream discharge and temperature shape the life-history of Lahontan cutthroat trout at multiple stages and can inform management actions to offset warming temperatures and facilitate life-history diversity and population resilience.

California, Nevada↗

USGS2024

Elemental composition and potential toxicity of the riverine macrophyte Podostemum ceratophyllum Michx. reflects land use in eastern North America

Land use influences surface water quality, often alleviating stoichiometric constraints on primary production and altering biogeochemical cycling. However, land use effects on nutrient content and potential trace metal accumulation in aquatic plants remain unclear, and high concentrations of metals and altered nutrient ratios could impact the health of herbivores and detritivores. We tested for land use effects on nutrient and trace metal accumulation in a widespread riverine macrophyte, Podostemum ceratophyllum , collected from 91 locations from Georgia to Maine, USA in 2014–2016. We quantified carbon (C), nitrogen (N), phosphorus (P), their molar and mass ratios, N and C stable isotopes, and 17 additional elements in dried plants collected from each location to estimate relationships between plant tissue content and watershed land use, which we quantified as agriculture, forest, and development. Decreasing forest cover was correlated with increasing δ 15 N, Mg, Mn, and P in Podostemum tissue. Increasing urban development was correlated with increasing δ 15 N, Mg and P, while increasing agriculture was correlated with a decrease in C: P and the concentrations of multiple metals, along with increases in P, Mg and δ 15 N. Decreases in ratios of N: P and C:P with increasing agriculture and urban development in the watershed indicate more rapid P storage relative to C and N in plant tissue, and increased resource quality of the plant to consumers in these watersheds. We also observed potentially toxic dietary concentrations of some trace metals (B, Cd, Tl, Zn) in plant tissue which could be related to the plant's natural herbivory defense system or to land use. We conclude that land use influences the elemental composition of P. ceratophyllum , and potentially the quality and toxicity of the plant to herbivores and detritivores in eastern North American rivers.

Science of the Total Environment↗

USGS2021

An approach for decomposing river water-quality trends into different flow classes

A number of statistical approaches have been developed to quantify the overall trend in river water quality, but most approaches are not intended for reporting separate trends for different flow conditions. We propose an approach called FN 2Q , which is an extension of the flow-normalization (FN) procedure of the well-established WRTDS (“Weighted Regressions on Time, Discharge, and Season”) method. The FN 2Q approach provides a daily time series of low-flow and high-flow FN flux estimates that represent the lower and upper half of daily riverflow observations that occurred on each calendar day across the period of record. These daily estimates can be summarized into any time period of interest (e.g., monthly, seasonal, or annual) for quantifying trends. The proposed approach is illustrated with an application to a record of total nitrogen concentration (632 samples) collected between 1985 and 2018 from the South Fork Shenandoah River at Front Royal, Virginia (USA). Results show that the overall FN flux of total nitrogen has declined in the period of 1985–2018, which is mainly attributable to FN flux decline in the low-flow class. Furthermore, the decline in the low-flow class was highly correlated with wastewater effluent loads, indicating that the upgrades of treatment technology at wastewater treatment facilities have likely led to water-quality improvement under low-flow conditions. The high-flow FN flux showed a spike around 2007, which was likely caused by increased delivery of particulate nitrogen associated with sediment transport. The case study demonstrates the utility of the FN 2Q approach toward not only characterizing the changes in river water quality but also guiding the direction of additional analysis for capturing the underlying drivers. The FN 2Q approach (and the published code) can easily be applied to widely available river monitoring records to quantify water-quality trends under different flow conditions to enhance understanding of river water-quality dynamics.

Virginia↗

USGS2016

Spatial and temporal patterns of mercury concentrations in freshwater fish across the Western United States and Canada

Methylmercury contamination of fish is a global threat to environmental health. Mercury (Hg) monitoring programs are valuable for generating data that can be compiled for spatially broad syntheses to identify emergent ecosystem properties that influence fish Hg bioaccumulation. Fish total Hg (THg) concentrations were evaluated across the Western United States (US) and Canada, a region defined by extreme gradients in habitat structure and water management. A database was compiled with THg concentrations in 96,310 fish that comprised 206 species from 4262 locations, and used to evaluate the spatial distribution of fish THg across the region and effects of species, foraging guilds, habitats, and ecoregions. Areas of elevated THg exposure were identified by developing a relativized estimate of fish mercury concentrations at a watershed scale that accounted for the variability associated with fish species, fish size, and site effects. THg concentrations in fish muscle ranged between 0.001 and 28.4 (μg/g wet weight (ww)) with a geometric mean of 0.17. Overall, 30% of individual fish samples and 17% of means by location exceeded the 0.30 μg/g ww US EPA fish tissue criterion. Fish THg concentrations differed among habitat types, with riverine habitats consistently higher than lacustrine habitats. Importantly, fish THg concentrations were not correlated with sediment THg concentrations at a watershed scale, but were weakly correlated with sediment MeHg concentrations, suggesting that factors influencing MeHg production may be more important than inorganic Hg loading for determining fish MeHg exposure. There was large heterogeneity in fish THg concentrations across the landscape; THg concentrations were generally higher in semi-arid and arid regions such as the Great Basin and Desert Southwest, than in temperate forests. Results suggest that fish mercury exposure is widespread throughout Western US and Canada, and that species, habitat type, and region play an important role in influencing ecological risk of mercury in aquatic ecosystems.

Science of the Total Environment↗

USGS2023

The influence of large woody debris on post-wildfire debris flow sediment storage

Debris flows transport large quantities of water and granular material, such as sediment and wood, and this mixture can have devastating impacts on life and infrastructure. The proportion of large woody debris (LWD) incorporated into debris flows can be enhanced in forested areas recently burned by wildfire, because wood recruitment into channels accelerates in burned forests. In this study, we explored how LWD transported in post-wildfire debris flows acts to retain sediment within small headwater channels. We found that the largest amount of debris flow sediment is stored in channel reaches with a ratio of LWD length to channel width between 0.25 and 1. Moreover, the largest debris flow deposits forced by LWD were found in low-angle channel reaches. We examined two models for debris flow volume estimation: (1) the current volume prediction model used in U.S. Geological Survey debris flow hazard assessments, and (2) a regional model developed to predict the sediment yield associated with debris-laden flows. We found that the regional model better matched the magnitude of the observed sediment at the terminal fan, suggesting the utility of regionally calibrated parameters for debris flow volume prediction. However, large wood created sediment storage upstream of the terminal fan, and this volume was of the same magnitude as the total volume at the terminal fans. Finally, we demonstrate a method for estimating debris flow velocity based on estimates of the critical velocity required to break wood, which can be used in future field studies to estimate minimum debris flow velocity values.

New Mexico↗