Search NASA⌕ Search

SEARCH · Search NASA

Results for “Runoff”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Role of underlying surface, rainstorm and antecedent wetness condition on flood responses in small and medium sized watersheds in the Yangtze River Delta region, China

Floods are becoming more frequent and are hard to control due to the shortage of water conservancy projects in small and medium sized basins, especially in developing regions. Understanding the hydrologic responses and estimated flood characteristics to storm events can help to predict flood disasters and form better effective mitigation and adaptation strategies. Here, using observations of several representative watersheds and Artificial Neural Networks (ANN) and Principal Component Regression (PCR) models, we conducted data-driven analyses to examine the flood responding characteristics. We quantified the relative contributions of the influencing factors to the variation of each flood characteristic. Statistical analysis of the observations shows that the drainage area plays a key role in determining the distribution of lag time and peak discharge. Rainstorm variability has direct influence on floods, and typhoon-induced rainstorms with high total rainfall and rainfall intensity generate higher lag time, flood peak, unit discharge and runoff depth, but lower runoff coefficient. The ANN and PCR models accurately predicted the variations of flood features using the driving factors including physical geographical characteristics, rainstorm features, and antecedent condition. Physical geographic characteristics are key influential factors of lag time, flood peak and runoff coefficient, while the rainstorm features control the magnitude of unit discharge and runoff depth. These results indicate that floods are mainly affected by rainstorm features and physical geographic characteristics in the Yangtze River delta, and it might become more damaging with the increasing rainfall extremes and sprawling impervious surfaces in a changing environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Coupled effects of future rainfall and land use on urban stormwater drainage system in Tampa, Florida (USA)

Future rainfall and land use are two important factors for flood management since they both directly and indirectly affect the functionality and performance of urban stormwater drainage infrastructures. Although some studies have already examined the coupled effects of future rainfall and land use conditions, they have not concluded how such effects on both surface runoff and flood areas would change with regional surface elevation variation and infrastructure conditions. This paper analyzed the coupled effects of future rainfall and land use on urban drainage systems in terms of surface runoff quantity and flood area changes using EPA SWMM. Future downscaled and bias-corrected precipitation projections for 2040–2060 and 2080–2099 from “Downscaled CMIP3 and CMIP5 Climate and Hydrology Projections” Archive and future land use for 2050 s and 2090 s from EPA ICLUS project were used for the City of Tampa in Florida (USA). It was shown that potential runoff volume and flood area changes ranged from –40% to 160%, and –40% to 400% due to rainfall change, while potential changes induced by land use change ranged from 0% to 3.5%, and 0% to 18% respectively. Additionally, this study examined the coupled effects of future rainfall and land use changes on flood area variations, considering the capacity of drainage infrastructure and elevation variation of study site. Results showed that flood area variations are not solely influenced by runoff quantity changes due to the two external drivers, but also by other factors. Specifically, in low-relief areas with inadequate drainage infrastructure and less elevation variation, these two drivers have no additive effects on the percentage of flood area changes. In contrast, in areas with adequate drainage infrastructure and greater surface elevation variation, there are additive effects from rainfall and land use changes on flood area expansions. These findings have important implications for stakeholders involved in city planning and flooding control, particularly for urban areas facing future precipitation and urbanization challenges.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Using SWMM for emergency response planning: A case study evaluating biological agent transport under various rainfall scenarios and urban surfaces

To assist in emergency preparedness for a biological agent terrorist attack or accidental pathogen release, potential contaminant levels and migration pathways of spores spread by urban stormwater were evaluated using a Storm Water Management Model (SWMM) of U.S. Coast Guard Base Elizabeth City, North Carolina. The high temporal-spatial resolution SWMM model was built using spore concentrations in stormwater runoff from asphalt, grass, and concrete collected from a point-scale field study. The subsequent modeled contamination scenarios included a notional plume release and point releases mimicking the field study under three rainfall conditions. The rainfall scenarios included a 6-hour natural rainfall event on Dec. 8, 2021 and two design storms (2-year and 100-year events). The observed spore concentrations from asphalt and concrete from the actual field experiment were applied to calibrate the washoff parameters in the SWMM model, using an exponential washoff function. The calibrated washoff coefficient (c 1 ) and exponent (c 2 ) were 0.01 and 1.00 for asphalt, 0.05 and 1.45 for grass, and 2.45 and 1.00 for concrete, respectively. The calibrated SWMM model simulated spore concentrations in runoff at times and magnitudes similar to the field study data. In the point release modeled scenario, the concrete surface generated 55.6% higher average spore concentrations than asphalt. Similarly, in the field experiment, a 175% (p < 0.05) higher average spore concentration in surface runoff was observed from concrete than from asphalt. Here, this study demonstrates how SWMM may be used to evaluate spore washoff from urban surfaces under different precipitation amounts, intensities, and durations, and how visualized spatial migration pathways in stormwater runoff may be used for emergency planning and remediation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Future Warming and Intensification of Precipitation Extremes: A “Double Whammy” Leading to Increasing Flood Risk in California

This study focuses on quantifying future anthropogenic changes in surface runoff associated with extreme precipitation in California's Sierra Nevada. The method involves driving a land surface model with output from a high resolution regional atmospheric simulation of the most extreme atmospheric rivers (ARs). AR events were selected from an ensemble of global climate model simulations of historical and late 21st century climate under the “high-emission” RCP8.5 scenario. Average precipitation during the future ARs increases by ~25% but a much lower proportion falls as snow. The resulting future runoff increase is dramatic—nearly 50%, reflecting both the precipitation increase and simultaneous conversion of snow to rain. The “double whammy” impact on runoff is largest in the 2,000–2,500 m elevation band, where the snowfall loss and precipitation increase are both especially large. Overall, this huge increase in runoff during the most extreme AR events could present major flood control challenges for the region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Validation of the Community Land Model Version 5 over the Contiguous United States (CONUS) using in situ and remote sensing data sets

The Community Land Model (CLM) is an effective tool to simulate the biophysical and biogeochemical processes and their interactions with the atmosphere. Although CLM Version 5 (CLM5) constitutes various updates in these processes, its performance in simulating energy, water and carbon cycles over the Contiguous United States (CONUS) at scales which land surface changes and hydrometeorological and hydroclimatological applications are more locally relevant is yet to be assessed. In this study, we conducted three simulations at 0.125? during 1979-2018 over the CONUS using different configurations of CLM, namely CLM5-biogeochemistry (CLM5BGC), CLM4.5BGC, and CLM5-satellite phenology (CLM5SP). We validated and compared their simulations against multiple remote-sensed and in-situ datasets. Overall, the parametric and structural updates (e.g., carbon cost for nitrogen uptake, variable soil thickness, dry surface layer) in CLM5 improve its ability in capturing terrestrial biogeochemical dynamics. The low evapotranspiration in CLM5BGC is associated with biases in simulating vegetation phenological characteristics rather than soil water limitations. The mismatch between CLM5BGC-simulated peak leaf area index and reference data can be attributed to CLM5BGC's inability in simulating phenology of trees and grasses. The differences between CLM-simulated irrigation and reference estimates can be attributed to differences between processes represented in models and in reality, and uncertainties in input and validation datasets. Evaluation against observations at small catchments suggest that hydrologic parameters needed to be calibrated to improve simulations of runoff, especially subsurface runoff. Additional efforts are needed to incorporate spatially-distributed plant phenology and physiology parameters and regional-specific agricultural management practices (e.g., planting, harvest).

Cheng, Yanyan↗
DOE OSTI2022

Causes of Missing Snowmelt Following Drought

Abstract Water management in snowy mountainous regions hinges on forecasting snowmelt runoff. However, droughts are altering snowpack‐runoff relationships with ongoing debate about the driving mechanisms. For example, in 2021 in California, less than half of predicted streamflow arrived. Mechanisms proposed for this “missing” streamflow included changes in evapotranspiration ( ET ), rainfall, and subsurface moisture conditions. Here, we demonstrate that ET in drought years generates dry subsurface conditions that reduce runoff in subsequent years. A model including this legacy of depleted moisture storage reduced median error in 2021 forecasts from 60% to 20% at 15 minimally disturbed basins and from 18% to 2% at 6 water supply basins in the Sierra Nevada (basins range in area from 5 to 23,051 km 2 and mean annual precipitation from 814 to 1,549 mm). Our findings indicate that the relationship between snowpack and runoff will evolve as plant ecosystems respond to climate change and alter subsurface water storage dynamics.

Lapides, Dana A.↗
DOE OSTI2023

Hydrologic Connectivity and Patch‐To‐Hillslope Scale Relations in Dryland Ecosystems

Abstract In drylands, runoff during storms redistributes water and nutrients from bare soil areas to vegetated patches, subsidizing vegetation with additional resources. The extent of this redistribution depends on the interplay between surface roughness and permeability; greater permeability in vegetated patches promotes run‐on to vegetation, but greater surface roughness diverts runoff, producing tortuous flow paths that bypass vegetation. Here, this interplay is examined in virtual experiments using the 2D Saint Venant Equations to measure runoff connectivity. Flowpaths are delineated using tracers advected by the flow. Distances between tracer sources and sinks along flowpaths measure hydrologic connectivity at two lengthscales: connectivity to the hillslope outlet and within‐slope source‐sink connectivity. Differences between these connectivity lengthscales indicate how flow may “by‐pass” vegetated patches within hillslopes. At the hillslope scale, a derived power‐law relation between the runoff coefficient and outlet connectivity describes hillslope water losses, providing a foundation for identifying landscapes likely to shed water.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Prolonged Drought in a Northern California Coastal Region Suppresses Wildfire Impacts on Hydrology

Wildfires naturally occur in many landscapes, however they are undergoing rapid regime shifts. Despite the emphasis in the literature on the most severe hydrological responses to wildfire, there remains a knowledge gap on the thresholds of wildfire (i.e., burned area/drainage area ratio, BAR) required to initiate hydrological responses. We investigated hydrological changes in the Russian River Watershed (RRW) in California, a coastal, Mediterranean, drought-prone, wildfire-adapted ecosystem, following ten wildfires that burned 30% of the watershed. Our findings suggest that sub-watersheds of the RRW have not burned beyond an intrinsic, unknown, threshold required to initiate change. Using paired watersheds, we examined spatiotemporal patterns of pre-and-post wildfire hydrology with a rainfall-runoff hydrological model. Even though these successive wildfires burned 1%–50% of each sub-watershed (1%–30% at moderate/high severity), we found little evidence of wildfire-related shifts in hydrology. As a function of BAR, wildfire imposed limited effects on runoff ratios (runoff/precipitation) and runoff residuals (observations—model simulations). Our findings that post-wildfire runoff enhancements asymptote beyond 30% burn indicate that when a watershed is burned beyond a certain threshold, the magnitude of the hydrologic response no longer increases. Drought and storm conditions explained much of the variability observed in streamflow, whereas wildfire explained only moderate variability in streamflow even when wildfire accounted for >45% BAR. While the BAR in the RRW was sufficiently beyond previously reported minimum disturbance thresholds (>20% burned forest), the lack of hydrological response is attributed to buffering effects of wildfire adaptation and drought factors that are unique to Mediterranean ecoregions.

Newcomer, Michelle E.↗
DOE OSTI2023

Mass‐Conserving Downscaling of Climate Model Precipitation Over Mountainous Terrain for Water Resource Applications

Abstract A mass‐conserving method to downscale precipitation from global climate models (GCMs) using sub‐grid‐scale topography and modeled 700‐mb wind direction is presented. Runoff is simulated using a stand‐alone hydrological model, with this and several other methods as inputs, and compared to runoff simulated using historical observations over the western contiguous United States. Results suggest the mitigation of grid‐scale biases is more critical than downscaling for some regions with large wet biases (e.g., the Great Basin and Upper Colorado). In other regions (e.g., the Pacific Northwest) the new method produces more realistic sub‐grid‐scale variability in runoff compared to unadjusted GCM output and a simpler downscaling method. The presented method also brings the runoff centroid timing closer to that simulated with observations for all subregions examined.

Geology↗
DOE OSTI2025

Impact of El Niño‐Southern Oscillation and Madden‐Julian Oscillation on the US Puget Sound Regional Hydroclimate

El Niño-Southern Oscillation (ENSO) and Madden-Julian Oscillation (MJO) are two major modes of climate variability with global hydroclimate impacts. However, their impacts often depend on the local climate and geography, resulting in large regional differences. In this study, we examined the connection of ENSO and MJO to the hydroclimate conditions and extremes in the Puget Sound (PS) basin located in the US Pacific Northwest coast. The results indicate that ENSO significantly modulates the cold season temperature and temperature-mediated hydrologic processes. El Niño cold seasons feature less snow accumulation and intensified surface runoff, even if the precipitation amount is similar to La Niña cold seasons. Therefore, El Niño causes more snow drought (in the form of compound dry and warm snow drought) and shifts the surface runoff seasonality by reducing runoff in the subsequent warm season. MJO phases 6–7 trigger more extreme precipitation, temperature, snowmelt, and runoff in the PS region at 0–9–day lags, and such connections are robust regardless of how the ENSO signals are removed. Meanwhile, MJO modulates large-scale extreme weather systems (e.g., atmospheric rivers) with significant enhancement during phases 6–7. ENSO impacts have intensified in the 2001–2020 period, whereas MJO impacts showed some phase shift in this period. This study reveals ENSO and MJO phases 6–7 as useful predictors of the PS hydroclimate anomalies/extremes at seasonal and daily scales, respectively. Utilizing these findings holds the potential to improve regional water resources prediction and management.

ENSO↗
DOE OSTI2023

The 1996 Mid-Atlantic Winter Flood: Exploring Climate Risk through a Storyline Approach

Abstract This article explores the application of thermodynamic perturbations to a historical midlatitude, wintertime, rain-on-snow flood event to evaluate how similar events may evolve under different climate forcings. In particular, we generate a hindcast of the 1996 Mid-Atlantic flood using an ensemble of 14-km variable-resolution simulations completed with the U.S. Department of Energy’s global Energy Exascale Earth System Model (E3SM). We show that the event is skillfully reproduced over the Susquehanna River Basin (SRB) by E3SM when benchmarked against in situ observational data and high-resolution reanalyses. In addition, we perform five counterfactual experiments to simulate the flood under preindustrial conditions and four different levels of warming as projected by the Community Earth System Model Large Ensemble. We find a nonlinear response in simulated surface runoff and streamflow as a function of atmospheric warming. This is attributed to changing contributions of liquid water input from a shallower initial snowpack (decreased snowmelt), increased surface temperatures and rainfall rates, and increased soil water storage. Flooding associated with this event peaks from around +1 to +2 K of global average surface warming and decreases with additional warming beyond this. There are noticeable timing shifts in peak runoff and streamflow associated with changes in the flashiness of the event. This work highlights the utility of using storyline approaches for communicating climate risk and demonstrates the potential nonlinearities associated with hydrologic extremes in areas that experience ephemeral snowpack, such as the SRB.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

Dayflow-PR: High-Resolution Streamflow Reanalysis for Puerto Rico, Version 1.0

This dataset presents a high-resolution historical streamflow reanalysis for NHDPlusV2 stream reaches across Puerto Rico (PR) spanning 1950 - 2019. The reanalysis is generated using the calibrated VIC-RAPID hydrologic modeling framework at the Hydrologic Unit Code Sub-basin (HUC08) scale, forced with sub-daily and daily meteorological forcings from Daymet. Runoff is simulated on 1- and 6-km grids, and the resulting total runoff is routed through the NHDPlusV2 river network using the RAPID routing model to produce Naturalized Streamflow Reanalysis. Where complete observational records are available over 1980 - 2019, streamflows are assimilated (substituted) and subsequently routed downstream through the river network to produce Assimilated Streamflow Reanalysis. The dataset includes streamflow outputs from eight distinct hydrologic modeling configurations along with key performanc evaluation metrics at daily and monthly scales, supporting a wide range of water resource applications. This dataset is derived to support the Non-Powered Dam Assessment, as well as 9505 Secure Water Assessment projects for the US Department of Energy (DOE) Water Power Technologies Office (WPTO). For further details, refer to Ghimire et al. (2023), Kao et al. (2024), and Ghimire et al. (2025).

13 HYDRO ENERGY↗
DOE OSTI2026

NGEE Arctic Integrated Modeling (IM2): Improved subgrid hillslope hydrologic connectivity

This data product represents the integration of new code capability for arctic tundra hillslope hydrologic processes into the Energy Exascale Earth System Model (E3SM), through the E3SM Land Model (ELM) component. This code integration is the result of collaborative effort between the NGEE Arctic project and the E3SM project. The current ELM represents water movement primarily through vertical processes, such as precipitation, canopy interception, evaporation, infiltration, and soil water movement. Lateral water movement—such as surface runoff, subsurface flow, and river transport—plays a significant role in the hydrological cycle, especially in regions with varied topography. While E3SM includes a runoff routing component representing water transport in the river network, the lateral transport of water at the subgrid scale within the land model has previously not been taken into account. With the recent development of topographic units within the ELM subgrid data structure, there is an opportunity to simulate hillslope hydrologic connectivity by introducing water transport along topographic gradients. We expect that more realistic representation of hillslope hydrologic processes will lead to improved predictions of both soil water content and river network flows. Lateral transport of water at and near the surface is represented as a sub-grid process in this new code development. Water is tracked as it moves from higher to lower elevations within a gridcell. This capability uses the nested hierarchical sub-grid scheme within ELM to connect water fluxes from sub-grid elements with higher elevation to those with lower elevation. This data record consists of a single document (pdf format) that describes the theoretical basis for the hillslope hydrology processes added to ELM, and describes the modifications made to the ELM code. The Methods section of this metadata record includes a link to the public E3SM code repository where the exact code modifications as integrated in E3SM can be accessed. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

Thornton, Peter E [ORNL] (ORCID:0000000247595158)↗
DOE OSTI2026

VIC-Global Parameter Dataset Sensitivity with the Variable Infiltration Capacity Model: Evaluating the importance of dynamic land surface parameters when using the VIC-Global parameter dataset

Accurate prediction of runoff is essential to water resources management, flood risk assessment, and ecosystem protection. However, many hydrological models still have relatively substantial limitations when representing the influence of land use and land cover (LULC) on runoff generation and routing. Changes in LULC, such as deforestation, urban expansion, agricultural intensification, and wetland loss, have been shown to alter the water balance at the land surface through fundamental hydrologic processes (e.g., interception, infiltration, evapotranspiration, and soil storage). However, it remains an open question what the exact magnitude and timing of these impacts are for the spatial and temporal scales commonly used in engineering applications. In this analysis we focus on one aspect of recent LULC change for assessing human impacts, which is urbanization. Specifically we seek to determine the impacts of urbanization on the magnitude and timing of surface runoff and baseflow in HUC-12 basins in Clark County, Nevada which has experienced rapid urbanization. We use the Variable Infiltration Capacity (VIC) hydrology model with a widely used off-the-shelf dataset of land surface parameters, VIC-Global, both of which have been commonly used in the past for water and energy balance modeling for large scale hydrologic studies. We examine two scenarios where the first scenario removes all urbanized land cover and parameterizes those areas of the basins as barren or open shrubland. The second scenario tests the opposite case where all areas of the basins are classified as urban regardless of their present classification. The results from the VIC model show there is a low sensitivity for daily surface runoff between scenarios. The daily baseflow values indicate similar low sensitivity to the classification change during specific periods, but then have substantial differences during other period when large precipitation events are occurring. This is likely due to the assumed parameter values for the urban land cover classification made by the VIC-Global dataset. Using a static land cover parameterization is reasonable for large domain hydrology models that are being used for near-term planning horizons (<30 years). However, longer planning horizons where feedbacks between the atmosphere and land surface are important, especially in transient climate situations, considerations for how to update land surface parameters should be incorporated.

42 ENGINEERING↗
DOE OSTI2024

Disentangling the hydrological and hydraulic controls on streamflow variability in Energy Exascale Earth System Model (E3SM) V2 – a case study in the Pantanal region

Abstract. Streamflow variability plays a crucial role in shaping the dynamics and sustainability of Earth's ecosystems, which can be simulated and projected by a river routing model coupled with a land surface model. However, the simulation of streamflow at large scales is subject to considerable uncertainties, primarily arising from two related processes: runoff generation (hydrological process) and river routing (hydraulic process). While both processes have impacts on streamflow variability, previous studies only calibrated one of the two processes to reduce biases in the simulated streamflow. Calibration focusing only on one process can result in unrealistic parameter values to compensate for the bias resulting from the other process; thus other water-related variables remain poorly simulated. In this study, we performed several experiments with the land and river components of the Energy Exascale Earth System Model (E3SM) over the Pantanal region to disentangle the hydrological and hydraulic controls on streamflow variability in coupled land–river simulations. Our results show that the generation of subsurface runoff is the most important factor for streamflow variability contributed by the runoff generation process, while floodplain storage effect and main-channel roughness have significant impacts on streamflow variability through the river routing process. We further propose a two-step procedure to robustly calibrate the two processes together. The impacts of runoff generation and river routing on streamflow are appropriately addressed with the two-step calibration, which may be adopted by developers of land surface and earth system models to improve the modeling of streamflow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Hydrological analysis in watersheds with a variable-resolution global climate model (VR-CESM)

Traditionally, watershed-scale hydrology is simulated by distributed hydrological models with offline meteorological forcing data, or by regional regional climate models that link atmospheric and land hydrology interactions. Global climate models (GCMs) are rarely used to study watershed-scale hydrology due to the relatively coarse grid resolution, computationally expensive downscaling, and simplified physical processes. Recently, however, watershed-scale hydrology analysis has become possible in GCMs due to the development of variable-resolution GCMs that dynamically couple the hydrological processes between atmosphere and land systems at fine resolutions in selected regions and coarse resolution elsewhere. In this study, we used the variable-resolution Community Earth System Model (VR-CESM) with refined-resolution (14 km) in the western U.S. and eastern China to evaluate watershed-scale hydrology. Furthermore, we compared the historical VR-CESM outputs with gauge measurements and other hydrological models (e.g., National Water Model in the U.S.) and calibrated the subsurface runoff capacities in four mountainous watersheds. An RCP8.5 projection from 2007 to 2050 is used to estimate the impact of changing precipitation and snow climatology on watershed hydrology. We also analyzed the long-term runoff variability and the possibility of extreme runoff events as simulated by the VR-CESM. Although calibration is not possible in larger-scale watersheds, VR-CESM simulates the long-term annual variability of watersheds and provides insights on climate change impacts on hydrology. We conclude that refined-resolution VR-CESM can be used for watershed-scale hydrology analysis to understand water resources and water balance, in addition to traditional watershed-scale hydrological models. It enables hydrological analysis at multiple watersheds in one simulation and can help understand the two-way dynamics between land surface hydrology and atmospheric processes, and is especially practical for projecting climate change impacts. However, it is challenging to apply VR-CESM for hydrologic analysis in regulated watersheds as human factors (e.g., pumping, irrigation, water diversion) have not been fully addressed in VR-CESM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Integrating Contaminant Source Indicators, Water Quality Measures, and Ecotoxicity to Characterize Contaminant Mixtures and Per- and Polyfluoroalkyl Substance (PFAS) Variability in an Urban Watershed

Thousands of chemical contaminants threaten watersheds but are time and cost prohibitive to monitor. Identifying their sources, transport, and ecological risk is limited in heterogeneous urban watersheds. We present an integrative watershed approach using source-specific indicator compounds, common water quality measures, and ecotoxicity assays to examine the distribution of contaminant mixtures in an urbanized watershed. Indicator compound concentrations were temporally and spatially distributed for treated/untreated sewage (sucralose, artificial sweetener), road runoff (diphenyl-guanidine [DPG] and 6PPD-quinone [6PPD-Q], automobile tire additives), and lawncare runoff (aminomethanephosphonic acid (AMPA), major degradant of the herbicide glyphosate). Sucralose was predominately sourced from treated wastewater; measurable concentrations in tributaries indicated raw sewage inputs. DPG and 6PPD-Q concentrations correlated to road density during base flow and were elevated during stormflow. AMPA was measurable spring through fall, especially where lawns were dense. When specific sources dominated flow, water quality measures correlated with wastewater (sulfate, potassium, chloride, and sodium) and road runoff (chromium and lead) indicators. The limited behavioral toxicity observed in exposed zebrafish (Danio rerio) (18%) was not well explained by source-indicators. PFAS concentrations were highly variable spatially but not well explained by our source-specific indicator compounds. Here, more costly compound-specific monitoring may be necessary when multiple sources exist or when unexpected toxicity trends occur.

computer simulations↗