Search NASA⌕ Search

SEARCH · Search NASA

Results for “model determination”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Improving radiograph analysis throughput through transfer learning and object detection

SIGN Fracture Care International partners with surgeons in low-resource hospitals worldwide to provide access to effective orthopedic care by donating educational materials and innovatively designed surgical implants. Over two decades, SIGN’s Online Surgical Database (SOSD) has grown to contain over 500,000 medical images, with radiographs holding the majority share. One challenge in working with hospitals worldwide is that both the radiographs uploaded to the SOSD and the data entry accompanying the uploads vary in quality. To improve the accuracy of data in the SOSD, we trained a model to detect surgical implants in radiographs. We first developed a tool to automatically detect radiographs, then trained an object detection model to determine the number and placement of surgical implants visible in the radiograph. Active learning was used to generate a training set containing 2,510 radiographs with screws, nails, and plates labeled by bounding boxes. Training a model to simultaneously recognize all three classes of implants gave a low average precision (AP) for the plate class, likely due to the low number of plate instances in our training set and the large variety of surgical plates used by SIGN-partnered surgeons. Applying standard image augmentation techniques to increase the plate count in our training set did not appreciably increase the AP of plate detection. To improve plate detection, we redrew the bounding boxes to account for correlations between the screw and plate classes. Training one model to detect nails and screws and a separate model to detect plates increased the AP of plate detection by 78.8 percentage points. The AP of each class was 80.7% for screws, 93.6% for nails, and 92.6% for plates; meanwhile, the sensitivity was 92% for screws, 86% for nails, and 81% for plates. We show that object detection methods can be used to detect surgical implants in radiographs of varying quality; however, the detection ability is dependent on the type of implant, and some implants, in our case plates, must be treated differently than others. Such tools can improve the throughput of radiograph analysis, assisting physicians and surgeons with the treatment of bone fractures.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

A Bayesian approach to evaluation of soil biogeochemical models

Abstract. To make predictions about the carbon cycling consequences of rising global surface temperatures, Earth system scientists rely on mathematical soil biogeochemical models (SBMs). However, it is not clear which models have better predictive accuracy, and a rigorous quantitative approach for comparing and validating the predictions has yet to be established. In this study, we present a Bayesian approach to SBM comparison that can be incorporated into a statistical model selection framework. We compared the fits of linear and nonlinear SBMs to soil respiration data compiled in a recent meta-analysis of soil warming field experiments. Fit quality was quantified using Bayesian goodness-of-fit metrics, including the widely applicable information criterion (WAIC) and leave-one-out cross validation (LOO). We found that the linear model generally outperformed the nonlinear model at fitting the meta-analysis data set. Both WAIC and LOO computed higher overfitting risk and effective numbers of parameters for the nonlinear model compared to the linear model, conditional on the data set. Goodness of fit for both models generally improved when they were initialized with lower and more realistic steady-state soil organic carbon densities. Still, testing whether linear models offer definitively superior predictive performance over nonlinear models on a global scale will require comparisons with additional site-specific data sets of suitable size and dimensionality. Such comparisons can build upon the approach defined in this study to make more rigorous statistical determinations about model accuracy while leveraging emerging data sets, such as those from long-term ecological research experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

From noise to information: The transfer function formalism for uncertainty quantification in reconstructing the nuclear density

The neutron distribution of neutron-rich nuclei provides critical information on the structure of finite nuclei and neutron stars. Parity violating experiments—such as PREX and CREX—provide a clean and largely model-independent determination of neutron densities. Such experiments, however, are challenging and expensive, which is why sound statistical arguments are required to maximize the information gained. We introduce a new framework, the “transfer function formalism,” aimed at uncertainty quantification, model selection, and experimental design in the context of neutron densities. The transfer functions (TFs) are built analytically by expressing the linear response of the objective function (e.g., χ2) to small perturbations of the data. Using the TF formalism, we are able to analyze the expected overall uncertainty—quantified in terms of bias and variance—of the mean square radius and interior density of 48 Ca and 208 Pb. Using relativistic mean field models as a proxy for the weak-charge density—and assuming that a total of five measurements could be performed on the weak form factor of 48 Ca and 208 Pb—we identify the optimal models and experimental locations that minimize the uncertainty in the extraction of the radius and interior density. We also explore the use of the TF formalism to understand the influence of prior distributions for the model parameters, as well as the optimization of model hyperparameters not constrained by the data. Here,we establish how the choice of experimental locations and the model that is used can have a significant impact on the final uncertainties of the extracted quantities of interest. For challenging experiments such as CREX and PREX, a proper quantification of such uncertainties is critical. We have demonstrated how the TF formalism provides several advantages for this type of analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

230 Th/ 234 U Model Age of International Atomic Energy Agency Sample 27107-14-04 Measured at Los Alamos National Laboratory

This report details the uranium age-dating of International Atomic Energy Agency (IAEA) bulk environmental cotton swipe sample 27107-14-04 at Los Alamos National Laboratory (LANL). The 230 Th/ 234 U radiochronometer was used to determine sample model age requiring the following analyses: U assay, U isotope composition, Th assay and Th isotope composition. Sample 27107-14-04 was received at LANL on November 3, 2025, along with cotton swipe blank 30124-24-20.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Optimizing green infrastructure placement under precipitation uncertainty

Increased urbanization, infrastructure degradation, and climate change threaten to overwhelm stormwater systems across the nation, rendering them ineffective. Green Infrastructure (GI) practices are low cost, low regret strategies that can contribute to urban runoff management. However, questions remain as to how to best distribute GI practices through urban watersheds given precipitation uncertainty and the variable hydrological responses to them. We develop stochastic programming models to determine the optimal placement of GI practices across a set of candidate locations in a watershed to minimize the total expected runoff under medium-term precipitation uncertainties. Specifically, we first develop a two-stage stochastic programming model. Next, we reformulate this model using perturbed parameters to reduce the requisite computational time and extend it to multi-stage. In addition, we introduce constraints that allow for incorporating sub-catchment-level runoff reduction considerations. Furthermore, we account for hydrological connectivity in the watershed using an underlying acyclic connectivity graph of sub-catchments and incorporate various practical considerations into the models. In addition, we develop a systemic approach to downscale the existing daily precipitation projections into hourly units and efficiently estimate the corresponding hydrological responses. These advancements are brought together in a case study for an urban watershed in a mid-sized city in the U.S., where we perform sensitivity analyses, evaluate the importance of the considered constraints, and provide insights.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

MSR Salt Spill Accident Testing Using Eutectic NaCl-UCl 3

Assessing the potential consequences of identified accident scenarios is an essential part of the licensing process for a new nuclear reactor and is achieved by using accident progression models. A likely accident scenario that will be evaluated by developers of molten salt reactors (MSRs) is a spill of radionuclide-bearing fuel salt onto the reactor containment floor (i.e., a salt spill accident). Models to determine the consequences of a molten salt spill accident need to be parameterized and validated using experimental data, but little experimental effort has been dedicated to fill these data gaps to date. Specifically, data is needed to enable model development for individual processes (e.g., spreading, heat transfer, corrosion, radionuclide vaporization, and aerosol generation) that quantifies the sensitivities towards the initial conditions of the spill, the ambient environment, and the features of the containment. In addition, integrated experiments that simulate salt spill accidents will need to be conducted to provide insight into coupled processes and data for model validation, but these experiments will require the use of proven methods to quantify the processes under evaluation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Data-Driven Simulation-Based Planning for Electric Airport Shuttle Systems: A Real-World Case Study

Many airports are adopting battery electric buses in their shuttle fleets due to concerns over air quality and regulations. This study proposes a simulation-based optimization modeling framework to help airport shuttle operators effectively deploy electric buses. We evaluated a planned airport electric shuttle system with an event-driven simulator. Empirical data collected from existing systems were used to drive the simulations. We then proposed a simulation-based optimization model to determine the battery capacity, charging power, and number of chargers so that predefined objective(s) (e.g., minimizing total capital cost, minimizing emissions) are optimized. Compared to existing studies, the primary contribution of the proposed method is that it can model the real-world stochastic nature of operations in an electric bus system with much higher fidelity. To demonstrate the proposed modeling framework, we study a real-world shuttle system at the Dallas-Fort Worth International Airport, and present extensive numerical studies. When considering partial fleet electrification, the model can provide a set of Pareto optimal solutions. When considering full fleet electrification, the optimal solution requires a 50-kWh battery capacity and four 210-kW chargers, resulting in a total capital cost of $26,744,000. The results demonstrate that the proposed modeling framework can effectively optimize the planning of electric airport shuttle systems with partial or full fleet electrification.

ADVANCED PROPULSION SYSTEMS,ENERGY PLANNING, POLIC↗
DOE OSTI2023

Knowledge-informed deep learning for hydrological model calibration: an application to Coal Creek Watershed in Colorado

Abstract. Deep learning (DL)-assisted inverse mapping has shown promise in hydrological model calibration by directly estimating parameters from observations. However, the increasing computational demand for running the state-of-the-art hydrological model limits sufficient ensemble runs for its calibration. In this work, we present a novel knowledge-informed deep learning method that can efficiently conduct the calibration using a few hundred realizations. The method involves two steps. First, we determine decisive model parameters from a complete parameter set based on the mutual information (MI) between model responses and each parameter computed by a limited number of realizations (∼50). Second, we perform more ensemble runs (e.g., several hundred) to generate the training sets for the inverse mapping, which selects informative model responses for estimating each parameter using MI-based parameter sensitivity. We applied this new DL-based method to calibrate a process-based integrated hydrological model, the Advanced Terrestrial Simulator (ATS), at Coal Creek Watershed, CO. The calibration is performed against observed stream discharge (Q) and remotely sensed evapotranspiration (ET) from the water year 2017 to 2019. Preliminary MI analysis on 50 realizations resulted in a down-selection of 7 out of 14 ATS model parameters. Then, we performed a complete MI analysis on 396 realizations and constructed the inverse mapping from informative responses to each of the selected parameters using a deep neural network. Compared with calibration using observations covering all time steps, the new inverse mapping improves parameter estimations, thus enhancing the performance of ATS forward model runs. The Nash–Sutcliffe efficiency (NSE) of streamflow predictions increases from 0.53 to 0.8 when calibrating against Q alone. Using ET observations, on the other hand, does not show much improvement on the performance of ATS modeling mainly due to both the uncertainty of the remotely sensed product and the insufficient coverage of the model ET ensemble in capturing the observation. By using observed Q only, we further performed a multiyear analysis and show that Q is best simulated (NSE > 0.8) by including in the calibration the dry-year flow dynamics that show more sensitivity to subsurface characteristics than the other wet years. Moreover, when continuing the forward runs till the end of 2021, the calibrated models show similar simulation performances during this evaluation period as the calibration period, demonstrating the ability of the estimated parameters in capturing climate sensitivity. Our success highlights the importance of leveraging data-driven knowledge in DL-assisted hydrological model calibration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Understanding and control of Zener pinning via phase field and ensemble learning

Zener pinning refers to the dispersion of fine particles which influences grain size distribution via movement of grain boundaries in a polycrystalline material. Grain size distribution in polycrystals has a significant impact on their properties including physical, chemical, mechanical, and optical to name a few. We explore the use of Phase-field modeling and machine-learning techniques to understand and improve the control of grain size distribution via Zener pinning in polycrystalline materials. We develop a machine learning model that determines the relative importance of various parameters to exercise microstructure control via Zener pinning. Our workflow combines high-throughput phase-field simulations and machine learning to address the computational bottlenecks associated with large-scale simulations as well as identify features necessary for microstructure control in polycrystals. A random forest (RF) regression model was developed to predict grain sizes based on five Phase-field model parameters, achieving an average prediction error of 0.72 nm for the training data and 1.44 nm for the test data. The importance of the input parameters is analyzed using the SHapley Additive exPlanations (SHAP) approach which reveals that diffusivity, volume fraction, and particle diameter are the most important parameters in determining the final grain size. These findings will allow us to select the best second-phase particles, optimize grain size distributions and thus design microstructures with the desired properties. The developed method is a highly versatile and generalizable approach that can be used to assess the combined effects of individual features in the presence of multiple variables.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hanford Waste Treatment Plant Effluent Management Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The Hanford Tank Waste Treatment and Immobilization Plant Effluent Management Facility (EMF) stack monitor location was qualified using a combination of scale model stacks to mitigate the risk of identifying that the sampling location does not meet the qualification criteria on the full-scale stack. The LV-S1 scale model stack was used as a baseline, augmented by the LB-S1 and LV-S2 scale model stacks to address the Direct Feed Low Activity Waste Effluent Management Facility Vessel Vent Process (DVP) injection into the main Active Confinement Ventilation (ACV) system duct. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale model and its sampling locations were geometrically similar to the actual stack, and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. The LV-S1, LB-S1, and LV-S2 scale model stack tests have met the criteria of the ANSI/HPS N13.1-1999 standard to demonstrate the stack sampling locations are well mixed. Verification tests of the EMF stack were performed at normal operating conditions. The minimum 1/6 DV value and the maximum 6 DV value from the scale model testing determine the range of stack flow rates for which the full-scale stack may be operated while remaining in compliance with the stack verification criterion. A practical range for the full-scale stack qualification uses the average DV through 6 DV from the scale model tests to compute the corresponding flow rates. Table S1 lists the operating flow rate along with the average and maximum qualified stack flow rate based on the LV-S1 scale model DV values. The operating flow is below the maximum qualified stack flow, which means that the scale model test results are acceptable for stack qualification. The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be ≤20°. Second, the velocity uniformity at the full-scale stack must be ≤20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack test at the EMF. Flow angle results were <5°; all flow angle results were within the ≤20° criterion. The velocity uniformity results for each test condition ranged between 2.2% COV and 4.3% COV, all of which were within the range of the target % COV values from the scale model tests on the LV-S1, LB-S1, and LV S2 scale models. Based on these stack verification test results, the EMF filtered exhaust stack sampling location meets the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all planned fan operating configurations. This includes each combination of ACV fans with DVP exhausters. Further changes to the system configuration or operating conditions that are outside the qualified flow rates described in this report may require additional tests or analyses to determine compliance with the standard.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Hanford Waste Treatment Plant LAB Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The Waste Treatment Plant laboratory (LAB) facility stack monitor locations were qualified using scale model stacks to mitigate the risk of identifying that sampling locations do not meet the qualification criteria on the full-scale stack. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale model and its sampling location were geometrically similar to the actual stack, and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the LAB stacks were performed at normal operating conditions. The minimum 1/6 DV value, along with the maximum 6 DV value from the scale model testing, determines the range of conditions for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. A practical range for the full-scale stack qualification uses the average DV through 6 DV from the scale model tests to compute the corresponding flow rates. Table S1 lists the operating flow rates along with the average and maximum qualified stack flow rates for each of the LAB facility stacks. For each stack, the operating flow is below the maximum qualified stack flow, which means that the scale model test results are acceptable for stack qualification. The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be ≤20°. Second, the velocity uniformity at the full-scale stack must be ≤20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack tests at the LAB facility. Flow angle results were primarily less than 10°, except for the LB-C2 Fan A results, which were an average of 13.7°; all flow angle results were within the ≤20° criterion. The velocity uniformity results for each test condition averaged between 1.5 and 4.1% COV, which were all within the range of the target percent coefficient of variation values from the scale model tests. Based on these stack verification test results, the three LAB filtered exhaust stack sampling locations meet the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all fan operating configurations. This includes single-fan as well as dual-fan operations for LB-C2, each of the dual-fan operating conditions for LB-S1, and each single-fan operating condition for LB-S2. Further changes to the system configuration or operating conditions that are outside the bounds described in this report may require additional tests or analyses to determine compliance with the standard.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Hanford Waste Treatment Plant LAB Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The WTP LAB stack monitor locations were qualified using scale model stacks to mitigate the risk of identifying that sampling locations do not meet the qualification criteria on the full-scale stack. As required by the ANSI/HPS N13.1-1999 standard, the scale model and its sampling location were geometrically similar to the actual stack and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the DV of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the LAB stacks were performed at normal operating conditions. The maximum 6 DV value from the scale model testing determines the maximum conditions for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be = 20°. Second, the velocity uniformity at the full-scale stack must be = 20% COV. Finally, the velocity uniformity results for the actual and scale model stacks must agree within 5% COV. In general, these criteria were met through the full-scale stack tests at the LAB facility. Some specific items for each scale model and full-scale stack comparison should be noted in assessing the validity of the verification tests. The LB-C2 scale model stacks were performed with single fan operations at then-minimum flow conditions. These conditions resulted in DV values that were too low to meet the DV range criterion, and therefore are unable to be used in qualifying single fan operations on the full-scale stack. However, the single fan operations demonstrated COV values for velocity uniformity as well as flow angles that were comparable to the dual-fan operations at high flows. Overall, the test conditions and test results were within the range of acceptable values based on current design flow rates. Single fan operations are expected to only occur infrequently for maintenance needs and is therefore not a planned operating condition at this time. The LB-S1 scale model stacks were performed at three fan combinations, and in each combination, at least one test was performed at the then-minimum flow conditions as well as at then-maximum flow conditions. While the then-minimum flow conditions result in a DV range that is lower than the verification test DV value, and therefore does not meet the criterion, these tests are un-necessary for the verification test acceptance. The normal and maximum flow conditions from the scale model stack tests meet the DV range criterion, and the velocity uniformity test results compare favorably with the full-scale stack results. The LB-S2 scale model stacks were performed at nominally the same maximum flow condition; however, the 6 DV value for Fan A operation was slightly lower than the full-scale stack test DV. The COV values from the scale model stacks were comparable between the Fan A and Fan B results, and the overall test conditions are within the range of acceptable values. The stack verification is therefore considered acceptable for both Fan A and Fan B operations. The verification tests were performed at flows that were appreciably higher than the design conditions, and further elevated flow rates would be beyond the range of acceptable DV. Based on these stack verification test results, the three LAB filtered exhaust stacks meet the qualification criteria provided in the ANSI/HPS N13.1-1999 standard. Further changes to the system configuration or operating conditions that are outside the bounds described in this and the scale model test reports (Glissmeyer, Flaherty, and Piepel (2001), Glissmeyer and Geeting (2013)) may require additional tests and additional analyses to determine compliance with the standard.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Black-Box Optimization for Design of Concentrating Solar Power and Photovoltaic Hybrid Systems with Optimal Dispatch Decisions

The hybridization of concentrating solar power (CSP) and photovoltaics (PV) can enable dispatchable renewable electricity generation at a lower price than current stand-alone CSP systems. However, designing a CSP-PV hybrid system can be challenging because of the many degrees of freedom in design that affect the internal and external system interactions and trade-offs. We develop a methodology to determine optimal designs for CSP-PV hybrids by implementing NLopt's derivative-free, or “black-box,” algorithms around pre-existing CSP-PV hybrid simulation software that utilizes the National Renewable Energy Laboratory’s System Advisor Model (SAM); we then employ a dispatch optimization model to determine operational decisions that maximize a plant’s profits. We present optimal designs for CSP-PV hybrid systems dispatching against four time-of-delivery (ToD) pricing structures. NLopt’s algorithms can improve the base case design’s power purchase agreement (PPA) price by 15% to 21%, depending on the ToD pricing structure. In addition, we present the resulting optimal CSP-PV hybrid design’s annual performance metrics, which tend to have capacity factors between 50% and 62%, but are able to generate electricity during the year’s highest-valued periods about 90% of the time. Lastly, we investigate the trade-offs between capacity factor and PPA price using Pareto fronts and demonstrate that, for some ToD pricing structures, the system capacity factor can increase by 20% but at the expense of a 2% increase in PPA price.

black-box↗
DOE OSTI2020

Using Fourier-Plane Imaging Microscopy for Determining Transition-Dipole-Moment Orientations in Organic Light-Emitting Devices

Here, we use Fourier-plane imaging microscopy (FIM) to determine the transition-dipole-moment orientation in doped organic emissive thin films. The use of FIM enables precise, sensitive, and rapid measurement of dipole orientation in the emission layer of an organic light-emitting device (OLED). An optical model of a stratified birefringent multilayer is introduced for interpreting results obtained by FIM. Using the model, we determine the average orientation of transition-dipole-moment vectors of three phosphorescent dopant emitters. The dipole alignment measured by FIM quantitatively explains the difference in OLED efficiencies using these archetype dopant molecules. FIM provides a nondestructive tool to measure and ultimately improve the outcoupling efficiency of OLEDs and other light-emitting devices.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

The Optical and Infrared Are Connected

Galaxies are often modeled as composites of separable components with distinct spectral signatures, implying that different wavelength ranges are only weakly correlated. They are not. We present a data-driven model that exploits subtle correlations between physical processes to accurately predict infrared (IR) Wide-field Infrared Survey Explorer (WISE) photometry from a neural summary of optical Sloan Digital Sky Survey spectra. The model achieves accuracies of $χ^{2}_{N} ≈ 1$ for all photometric bands in WISE, as well as good colors. We are able to tightly constrain typically IR-derived properties, e.g., the bolometric luminosities of active galactic nuclei (AGN) and dust parameters such as q PAH . We also test whether current spectral energy distribution (SED) fitting methods reproduce such panchromatic relations, but find their predictions biased and overconfident, likely due to model misspecification, with correlated biases in star-formation rates (SFRs) and AGN luminosities being most evident. To help improve SED models, we determine which features of the optical spectrum are responsible for our improved predictions, and identify several lines (Ca II , Sr II , Fe I , [O II ], and Hα), which point to the complex chronology of star formation and chemical enrichment being incorrectly modeled.

Jespersen, Christian Kragh [Princeton Univ., NJ (U↗
DOE OSTI2023

Analytical population dynamics underlying harmful algal blooms triggered by prey avoidance

A modified version of the NPZD ecosystem model is used to analytically examine the effects of predation avoidance, a possible mechanism for triggering harmful algal blooms (HAB). To resolve HAB development caused by predation avoidance, an additional phytoplankton functional group is considered, one that has slower nutrient uptake and better predation avoidance characteristics than the non-harmful phytoplankton group used in traditional NPZD models. Because the two phytoplankton groups (one non-harmful and one HAB) compete for only one resource within the same system, steady state (equilibrium) conditions cannot occur without the presence of zooplankton; only the non-harmful phytoplankton group, which defeats the HAB group in the resource competition, can survive in the equilibrium. The presence of sufficient zooplankton effectively acts to replenish the nutrient pool by consuming the non-harmful phytoplankton. When this occurs, two equilibrium states are found: one with both phytoplankton groups coexisting, and one that only includes the HAB group. Further, the condition required for equilibrium is that the total nitrogen within the system should be larger than a threshold determined by model coefficients. The threshold and feasibility of the equilibrium are sensitive to the relative HAB predation avoidance coefficient. If the coefficient is larger than the ratio of net growth rates between the HAB and non-harmful phytoplankton group, the threshold becomes infinite, and an equilibrium is not feasible. The time scale for the system to reach an equilibrium state that includes a HAB group is determined asymptotically. The dependence of a threshold condition as a controlling factor may explain the regime shift of dominant species causing HABs. The ecosystem model is fully implemented into the Regional Ocean Modeling System and applied to an idealized coastal embayment (with depths and geometry taken from San Francisco Bay) to show numerically the dominance of prey avoidance dynamics in a natural shallow water environment that includes advection and diffusion. The analytical results improve strategies for HAB modeling and provide guidance for setting model coefficients necessary to resolve a HAB event.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Geochemistry and Multiomics Data Differentiate Streams in Pennsylvania Based on Unconventional Oil and Gas Activity

Unconventional oil and gas (UOG) extraction is increasing exponentially around the world, as new technological advances have provided cost-effective methods to extract hard-to-reach hydrocarbons. While UOG has increased the energy output of some countries, past research indicates potential impacts in nearby stream ecosystems as measured by geochemical and microbial markers. Here, we utilized a robust data set that combines 16S rRNA gene amplicon sequencing (DNA), metatranscriptomics (RNA), geochemistry, and trace element analyses to establish the impact of UOG activity in 21 sites in northern Pennsylvania. These data were also used to design predictive machine learning models to determine the UOG impact on streams. We identified multiple biomarkers of UOG activity and contributors of antimicrobial resistance within the order Burkholderiales. Furthermore, we identified expressed antimicrobial resistance genes, land coverage, geochemistry, and specific microbes as strong predictors of UOG status. Of the predictive models constructed (n = 30), 15 had accuracies higher than expected by chance and area under the curve values above 0.70. The supervised random forest models with the highest accuracy were constructed with 16S rRNA gene profiles, metatranscriptomics active microbial composition, metatranscriptomics active antimicrobial resistance genes, land coverage, and geochemistry (n = 23). The models identified the most important features within those data sets for classifying UOG status. These findings identified specific shifts in gene presence and expression, as well as geochemical measures, that can be used to build robust models to identify impacts of UOG development.

16S rRNA↗
DOE OSTI2023

Synchronization of Alternative Models in a Supermodel and the Learning of Critical Behavior

“Supermodeling” climate by allowing different models to assimilate data from one another in run time has been shown to give results superior to those of any one model and superior to any weighted average of model outputs. The only free parameters, connection strengths between corresponding variables in each pair of models, are determined using some form of machine learning. It is demonstrated that supermodeling succeeds because near critical states, interscale interactions are important but unresolved processes cannot be effectively represented diagnostically in any single parameterization scheme. In two examples, a pair of toy quasigeostrophic (QG) channel models of the midlatitudes and a pair of ECHAM5 models of the tropical Pacific atmosphere with a common ocean, supermodels dynamically combine parameterization schemes so as to capture criticality, associated critical structures, and the supporting scale interactions. The QG supermodeling scheme extends a previous configuration in which two such models synchronize with intermodel connections only between medium-scale components of the flow; here the connections are trained against a third “real” model. Intermittent blocking patterns characterize the critical behavior thus obtained, even where such patterns are missing in the constituent models. In the ECHAM-based climate supermodel, the corresponding critical structure is the single ITCZ pattern, a pattern that occurs in neither of the constituent models. In conclusion, for supermodels of both types, power spectra indicate enhanced interscale interactions in frequency or energy ranges of physical interest, in agreement with observed data, and supporting a generalized form of the self-organized criticality hypothesis.

54 ENVIRONMENTAL SCIENCES↗