Search NASA⌕ Search

SEARCH · Search NASA

Results for “Surface 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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The implications of future climate change on the blue water footprint of hydropower in the contiguous US

As the largest renewable energy source, hydropower is essential to the sustainability of the global energy market. However, a considerable amount of water can be lost in the form of evaporation from the associated multipurpose reservoirs, and hence enlarge the blue water footprint (BWF) of hydropower in a warming climate. To facilitate the sustainable management of both water and energy resources under the impact of climate change in the contiguous United States (CONUS), the BWF values of 143 major multipurpose reservoirs were evaluated during the historical period (1985–2014) and two future periods (2020–2049 and 2070–2099). The historical reservoir evaporation loss was calculated using the Landsat-based reservoir surface area and a new evaporation rate algorithm that considers the heat storage effect. Future projections of runoff availability, hydropower generation, and reservoir evaporation were estimated based on the downscaled climate model ensemble from phase 5 of the Coupled Model Intercomparison Project (CMIP5). It was found that the BWF for the CONUS is highly spatially heterogeneous, with an average value of 26.2 m 3 /MWh in the historical period. In the future, the BWF values are projected to increase under both Representative Concentration Pathway (RCP) 4.5 and 8.5 scenarios. This is especially noticeable under RCP 8.5, which has an average BWF value of 30.2 m 3 /MWh for 2070–2099 (increasing by 15.3% from 26.2 m 3 /MWh). The uncertainty ranges increase even more, from 3.4 m 3 /MWh during 2020–2049 to 5.7 m 3 /MWh during 2070–2099. These findings can benefit water and energy resources management in identifying suitable environmental, economic, operational, and investment strategies for multipurpose reservoirs in a changing environment.

13 HYDRO ENERGY↗
DOE OSTI2026

Public water supply infrastructure extensification and diversification in surface waters is insufficient to meet future demands in Texas

The data were developed to evaluate the capacity of existing and potential new surface water supply infrastructure to meet projected public water demands across districts in Texas under multiple future socioeconomic and climate scenarios. The database integrates hydrologic, water quality, infrastructure, energy, cost, demographic, and demand-projection information for candidate surface water supply locations. Candidate sites include stream reaches, waterbodies, reservoir surplus locations, and potential new reservoir sites. Water availability is characterized using historical and projected flow conditions, while site suitability is evaluated using five indicators: Water Availability Index (WAI), Water Quality Index (WQI), Energy Requirement Index (ERI), Water Treatment Cost (WTC), and Water Infrastructure Cost (WIC). The datasets include statewide candidate-site information, district-level demand projections under Shared Socioeconomic Pathways (SSPs), runoff-based allocation constraints, climate-stress metrics, and optimization outputs evaluating alternative infrastructure planning strategies. Optimization results compare Business-as-Usual (BAU) and All Surface Water (AllSW) demand-management approaches under both scaled and fixed cost-cap strategies. Associated validation datasets provide district-level feasibility assessments, infrastructure selection outcomes, cost-cap utilization, demand satisfaction metrics, and constraint diagnostics. Additional datasets quantify projected changes in storage and flow conditions as well as water availability stress for both existing and newly selected intake locations under the SSP5 scenario for mid-century and late-century climate conditions. Together, these datasets support assessment of the extent to which surface-water infrastructure expansion and diversification strategies can satisfy future public water demands while accounting for hydrologic, economic, and planning constraints across Texas. Dataset(s) Description Dataset_preoptimization.xlsx Comprehensive pre-optimization dataset containing candidate water-supply sites and associated hydrologic, water-quality, infrastructure, climate, demographic, runoff, and demand-projection variables used as inputs to the optimization analyses. Includes variable descriptions and the full statewide candidate-site database. District_level_site_selection.zip - Compressed archive containing all SSP-specific district-level optimization result files MESIO_ssp1_results.xlsx District-level site selection results for SSP1 (MESIO). Includes variable descriptions, BAU and AllSW site-selection results under scaled and fixed cost strategies, and district-level validation diagnostics. MESID_ssp2_results.xlsx District-level site selection results for SSP2 (MESID). Includes variable descriptions, BAU and AllSW site-selection results under scaled and fixed cost strategies, and district-level validation diagnostics. LCMRD_ssp3_results.xlsx District-level site selection results for SSP3 (LCMRD). Includes variable descriptions, BAU and AllSW site-selection results under scaled and fixed cost strategies, and district-level validation diagnostics. IRDev-Low_ssp4l_results.xlsx District-level site selection results for SSP4-Low (IRDev-Low). Includes variable descriptions, BAU and AllSW site-selection results under scaled and fixed cost strategies, and district-level validation diagnostics. IRDev-High_ssp4h_results.xlsx District-level site selection results for SSP4-High (IRDev-High). Includes variable descriptions, BAU and AllSW site-selection results under scaled and fixed cost strategies, and district-level validation diagnostics. RSIM_ssp5_results.xlsx District-level site selection results for SSP5 (RSIM). Includes variable descriptions, BAU and AllSW site-selection results under scaled and fixed cost strategies, and district-level validation diagnostics. tx_hydrological_stress.xlsx Hydrological stress dataset for existing and newly selected intake locations. Includes projected mid-century and late-century changes, gain/loss classifications, planning strategy information, and accompanying variable descriptions. Also includes water-stress metrics derived from historical and projected low-flow conditions.

Okoye, Perpetua I. (ORCID:0000000215545033)↗
DOE OSTI2025

Computational Modeling of Water Infiltration From Ditches Into Roadway Embankments

Ditches are built along roadways to convey the stormwater runoff. They collect the runoff from the pavement and adjacent drainage areas. Water can move freely in a ditch away from the road surface, thus making it safer for the vehicular traffic. South Carolina DOT expressed interest in studying infiltration of water from ditches into roadway embankments. As part of the current design criteria, it is recommended that the road subgrades be 1.0 foot above the design event water surface elevation in the ditch. The distance between the water surface and bottom of the subbase is referred to as “freeboard”. The purpose of this requirement is to lower the risk of pavement damage by moving the flow away from the pavement, containing the captured flow, and conveying flow received from adjacent properties to outfalls. With changes in roadway design standards, larger shoulders and flatter cross-slopes have moved ditches further away from the pavement and the subbase. These changes could allow for a lower freeboard standard without reducing safety margins due to the increased lateral distance from the subbase to the ditch. The objective of this study is to evaluate whether the current SCDOT criterion requiring the roadway subgrade to be 1 ft above the design-event ditch water surface can be relaxed under selected geometric and soil conditions, particularly where larger shoulder widths and flatter side slopes increase horizontal separation between ditch water and pavement structure.

42 ENGINEERING↗
DOE OSTI2022

Storm and Annual Time Scale Hydrological Data for the Russian River Watershed 1996-2022

These files contains observed and simulated hydrological data (discharge, precipitation) that have been aggregated to the storm-event scale and to the annual time scale for the Russian River Watershed (RRW), California. Observed data are obtained from 12 different USGS hydrological stations located throughout the watershed, and discharge is simulated using the GR5H hourly hydrological model. Precipitation data are obtained from the NASA NLDAS NOAH Community Land Surface Model. All simulated and observed discharge are aggregated to the storm and annual time scales. We investigated changes in watershed hydrological conditions in the Russian River Watershed, a Mediterranean, drought prone, wildfire-adapted ecosystem, following eleven wildfires that occurred from 2017-2020. We ask two research questions: 1) How do interacting wildfire events, drought, and atmospheric rivers impact hydrological conditions of the watershed and in particular streamflow?, and 2) What percentage of new and/or cumulative wildfire disturbance is required to initiate hydrological change? We hypothesize that sub-watersheds of the RRW Mediterranean ecoregion have not burned beyond an intrinsic, and still unknown, threshold required to initiate change. Using a series of paired burned/unburned catchments nested within the larger watershed, we examined temporal and spatial patterns of pre-and-post wildfire water hydrology using a rainfall runoff hydrological model compared with data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Projected Changes to Hydroclimate Seasonality in the Continental United States

Abstract Future changes to the hydrological cycle are projected in a warming world, and any shifts in drought risk may prove extremely consequential for natural and human systems. In addition to long‐term moistening, drying, or warming trends, perturbations to the annual cycle of regional hydroclimate variables may also have substantial impacts. We analyze projected changes in several hydroclimate variables across the continental United States, along with shifts in the amplitude and phase of their annual cycles. We find that even in regions where no robust change in the annual mean is expected, coherent changes to the annual cycle are projected. In particular, we identify robust regional phase shifts toward earlier arrival of peak evaporation in the northern regions, and peak runoff and total soil moisture in the western regions. Changes in the amplitude of the annual cycle of total and surface soil moisture are also projected, and reflect changes to the annual cycle in surface water supply and demand. Whether changes become detectable above the background noise of internal variability depends strongly on the future scenario considered, and significant changes to the annual cycle are largely avoided in the lowest‐forcing scenario.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Comparison of potential drinking water source contamination across one hundred U.S. cities

Abstract Drinking water supplies of cities are exposed to potential contamination arising from land use and other anthropogenic activities in local and distal source watersheds. Because water quality sampling surveys are often piecemeal, regionally inconsistent, and incomplete with respect to unregulated contaminants, the United States lacks a detailed comparison of potential source water contamination across all of its large cities. Here we combine national-scale geospatial datasets with hydrologic simulations to compute two metrics representing potential contamination of water supplies from point and nonpoint sources for over a hundred U.S. cities. We reveal enormous diversity in anthropogenic activities across watersheds with corresponding disparities in the potential contamination of drinking water supplies to cities. Approximately 5% of large cities rely on water that is composed primarily of runoff from non-pristine lands (e.g., agriculture, residential, industrial), while four-fifths of all large cities that withdraw surface water are exposed to treated wastewater in their supplies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Machine Learning Calibration of Groundwater Table Depth in ELM: Impact on Land Surface Hydrology and Land‐Atmosphere Fluxes

Accurate representation of groundwater table depth (GWTD) is crucial for simulating hydrological cycling in Earth system models (ESM). Nevertheless, there is a notable gap in the literature regarding the validation of GWTD simulations in ESMs and their subsequent impact on downstream hydrological components. This study explores the calibration of parameterization of global GWTD using machine learning within the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Despite achieving significant gains in simulating GWTD through calibration, offline ELM simulations unexpectedly show that these improvements do not translate to substantial enhancements in model performance for other key hydrological variables, including soil moisture (SM), runoff, groundwater contribution to runoff or base flow index (BFI), and evapotranspiration and its partitioning. The performance in SM and runoff was even degraded in some regions, while BFI was mostly overestimated. Although there is significant improvement in GWTD within the critical range of 1–5 m, where groundwater traditionally influences land surface energy fluxes, these improvements occurred mostly in humid areas where the impact of GWTD on surface processes is minimal. Although the impacts of model calibration are generally small in offline ELM simulations, coupled land-atmosphere simulations exhibit much stronger responses to GWTD calibration, highlighting the role of land-atmosphere feedbacks in Earth system modeling. These findings underscore the need for integrated calibration strategies that simultaneously optimize multiple hydrological variables. However, if a single-variable approach is necessary, it is crucial to establish clear priorities for calibration, identifying the most critical variables that have the greatest impact on overall model performance.

Fang, Yilin [Pacific Northwest National Laboratory↗
DOE OSTI2024

Wildland–urban interface wildfire increases metal contributions to stormwater runoff in Paradise, California

The 2018 Camp Fire was a large late-year (November) wildfire that produced an urban firestorm in the Town of Paradise, California, USA, and destroyed more than 18 000 structures. Runoff from burned wildland areas is known to contain ash, which can transport contaminants including metals into nearby watersheds. However, due to historically infrequent occurrences, the effect of wildland-urban interface (WUI) fires, such as the Camp Fire, on surface water quality has not been well-characterized. Therefore, this study investigated the effects of widespread urban burning on surface water quality in major watersheds of the Camp Fire area. Between November 2018 and May 2019, 140 surface water samples were collected, including baseflow and stormflow, from burned and unburned watersheds with varying extent of urban development. Samples were analyzed for total and filter-passing metals, dissolved organic carbon, major anions, and total suspended solids. Ash and debris from the Camp Fire contributed metals to downstream watersheds via runoff throughout the storm season. Increases in concentration up to 200-fold were found for metals Cr, Cu, Ni, Pb, and Zn in burned watersheds compared to pre-fire values. Total concentrations of Al, Cd, Cu, Pb, and Zn exceeded EPA aquatic habitat acute criteria by up to 16-fold for up to five months after the fire. To assess possible transport mechanisms and bioavailability, a subset of 18 samples was analyzed using four filters with nominal pore sizes ranging from 0.22 to 1.2 μm to determine the particulate size distribution of metals. Trace and major metals (Al, Ba, Co, Cr, Cu, Fe, Hg, Mn, Ni, Pb, and Zn) were found mostly associated with larger grain sizes (>0.45 μm), and some metals (Al, Cr, Fe, and Pb) also included a substantial colloidal phase (0.22 to 0.45 μm). This study suggests that fires in the wildland–urban interface increase metal concentrations, mainly through particulate driven transport. The metals with the largest increases are likely from anthropogenic disaster materials, though biomass ash also is a major contributor to water quality. In conclusion, the increase in metals following WUI burning may have adverse ecological impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Permafrost Thaw, Uneven Subsidence and Projected Drying of Ice-wedge Polygon Tundra: Modeling Archive

This dataset is a model archive of the paper Permafrost Thaw, Uneven Subsidence and Projected Drying of Ice-wedge Polygon Tundra (in prep) to support a modeling study investigating how projected increases in Arctic temperature and precipitation will jointly influence hydrologic conditions in ice-rich tundra landscapes. With this dataset, this study is to address the research question: Will Arctic tundra landscapes become wetter or drier with increasing precipitation and temperature in the future when thaw-induced ground subsidence and associated microtopographic evolution are represented? The simulations focus on ice-wedge polygon tundra, a widespread form of ice-rich permafrost terrain that is highly sensitive to thaw-driven landscape change. This dataset contains model input and output data for four study watersheds in Alaska: Anaktuvuk, Utqiagvik (formerly Barrow), Brooks Foothills, and Prudhoe Bay. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.5), a physics-rich integrated surface–subsurface hydrologic model. For each watershed, ten modeling cases were performed representing two landscape evolution conditions (with subsidence and without subsidence) combined with five climate forcing scenarios derived from Shared Socioeconomic Pathways (SSP5, SSP5 with precipitation trend, SSP2, SSP2 with precipitation trend, and SSP2 with double precipitation trend). Particularly, for each watershed under the forcing SSP2 with precipitation trend, there are two additional simulations considering spatially heterogeneous subsidence distributions: one assumes randomly distributed scaling and the other includes elevation dependent distribution scaling. These simulations span 1980–2099 and include spin-up runs (1980–2009) followed by transient projections (2010–2099). To facilitate reproducibility of simulations, all datasets are organized by watershed. For each study watershed, the dataset contains: (1) Pre-partitioned mesh files for 32-core modeling (.par.32.XX), located in EACH_WATERSHED/mesh/basin; and also a non-partitioned mesh file (.exo) located in EACH_WATERSHED/mesh; (2) Climate forcings corresponding to the five SSP scenarios (.h5), located in EACH_WATERSHED/data; (3) Final states (.h5) from column spin-up modeling used to initialize historical watershed-scale spin-up runs from 1980 to 2009, located in EACH_WATERSHED/PreSpinupHistorical; (4) Final states (.h5) of historical watershed-scale spin-up runs from 1980 to 2009 used to initialize projection runs, located in EACH_WATERSHED/Spinup_daymetERA5; (5) ATS modeling input files (.xml), located in EACH_WATERSHED/EACH_SIMULATION_SCENARIO/inputfiles; (6) ATS modeling output files (.dat), located in in EACH_WATERSHED/EACH_SIMULATION_SCENARIO/combined_obs; (7) For the Brooks Foothills watershed, additional spatial model outputs are provided (.h5) for selected years (2033 and 2093) used to generate spatial figures in this study, located in Brooksfoothills/EACH_SIMULATION_SCENARIO/results-WITH/WITHOUT_SUBSIDENCE-year2033/2093. All data files with suffix .h5 can be accessible through Python h5py, and all data files with suffix of .dat can be imported by Python pandas. Mesh file with .exo can be visualized through Paraview or read by Python netCDF. The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research. Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

EARTH SCIENCE > ATMOSPHERE > PRECIPITATION↗
DOE OSTI2020

Assessment of Spatiotemporal Variability of Evapotranspiration and Its Governing Factors in a Mountainous Watershed at the East River, Colorado, Water: Dataset

This data package contains simulation results of daily evapotranspiration for 17 locations of meteorological stations with different precipitation, air temperature and vegetation within the East River Watershed, Colorado (CLMResult.csv). Daily precipitation data is included for reference. These are the main results of the paper Tran, 2019. Meteorological stations names, latitude and longitude coordinates of meteorological stations, types of lands from the National Land Cover Database, and a comparison with the community land model (CLM) data are also included (Met_LC.csv). The purpose of this data is to analyze the spatio-temporal variations of evapotranspiration at the East River Watershed and understand the factors that control these variations. More information regarding the study, methods, and findings can be found in the Related Reference paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Mountain Basin Controls on the Snow-to-Streamflow Signal: An AIC-Weighted Multiple Linear Regression Framework

A regression-based analysis quantifies how basin characteristics modulate the snow-to-streamflow signal. First, we use the ERA5-Land reanalysis gridded product (European Centre for Medium Range Weather Forecasts reanalysis 5 -Land component) for 4,655 hydrologic unit code - 10 (HUC10) mountain basins across the western United States (US) for water years 1987–2024. Linear regressions are performed for peak snow water equivalent (SWE) and annual streamflow for each mountain basin. Models use ordinary least squares in Python’s statsmodels package. After which, an Akaike Information Criterion (AIC)–weighted ensemble multiple linear regression (MLR) framework with 47 watershed traits is used to predict the linear regression coefficient of determination (r-squared) defining the ability of peak SWE to predict annual streamflow across all mountain basin. Predictor sets are constrained to avoid multicollinearity by excluding models with variance inflation factors (VIF) greater than 5. Mountain basin traits included in the MLR include seasonal climate, topography, vegetation type and structure, and bedrock geology. Accepted models are considered if their AIC is within 2.0 of the model with the minimum AIC, or best model. To compare predictor influence across acceptable models, we computed standardized regression coefficients. To evaluate structural redundancy among models, we constructed binary inclusion vectors for each acceptable model, denoting whether a predictor was present (1) or absent (0). Core predictor variables are defined as occurring in at least 67% of the acceptable models. For this regional analysis, only one model was found acceptable, with higher snow-to-streamflow translation (higher r-squared) occurring in colder mountain basins with higher relative winter precipitation, more snow accumulation and a lower fraction of annual precipitation that falls in the spring and summer. The second component of the data package uses previously published, high-resolution output from an integrated hydrological model of the East River watershed using the U.S. Geological Survey Groundwater and Surface water Flow model (GSFLOW, doi:10.15485/1998576). East River MLR expands upon the approach described above to explore the response of five streamflow metrics—annual streamflow, runoff efficiency, 7-day minimum flow, low-flow duration, and non-perennial stream fraction to snow system indicators including peak SWE, snow-covered area, snow disappearance date, and the fraction of basin area characterized by low-to-no snow, as well as seasonal precipitation and temperature, and annual hydrologic variables representing soil moisture, evapotranspiration (ET), the partitioning of incoming precipitation to evapotranspiration (ET/P), groundwater storage, and groundwater inflow to streams. MLR was done on all water years (P0: 1987-2024) and for each period as determined in the split analysis using pooled regression techniques (P1: 1987-2011 and P2: 2012-2024) to evaluate shifting predictor variable emphasis on streamflow generation. Results indicate that since 2012, peak SWE has lost statistical strength in its prediction of annual streamflow and runoff efficiency, and the indirect influence of spring temperature has emerged as critically important. Low-flow metrics remain largely influenced by soil moisture, vegetation water use and groundwater inflows with summer precipitation becoming a direct influence on minimum summer flow. Together, these data and Python-based analysis tools provide a framework for identifying the key watershed characteristics that control how streamflow responds to snow from year to year. The package also helps quantify uncertainty in statistical models and assess how snow–streamflow relationships vary across regions and over time. This dataset contains comma-separated values files (.csv), text files (.txt), python code files (.py), figure files (.png), and shapefiles (.cpg, .dbf, .prj, .sbn, .sbx, .shp, .xml). Further details on file contents and MLR execution can be found in the readme file and the FLMD files. Work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Impacts of Topography-Based Subgrid Scheme and Downscaling of Atmospheric Forcing on Modeling Land Surface Processes in the Conterminous US

The effects of small-scale topography-induced land surface heterogeneity are not well represented in current Earth System Models (ESMs). In this study, a new topography-based subgrid structure referred to as topographic units (TGU) designed to better capture subgrid topographic effects, and methods to downscale atmospheric forcing to the land TGUs have been implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Effects of the subgrid scheme and downscaling methods on ELM simulated land surface processes are evaluated over the conterminous United States (CONUS). For this purpose, ELM simulations are performed using two configurations without (NoD ELM) and with (D ELM) downscaling, both using TGUs derived for the 0.5-degree grids and the same land surface parameters. Simulations using the two ELM configurations are compared over the CONUS domain, regional levels, and at observational sites (e.g., SNOTEL). The CONUS-level results suggest that D ELM simulates more snowfall and snow water equivalent (SWE), higher runoff, and less ET during spring and summer. Regional-level results suggest more pronounced impacts of downscaling over regions dominated by higher elevation TGUs and regions with maximum precipitation occurring during cool seasons. Results at the SNOTEL sites suggest that D ELM has superior capability of reproducing the observed SWE at 83% of the sites, with more pronounced performance over topographically heterogeneous TGUs with their maximum precipitation occurring during cool seasons. The results highlight the importance of improving representation of small-scale surface heterogeneity in ESMs and motivate future research to understand their effects on land-atmosphere interactions, streamflow, and water resources management over mountainous regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Nitrate loading projection is sensitive to freeze-thaw cycle representation

Climate change can have substantial impacts on nitrogen runoff, which is a major cause of eutrophication, harmful algal blooms, and hypoxia in freshwaters and coastal regions. We examined responses of nitrate loading to climate change in the Upper Mississippi River Basin (UMRB) with an enhanced Soil and Water Assessment Tool with physically based Freeze-Thaw cycle representation (SWAT-FT), as compared with the original SWAT model that employs an empirical equation. Driven by future climate projections from five General Circulation Models (GCMs) from 1960 to 2099 under the Representative Concentrations Pathways (RCP) 8.5 scenario, we analyzed changes in riverine nitrate loadings, as well as terrestrial surface and subsurface contributions of the UMRB in the 21st century relative to the baseline period of 1960–1999. By the end of the 21st century, the original SWAT model predicted about a 50% increase in riverine nitrate loadings which is nearly twice as much as that estimated by SWAT-FT (ca. 25%). Such a large difference in projected nitrate changes can potentially mislead mitigation strategies that aim to reduce nitrogen runoff from the UMRB. Further analysis shows that the difference between the original SWAT model and SWAT-FT led to substantial discrepancies in the spatial distribution of surface and subsurface nitrate loadings in the UMRB. In general, SWAT-FT predicted more nitrate leaching for northwestern parts of the UMRB which are more sensitive to freeze-thaw cycle, mainly because SWAT-FT simulated less frequent frozen soils. This study highlights the importance of using physically based freeze-thaw cycle representation in water quality modeling. Design of future nitrogen runoff reduction strategies should include careful assessment of effects that land management has on the freeze-thaw cycles to provide reliable projection of water quality under climate change.

Wang, Qianfeng↗
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 OSTI2022

Characterizing Drought Behavior in the Colorado River Basin Using Unsupervised Machine Learning

Drought is a pressing issue for the Colorado River Basin (CRB) due to the social and economic value of water resources in the region and the significant uncertainty of future drought under climate change. Here, we use climate simulations from various Earth System Models (ESMs) to force the Variable Infiltration Capacity hydrologic model and project multiple drought indicators for the sub-watersheds within the CRB. We apply an unsupervised machine learning (ML) based on Non-Negative Matrix Factorization using K-means clustering (NMFk) to synthesize the simulated historical, future, and change in drought indicators. The unsupervised ML approach can identify sub-watersheds where key changes to drought indicator behavior occur, including shifts in snowpack, snowmelt timing, precipitation, and evapotranspiration. While changes in future precipitation vary across ESMs, the results indicate that the Upper CRB will experience increasing evaporative demand and surface-water scarcity, with some locations experiencing a shift from a radiation-limited to a water-limited evaporation regime in the summer. Large shifts in peak runoff are observed in snowmelt-dominant sub-watersheds, with complete disappearance of the snowmelt signal for some sub-watersheds. The work demonstrates the utility of the NMFk algorithm to efficiently identify behavioral changes of drought indicators across space and time and to quickly analyze and interpret hydro climate model results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Science-integrated Artificial-intelligence for Flooding and precipitation Extremes (SAFE)

A grand challenge in hydrologic science is to understand why signals of climate change and variability, which are often visible in precipitation extremes at aggregate scales, are not consistently observed in the case of extreme flooding. However, a solution to this challenge may prove elusive unless the water cycle is viewed in an integrative manner. Thus, for riverine flooding, while Hortonian (infiltration excess) runoff may have stronger correlation with precipitation extremes and hence perhaps to warming trends or climate oscillators, Dunne (saturation excess) runoff may have a more complex relationships with time series of precipitation and with evaporation and transpiration, but rain-on-snow and snowmelt events may depend on land-surface and atmospheric temperatures. Atmospheric rivers and tropical cyclones lead to precipitation or flooding and are impacted by climate. Flooding assessments need to consider long-term baselines, evolving risk factors, coupled natural-human systems, and novel adaptation such as nature-inspired design.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Water balance model (WBM) v.1.0.0: a scalable gridded global hydrologic model with water-tracking functionality

Abstract. This paper describes the University of New Hampshire Water Balance Model, WBM, a process-based gridded global hydrologic model that simulates the land surface components of the global water cycle and includes water extraction for use in agriculture and domestic sectors. The WBM was first published in 1989; here, we describe the first fully open-source WBM version (v.1.0.0). Earlier descriptions of WBM methods provide the foundation for the most recent model version that is detailed here. We present an overview of the model functionality, utility, and evaluation of simulated global river discharge and irrigation water use. This new version adds a novel suite of water source tracking modules that enable the analysis of flow-path histories on water supply. A key feature of WBM v.1.0.0 is the ability to identify the partitioning of sources for each stock or flux within the model. Three different categories of tracking are available: (1) primary inputs of water to the surface of the terrestrial hydrologic cycle (liquid precipitation, snowmelt, glacier melt, and unsustainable groundwater); (2) water that has been extracted for human use and returned to the terrestrial hydrologic system; and (3) runoff originating from user-defined spatial land units. Such component tracking provides a more fully transparent model in that users can identify the underlying mechanisms generating the simulated behavior. We find that WBM v.1.0.0 simulates global river discharge and irrigation water withdrawals well, even with default parameter settings, and for the first time, we are able to show how the simulation arrives at these fluxes by using the novel tracking functions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Spatial heterogeneity effects on land surface modeling of water and energy partitioning

Abstract. Understanding the influence of land surface heterogeneity on surface water and energy fluxes is crucial for modeling earth system variability and change. This study investigates the effects of four dominant heterogeneity sources on land surface modeling, including atmospheric forcing (ATM), soil properties (SOIL), land use and land cover (LULC), and topography (TOPO). Our analysis focused on their impacts on the partitioning of precipitation (P) into evapotranspiration (ET) and runoff (R), partitioning of net radiation into sensible heat and latent heat, and corresponding water and energy fluxes. An initial set of 16 experiments were performed over the continental US (CONUS) using the E3SM land model (ELMv1) with different combinations of heterogeneous and homogeneous datasets. The Sobol' total and first-order sensitivity indices were utilized to quantify the relative importance of the four heterogeneity sources. Sobol' total sensitivity index measures the total heterogeneity effects induced by a given heterogeneity source, consisting of the contribution from its own heterogeneity (i.e., the first-order index) and its interactions with other heterogeneity sources. ATM and LULC are the most dominant heterogeneity sources in determining spatial variability of water and energy partitioning, mainly contributed by their own heterogeneity and slightly contributed by their interactions with other heterogeneity sources. Their heterogeneity effects are complementary, both spatially and temporally. The overall impacts of SOIL and TOPO are negligible, except TOPO dominates the spatial variability of R/P across the transitional climate zone between the arid western and humid eastern CONUS. Accounting for more heterogeneity sources improves the simulated spatial variability of water and energy fluxes when compared with ERA5-Land reanalysis dataset. An additional set of 13 experiments identified the most critical components within each heterogeneity source, which are precipitation, temperature, and longwave radiation for ATM, soil texture, and soil color for SOIL and maximum fractional saturated area parameter for TOPO.

54 ENVIRONMENTAL SCIENCES↗