A Platform for the Next Generation of Data-driven, River and Watershed-centric Modeling: Veins of the Earth
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
In 1991 a dam-break study was conducted for the high hazard dams located at L Lake and PAR Pond on the Savannah River Site. Two scenarios were considered, over topping from a Probable Maximum Flood (PMF), and a fair weather dam-break for either or both dams. Unfortunately, no inundation map was developed from the study. The purpose of this project was to redo the original dam-break study with improved data and methodology to generate Inundation maps to assist with emergency response and evacuation plans. The Hydrologic Engineering Center's River Analysis System (HEC-RAS) is a free to download river analysis modeling program developed by the US Army Corps of Engineers capable of 1D and 2D hydraulic calculations. Version 5.0.7 (released March 2019) was used for this project. Digital elevation models for the area were retrieved from the US Geological Survey database and converted to a .hdf file within the program. Then the 2D flood area was identified from the contours. Both L Lake and PAR Pond were inputted as 1D storage areas because DEM data does not contain elevation values under water bodies. An elevation vs volume curve was available for both storage areas. Initial elevations were set for both scenarios. Both dams are earthen dams. The Steel Creek dam at L Lake has 6 ft diameter conduit with an upper and lower sluice gate. PAR pond dam consist of a weir and sluice gate connected to an 8 x 8 ft channel. Both outlets were modeled with a pool elevation vs discharge curve. The steel creek dam sluice gates were assumed to be fully open in all cases. As in the previous study, the dams were set to breach when they were overtopped by 1.5 ft during PMF conditions (Figure 1). A fair weather breach was set to be due to a piping failure (Figure 1). In the dual dam break during fair weather conditions the PAR pond dam fails 3 hours after Steel Creek to achieve maximum flooding in the down stream reaches. Simulation was run 3 days for each case and with a 1 minute computational interval. Maximum flooding occurs under PMF conditions with the failure of both dams. PAR pond dam fails first 16 hours and 32 minutes after the start of the simulation with the Steel creek dam failing 6 minutes later. In all cases, the bridges and roads spanning Steel Creek and Lower Three Runs will be inundated and potentially washed away. The Burtons Ferry Highway south of the storage areas will be partially flooded during PMF failure, dual fair weather failure, and PAR pond failure under fair weather conditions.
Simultaneous flooding&heat and droughts&heat events can potentially destabilize hydro-meteorological conditions to deteriorate the water quality of Neches River. Machine learning (ML) models utilizing wireless sensor measurements have been applied to predict water quality and optimize various water management strategies. This study aims to develop ML models to predict dissolved oxygen (DO) prediction under various hydro-meteorological conditions and enhance water management decision-making. Wireless sensor measurements of DO, water temperature, sample depth, conductivity, turbidity, and pH, along with discharge from the United States Geological Survey stations, are collected for model inputs at the Pine Island Bayou C749 station (PIB-C749) and Neches River Saltwater Barrier (SWB). Multilayer perceptron neural networks, recurrent neural networks, long short-term memory (LSTM), and bidirectional LSTM (BiLSTM) with and without attention mechanism (AT) are tested to determine the best model, which is applied the rolling forecast method to predict 14-day DO. Traditional and recurrent transfer learning (TL and RTL) methods are adopted to overcome insufficient data at the SWB. The input feature importance analysis using the integrated gradients (IG) algorithm is applied to determine dominant inputs. The results show LSTM-based models are capable handling long sequential data. AT-BiLSTM and RTL-LSTM demonstrate the best performance at the PIB-C749 (RMSE=0.054) and the SWB (RMSE=0.028), respectively. TL and RTL methods significantly improve model performance at the SWB. DO, temperature, and pH show higher importance, consistent with hydrodynamics and water chemistry. Both best models are applied to predict 14-day DO and demonstrate reasonable performance for decision-making. Hydro-meteorological conditions of 2017 flood and 2012 drought events are simulated and reveal that possible hypoxia occurs after flooding due to increasing temperature and turbidity, and DO concentration decreases significantly under heat and drought conditions. In conclusion, LSTM-based models utilizing wireless sensor data can be a timely and effective approach to make appropriate decisions on water resource management.
Recently, rainfall-runoff simulations in small headwater basins have been improved by methodological advances such as deep neural networks (NNs) and hybrid physics-NN models—particularly, a genre called differentiable modeling that intermingles NNs with physics to learn relationships between variables. However, hydrologic routing simulations, necessary for simulating floods in stem rivers downstream of large heterogeneous basins, had not yet benefited from these advances and it was unclear if the routing process could be improved via coupled NNs. We present a novel differentiable routing method (δMC-Juniata-hydroDL2) that mimics the classical Muskingum-Cunge routing model over a river network but embeds an NN to infer parameterizations for Manning's roughness (n) and channel geometries from raw reach-scale attributes like catchment areas and sinuosity. The NN was trained solely on downstream hydrographs. Synthetic experiments show that while the channel geometry parameter was unidentifiable, n can be identified with moderate precision. With real-world data, the trained differentiable routing model produced more accurate long-term routing results for both the training gage and untrained inner gages for larger subbasins (>2,000 km2) than either a machine learning model assuming homogeneity, or simply using the sum of runoff from subbasins. The n parameterization trained on short periods gave high performance in other periods, despite significant errors in runoff inputs. The learned n pattern was consistent with literature expectations, demonstrating the framework's potential for knowledge discovery, but the absolute values can vary depending on training periods. The trained n parameterization can be coupled with traditional models to improve national-scale hydrologic flood simulations.
The hydrologic flows across the river–aquifer interface play an important role in groundwater dynamics and biogeochemical reactions within the subsurface; however, little is known about the effects of river–aquifer interactions on land surface processes. In this study, we developed a fully coupled three-dimensional (3D) land surface and subsurface model at a high resolution (~1 km) that accounts for high-frequency hydrologic exchange flow conditions to investigate how river–aquifer interactions modulate surface water budgets in the Upper Columbia-Priest Rapids watershed, a typical semiarid watershed located in the northwestern United States where river stage fluctuates in response to reservoir releases changing. Our results show that the spatiotemporal dynamics of river–aquifer interactions are highly heterogeneous, driven mainly by river-stage fluctuations. Adding 6.64 × 10 6 m 3 year –1 of water over the watershed from the river to groundwater owing to the lateral flow, river–aquifer interactions led to an increase in soil evaporation and transpiration supplied by higher soil moisture content, particularly in deeper subsurface. In a hypothetic future scenarios where a 5-m rise in river stage was assumed, the hydrologic flow exchange rates were intensified, resulting in higher surface water over the entire watershed. Overall, lateral flow induced by river–aquifer exchanges leads to an increase in evapotranspiration of ~75% in the historical period and of ~83% in the hypothetical future scenario. Finally, our study demonstrates the potential of coupled model as an effective tool for understanding river–aquifer–land surface interactions, and indicates that river–aquifer interactions fundamentally alter the water balance of the riparian zone for the semiarid watershed and will likely become more frequent and intense in the future under the effects of climate change.
A recently introduced multiscale model for representing the effects of small-scale hyporheic-zone biogeochemical processes is extended from the reach scale to river network scales. The model uses advection-dispersion-reaction equations for the channel network and one-dimensional advection-reaction subgrid models for the hyporheic zone. We summarize implementation in the integrated surface/subsurface hydrology modeling system Amanzi-ATS. The extension uses topologically defined meshes to represent stream/river networks and associates a hyporheic subgrid model with each channel grid cell. Biogeochemical reaction modeling capability residing in community software is accessed through an application programming interface. The implementation is verified against independent numerical solutions on a single reach. Mesh convergence studies show that commonly used semi-distributed representations can introduce significant spatial discretization error. Denitrification of farm runoff in a subbasin of the Portage River Basin in Ohio, USA is used to demonstrate the general-purpose reactive transport capability.
A global atmospheric model is used to calculate the annual river runoff for 33 of the world's major rivers for the present climate and for a doubled CO2 climate. The model has a horizontal resolution of 4 x 5 deg, but the runoff from each model grid box is quartered and added to the appropriate river drainage basin on a 2 x 2.5 deg resolution. The computed runoff depends on the model's precipitation, evapotranspiration, and soil moisture storage. For the doubled CO2 climate, the runoff increased for 25 of the 33 rivers, and in most cases the increases coincide with increased rainfall within the drainage basins. There were runoff increases in all rivers in high northern latitudes, with a maximum increase of 47 percent. At low latitudes there were both increases and decreases ranging from a 96 increase to a 43 percent decrease. The effect of the simplified model assumptions of land-atmosphere interactions on the results is discussed.
This data package is associated with the publication “Spatial microbial respiration variations in the hyporheic zones within the Columbia River Basin” published in the Journal of Geophysical Research: Biogeosciences (Son et al. 2022) available at doi: 10.1029/2021JG006654. This data package includes the key model inputs/outputs of the river corridor model for the Columbia River Basin (CRB) and the model source codes, which were used in the manuscript. The model is a carbon-nitrogen-coupled river corridor model (RCM), and the model is used to quantify hyporheic zone (HZ) aerobic and anaerobic respiration at the NHDPLUS stream reach scales. The RCM used in this study combines empirical substrate models derived from observations and three microbially driven reactions to compute respiration of the HZ for each National Hydrography Dataset (NHD) reach within the CRB. The reactions in HZs of each NHD reach include anaerobic respiration and two-step anaerobic respiration via denitrification. Our HZ respiration estimates are limited to the lotic (or flowing) stream/river systems, and do not account for the respiration process in water column. Note that the RCM only simulates the HZ’s contribution to the dissolved carbon dioxide (CO2) concentrations in the streams, and the CO2 emissions to the atmosphere are not modelled. The model computes at hourly timesteps because of the fast reaction rates. The key input data of the model are exchange flux, residence time, and stream solute (dissolved organic carbon (DOC), dissolved oxygen (DO), and nitrate concentrations). These inputs are constant over time and represent long-term averaged values.This modeling framework successfully quantified HZ respiration components over multiple scales. It revealed key mechanisms driving the spatial variation of HZ aerobic and anaerobic respiration in reaches with varying hydrologic and substrate conditions. Thus, this modeling study offers a testing hypothesis in different river system (e.g., climate and biomes) for the HZ respiration processes, and can be used as a sampling design tool for large-scale HZ experimental studies.This dataset contains five folders: (1) model_inputs, (2) model_outputs, (3) Rscripts, (4) figures, and (5) model_codes. It also contains a readme, file level metadata (FLMD), and data dictionary (dd). Please see the FLMD for a list of all the files contained in this data package and descriptions for each. The model_inputs folder contains the model inputs used to drive the model simulations. The model_outputs folder contains key model output files from the river corridor model. The Rscripts folder contains the Rscripts for pre- and post- processing model results. The figures folder contains the raw figures associated with the manuscript. The model_codes folder includes key model source codes/input files. All files are .jpg, .jpeg, .out, .e, .od, .dat, .sub, .F90, .0, .R, .sbx, .cpg, .sbn, .shx, .shp, .dbf, .prj, .tfw, .tif, .xml, .pdf, or .csv.
In this work, we have used the Soil & Water Assessment Tool (SWAT) to examine streamflow variability of the Lower Mekong River Basin (LMRB) associated with changes in the Upper Mekong River Basin (UMRB) inflows. Two hypothetical experiments were formulated and evaluated for the LMRB, where we conducted runoff simulations with multiple inflow changes that include upstream runoff yield increase and decrease scenarios. Streamflow variability of the LMRB was quantified by two streamflow metrics that explain flow variability and predictability, and high flow disturbance. The model experiments were performed for the Lower Mekong River Basin with identical climate, soil, and other watershed characteristics data. Remote sensing precipitation (Tropical Rainfall Measurement Mission, TRMM, and Global Precipitation Measurement mission, GPM), meteorological data as well as spatial data that include a digital elevation model, newly developed soil information (Harmonized World Soil Database, HWSD), and land use and land cover were processed as input to the LMRB model simulations. Observed daily streamflow data along the Lower Mekong River from Chiang Sean, Thailand to Kratie, Cambodia were used for calibration and validation. Our work results suggest that the Lower Mekong River streamflow is highly variable and has a low predictability (Colwell index of about 32%). We found that releasing more water from upstream Mekong during rainfall months by 30% would result in a reduction in the Lower Mekong streamflow predictability by about 21%. This reduction in predictability is mainly attributed to a decrease in the Contingency index. Our work shows that the ability to predict floods/droughts at the Lower Mekong River would be reduced if there is any anticipated change (i.e., increase/decrease) from UMRB releases. Our results also show that releasing more flows from the upstream Mekong would also affect flood duration and the frequency of flood occurrences downstream. The results of this work thus help to quantify the sensitivity of streamflow variability at the Lower Mekong River Basin to upstream anthropogenic changes.
Air quality dispersion modeling is performed for the Savannah River Site (SRS) to demonstrate compliance with applicable regulations. The AMS/EPA Regulatory Model (AERMOD) modeling system is an EPA recommended model for air quality applications with a data preprocessor (AERMET) to incorporate meteorological data collected on site. AERMET parameterizes or calculates meteorological variables that are not directly measured onsite. One of the parameters estimated by AERMET is the atmospheric mixing height. While the mixing height is not currently a measurement input into AERMET, SRS has the capability to measure the local mixing height. The Savannah River National Laboratory (SRNL) operates a Vaisala CL31 Lidar Ceilometer which estimates mixing height from aerosol backscatter. This study compares the parameterized mixing height from AERMET to the ceilometer estimated mixing height for the current regulatory period at SRS incorporating data from 2015-2019. Results from this study showed the average daily minimum values (morning) from AERMET were an order of magnitude lower than the commonly used Holzworth (1972) method and the ceilometer estimated mixing heights. Additionally, on average, the ceilometer exhibited a daily maximum mixing height value that occurred 1-3 hours later than the AERMET estimated maximum. This difference is likely due to the nighttime atmospheric mixing height assumptions and calculations used by AERMET. The AERMET algorithm cuts off mixing height growth at sunset while the ceilometer data show ongoing evening convection typical of the southeastern United States. These results suggest that the AERMET parametrization scheme assumptions may not be representative of a forested landscape and evening convection which could account for more mixing overnight. The results obtained in this study are significant for air dispersion modeling applications for regulatory purposes and worker safety. Mixing height can impact model estimated pollutant concentrations. A greater mixing height will provide more volume for pollutant dispersion. This report documents efforts to quantify the dependence of mixing height inputs toward a conservative estimated pollutant concentration.
Riverine dissolved organic carbon (DOC), primarily sourced from soil organic carbon (SOC), plays a crucial role in regional and global carbon cycles. However, the complexities of the underlying mechanisms and limited observations present significant challenges for predictive understanding of DOC at regional or larger scales. Recently, we developed a machine learning‐based (ML) map of DOC transformation rates, bridging the gap between SOC and DOC leaching flux and simplifying terrestrial DOC representation. Building on this advancement, we introduce ELM‐MOSART‐DOC, a DOC module integrated into the riverine component of the Energy Exascale Earth System Model (E3SM)—the Model for Scale Adaptive River Transport (MOSART). ELM‐MOSART‐DOC simulates DOC transport and transformation across both headwater streams and river networks, including those managed. Model validation demonstrates the ability of ELM‐MOSART‐DOC to accurately capture long‐term average DOC concentrations, with Kling‐Gupta Efficiency (KGE) scores of 0.58 and 0.76 at large and local stations, respectively. We further assess the impact of reservoirs through different simulation schemes, revealing that reservoirs significantly alter DOC fluxes by regulating streamflow patterns and promoting DOC mineralization. Model simulations indicate that reservoirs reduce total DOC flux from the Mississippi River into the ocean by 7.5%, with the long‐term average annual export decreasing from 3.34 to 3.14 teragrams (Tg) per year. ELM‐MOSART‐DOC integrates process‐based modeling with ML parameterization to enhance the predictive understanding of riverine biogeochemical processes. This approach reduces uncertainties in modeling regional and global carbon cycle ESMs and provides new insights into carbon cycling and its implications for global environmental change.
Streamflow estimation across areas facing water scarcity is important considering the evolving climatic conditions. In this study, the impact of assimilating Soil Moisture Active Passive (SMAP) soil moisture retrievals on the modeled streamflow across three main river basins in South Asia is explored. Model estimated runoff was hydraulically routed using the HYMAP routing model within the NASA Land Information System (LIS) framework. Streamflow estimates via soil moisture assimilation highlighted the potential improvements in streamflow across low-flow river tributaries in irrigated regions. Limited relative change was noted along high-flow river tributaries. Modeled streamflow tended to underestimate the flow magnitude at upstream stations and overestimated the streamflow at downstream stations within the Indus basin due to missing physics related to reservoir operation. Improvement of the modeling systems representativeness of the ground conditions by the inclusion of water management information would potentially improve the current streamflow modeling capability across South Asia.
This package contains all reactive transport simulation outputs used in "Modeling the impact of riparian hollows on river corridor nitrogen exports" (Rogers, D.B., et al., 2021). The simulation outputs reflect hydrological and biogeochemical behavior of a single, generalized riparian hollow located in the intra-meander region of a sub-alpine floodplain under transient hydrological conditions. The context for, purpose of, and conclusions drawn from this model are detailed in Rogers, D.B., et al., 2021.The .zip file contains simulation outputs for each of the five time periods analyzed (see metadata for specific dates), with each folder representing one continuous simulation period. Within each folder are the raw simulation outputs, which are comprised of hundreds of files. Files ending in ".vtk" are best visualized and analyzed using ParaView (https://www.paraview.org). All other files are best viewed or analyzed using R. For a detailed description of the meaning of each file, refer to the ".fls" file located in each folder.
We state that the goal of run-of-river hydropower is to produce cost-competitive renewable electricity with minimal disruption of the natural riverine ecosystem. Modeling and feasibility analysis of alternative design options are crucial for developing new run-of-river hydropower projects. Our review shows that existing run-of-river hydropower design models focus on maximizing economic potential at high-head diversion schemes with limited consideration of environmental outcomes. Since nearly three-quarters of new hydropower potential in the United States is found at low-head sites and environmental performance standards are imperative to project success, new models are needed to address the multi-dimensional design challenges at these sites. To aid in formulating holistic models, we synthesize the performance objectives and design variables related to early-stage run-of-river facility design. The objectives span six potential impact areas, including hydrologic alteration, sediment continuity, water quality, aquatic species passage, social, and economic. Based on these reviews, we identify three key areas to enhance the capabilities of run-of-river hydropower design models. These are 1) expanded model formulations, 2) assessment of barrier effects, and 3) explicit environmental objectives. The resulting modeling improvements would accelerate the identification of run-of-river hydropower designs that minimize environmental impacts, promote economic competitiveness, and incorporate the value of non-power benefits.
Predictive modeling of watershed to river basin scale biogeochemistry through the use of AI techniques and AI-derived model components and the use of AI to design a prediction system comprised of a hierarchy of multi-fidelity models, including AI driven model/component/parameterization selection.
As the world's population continues to increase, additional stress is placed on water resources. This stress, coupled with future uncertainties regarding climate change, makes arid and semi-arid regions particularly vulnerable. One example is the Aral Sea where the freshwater inflow, which is dominated by snowmelt runoff, has decreased significantly since the expansion of intensive irrigation in the 1960s. The purpose of this paper is to use a river routing scheme from a global climate model to examine the flow of the Amu Dar'ya River into the Aral Sea. The river routing scheme is modified to include groundwater flow, flooding, and evaporative losses in the river's wetlands and floodplain, and anthropogenic withdrawals for irrigation. A set of scenarios is designed to test the sensitivity of river flow to the inclusion of these modifications into the river routing scheme. When riverine wetlands and floodplains are present, the river flow is reduced significantly and is similar to the observed flow. In addition the model results show that it is essential to incorporate human diversions to accurately represent the inflow to the Aral Sea, and they also indicate potential management strategies that might be appropriate to maintain a balance between inflow to the Sea and upstream diversions for irrigation.
This data package is associated with the publication “Combined effects of stream hydrology and land use on basin‐scale hyporheic zone denitrification in the Columbia River Basin”, published in Water Resource Research (Son et al.2022) available at https://doi.org/10.1029/2021WR031131. This data package includes the key model inputs/outputs of the river corridor model for the Columbia River Basin (CRB) and the model source codes used in the manuscript. The model is a carbon-nitrogen-coupled river corridor model (RCM), and the model is used to quantify hyporheic zone (HZ) denitrification at the NHDPLUS stream reach scales. The RCM used in this study combines empirical substrate models derived from observations and three microbially driven reactions, including two-step denitrification and aerobic respiration, are considered within the HZ. The key input data of the model are exchange flux, residence time, and stream solute (dissolved organic carbon (DOC), dissolved oxygen (DO), and nitrate concentrations). These inputs are constant over time and represent long-term averaged values. This study uses the RCM to explore the spatial patterns of HZ denitrification across reaches with different sizes and land use in the CRB. Our main objective is to use the RCM as a virtual reality model, and the machine-learning models as surrogates that encapsulate the complexities of the physics-based model while identifying the importance of different variables that are not evident in the model conceptualization. We do not include a direct comparison of the modeled HZ denitrification and measurements; however, the RCM can capture the overall spatial patterns of the HZ denitrification because the model inputs and its reaction networks are based on well-established theory and a physical-based model. The combination of the model-based predictions and a machine-learning approach (e.g., random forest) is used to improve our understanding of what variables of the model are associated with spatial patterns of the modeled denitrification across reaches with different sizes and land uses, and to develop a proxy model using measurable variables to reproduce the simulated patterns.This dataset contains five folders: (1) model_inputs, (2) model_outputs, (3) Rscripts, (4) figures, and (5) model_codes. It also contains a readme, file level metadata (FLMD), and data dictionary (dd). Please see the FLMD for a list of all the files contained in this data package and descriptions for each. The model_inputs folder contains the model inputs used to drive the model simulations. The model_outputs folder contains key model output files from the river corridor model. The Rscripts folder contains the Rscripts for pre- and post- processing model results. The figures folder contains the raw figures associated with the manuscript. The model_codes folder includes key model source codes/input files. All files are .jpg, .jpeg, .out, .e, .od, .dat, .sub, .F90, .0, .R, .sbx, .cpg, .sbn, .shx, .shp, .dbf, .prj, .tfw, .tif, .xml, .pdf, or .csv.