Engineering Papers⌕ Search

Engineering topics

Painter, Scott

Publications and source records attributed to Painter, Scott.

Reviews and syntheses: Variable inundation across Earth's terrestrial ecosystems

Abstract. The structure, function, and dynamics of Earth's terrestrial ecosystems are profoundly influenced by how often (frequency) and how long (duration) they are inundated with water. A diverse array of natural and human-engineered systems experience temporally variable inundation whereby they fluctuate between inundated and non-inundated states. Variable inundation spans extreme events to predictable sub-daily cycles. Variably inundated ecosystems (VIEs) include hillslopes, non-perennial streams, wetlands, floodplains, temporary ponds, tidal systems, storm-impacted coastal zones, and human-engineered systems. VIEs are diverse in terms of inundation regimes, water chemistry and flow velocity, soil and sediment properties, vegetation, and many other properties. The spatial and temporal scales of variable inundation are vast, ranging from sub-meter to whole landscapes and from sub-hourly to multi-decadal. The broad range of system types and scales makes it challenging to predict the hydrology, biogeochemistry, ecology, and physical evolution of VIEs. Despite all experiencing the loss and gain of an overlying water column, VIEs are rarely considered together in conceptual, theoretical, modeling, or measurement frameworks and approaches. Studying VIEs together has the potential to generate mechanistic understanding that is transferable across a much broader range of environmental conditions, relative to knowledge generated by studying any one VIE type. We postulate that enhanced transferability will be important for predicting changes in VIE function in response to global change. Here we aim to catalyze cross-VIE science that studies drivers and impacts of variable inundation across Earth's VIEs. To this end, we complement expert mini-reviews of eight major VIE systems with overviews of VIE-relevant methods and challenges associated with scale. We conclude with perspectives on how cross-VIE science can derive transferable understanding via unifying conceptual models in which the impacts of variable inundation are studied across multi-dimensional environmental space.

54 ENVIRONMENTAL SCIENCES↗

Modeling the Effects of Artificial Drainage on Agriculture-dominated Watersheds using a Fully Distributed Integrated Hydrology Model: Datasets, code, models files

This model-data archive is for a modeling study focused on the effects of artificial drainage on agriculture-dominated watersheds using a fully distributed integrated hydrology model. This study introduced a novel strategy to represent the effects of tile drains and agricultural ditches in a fully distributed basin-scale integrated hydrology model. The numerical experiments in this study reveal the effects of surface and subsurface drainage on various hydrological states and fluxes. The details about the datasets can be found in the "readme" document attached with the dataset.

Agricultural Watersheds↗

Model data for numerical evaluation of photosensitive tracers as a strategy for separating surface and subsurface transient storage in streams

This model-data archive pertains to a study aimed at the numerical evaluation of the photosensitive tracers as a potential strategy for separating the effects of surface and hyporheic storage zones (SSZs and HSZs, respectively). Separating effects of SSZ and HSZ are important for accurately representing stream function as HSZs and SSZs expose solutes to significantly different biogeochemical conditions. We perform numerical experiments using a multiscale reactive transport model for stream corridors implemented in ATS code, which allows for representing multiple storage zones with their respective travel time distributions and biogeochemistry. For each of the numerical experiment, we provide in this model-data archive python wrapper to drive ATS (forward_model.py), synthetic observation (synthetic_btc.py, BTC_observed.csv), forward model files including ATS input (multiscale_transport.tpl), PFLOTRAN inputs (reactions_channel.tpl, reactions_hz.tpl, reactions_sz.tpl), MCMC (mcmc_run.py), predictive uncertainty (pred_uncert.py, BTCs_simulated_day.csv, BTCs_simulated_night.csv), MCMC outputs (tracer_test_GR.npy, tracer_test_logps.npy, tracer_test_parameters.npy) and Jupyter notebook for post-processing and visualization (post-processing.ipynb). For the denitrification application, ATS input file (denitrification_multisubgrid.tpl), PFLOTRAN input files (denitrification_channel.in, denitrification_hz.in), predictive uncertainty (pred_uncert_denitrification.py, BTCs_simulated_DO.csv, BTCs_simulated_DOC.csv, BTCs_simulated_Nitrate.csv).

54 ENVIRONMENTAL SCIENCES↗

Advanced Terrestrial Simulator (ATS) evaluation dataset at 7 catchments across the continental United States

This dataset comprises of the input files and other files required for Advanced Terrestrial Simulator (ATS) simulations at 7 catchments across the continental United States. ATS is an integrated surface-subsurface hydrology model. We include Jupyter notebooks (within scripts folder) for individual catchments showing information (including data sources, river network, soil, geology, landuse types etc.) on preparing the machine readable input files. ATS observation output files are provided in the output folder. Figures and analyses (.xlsx sheets) are also provided. The catchments include, Taylor River Upstream (Colorado); (b) Cossatot River (Arkansas); (c) Panther Creek (Alabama); (d) Little Tennessee River (North Carolina and Georgia); (e) Mayo River (Virginia); (f) Flat Brook (New Jersey); (g) Neversink River headwaters (New York). Readme files are provided inside the directories providing more details. Files types include: .xml, .h5, .xlsx, .png, .ipynb, .py, .nc, .txt. All of the files types can be accessed by open source software, details on software requirements are following: .xml (any text editors including notepad and textedit), .h5 (in python using hdf libraries), .xlsx (WPS Office Spreadsheets, OpenOffice Calc, LibreOffice Calc, Microsoft Office etc.), .png (any image viewer), .ipynb (Jupyter notebook), .py (any text editors including notepad and textedit), .nc (using python or other open source software).

54 ENVIRONMENTAL SCIENCES↗

Accurate and Timely Forecasts of Geologic Carbon Storage using Machine Learning Methods

Carbon capture and storage is one strategy to reduce greenhouse gas emissions. One approach to storing the captured CO2 is to inject it into deep saline aquifers. However, dynamics of the injected CO2 plume is uncertain and the potential for leakage back to the atmosphere must be assessed. Thus, accurate and timely forecasts of CO2 storage via real-time measurements integration becomes very crucial. This study proposes a learning-based, inverse-free prediction method that can accurately and rapidly forecast CO2 movement and distribution with uncertainty quantification based on limited simulation and observation data. The machine learning techniques include dimension reduction, multivariate data analysis, and Bayesian learning. The outcome is expected to provide CO2 storage site operators with an effective tool for real-time decision making.

Lu, Dan↗

Toward more mechanistic representations of biogeochemical processes in river networks: Implementation and demonstration of a multiscale model

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.

54 ENVIRONMENTAL SCIENCES↗

Toward more mechanistic representations of biogeochemical processes in stream and river networks: Implementation and demonstration of a multiscale model : Modeling Archive

This Modeling Archive supports the ORNL-SFA and IDEAS-Waterhsed publication submitted to Environmental Modeling and Software. A recently introduced multiscale model for representing the combined effect of hyporheic exchange flows and small-scale hyporheic-zone biogeochemical processes is extended from the reach scale to river network scales. The ADELS (Advection Dispersion Equation with Lagrangian Subgrid) model uses a one-dimensional advection-dispersion-reaction equation for the channel network and couples that equation at each channel location to a one-dimensional advection-reaction subgrid model representing an ensemble of streamlines that are diverted into the hyporheic zone before returning to the channel. The subgrid model is written in a computationally advantageous Lagrangian form with hyporheic age replacing the hyporheic travel distance. In the paper, we summarized implementation in the integrated surface/subsurface hydrology modeling system Advanced Terrestrial Simulator (ATS).

East Fork Poplar Creek↗

Streamflow simulation in data-scarce basins using Bayesian and physics-informed machine learning models

Hydrologic predictions at rural watersheds are important but also challenging due to data shortage. Long short-term memory (LSTM) networks are a promising machine learning approach and have demonstrated good performance in streamflow predictions. However, due to its data-hungry nature, most LSTM applications focus on well-monitored catchments with abundant and high-quality observations. In this work, we investigate predictive capabilities of LSTM in poorly monitored watersheds with short observation records. To address three main challenges of LSTM applications in data-scarce locations, i.e., overfitting, uncertainty quantification (UQ), and out-of-distribution prediction, we evaluate different regularization techniques to prevent overfitting, apply a Bayesian LSTM for UQ, and introduce a physics-informed hybrid LSTM to enhance out-of-distribution prediction. Through case studies in two diverse sets of catchments with and without snow influence, we demonstrate that 1) when hydrologic variability in the prediction period is similar to the calibration period, LSTM models can reasonably predict daily streamflow with Nash–Sutcliffe efficiency above 0.8, even with only 2 years of calibration data; 2) when the hydrologic variability in the prediction and calibration periods is dramatically different, LSTM alone does not predict well, but the hybrid model can improve the out-of-distribution prediction with acceptable generalization accuracy; 3) L2 norm penalty and dropout can mitigate overfitting, and Bayesian and hybrid LSTM have no overfitting; and 4) Bayesian LSTM provides useful uncertainty information to improve prediction understanding and credibility. In conclusion, these insights have vital implications for streamflow simulation in watersheds where data quality and availability are a critical issue.

54 ENVIRONMENTAL SCIENCES↗

ORNL Input to GDSA Repository Systems Analysis FY21

This document satisfies the M3 milestone M3SF-21OR010304082 titled “ORNL Input to GDSA Repository Systems.” This document describes the current status of Oak Ridge National Laboratory’s (ORNL’s) efforts related to analyses of dual-purpose canister (DPC) disposal in unsaturated alluvium with a focus on thermal hydrological constraints on criticality timing and power output. This analysis updates previous analyses of thermal hydrological constraints on timing and power output of a potential criticality event in Dual-Purpose Canisters (DPCs) in a hypothetical repository in unsaturated alluvium using a more realistic representation of heat transport inside dry DPCs. PFLOTRAN was used for the simulations of multiphase thermal hydrology near a single DPC. The scenario considers a DPC failure at 9000 years, allowing water to enter and eventually fill the DPC. Once the DPC is filled to a level that could support a criticality event, different values for criticality power output were added to the decay heat. The main objective is to bound the power output that could be produced by a criticality event without driving water out of the package. For the conditions analyzed here, following a package breach, the alluvial formation could supply enough water to allow enough accumulation in the DPC to support a criticality event. However, the power output that would be generated is limited to modest values by loss of water moderator due to evaporation and vapor diffusion. In the reference case scenario, the DPC would not start to fill with water until about 16,000 years post closure and would not fill to a level that allows a criticality event until approximately 25,000 years post closure. The long-term average power output that could be sustained without driving off the water and terminating the criticality event is limited to approximately 100 W. Sensitivity to assumed conditions and parameters in the reference case, especially the deep percolation rate, need to be addressed but could not be undertaken because of numerical failures of the PFLOTRAN code running in general mode in the dry conditions of an unsaturated alluvium repository

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

On the Representation of Hyporheic Exchange in Models for Reactive Transport in Stream and River Corridors

Efforts to include more detailed representations of biogeochemical processes in basin-scale water quality simulation tools face the challenge of how to tractably represent mass exchange between the flowing channels of streams and rivers and biogeochemical hotspots in the hyporheic zones. Multiscale models that use relatively coarse representations of the channel network with subgrid models for mass exchange and reactions in the hyporheic zone have started to emerge to address that challenge. Two such multiscale models are considered here, one based on a stochastic Lagrangian travel time representation of advective pumping and one on multirate diffusive exchange. The two models are formally equivalent to well-established integrodifferential representations for transport of non-reacting tracers in steady stream flow, which have been very successful in reproducing stream tracer tests. Despite that equivalence, the two models are based on very different model structures and produce significantly different results in reactive transport. In a simple denitrification example, denitrification is two to three times greater for the advection-based model because the multirate diffusive model has direct connections between the stream channel and transient storage zones and an assumption of mixing in the transient storage zones that prevent oxygen levels from dropping to the point where denitrification can progress uninhibited. By contrast, the advection-based model produces distinct redox zonation, allowing for denitrification to proceed uninhibited on part of the hyporheic flowpaths. These results demonstrate that conservative tracer tests alone are inadequate for constraining representation of mass transfer in models for reactive transport in streams and rivers.

54 ENVIRONMENTAL SCIENCES↗

Recovering the effects of subgrid heterogeneity in simulations of radionuclide transport through fractured media

Groundwater flow and contaminant transport through fractured media can be simulated using Discrete Fracture Network (DFN) models which provide a natural description of structural heterogeneity. However, this approach is computationally expensive, with the large number of intersecting fractures necessitated by many real-world applications requiring modeling simplifications to be made for calculations to be tractable. Upscaling methods commonly used for this purpose can result in some loss of local-scale variability in the groundwater flow velocity field, resulting in underestimation of particle travel times, transport resistance and retention in transport calculations. In this paper, a transport downscaling algorithm to recover the transport effects of heterogeneity is tested on a synthetic Brittle Fault Zone model, motivated by the problem of large safety assessment calculations for geological repositories of spent nuclear fuel. We show that the variability in the local-scale velocity field which is lost by upscaling can be recovered by sampling from a library of DFN transport paths, accurately reproducing DFN transport statistic distributions and radionuclide breakthrough curves in an upscaled model.

54 ENVIRONMENTAL SCIENCES↗