SEARCH · Engineering Papers

Results for “Variability analysis”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2026

CMIP7 data request: land and land ice priorities and opportunities

The Land and Land Ice Theme in the Coupled Model Intercomparison Project Phase 7 (CMIP7) represents the current understanding of physical processes in land surface ecosystems, hydrology, cryosphere, and their physical interactions with other Earth system components. Simulations from Earth system models (ESMs) could provide crucial information for assessing planetary safety, such as critical tipping elements, and be used to inform climate risks for improving climate impact assessments and policy decisions. This paper presents a collaborative effort to identify scientific opportunities in the Land and Land Ice Theme of the CMIP7 Data Request. The proposed opportunities build upon advances in ESMs, including new freshwater system and land ice processes being included in CMIP7, as well as the scientific community's demand for high-frequency and sub-grid-scale land surface outputs. In total, 25 variable groups that contain 716 variables have been identified to be potentially available to the broad scientific audience for performing analysis in land–atmosphere coupling, hydrological processes and freshwater systems, glacier and ice sheet mass balance and their influence on the sea levels, land use, and plant phenology. Key reflections from this data request effort include advocacy for closer engagement between the user community and modeling groups, reduction in the technical barriers to tracking existing parameters and defining new variables, and more streamlined variable management. These will be essential to enhance the usability and reliability of CMIP7 outputs for climate and Earth system research and applications to a broad audience that relies on the CMIP7 endeavor.

Li, Yue [Univ. of California, Los Angeles, CA (Uni↗

DOE OSTI2025

Symmetry Determining Equations of the Rankine-Hugoniot Equations for Variable Velocity Shock Waves

The “constant velocity piston” problem (Fig. 1), also known as the “piston problem,” is a standard model for a one dimensional, in our case linear, symmetric shock wave moving through an inviscid, perfect gas. The model can be divided into two regions - a perturbed section on the left and an unperturbed section on the right - by a moving shock wave moving left to right. Both the perturbed and unperturbed sections, i.e. the shocked and unshocked regions, respectively, obey the Eulerian conservation equations; however, at the exact location of the shock, there is a mathematical discontinuity not satisfied by the Euler equations. To ensure continuity and conservation of certain quantities when crossing between the unshocked and shocked regions, we evoke a series of equations derived from the Eulerian conservation equations, called the Rankine-Hugoniot equations, or “jump” equations as it is often referred to in the literature on the topic. The classical constant-velocity piston problem assumes the piston features a constant driving velocity (among many other willing suspensions of belief required in the pursuit of a first principles equation model); consequent to this assumption is a constant-velocity shock and a constant-velocity shocked flow state. However, using Lie Group Theory (LGT), also known as symmetry analysis, we can attempt to reinterpret the model with a shock wave of variable velocity in time and space. An extension of the model in this way opens up the possibility for obtaining new analytical solutions to the piston problem for certain shock velocity models. In this report, we use LGT to derive the symmetry determining equations (SDEs), whose solutions are Lie groups, which permit analytical solutions. In the future, we can then use the SDEs to define constraint equations on the shock velocity model and what the successive solutions to the Euler equations might be based off such constraints. This report is structured as follows: Section 2 provides a brief derivation of the Rankine-Hugoniot (“jump”) equations; Section 3 gives an overview of Lie group theory; Section 4 derives the SDEs of the jump equations; Section 5 derives the Euler conservation equations for fluids; and Section 6 presents concluding remarks and opportunities for future studies.

42 ENGINEERING↗

DOE OSTI2026

Mountain Basin Controls on the Snow-to-Streamflow Signal: An AIC-Weighted Multiple Linear Regression Framework

A regression-based analysis quantifies how basin characteristics modulate the snow-to-streamflow signal. First, we use the ERA5-Land reanalysis gridded product (European Centre for Medium Range Weather Forecasts reanalysis 5 -Land component) for 4,655 hydrologic unit code - 10 (HUC10) mountain basins across the western United States (US) for water years 1987–2024. Linear regressions are performed for peak snow water equivalent (SWE) and annual streamflow for each mountain basin. Models use ordinary least squares in Python’s statsmodels package. After which, an Akaike Information Criterion (AIC)–weighted ensemble multiple linear regression (MLR) framework with 47 watershed traits is used to predict the linear regression coefficient of determination (r-squared) defining the ability of peak SWE to predict annual streamflow across all mountain basin. Predictor sets are constrained to avoid multicollinearity by excluding models with variance inflation factors (VIF) greater than 5. Mountain basin traits included in the MLR include seasonal climate, topography, vegetation type and structure, and bedrock geology. Accepted models are considered if their AIC is within 2.0 of the model with the minimum AIC, or best model. To compare predictor influence across acceptable models, we computed standardized regression coefficients. To evaluate structural redundancy among models, we constructed binary inclusion vectors for each acceptable model, denoting whether a predictor was present (1) or absent (0). Core predictor variables are defined as occurring in at least 67% of the acceptable models. For this regional analysis, only one model was found acceptable, with higher snow-to-streamflow translation (higher r-squared) occurring in colder mountain basins with higher relative winter precipitation, more snow accumulation and a lower fraction of annual precipitation that falls in the spring and summer. The second component of the data package uses previously published, high-resolution output from an integrated hydrological model of the East River watershed using the U.S. Geological Survey Groundwater and Surface water Flow model (GSFLOW, doi:10.15485/1998576). East River MLR expands upon the approach described above to explore the response of five streamflow metrics—annual streamflow, runoff efficiency, 7-day minimum flow, low-flow duration, and non-perennial stream fraction to snow system indicators including peak SWE, snow-covered area, snow disappearance date, and the fraction of basin area characterized by low-to-no snow, as well as seasonal precipitation and temperature, and annual hydrologic variables representing soil moisture, evapotranspiration (ET), the partitioning of incoming precipitation to evapotranspiration (ET/P), groundwater storage, and groundwater inflow to streams. MLR was done on all water years (P0: 1987-2024) and for each period as determined in the split analysis using pooled regression techniques (P1: 1987-2011 and P2: 2012-2024) to evaluate shifting predictor variable emphasis on streamflow generation. Results indicate that since 2012, peak SWE has lost statistical strength in its prediction of annual streamflow and runoff efficiency, and the indirect influence of spring temperature has emerged as critically important. Low-flow metrics remain largely influenced by soil moisture, vegetation water use and groundwater inflows with summer precipitation becoming a direct influence on minimum summer flow. Together, these data and Python-based analysis tools provide a framework for identifying the key watershed characteristics that control how streamflow responds to snow from year to year. The package also helps quantify uncertainty in statistical models and assess how snow–streamflow relationships vary across regions and over time. This dataset contains comma-separated values files (.csv), text files (.txt), python code files (.py), figure files (.png), and shapefiles (.cpg, .dbf, .prj, .sbn, .sbx, .shp, .xml). Further details on file contents and MLR execution can be found in the readme file and the FLMD files. Work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Exploring the Role of Cryptic Nitrogen Fixers in Terrestrial Ecosystems: A Frontier in Nitrogen Cycling Research

Biological nitrogen fixation represents the largest natural flux of new nitrogen (N) into terrestrial ecosystems, providing a critical N source to support net primary productivity of both natural and agricultural systems. When they are common, symbiotic associations between plants and bacteria can add more than 100 kg N ha –1 y –1 to ecosystems. Yet, these associations are uncommon in many terrestrial ecosystems. In most cases, N inputs derive from more cryptic sources, including mutualistic and/or free-living microorganisms in soil, plant litter, decomposing roots and wood, lichens, insects, and mosses, among others. Unfortunately, large gaps remain in the understanding of cryptic N fixation. Here, we conducted a literature review to explore rates, patterns, and controls of cryptic N fixation in both unmanaged and agricultural ecosystems. Our analysis indicates that, as is common with N fixation, rates are highly variable across most cryptic niches, with N inputs in any particular cryptic niche ranging from near zero to more than 20 kg ha –1 y –1 . Such large variation underscores the need for more comprehensive measurements of N fixation by organisms not in symbiotic relationships with vascular plants in terrestrial ecosystems, as well as identifying the factors that govern cryptic N fixation rates. We highlight several challenges, opportunities, and priorities in this important research area, and we propose a conceptual model that posits an interacting hierarchy of biophysical and biogeochemical controls over N fixation that should generate valuable new hypotheses and research.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Data-Driven Insights into the Structural Essence of Plasticity in High-Entropy Alloys

The heterogeneous mechanical response of a crystalline alloy with multiple principal elements was investigated using molecular dynamics simulations. The local configuration of the alloy in its quiescent state was characterized by the variables derived from the gyration tensor and the atomic electronegativity. A multivariate analysis identified the geometric and chemical factors that influenced the atomic packing variations. Further, upon straining, the non-affine displacement exhibited spatial heterogeneity. A statistical correlation was established between the local yield events and the specific features of the local configuration. Our findings, validated by the performance metrics analysis, provided a structural criterion for the instability mechanisms in high-entropy alloys (HEAs) and enhanced the understanding of their plasticity.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Investigating capacity credit sensitivity to reliability metrics and computational methodologies

Assigning capacity value to renewable energy sources (RES) is a challenge faced in planning their integration with the grid. The difficulties stem from the natural characteristics of variability and intermittency of wind and solar sources. The capacity credit (CC) analysis evaluates the system’s actual power output compared with a constant capacity generator, i.e., conventional generator and determines an effective capacity to use for planning and operation. Herein this paper presents different factors that could affect the CC of a system. Two methods are proposed to determine the CC, namely equivalent firm capacity (EFC) and effective load carrying capability (ELCC). Since these methods are based on satisfying reliability criteria, daily loss of load expectation (LOLE), hourly loss of load (LOLH), and expected energy not served (EENS) have been employed as indices. To obtain the CC value, both methods apply two techniques: traditional and optimization. Genetic algorithm (GA) is the optimization approach used in this paper. Then, this work compares the two techniques and shows the superior performance of the optimization approach. Two hybrid systems, stand-alone (SA) and grid-connected (GC) modes, are proposed and used as case studies. The hybrid systems consist of photovoltaic (PV), wind turbine (WT), and battery energy storage system (BESS). In this work, three different scenarios are used to compare capacity credit: system as a whole, only wind, and no batteries. Finally, sensitivity analysis is carried out to examine the impact of varying the wind speed, solar irradiation, and load. It is demonstrated that the choice of reliability index plays an important role in determining the capacity credit and it is shown that EENS is a more comprehensive and consistent index of reliability.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

The root-zone soil moisture spectrum in a mediterranean ecosystem

Storage of water within soil pores of the root zone introduce memory effects in the dynamics of soil moisture that are considerably longer than the integral timescale of many atmospheric processes. Thus, hydro-climatic states can be “sustained” through land-surface heat and water vapor fluxes primarily because they can “feed off” on this long-term soil moisture memory. Root-zone soil moisture memory is only but one feature characterizing the spectrum of soil moisture dynamics, which is analyzed here using a combination of long-term measurements and models. In particular, the spectrum of root-zone soil moisture content in a Mediterranean ecosystem is examined using 14-years of half-hourly measurements. A distinguishing hydro-climatic feature in such ecosystems is that sources (mainly rainfall) and sinks (mainly evapotranspiration) of soil moisture are roughly out of phase with each other. For over 4 decades of time scales and 7 decades of energy, the canonical shape of the measured soil moisture spectrum is shown to be approximately Lorentzian determined by the soil moisture variance and its memory but with two exceptions: the occurrences of a peak at diurnal-to-daily time scales and a weaker peak at near annual time scales. Model calculations and spectral analysis demonstrate that diurnal and seasonal variations in hydroclimate forcing responsible for variability in evapotranspiration had minor impact on the normalized shape of the soil moisture spectrum. However, their impact was captured by adjustments in the temporal variance. These findings indicate that precipitation and not evapotranspiration variability dominates the multi-scaling properties of soil moisture variability consistent with prior climate model simulations. Furthermore, the soil moisture memory inferred by the annual peak of soil moisture (340 d) is consistent with climate model simulations, while the memory evaluated from the loss function of a linearized mass balance approach leads to a smaller value (50 d), highlighting the effect of weak non-stationarity on soil moisture variability.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Direct sensitivity analysis on the parameterization of crystal plasticity models

Various methods for calibrating crystal plasticity finite element (CPFE) models lead to non-unique input parameter values, which subsequently introduce uncertainty in the predicted mechanical response. Sensitivity analysis (SA) conducted on crystal plasticity models is used to identify how variability in these parameters contribute to output uncertainty. Traditional SA on CPFE parameters uses simplified surrogate models to save computational time. However, the accuracy of the surrogate models depends on the quantity of training data used, and any modeling error can propagate into the SA results, potentially affecting their reliability. In this work, the elementary effects test (EET) method, a global SA technique using direct CPFE simulations was employed, and the results obtained were compared with the First Order Second Moment (FOSM) method. ExaConstit, an open-source GPU-enabled CPFE code, was used to perform the simulations and direct SA. The EET method was accurately able to capture the non-linear effects of all the input parameters on the output and is a valuable approach for reliably attributing parameter sensitivities in CPFE models. Based on the results, efficient strategies to perform future parameter calibration and SA are discussed. Additionally, the SA trends observed in different single crystal orientations closely mirrored the activity of the slip systems.

Elementary Effects Test↗

DOE OSTI2020

Probing Conformational Evolution and Associated Dynamics of Mg(N(SO 2 CF 3 ) 2 ) 2 ·Dimethoxyethane Adduct Using Solid-State 19 F and 1 H NMR

Bis(trifluoromethanesulfonimide) or TFSI is widely used as a counter anion in electrolyte design due to its structural flexibility and chemical stability. Here we studied the conformational variations and associated dynamics of TFSI in adduct of Mg(TFSI) 2 with dimethoxyethane (DME), a solvate crystalline material using solid-state 1 H and 19 F NMR. TFSI molecular motion in this solvate structure falls within the timescale of the 19 F NMR experiment, yielding spectroscopic signatures for unique TFSI conformers under the coordination environment of Mg 2+ cation. Within the temperature range of -5 to 82 °C, we observe nine distinct TFSI sites in both crystalline and disordered regions using 19 F NMR, reflecting complexity of structural and dynamics of TFS1 anions within solvate structure. The four distinguishable sites in the disordered region for the two CF 3 groups of the same TFSI molecule are identified using chemical shift analysis. The exchange rate constants from site to site are calculated through variable temperature 19 F NMR and two-dimensional (2D) exchange spectroscopy (EXSY) experiments, along with respective activation enthalpies using Eyring's formulation. The flip rate of CF 3 around the S-C bond is estimated as ~ 15 s -1 at 8 °C with ΔH ≠ ~ 22 kJ/mol, but the rotation of the entire TFSI is 4.8 s -1 at 8 °C with a significantly greater ΔH ≠ = 98 ± 10 kJ/mol. Furthermore, the slow conversion of trans to cis conformers at a lower temperature (T ≤ 1 °C) in the crystalline region is monitored, with a conversion rate of ~ 2 x 10 -5 s -1 at -5 °C. Density functional theory (DFT)-based calculations were performed to support further the assignment of experimental chemical shifts, and the activation energy E a = 21.1 kJ/mol obtained for the cis to trans transition is consistent with experimental values. The combined set of 19 F and 1 H under both one-dimensional (1D) and 2D NMR methods demonstrated here can be further used for examining electrode-electrolyte interfaces to probe the motions of various constituents that can enable detailed studies of interfacial processes and dynamics. Ultimately, such studies will aid in the design and discovery of interfacial constructs in which directed defect chemistry, chemical moiety distribution, and nanostructure are employed to drive efficient charge transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Resolving the Valence of Iron Oxides by Resonant Photoemission Spectroscopy

Precisely determining the oxidation states of metal cations within variable-valence transition metal oxides remains a significant challenge, yet it is crucial for understanding and predicting the properties of these technologically important materials. Iron oxides, in particular, exhibit a remarkable diversity of electronic structures due to the variable valence states of iron (Fe 2+ and Fe 3+ ). A quantitative analysis using conventional X-ray photoelectron spectroscopy (XPS) is challenging because of the strong overlap of the Fe 2p XPS peaks from different oxidation states. Here, in this study, we show how this problem can be resolved using Resonant Photoemission Spectroscopy (ResPES), which unambiguously distinguishes Fe oxidation states and spectroscopically estimates the composition ratio of Fe cation valence states in the complex Fe oxides. We demonstrate this in the model case of a FeO 2 monolayer film on Pt(111), showing that the FeO 2 film consists of an equal mixture of Fe 2+ and Fe 3+ cations, yielding an average valence of +2.5, contrary to the +3 valence proposed based on density functional theory (DFT).

Chen, Hao [Lawrence Berkeley National Laboratory (↗

DOE OSTI2020

Interactions Between the Amazonian Rainforest and Cumuli Clouds: A Large‐Eddy Simulation, High‐Resolution ECMWF, and Observational Intercomparison Study

Abstract The explicit coupling at meter and second scales of vegetation's responses to the atmospheric‐boundary layer dynamics drives a dynamic heterogeneity that influences canopy‐top fluxes and cloud formation. Focusing on a representative day during the Amazonian dry season, we investigate the diurnal cycle of energy, moisture and carbon dioxide at the canopy top, and the transition from clear to cloudy conditions. To this end, we compare results from a large‐eddy simulation technique, a high‐resolution global weather model, and a complete observational data set collected during the GoAmazon14/15 campaign. The overall model‐observation comparisons of radiation and canopy‐top fluxes, turbulence, and cloud dynamics are very satisfactory, with all the modeled variables lying within the standard deviation of the monthly aggregated observations. Our analysis indicates that the timing of the change in the daylight carbon exchange, from a sink to a source, remains uncertain and is probably related to the stomata closure caused by the increase in vapor pressure deficit during the afternoon. We demonstrate quantitatively that heat and moisture transport from the subcloud layer into the cloud layer are misrepresented by the global model, yielding low values of specific humidity and thermal instability above the cloud base. Finally, the numerical simulations and observational data are adequate settings for benchmarking more comprehensive studies of plant responses, microphysics, and radiation.

Vilà‐Guerau de Arellano, J.↗

DOE OSTI2022

Sulfur Biogeochemical Cycling and Redox Dynamics in a Shale-Dominated Mountainous Watershed

Sulfur (S) is an essential macronutrient and important component of the earth’s crust, and its cycling has critical impacts on trace metal mobility, water quality, and human health. Pyrite weathering is the primary pathway by which sulfur enters surface waters. However, biogeochemical cycling of sulfur in soils and the river corridor mediates sulfate exports. In this study, we identified the major forms of sulfur across multiple compartments and scales in a pristine mountainous watershed, including shale bedrock weathering profiles, hillslope soils, and alluvial floodplain sediments, in order to provide insight into biogeochemical sulfur cycling in a hydrologically variable alpine system. X-ray absorption near-edge spectroscopy (XANES) analysis of shale weathering profiles showed clear evidence of pyrite oxidation to sulfate, with large accumulations of intermediate S(0) (20%–53%). Micro-scale XANES showed evidence of reprecipitation of pyrite at fracture surfaces within the permanently saturated zone. Organic sulfur dominated S speciation in shallow hillslope soil and floodplain sediment, with little evidence of reduced inorganic S. However, mackinawite formation, representing active sulfate reduction, was observed in saturated oxbow sediments and saturated weathered shale underlying floodplain sediments. Further evidence of sulfate reduction from aqueous sulfur isotopic analysis was observed in shallow groundwater transects across an Fe-reducing meander, whereas increases in pore water sulfate concentrations implied sulfur oxidation at other locations. Furthermore, the data present an integrated picture of sulfur cycling in a shale-dominated watershed, where riverine sulfate exports are mediated by biological cycling, particularly in redox-stratified and temporally dynamic hyporheic zone sediments.

58 GEOSCIENCES↗

DOE OSTI2023

The Marginal Ice Zone as a dominant source region of atmospheric mercury during central Arctic summertime

Atmospheric gaseous elemental mercury (GEM) concentrations in the Arctic exhibit a clear summertime maximum, while the origin of this peak is still a matter of debate in the community. Based on summertime observations during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition and a modeling approach, we further investigate the sources of atmospheric Hg in the central Arctic. Simulations with a generalized additive model (GAM) show that long-range transport of anthropogenic and terrestrial Hg from lower latitudes is a minor contribution (~2%), and more than 50% of the explained GEM variability is caused by oceanic evasion. A potential source contribution function (PSCF) analysis further shows that oceanic evasion is not significant throughout the ice-covered central Arctic Ocean but mainly occurs in the Marginal Ice Zone (MIZ) due to the specific environmental conditions in that region. Our results suggest that this regional process could be the leading contributor to the observed summertime GEM maximum. In the context of rapid Arctic warming and the observed increase in width of the MIZ, oceanic Hg evasion may become more significant and strengthen the role of the central Arctic Ocean as a summertime source of atmospheric Hg.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Long-term socioeconomic trends and climate variability as drivers of virtual water scarcity in China

Water scarcity can have far-reaching sectoral impacts beyond its physical location through the propagation of virtual water flows. Socioeconomic and hydroclimatic changes affect local and virtual water scarcity by altering availability and demand. Yet most studies of this phenomenon focus on volumetric footprints, and the few on water scarcity risk have not examined hydroclimatic variability beyond long-term trends. In this study, we ask how gross domestic product (GDP) and population changes, long-term meteorological trends, sea surface temperature (SST) patterns, and interannual meteorological variability affect water scarcity in China, both locally (through the local water scarcity risk index, LWSR) and remotely (through the virtual water scarcity risk index, VWSR). Counterfactual scenarios were compared in a regression-and-simulation framework, with the socioeconomic and meteorological drivers varying over 1923–2019 and the multi-regional input–output structure staying fixed at 2017. Relative to a 5 year baseline centered on 2017, GDP and population changes have induced a cumulative 17%–50% increase in LWSR and 13%–21% increase in VWSR, outweighing the effect of long-term meteorological trends. phase change in one of two examined SST patterns induce 4%–13% differences in LWSR and 1%–4% differences in VWSR. Interannual meteorological variability induces 10%–20% standard deviations in LWSR and 3%–7% in VWSR. The findings highlight the importance of using longer time series to accurately assess local and virtual water scarcity situations. Water scarcity management should prioritize socioeconomic factors when planning at century-long timescales and prioritize hydroclimatic factors at multidecadal or shorter timescales. water managers should consider interannual variabilities in LWSR and VWSR and plan for potential occurrences of extreme conditions.

climate variability↗

DOE OSTI2021

Anomaly detection in the Zwicky Transient Facility DR3

We present results from applying the SNAD anomaly detection pipeline to the third public data release of the Zwicky Transient Facility (ZTF DR3). The pipeline is composed of three stages: feature extraction, search of outliers with machine learning algorithms, and anomaly identification with followup by human experts. Our analysis concentrates in three ZTF fields, comprising more than 2.25 million objects. A set of four automatic learning algorithms was used to identify 277 outliers, which were subsequently scrutinized by an expert. From these, 188 (68 per cent) were found to be bogus light curves – including effects from the image subtraction pipeline as well as overlapping between a star and a known asteroid, 66 (24 per cent) were previously reported sources whereas 23 (8 per cent) correspond to non-catalogued objects, with the two latter cases of potential scientific interest (e.g. one spectroscopically confirmed RS Canum Venaticorum star, four supernovae candidates, one red dwarf flare). Moreover, using results from the expert analysis, we were able to identify a simple bi-dimensional relation that can be used to aid filtering potentially bogus light curves in future studies. We provide a complete list of objects with potential scientific application so they can be further scrutinised by the community. These results confirm the importance of combining automatic machine learning algorithms with domain knowledge in the construction of recommendation systems for astronomy. Our code is publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

eDNAjoint: An R package for interpreting paired or semi‐paired environmental DNA and traditional survey data in a Bayesian framework

Abstract Environmental DNA (eDNA) sampling is increasingly used in surveys of species distribution as a potentially sensitive and efficient monitoring method. Yet access to modelling tools designed specifically for interpreting this new data type lags behind its ubiquity. While occupancy modelling software has dominated the analytical landscape for eDNA data analysis of single species, this type of model may not always be the most appropriate. The rate of eDNA detection often corresponds to species density, rather than just occupancy, and researchers often have access to observations from non‐genetic sampling methods at the same sites. To provide users access to a modelling framework designed to maximize the use of all available data, we developed an R package, eDNAjoint . The package provides an easy‐to‐use interface for fitting a ‘joint’ model that integrates data from paired or semi‐paired eDNA and traditional surveys in a Bayesian framework. The model can be used to estimate parameters like the probability of a false positive eDNA detection and mean catch rate at a site, and the package allows access to multiple model variations and Bayesian prior customization. Additional functionality can be used for model selection, summarising posteriors and comparing the relative sensitivities of the two survey methods. We demonstrate the use of eDNAjoint by fitting a variation of the model with site‐level covariates that scale the sensitivity of eDNA sampling relative to traditional sampling. The example workflow uses binary eDNA and seine count data for the endangered tidewater goby ( Eucyclogobius newberryi ) from a study by Schmelzle and Kinziger (2016). This use case includes a prior sensitivity analysis and an evaluation of the relationship between detection rates and environmental variables. eDNAjoint has the potential to greatly increase the range of users who will be able to rigorously analyse eDNA and traditional survey data in a Bayesian framework, understand if and how eDNA can improve monitoring practices, and gain confidence in the interpretability of eDNA data.

Keller, Abigail G. [Department of Environment Scie↗

DOE OSTI2023

maslov-group/FUN-PROSE

mRNA levels of all genes in a genome is a critical piece of information defining the overall state of the cell in a given environmental condition. Being able to reconstruct such condition-specific expression in fungal genomes is particularly important to metabolically engineer these organisms to produce desired chemicals in industrially scalable conditions. Most previous deep learning approaches focused on predicting the average expression levels of a gene based on its promoter sequence, ignoring its variation across different conditions. Here we present FUN-PROSE—a deep learning model trained to predict differential expression of individual genes across various conditions using their promoter sequences and expression levels of all transcription factors. We train and test our model on three fungal species and get the correlation between predicted and observed condition-specific gene expression as high as 0.85. We then interpret our model to extract promoter sequence motifs responsible for variable expression of individual genes. We also carried out input feature importance analysis to connect individual transcription factors to their gene targets. A sizeable fraction of both sequence motifs and TF-gene interactions learned by our model agree with previously known biological information, while the rest corresponds to either novel biological facts or indirect correlations.

Liu, Simon↗

DOE OSTI2023

FUN-PROSE: A deep learning approach to predict condition-specific gene expression in fungi

59 BASIC BIOLOGICAL SCIENCES↗