Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Water resources”

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 1,477 records · Page 82

Potential Chemical Impacts of Subsurface CO 2 : An Integrated Experimental and Numerical Assessment for a Case Study of the Ogallala Aquifer

Leakage from geologic CO 2 sequestration (GCS) reservoirs to overlying underground sources of drinking water (USDW) is a tangible risk. This study is an integrated assessment that combines column experiments and reactive transport simulations of sediments sampled from the Ogallala aquifer above an active commercial-scale GCS site (the Farnsworth Unit in northern Texas). Experimental and simulation results suggest that carbonate mineral (calcite and dolomite) dissolution is the most significant reaction following CO 2 intrusion, and is also the dominant source of trace metal release. Cation exchange is another key mechanism controlling trace metal release by cation interference. Most of the trace metals, including Ba, Sr, As, Pb, and Zn, show a short-term release and quickly drop to the baseline values, suggesting low risk to the overlying USDW quality. Other trace metals, such as Mn and U, exhibit a tangible increase of their concentrations in the beginning, and drop to a higher level compared to the baseline, which may become a potential concern of long-term USDW quality change with CO2introduction. This study provides a comprehensive example of combining laboratory experiments and simulations for assessment of CO 2 -sediment interactions with combined release mechanisms in shallow groundwater aquifers. Data presented here provides useful insights for quantitative risk assessment and effective public education regarding CO 2 geological sequestration, and trace metal reactive transport studies in shallow groundwater aquifers.

54 ENVIRONMENTAL SCIENCES↗

Hyporheic Oxygen Dynamics in the East River, Colorado: Insights From an In–Situ, High Frequency Time Series During Two Distinct Flow Seasons

Dissolved oxygen (DO) is critical for aquatic ecosystems, however, few studies have focused on the long-term DO dynamics in hyporheic zones, which are a function of both transport (hydrologic exchange between river and hyporheic zone) and uptake by biogeochemical reactions or respiration. We explore the dynamics of temperature and DO at 10, 20, and 35 cm depth beneath the East River, Colorado, from July–October 2017 (relatively normal water year) and April to October 2018 (comparatively low flow year), enabled by distinctive, in-situ, high frequency (Δt = 5 min) sensors that provided continuous time-series from the undisturbed study site over 14 months. We expect that hyporheic DO, which has a regular daily fluctuation pattern, is supplied by the surface water (at all times we estimate downwelling) and that diurnal hyporheic DO temporal patterns should be aligned with diurnal hyporheic temperature patterns. However, this was not found to be the case. Hyporheic DO becomes depleted briefly at 20 and 35 cm depths in 2017, and at all three hyporheic depths for extended periods in 2018. Whereas diurnal temperature fluctuations have consistent timings of maxima and minima, hyporheic DO rarely has as regular a pattern, and daily ranges are inconsistent. Rainfall events caused some of these changes to diurnal hyporheic DO patterns without repeatable effects. Antecedent snowpack conditions influence streamflow dynamics and therefore hyporheic DO dynamics in this alpine river. These results also point to the strong and variable influence of hyporheic microbial communities regulating hyporheic DO.

54 ENVIRONMENTAL SCIENCES↗

Method of Distributions for Two‐Phase Flow in Heterogeneous Porous Media

Abstract Multiscale heterogeneity and insufficient characterization data for a specific subsurface formation of interest render predictions of multi‐phase fluid flow in geologic formations highly uncertain. Quantification of the uncertainty propagation from the geomodel to the fluid‐flow response is typically done within a probabilistic framework. This task is computationally demanding due to, for example, the slow convergence of Monte Carlo simulations (MCS), especially when computing the tails of a distribution that are necessary for risk assessment and decision‐making under uncertainty. The frozen streamlines method (FROST) accelerates probabilistic predictions of immiscible two‐phase fluid flow problems; however, FROST relies on MCS to compute the travel‐time distribution, which is then used to perform the transport (phase saturation) computations. To alleviate this computational bottleneck, we replace MCS with a deterministic equation for the cumulative distribution function (CDF) of travel time. The resulting CDF‐FROST approach yields the CDF of the saturation field without resorting to sampling‐based strategies. Our numerical experiments demonstrate the high accuracy of CDF‐FROST in computing the CDFs of both saturation and travel time. For the same accuracy, it is about 5 and 10 times faster than FROST and MCS, respectively.

Yang, Hyung Jun↗

Testing Hidden Assumptions of Representativeness in Reach-Scale Studies of Hyporheic Exchange

Field studies of hyporheic exchange in mountain systems are often conducted using short study reaches and a limited number of observations. It is common practice to assume these study reaches represent hyporheic exchange at larger scales or different sites and to infer general relationships among potential causal mechanisms from the limited number of observations. However, these assumptions of representativeness are rarely tested. In this study, we develop numerical models from four segments of mountain streams in different geomorphologic settings and extract shorter reaches to test how representative exchange metrics are in shorter reaches compared to their reference segments. We also map the locations of the representative reaches to determine if a pattern exists based on location. Finally, we compare variance of these shorter within-site reaches to 29 additional reaches across the same basin to understand the impacts of inferring causal mechanisms, for example, the expectation that wide and narrow valley bottoms will yield different hyporheic exchange patterns. Our results show that the location and length strategy of the study reach must be considered before assuming an exchange metric to be representative of anything other than the exact segment studied. Further, it is necessary to quantify within and between site variations before making causal inferences based on observable characteristics, such as valley width or stream morphology. Our findings have implications for future field practices and how those practices are translated into models.

54 ENVIRONMENTAL SCIENCES↗

Combining Remotely Sensed Evapotranspiration and an Agroecosystem Model to Estimate Center‐Pivot Irrigation Water Use at High Spatio‐Temporal Resolution

Abstract Estimating irrigation water use accurately is critical for sustainable irrigation and studying terrestrial water cycle in irrigated croplands. However, irrigation is not monitored in most places, and current estimations of irrigation water use has coarse spatial and/or temporal resolutions. This study aims to estimate irrigation water use at the daily and field scale through the proposed model‐data fusion framework, which is achieved by particle filtering with two configurations (concurrent, CON, and sequential, SEQ) by assimilating satellite‐based evapotranspiration (ET) observations into an advanced agroecosystem model, ecosys . Two types of experiments using synthetic and real ET observations were conducted to study the efficacy of the proposed framework for estimating irrigation water use at the irrigated fields in eastern and western Nebraska, United States. The experiments using synthetic ET observations indicated that, for two major sources of uncertainties of ET difference between observations and model simulations, which are bias and noise, noise had larger impacts on degrading the estimation performance of irrigation water use than bias. For the experiments using real ET observations, monthly and annual estimations of irrigation water use matched well with farmer irrigation records, with Pearson correlation coefficient ( r ) around 0.80 and 0.50, respectively. Although detecting daily irrigation records was very challenging, our method still gave a good performance with RMSE, BIAS, and r around 2.90, 0.03, and 0.4 mm/d, respectively. Our proposed model‐data fusion framework for estimating irrigation water use at high spatio‐temporal resolution could contribute to regional water management, sustainable irrigation, and better tracking terrestrial water cycle.

Environmental Sciences & Ecology↗

Automated Symbolic Upscaling: 1. Model Generation for Extended Applicability Regimes

Abstract In porous media theory, upscaling techniques are fundamental to deriving rigorous Darcy‐scale models for flow and reactive transport in subsurface systems. Due to limitations in classical techniques, a number of ad hoc approaches have been proposed to address physical regimes in which reactive time scales are similar to, or faster than, diffusive time scales. In Part 1 of this two part series, we present a strategy for expanding the applicability of classical homogenization theory by generalizing the assumed closure form. We detail the implementation of this strategy on two reactive mass transport problems with moderately reactive physics. The strategy produces nontrivial homogenized models with emergent terms and effective parameters that couple reactive, diffusive, and advective transport. The differences in equation forms between the macroscopic and pore‐scale descriptions advise caution to further studies where the forms of macroscopic equations are assumed, as opposed to rigorously derived. Numerical validation is provided for each problem to show that the error estimates of homogenization theory are satisfied, and to justify the implemented strategy. In Part 2, the presented strategy is automated using symbolic computing to expedite its implementation.

Pietrzyk, Kyle↗

Evaporation and Transpiration From Multiple Proximal Forests and Wetlands

Abstract Climate change is intensifying the hydrologic cycle and altering ecosystem function, including water flux to the atmosphere through evapotranspiration (ET). ET is made up of evaporation (E) via non‐stomatal surfaces, and transpiration (T) through plant stomata which are impacted by global changes in different ways. E and T are difficult to measure independently at the ecosystem scale, especially across multiple sites that represent different land use and land management strategies. To address this gap in understanding, we applied flux variance similarity (FVS) to quantify how E and T differ across 13 different ecosystems measured using eddy covariance in a 10 × 10 km area from the CHEESEHEAD19 experiment in northern Wisconsin, USA. The study sites included eight forests with a large deciduous broadleaf component, three evergreen needleleaf forests, and two wetlands. Average T/ET for the study period averaged nearly 52% in forested sites and 45% in wetlands, with larger values after excluding periods following rain events when evaporation from canopy interception may be expected. A dominance analysis revealed that environmental variables explained on average 69% of the variance of half‐hourly T, which decreased from summer to autumn. Deciduous and evergreen forests showed similar E trajectories over time despite differences in vegetation phenology, and vapor pressure deficit explained some 13% of the variance E in wetlands but only 5% or less in forests. Retrieval of E and T within a dense network of flux towers lends confidence that FVS is a promising approach for comparing ecosystem hydrology across multiple sites to improve our process‐based understanding of ecosystem water fluxes.

Environmental Sciences & Ecology↗

Ensemble Estimation of Historical Evapotranspiration for the Conterminous U.S.

Abstract Evapotranspiration (ET) is the largest component of the water budget, accounting for the majority of the water available from precipitation. ET is challenging to quantify because of the uncertainties associated with the many ET equations currently in use, and because observations of ET are uncertain and sparse. In this study, we combine information provided by available ET data and equations to produce a new monthly data set for ET for the conterminous U.S. (CONUS). These maps are produced from 1895 to 2018 at an 800 m spatial scale, marking a finer resolution than currently available products over this time period. In our approach, the relative performance of a suite of ET equations is assessed using water balance, flux tower, and remotely sensed ET estimates. At the observation locations, we use error distributions to quantify relative weights for the equations and use these in a modified Bayesian model averaging weighted ensemble approach. The relative weights are spatially generalized using a random forest regression, which is applied to wall‐to‐wall explanatory variable maps to generate CONUS‐wide relative weight maps and ensemble estimates. We assess the performance of the ensemble using a reserved subset of the observations and compare this performance against other national‐scale map products for historical to modern ET. The ensemble ET maps are shown to provide an improved accuracy over the alternative comparison products. These ET maps could be useful for a variety of hydrologic modeling and assessment applications that benefit from a long record, such as the study of periods of water scarcity through time.

Environmental Sciences & Ecology↗

Hydrodynamic Modeling of Stratification and Mixing in a Shallow, Tropical Floodplain Lake

Abstract Floodplain lakes are widespread and ecologically important throughout tropical river systems, however data are rare that describe how temporal variations in hydrological, meteorological and optical conditions moderate stratification and mixing in these shallow lakes. Using time series measurements of meteorology and water‐column temperatures from 17 several day campaigns spanning two hydrological years in a representative Amazon floodplain lake, we calculated surface energy fluxes and thermal stratification, and applied and evaluated a 3‐dimensional hydrodynamic model. The model successfully simulated diel cycles in thermal structure characterized by buoyancy frequency, depth of the actively mixing layer, and other terms associated with the surface energy budget. Diurnal heating with strong stratification and nocturnal mixing were common; despite considerable heat loss at night, the strong stratification during the day meant that mixing only infrequently extended to the bottom at night. Simulations indicated that the diurnal thermocline up and downwelled creating lake‐wide differences in near‐surface temperatures and mixing depths. Infrequent full mixing creates conditions conducive to anoxia in these shallow lakes given their warm temperatures.

Environmental Sciences & Ecology↗

Evaluating the Effects of Precipitation and Evapotranspiration on Soil Moisture Variability Within CMIP5 Using SMAP and ERA5 Data

=The effects of precipitation (Pr) and evapotranspiration (ET) on surface soil moisture (SSM) play an essential role in the land-atmosphere system. Here we evaluate multimodel differences of these effects within the Coupled Model Intercomparison Project Phase 5 (CMIP5) compared to Soil Moisture Active Passive (SMAP) products and ECMWF Reanalysis v5 (ERA5) as references in a frequency domain. The variability of SSM, Pr, and ET within three frequency bands (1/7 ~ 1/30 days -1 , 1/30 ~ 1/90 days -1 , and 1/90 ~ 1/365 days -1 ) after normalization is quantified using Fourier transform. We analyze the impact of ET and Pr on SSM variability based on a transfer function assuming that these variables form a linear time-invariant (LTI) system. For the total effects of ET and Pr on SSM variability, the CMIP5 estimations are smaller than the reference data in the two higher frequency bands and are larger than the reference data in the lowest frequency band. Besides, the effects on SSM by Pr and ET are found to be different across the three frequency bands. In each frequency band, the variability of the factor that dominates SSM (i.e., Pr or ET) from CMIP5 is smaller than that from the references. This study identifies the spatiotemporal distribution of differences between CMIP5 models and references (SMAP and ERA5) in simulating ET and Pr effects on SSM within three frequency bands. This study provides insightful information on how soil moisture variability is affected by varying precipitation and evapotranspiration at different time scales within Earth System Models.

54 ENVIRONMENTAL SCIENCES↗

Toward Improved Regional Hydrological Model Performance Using State-Of-The-Science Data-Informed Soil Parameters

Accurate soil moisture and streamflow data are an aspirational need of many hydrologically relevant fields. Model simulated soil moisture and streamflow hold promise but models require validation prior to application. Calibration methods are commonly used to improve model fidelity but misrepresentation of the true dynamics remains a challenge. In this study, we leverage soil parameter estimates from the Soil Survey Geographic (SSURGO) database and the probability mapping of SSURGO (POLARIS) to improve the representation of hydrologic processes in the Weather Research and Forecasting Hydrological modeling system (WRF-Hydro) over a central California domain. Our results show WRF-Hydro soil moisture exhibits increased correlation coefficients ( r ), reduced biases, and increased Kling-Gupta Efficiencies (KGEs) across seven in situ soil moisture observing stations after updating the model's soil parameters according to POLARIS. Compared to four well-established soil moisture data sets including Soil Moisture Active Passive data and three Phase 2 North American Land Data Assimilation System land surface models, our POLARIS-adjusted WRF-Hydro simulations produce the highest mean KGE (0.69) across the seven stations. More importantly, WRF-Hydro streamflow fidelity also increases, especially in the case where the model domain is set up with SSURGO-informed total soil thickness. The magnitude and timing of peak flow events are better captured, r increases across nine United States Geological Survey stream gages, and the mean KGE across seven of the nine gages increases from 0.12 to 0.66. Our pre-calibration parameter estimate approach, which is transferable to other spatially distributed hydrological models, can substantially improve a model's performance, helping reduce calibration efforts and computational costs.

54 ENVIRONMENTAL SCIENCES↗

Automated Symbolic Upscaling: 2. Model Generation for Extended Applicability Regimes

Abstract In this second part of the two paper series, we detail an algorithmic procedure for systematically implementing the generalized closure form strategy presented in Part 1. This strategy extends the applicability of homogenized models with respect to classical homogenization theory, as demonstrated in Part 1 where upscaled models are rigorously derived in moderately reactive physical regimes. After encoding the algorithm into Symbolica, an automated upscaling framework, we upscale two reactive mass transport problems and numerically validate the resulting nonlinear homogenized models by showing the absolute error estimates predicted by homogenization theory are satisfied. In both problems, nontrivial closure forms and closure problems are automatically formulated using the encoded strategy with no human interaction, nor prior knowledge regarding the closure required for the systems. We hope these demonstrations spark further interest in automated analytical frameworks for multiscale modeling, as such capabilities are invaluable for generating rigorous multiscale models of complex phenomena in porous media.

Pietrzyk, Kyle↗

Understanding the Shift of Drivers of Soil Erosion and Sedimentation Based on Regional Process-Based Modeling in the Mississippi River Basin During the Past Century

Soil erosion and sedimentation problems remain a major water quality concern for making watershed management policies in the Mississippi River Basin (MRB). It is unclear whether the observed decreasing trend of stream suspended sediment loading to the mouth of the MRB over the last eight decades truly reflects a decline in upland soil erosion in this large basin. Here, for this work, we improved a distributed regional land surface model, the Dynamic Land Ecosystem Model, to evaluate how climate and land use changes have impacted soil erosion and sediment yield over the entire MRB during the past century. Model results indicate that total sediment yield significantly increased during 1980–2018, despite no significant increase in annual precipitation and runoff. The increased soil erosion and sediment yield are mainly driven by intensified extreme precipitation (EP). Spatially, we found notable intensified EP events in the cropland-dominated Midwest region, resulting in a substantial increase in soil erosion and sediment yield. Land use change played a critical role in determining sediment yield from the 1910s to the 1930s, thereafter, climate variability increasingly became the dominant driver of soil erosion, which peaked in the 2010s. This study highlights the increasing influences of extreme climate in affecting soil erosion and sedimentation, thus, water quality. Therefore, existing forest and cropland Best Management Practices should be revisited to confront the impacts of climate change on water quality in the MRB.

54 ENVIRONMENTAL SCIENCES↗

Automatic Time Step Control to Resolve Hydromechanically Driven Fault Reactivation, Spontaneous Nucleation, and Seismic Arrest

Abstract A physical understanding of the progression from flow‐driven (quasi‐static) poromechanical deformation to dynamic fault rupture is critical to the resilient operations of several engineering systems. These processes are bridged by a progression from fault reactivation to the spontaneous nucleation of unstable sliding. Toward addressing this challenge, novel automatic time step size control methods are developed to enable accurate and efficient simulation of these dynamics and transitions from the first principles. The controllers combine local models for discretization error and Coulomb failure conditions to automatically adjust the time step size across several orders of magnitude. The methods do not require additional empirical or theoretical input and can resolve the pre‐rupture, interseismic, and seismic periods to the allowed accuracy. The computational results reveal that the proposed methods automatically capture the onset of reactivation and nucleation for homogeneous and heterogeneous fields. Hydrodynamic and structural heterogeneity lead to disparate critical nucleation sizes compared to those predicted by theory. The results highlight its potential in predicting induced seismicity in realistic subsurface engineering systems and at practical scales.

Environmental Sciences & Ecology↗

Stress Testing California's Hydroclimatic Whiplash: Potential Challenges, Trade‐Offs and Adaptations in Water Management and Hydropower Generation

Abstract Inter‐annual precipitation in California is highly variable, and future projections indicate an increase in the intensity and frequency of hydroclimatic “whiplash.” Understanding the implications of these shocks on California's water system and its degree of resiliency is critical from a planning perspective. Therefore, we quantify the resilience of reservoir services provided by water and hydropower systems in four basins in the western Sierra Nevada. Using downscaled runoff from 10 climate model outputs, we generated 200 synthetic hydrologic whiplash sequences of alternating dry and wet years to represent a wide range of extremes and transitional conditions used as inputs to a water system simulation model. Sequences were derived from upper (wet) and lower (dry) quintiles of future streamflow projections (2030–2060). Results show that carryover storage was negatively affected in all basins, particularly in those with lower storage capacity. All basins experienced negative impacts on hydropower generation, with losses ranging from 5% to nearly 90%. Reservoir sizes and inflexible operating rules are a particular challenge for flood control, as in extremely wet years spillage averaged nearly the annual basins' total discharge. The reliability of environmental flows and agricultural deliveries varied depending on the basin, intensity, and duration of whiplash sequences. Overall, wet years temporarily rebound negative drought effects, and greater storage capacity results in higher reliability and resiliency, and lesser volatility in services. We highlight potential policy changes to improve flexibility, increase resilience, and better equip managers to face challenges posed by whiplash while meeting human and environmental needs.

Environmental Sciences & Ecology↗

Microtomographic Measurements of Total Air‐Water Interfacial Areas for Soils

Abstract Synchrotron X‐ray microtomography (XMT) was used to measure total air‐water interfacial areas (A aw ) as a function of water saturation (S w ) for several soils that comprise a range of physical and geochemical properties. Measurements were also conducted for glass beads and quartz sands for comparison. Apparent near‐linearA aw ‐S w relationships are observed for the three sands and the three sandy soils. In contrast, the measured interfacial areas for two soils that contain greater proportions of silt and clay are strongly nonlinear functions of water saturation. The greater degree of nonlinearity observed for these two soils is due to their much greater particle‐size distributions (i.e., uniformity coefficients) and their concomitant greater range in pore sizes. Interfacial areas determined with the thermodynamic method were used to benchmark the XMT measurements. XMT‐measured interfacial areas compare well to the thermodynamic‐determined values for the sands and sandy soils. In contrast, the XMT‐measured interfacial areas for the two soils with larger particle‐size distributions are not fully congruent with the thermodynamic‐determined values. Both of these soils have large fractions of pore space comprising nominal pore diameters smaller than the resolution of the XMT imaging. These results suggest that air‐water interfacial area may not always be fully characterized by standard XMT for soils with large particle‐size distributions.

Environmental Sciences & Ecology↗

Explicit Consideration of Plant Xylem Hydraulic Transport Improves the Simulation of Crop Response to Atmospheric Dryness in the U.S. Corn Belt

Abstract Atmospheric dryness (i.e., high vapor pressure deficit, VPD), together with soil moisture stress, limits plant photosynthesis and threatens ecosystem functioning. Regions where rainfall and soil moisture are relatively sufficient, such as the rainfed part of the U.S. Corn Belt, are especially prone to high VPD stress. With globally projected rising VPD under climate change, it is crucial to understand, simulate, and manage its negative impacts on agricultural ecosystems. However, most existing models simulating crop response to VPD are highly empirical and insufficient in capturing plant response to high VPD, and improved modeling approaches are urgently required. In this study, by leveraging recent advances in plant hydraulic theory, we demonstrate that the VPD constraints in the widely used coupled photosynthesis‐stomatal conductance models alone are inadequate to fully capture VPD stress effects. Incorporating plant xylem hydraulic transport significantly improves the simulation of transpiration under high VPD, even when soil moisture is sufficient. Our results indicate that the limited water transport capability from the plant root to the leaf stoma could be a major mechanism of plant response to high VPD stress. We then introduce a Demand‐side Hydraulic Limitation Factor (DHLF) that simplifies the xylem and the leaf segments of the plant hydraulic model to only one parameter yet captures the effect of plant hydraulic transport on transpiration response to high VPD with similar accuracy. We expect the improved understanding and modeling of crop response to high VPD to help contribute to better management and adaptation of agricultural systems in a changing climate.

Environmental Sciences & Ecology↗

Time‐And‐Space Averaging Applied to Intermittent Multiphase Flow Experiments

Abstract Various researchers have studied fluctuations in pore‐scale phase occupancy during multiphase flow in porous media using synchrotron‐based X‐ray microcomputed tomography (micro‐CT). However, the impact of these fluctuations on the concept of a representative volume is not yet fully understood. In this study, we performed spatial and temporal averaging of multiphase flow experiments visualized with synchrotron‐based micro‐CT, focusing on oil saturation as the key parameter to determine a representative time‐and‐space average. Our findings revealed that a saturation value representative of both time and space was achieved during fractional flow experiments in drainage mode with fractional flows of 0.8, 0.5, and 0.3. Furthermore, we computed a range of relative permeabilities on the basis of whether momentaneous saturation or time‐and‐space averaged saturation was utilized for direct simulation. Our results highlighted the importance of time‐and‐space averaging in determining a representative relative permeability and indicated that the temporal and spatial scales covered in a typical micro‐CT flow experiment were sufficient to obtain a representative saturation value for sandstone rock under intermittent flow conditions.

Environmental Sciences & Ecology↗