Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ecological modeling”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Predictability and empirical dynamics of fisheries time series in the North Pacific

Previous studies have documented a strong relationship between marine ecosystems and large-scale modes of sea surface height (SSH) and sea surface temperature (SST) variability in the North Pacific such as the Pacific Decadal Oscillation and the North Pacific Gyre Oscillation. In the central and western North Pacific along the Kuroshio-Oyashio Extension (KOE), the expression of these modes in SSH and SST is linked to the propagation of long oceanic Rossby waves, which extend the predictability of the climate system to ~3 years. Using a multivariate physical-biological linear inverse model (LIM) we explore the extent to which this physical predictability leads to multi-year prediction of dominant fishery indicators inferred from three datasets (i.e., estimated biomasses, landings, and catches). We find that despite the strong autocorrelation in the fish indicators, the LIM adds dynamical forecast skill beyond persistence up to 5-6 years. By performing a sensitivity analysis of the LIM forecast model, we find that two main factors are essential for extending the dynamical predictability of the fishery indicators beyond persistence. The first is the interaction of the fishery indicators with the SST/SSH of the North and tropical Pacific. The second is the empirical relationship among the fisheries time series. This latter component reflects stock-stock interactions as well as common technological and human socioeconomic factors that may influence multiple fisheries and are captured in the training of the LIM. These results suggest that empirical dynamical models and machine learning algorithms, such as the LIM, provide an alternative and promising approach for forecasting key ecological indicators beyond the skill of persistence.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Data and scripts associated with “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” (v3)

This data package is associated with the publication “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” submitted to Journal of Advances in Modeling Earth Systems (Butler et al. 2025). This study developed the Sequential Precipitation Input Tagging (SPIT) framework to tag input precipitation and estimate water transit times and hydrologic tracers. SPIT tags all precipitation events at regular intervals over an extended period (monthly tags over seven years) in a hydrologic model from 2016-2022. SPIT is applied at six National Ecological Observatory Network (NEON) sites across the continental United States to calculate transit time distributions (TTD) and derive from these mean transit times (MTT), fractions of young water (Fyw), and hydrologic tracer concentrations in stream water (δ18O) within a water-tagging enabled version of the Weather Research and Forecast (WT-WRF-Hydro) model with national water model (NWM) configurations. We go on to validate WT-WRF-Hydro estimates against Butler et al. (2023), who analyzed the same NEON sites using stable water isotope data to estimate water transit times. This new tracking method provides a detailed picture of water movement and helps improve predictions about water availability in the future. This data package was originally published in January 2025. It was updated May 2025 (v2; new and modified files) and October 2025 (v3; new and modified files). File and folder names were not revised to indicate changes. See the change history section in the readme for more details. This data package contains the data and scripts used to develop the SPIT framework WT-WRF-Hydro (Water Tagging Weather Research and Forecasting Hydrologic) model and is associated with the following GitHub repository: https://github.com/zbutler33/SPIT-Framework. This data package contains five parent folders: (1) “Manipulated_outputs”, (2) “Metadata”, (3) “Observed”, (4) “Outputs”, and (5) “Scripts”. Each of these parent folders contains additional subfolders and files. Please see the FLMD (“v*_Butler_2024_WT_WRF_Hydro_flmd.csv”) for a list of all the files contained in this data package and descriptions for each. See the data dictionary (“v*_Butler_2024_WT_WRF_Hydro_dd.csv”) for definitions and units of all of the tabular (files ending in “.csv” and ".tsv") column headers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Deep Learning for Fish Identification from Sonar Data: CRADA 481 [Abstract only]

To help solve the challenges of hydropower energy production related to the potential for eel injury and mortality from passage through hydropower turbines, we will develop a deep learning method for identifying migrating eels from imaging sonar. This project continues with a prior project conducted by the Pacific Northwest National Laboratory (PNNL) and the Electric Power Research Institute (EPRI) in FY2018-2019. The proposed method employs Convolution Neural Network (CNN), a powerful deep learning method for image classification, to distinguish between images of eels and non-eel moving objects. We propose to collect more laboratory data and add more existing field data to train a powerful deep learning model. In addition to eels and sticks as classified in previous studies, we will add images containing several non-eel fish species and macrophyte mats to the training data. A multi-class classification model will be developed to distinguish these objects. Object detection algorithm will be explored and developed to locate and identify multiple objects in each sonar frame. Motion analysis will be performed to track the movement of objects in sonar video clips. We will also improve the data conversion algorithm so that it can read in both DIDSON and ARIS (both are imaging sonars developed by Sound Metrics Corp) data files and convert them to images with comparably high resolution, regardless of the varying detection ranges in different environments. The developed algorithms will be packaged as a software with a graphic user interface. The software will be evaluated by external collaborators in the field. The developed framework can be generalized for automatic monitoring of fish passage and migration using other imaging sonars like ARIS and will benefit the design and operation of ecologically friendly hydroelectric projects. The developed wavelet and CNN model configuration parameters can potentially be transferred to lamprey detection in similar riverine environments.

13 HYDRO ENERGY↗
DOE OSTI2022

In Silico Prediction of the Toxicity of Nitroaromatic Compounds: Application of Ensemble Learning QSAR Approach

In this work, a dataset of more than 200 nitroaromatic compounds is used to develop Quantitative Structure–Activity Relationship (QSAR) models for the estimation of in vivo toxicity based on 50% lethal dose to rats (LD 50 ). An initial set of 4885 molecular descriptors was generated and applied to build Support Vector Regression (SVR) models. The best two SVR models, SVR_A and SVR_B, were selected to build an Ensemble Model by means of Multiple Linear Regression (MLR). The obtained Ensemble Model showed improved performance over the base SVR models in the training set (R 2 = 0.88), validation set (R 2 = 0.95), and true external test set (R 2 = 0.92). The models were also internally validated by 5-fold cross-validation and Y-scrambling experiments, showing that the models have high levels of goodness-of-fit, robustness and predictivity. The contribution of descriptors to the toxicity in the models was assessed using the Accumulated Local Effect (ALE) technique. The proposed approach provides an important tool to assess toxicity of nitroaromatic compounds, based on the ensemble QSAR model and the structural relationship to toxicity by analyzed contribution of the involved descriptors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Comparative Study of Enhanced Geothermal System Supply Curves Across CONUS from Two Temperature Models Using the Renewable Energy Potential Model (reV)

Enhanced geothermal systems (EGS) have had recent breakthroughs within the geothermal sector. These breakthroughs are reflected in the National Renewable Energy Laboratory (NREL) 2024 Annual Technology Baseline and will result in updated EGS supply curves (i.e., the available resource capacity relative to cost). Our research uses NREL's Renewable Energy Potential (reV) model to compare EGS supply curves across the conterminous United States (CONUS) for two different temperature models: the Stanford temperature model (STM) and the Southern Methodist University temperature model (SMU). The reV model provides the levelized cost of energy (LCOE) at a consistent resolution across CONUS, taking into consideration transmission costs and constraints as well as technical exclusions pertaining to sensitive cultural, ecological, or infrastructure locations. In addition to a countrywide analysis of both models, we also conducted a regional analysis of Texas. We observed the STM had, on average, lower temperatures across different depths, resulting in slightly higher mean and median LCOEs as compared to the SMU temperature model at the same depths. In the regional analysis for Texas, however, when we compared only the common points between the two temperature models, the STM had lower median and mean LCOEs compared to SMU due to higher temperatures at depths greater than 5 km.

enhanced geothermal systems↗
DOE OSTI2020

Local-scale Arctic tundra heterogeneity affects regional-scale carbon dynamics

In northern Alaska nearly 65% of the terrestrial surface is composed of polygonal ground, where geomorphic tundra landforms disproportionately influence carbon and nutrient cycling over fine spatial scales. Process-based biogeochemical models used for local to Pan-Arctic projections of ecological responses to climate change typically operate at coarse-scales (1km 2 –0.5°) at which fine-scale (<1km 2 ) tundra heterogeneity is often aggregated to the dominant land cover unit. Here, we evaluate the importance of tundra heterogeneity for representing soil carbon dynamics at fine to coarse spatial scales. We leveraged the legacy of data collected near Utqiagvik, Alaska between 1973 and 2016 for model initiation, parameterization, and validation. Simulation uncertainty increased with a reduced representation of tundra heterogeneity and coarsening of spatial scale. Hierarchical cluster analysis of an ensemble of 21 st -century simulations reveals that a minimum of two tundra landforms (dry and wet) and a maximum of 4km 2 spatial scale is necessary for minimizing uncertainties (<10%) in regional to Pan-Arctic modeling applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Insights on Simulating Summer Warming of the Great Lakes: Understanding the Behavior of a Newly Developed Coupled Lake-Atmosphere Modeling System

The Laurentian Great Lakes are the world's largest freshwater system and regulate the climate of the Great Lakes region, which has been increasingly experiencing climatic, hydrological, and ecological changes. An accurate mechanistic representation of the Great Lakes thermal structure in Regional Climate Models (RCMs) is paramount to studying the climate of this region. Currently, RCMs have primarily represented the Great Lakes through coupled one-dimensional (1D) column lake models; this approach works well for small inland lakes but is unable to resolve the realistic hydrodynamics of the Great Lakes and leads to inaccurate representations of lake surface temperature (LST) that influence regional climate and weather patterns. This work overcomes this limitation by developing a fully two-way coupled modeling system using the Weather Research and Forecasting model and a three-dimensional (3D) hydrodynamic model. The coupled model system resolves the interactive physical processes between the atmosphere, lake, and surrounding watersheds; and validated against a range of observational data. The model is then used to investigate the potential impacts of lake-atmosphere coupling on the simulated summer LST of Lake Superior. By evaluating the difference between our two-way coupled modeling system and our observation-driven modeling system, we find that coupled-lake atmosphere dynamics can lead to a higher LST during June-September through higher net surface heat flux entering the lake in June and July and a lower net surface heat flux entering the lake in August and September. The unstratified water in June distributes the entering surface heat flux throughout the water column leading to a minor LST increase, while the stratified waters of July create a conducive thermal structure for the water surface to warm rapidly under the higher incoming surface heat flux. This research provides insight into the coupled modeling system behavior, which is critical for enhancing our predictive understanding of the Great Lakes climate system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evidence for the existence and ecological relevance of fast-cycling mineral-associated organic matter

Longstanding theories and models classify mineral-associated organic matter as the large ( ~ 60%) but slow-cycling and persistent portion of soil organic matter. Strong physico-chemical interactions and diffusion limitations restrict the turnover of mineral-associated organic matter, allowing carbon and nitrogen bound therein to persist in soil for as long as centuries to millennia. However, mineral-associated organic matter is a chemically and functionally diverse pool with a substantial portion cycling at relatively fast (i.e., minutes to years) timescales. Despite a growing body of evidence for the heterogenous and multi-pool nature of mineral-associated organic matter, we lack consensus on how to conceptualize and directly quantify fast-cycling mineral-associated organic matter and its ecological significance. We demonstrate that the dynamic qualities of fast-cycling mineral-associated organic matter vary based on 1) the chemistry of the mineral particles and organic matter, 2) the complex set of interactions between organic matter and the mineral matrix, and 3) the presence and strength of destabilizing forces that lead to decomposition or loss of mineral-associated organic matter (i.e., plant-microbe interactions, agricultural intensification, and climate change). Finally, we discuss potential implications and research opportunities for how we measure, manage, and model the dynamic subfraction of this otherwise persistent pool of soil organic matter.

Biological and medical sciences↗
DOE OSTI2020

NEON AOP foliar trait maps, maps of model uncertainty estimates, and conifer map, East River, CO 2018

This data package contains mapped trait estimates and their uncertainties, and conifer map, for the National Ecological Observatory Network's Airborne Observation Platform survey data acquired over the Upper East River, Colorado in 2018. For full details, please see associated reference. in brief, trait models were developed independently for needle and non-needle leaf species using partial least squares regression (PLSR) using ground data from additional datasets: doi:10.15485/1618130, doi:10.15485/1618132, and doi:10.15485/1631278, merged with extracted spectral data from doi:10.15485/1618131. We separated vegetated pixels into needle and non-needle classes in order to generate a classification map based on the spectral differences between these leaf types (conifer.tif). We trained a deep learning model with custom architecture, detailed in Chadwick et al. In Press. The model performed with 0.998 true positive rate and 0.982 true negative rate, with ‘positives’ being non-needle identification. We then utilized PLSR to generate models of foliar traits for each leaf type. So that we could also map uncertainty in these predictions, we generated ten different models for needle and non-needle leaf species using different testing holdout sets of discrete sites. Each of these models was developed with a 100-fold cross validation procedure that utilized a 70% training set and 30% validation set with each fold, and then assessed based on the 10% of testing sites that were not included in that model’s development. The mean predicted value across the 10 models is used for the trait estimate in each pixel across the study area. The models are applied according to the leaf type designation in the conifer.tif map. The errors are the standard deviation across the 10 different models developed, with high error suggesting instability in model prediction and areas where values may not be reliable for ecological inference. These maps are only applied to areas with a NDVI > 0.5 to exclude non-vegetated areas. Shade masks could be applied to these data (doi:10.15485/1618131), but have not been for this data package. These data are also available on Google Earth Engine: https://code.earthengine.google.com/?asset=users/kdc/ER_NEON

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Hydropower operation in future power grid with various renewable power integration

Hydropower generation may play an increasingly important role in the power grid under increasing contribution of variable renewable sources such as wind and solar. An improved understanding of the changes to hydropower dispatch under future higher VRE grid conditions reveals research gap that should be informed power grid planning and reservoir water releases policies considering multiple other water uses and varying hydrologic condition. Here this study aims to understand the role of hydropower in a changing power grid by employing a production cost model, PLEXOS, across future power system scenarios, planning horizons, and regions. We explore optimized hydropower dispatch to understand its potential role in minimizing the system cost and renewable curtailment. We also examine the sensitivity of hydropower revenue under various grid scenarios of the Eastern U.S. and hydrology conditions. Results indicate hydropower generation follows net load and compensates for the variability of solar and wind generation. Although energy prices are lower during some periods in the future grid scenarios, there is a potential for higher revenue for hydropower by providing both energy and ancillary services during times of stress. Additionally, hydropower revenue is sensitive to hydrology in the SERC region, which we considered as an example. The feasibility of hydropower dispatching with higher ramps between low and high hourly-capacity factors, as indicated in the optimization model, requires further study to consider other water use and ecology constraints.

13 HYDRO ENERGY↗
DOE OSTI2024

Predicting responses to climate change using a joint species, spatially dependent physiologically guided abundance model

Abstract Predicting the effects of warming temperatures on the abundance and distribution of organisms under future climate scenarios often requires extrapolating species–environment correlations to climatic conditions not currently experienced by a species, which can result in unrealistic predictions. For poikilotherms, incorporating species' thermal physiology to inform extrapolations under novel thermal conditions can result in more realistic predictions. Furthermore, models that incorporate species and spatial dependencies may improve predictions by capturing correlations present in ecological data that are not accounted for by predictor variables. Here, we present a joint species, spatially dependent physiologically guided abundance (jsPGA) model for predicting multispecies responses to climate warming. The jsPGA model uses a basis function approach to capture both species and spatial dependencies. We apply the jsPGA model to predict the response of eight fish species to projected climate warming in thousands of lakes in Minnesota, USA. By the end of the century, the cold‐adapted species was predicted to have high probabilities of extirpation across its current range—with 10% of lakes currently inhabited by this species having an extirpation probability >0.90. The remaining species had varying levels of predicted changes in abundance, reflecting differences in their thermal physiology. Though the model did not identify many strong species dependencies, the variation in estimated spatial dependence across species suggested that accounting for both dependencies was important for predicting the abundance of these fishes. The jsPGA model provides a new tool for predicting changes in the abundance, distribution, and extirpation probability of poikilotherms under novel thermal conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Growth and survival rates of dispersing free embryos and settled larvae of pallid sturgeon ( Scaphirhynchus albus ) in the Missouri River, Montana and North Dakota

We released nearly 1.0 million 1-day post-hatch (dph) and 5-dph pallid sturgeon (Scaphirhynchus albus) free embryos in the Missouri River on 1 July 2019 and sequentially captured survivors at multiple sites through a 240-km river reach to quantify daily growth and survival rates during the early life stages. Genetic analysis was used to assign captured fish to released family lots and known ages. Growth rate was similar (0.74–0.75 mm day –1 ) between the 1- and 5-dph age groups during the 3–4-day dispersal period when water temperature averaged 16.8 °C. Daily survival rate was 0.64 during 1–4 dph for the original 1-dph age group and 0.80 during 5–7 dph for the original 5-dph age group. Total survival during free embryo dispersal (hatch to 9 dph) was estimated as 0.0437. The transition from dispersing as free embryos to settling as benthic larvae was verified for fish originally released as 5 dph. Growth of settled larvae was quantified with a Gompertz model through 75 dph (9 September; 112 mm) when water temperature was 18.8–21.0 °C in the rearing areas. Settled larvae had an estimated daily survival rate of 0.96, and estimated total survival during 9–75 dph was 0.0714. This study provides the first empirical survival estimates for pallid sturgeon early life stages in natural settings and is one of few studies reporting similar information for other sturgeon species. Applications of this work extend to pallid sturgeon restoration programs where population models are being developed to predict recruitment potential and population responses to river management alternatives.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Sea level extremes and compounding marine heatwaves in coastal Indonesia

Low-lying island nations like Indonesia are vulnerable to sea level Height EXtremes (HEXs). When compounded by marine heatwaves, HEXs have larger ecological and societal impact. Here we combine observations with model simulations, to investigate the HEXs and Compound Height-Heat Extremes (CHHEXs) along the Indian Ocean coast of Indonesia in recent decades. We find that anthropogenic sea level rise combined with decadal climate variability causes increased occurrence of HEXs during 2010–2017. Both HEXs and CHHEXs are driven by equatorial westerly and longshore northwesterly wind anomalies. For most HEXs, which occur during December-March, downwelling favorable northwest monsoon winds are enhanced but enhanced vertical mixing limits surface warming. For most CHHEXs, wind anomalies associated with a negative Indian Ocean Dipole (IOD) and co-occurring La Niña weaken the southeasterlies and cooling from coastal upwelling during May-June and November-December. Our findings emphasize the important interplay between anthropogenic warming and climate variability in affecting regional extremes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Structured interactions drive abrupt transitions in the spatial organization of microbial communities

Bacteria possess diverse mechanisms to regulate their motility in response to environmental and physiological signals, enabling them to navigate complex habitats and adapt their behavior. Some of these mechanisms are species specific and enable cells to modulate their movement based on the ecological identity of neighboring species. Here, we introduce a model in which bacteria interact via local signals that either enhance or suppress the motility of neighboring cells depending on species type. Through large-scale simulations and a coarse-grained stochastic model, we demonstrate the emergence of a sharp transition driven by nucleation processes: increasing the density of motility-suppressing interactions drives the system from a fully mixed, motile phase to a state characterized by large, stationary bacterial clusters. Remarkably, in systems with a large number of interacting species, this transition can be triggered solely by altering the structure of the motility-regulation interaction matrix while maintaining species and interaction densities constant. In particular, we find that heterogeneous and modular interactions promote the transition more readily than homogeneous random ones. These findings add a dimension to the theory of motility-induced phase separation and contribute to the ongoing effort to understand microbial interactions, suggesting that structured, nonrandom ones may be key to reproducing commonly observed spatial patterns in microbial communities.

bacterial communities↗
DOE OSTI2021

Data-Driven Forgetting and Discount Factors for Vehicle Speed Forecasting in Ecological Adaptive Cruise Control

Abstract This paper investigates temporal correlations in human driving behavior using real-world driving to improve speed forecasting accuracy. These correlations can point to a measurement weighting function with two parameters: a forgetting factor for past speed measurements that the vehicle itself drove with, and a discount factor for the speeds of vehicles ahead based on information from vehicle-to-vehicle communication. The developed weighting approach is applied to a vehicle speed predictor using polynomial regression, a prediction method well-known in the literature. The performance of the developed approach is then assessed in both real-world and simulated traffic scenarios for accuracy and robustness. The new weighting method is applied to an ecological adaptive cruise control system, and its influence is analyzed on the prediction accuracy and the performance of the ecological adaptive cruise control in an electric vehicle powertrain model. The results show that the new prediction method improves energy saving from the eco-driving by up to 4.7% compared to a baseline least-square-based polynomial regression. This is a 10% improvement over the constant speed/acceleration model, a conventional speed predictor.

Automation & Control Systems↗
DOE OSTI2023

Perspectives on biorefineries in microbial production of fuels and chemicals

Microbes drive our complex biosphere by regulating the global ecosystem through cycling elements and energy. Humankind has barely begun leveraging this biotransformation capacity to impact global economies and ecologies. Advances in genetic engineering, molecular analysis, metabolic flux modeling, microbial consortia/biome mapping and engineering, cell-free bioproduction, artificial intelligence/machine learning and the ever expanding -omics frontiers have set the stage for paradigm changes to how humankind produces, uses, transforms, and recycles carbon and energy through microbes. Harnessing this enormous potential could drive a global bioeconomy and manage carbon at a planetary level but requires understanding and application at a grand scale across a broad range of science and engineering disciplines. The penultimate manifestation of these advances is the “bio-refinery”, which is often referenced, but is a long way from being fully developed as a global carbon management platform. Broadening the feed stocks, processing operations, and product portfolio to a sequential cascade optimizing the conversion as a whole instead of limited outputs could greatly advance deployment and stability of a bioeconomy.

09 BIOMASS FUELS↗
DOE OSTI2023

Variability in modelled reservoir greenhouse gas emissions: comparison of select US hydropower reservoirs against global estimates

Abstract Greenhouse gas (GHG) emissions from reservoirs have most often been evaluated on a global extent through areal scaling or linear-regression models. These models typically rely on a limited number of characteristics such as age, size, and average temperature to estimate per reservoir or areal flux. Such approaches may not be sufficient for describing conditions at all types of reservoirs. Emissions from hydropower reservoirs have received increasing attention as industry and policy makers seek to better understand the role of hydropower in sustainable energy solutions. In the United States (US), hydropower reservoirs span a wide range of climate regions and have diverse design and operational characteristics compared to those most heavily represented in model literature (i.e., large, tropical reservoirs). It is not clear whether estimates based on measurements and modeling of other subsets of reservoirs describe the diverse types of hydropower reservoirs in the US. We applied the Greenhouse Gas from Reservoirs (G‐res) emissions model to 28 hydropower reservoirs located in a variety of ecological, hydrological, and climate settings that represent the range of sizes and types of facilities within the US hydropower fleet. The dominant pathways for resulting GHG emissions estimates in the case-study reservoirs were diffusion of carbon dioxide, followed by methane ebullition. Among these case-study reservoirs, total post-impoundment areal flux of carbon ranges from 84 to 767 mgCm −2 d −1 , which is less variable than what has been reported through measurements at other US and global reservoirs. The net GHG reservoir footprint was less variable and towards the lower end of the range observed from modeling larger global reservoirs, with a range of 138 to 1,052 g CO 2 eq m −2 y −1 , while the global study reported a range of 115 to 145,472 g CO 2 eq m −2 y −1 . High variation in emissions normalized with respect to area and generation highlights the need to be cautious when using area or generation in predicting or communicating emissions footprints for reservoirs relative to those of other energy sources, especially given that many of the hydropower reservoirs in the US serve multiple purposes beyond power generation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Using Complex Probability Amplitudes to Simulate Solute Transport in Composite Porous Media

Probability amplitudes are fundamental to quantum mechanics and offer robust descriptions of complicated systems, which have allowed physicists to explain behaviors inaccessible to classical physics. This article ponders how some of the same conceptual underpinnings of the mathematics used for modeling quantum systems might be applied to subsurface water resources problems and speculates how these tools could facilitate applications on quantum computers. A probability amplitude-based model for describing advective-dispersive transport in porous media using linear operators is investigated. The proposed complex valued model decomposes spreading into two “sub-continuum partial dispersion” coefficients, and this recovers classical spreading when the sum of these coefficients is the Fickian dispersion coefficient. However, the probability amplitudes have a manyto-one relationship to a probability distribution, so it embeds a level of heterogeneity into seemingly equivalent functions. Two propagators with different sub-continuum coefficients may have the same macroscopic behavior when either is considered in isolation, but when they act on the other the system’s behavior changes. Additionally, differences in the amplitudes cause a reduction in spreading as velocity correlations are disrupted, despite both propagators having identical dispersion coefficients, and this cannot be achieved using classical methods without changing the dispersion coefficient. The main point is that these amplitude-based models offer a way to embed information about the system into the propagators, instead of just “averaging it out” when making an upscaled model.

54 ENVIRONMENTAL SCIENCES↗