Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variables”

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 91 records · Page 5

Adaptive Discovery and Mixed-Variable Optimization of Next Generation Synthesizable Microelectronic Materials

Design of new microelectronic materials is characterized by several challenges such as high-dimensionality of the atomic structure-composition variable space, formidable cost of directly using high-fidelity simulations for design optimization, dispersity in literature-reported similar materials and synthesis methods, complex physical mechanisms, and mixed qualitative and quantitative design variables that lead to a disjointed design space. Even though machine learning (ML) techniques have been employed to expedite materials innovation, existing methods treat ML and design optimization as two separate processes, failing to resolve the fundamental challenges associated with high dimensionality and mixed-variable complexity. We have developed a ML enhanced mixed-variable material design optimization framework to efficiently extract useful information from existing data in literature and physics-based simulations to guide the autonomous search for optimal materials. Our proposed framework is composed of four computational modules: (1) a natural language processing (NLP) based virtual screening module, (2) classification based concept exploration module, (3) a density functional theory (DFT)-based high-fidelity evaluation model, and (4) a novel latent-variable Gaussian process (LVGP) ML model for mixed-variable problems with uncertainty quantification, which seamlessly integrates with Bayesian Optimization (BO) and achieves superb efficiency through embedded physics-based dimension reduction. Our approach is demonstrated and validated using the testbed of functional materials exhibiting metal-insulation transitions (MITs), with the targeted reversible resistivity changes (∼10^5) near room temperature. At the end of the 30-month project, we have developed a series of new ML techniques using NLP, conditional variational autoencoders, active learning, latent-variable Gaussian processes, integrated with Bayesian optimization. Our project has resulted in new predicted MITs compounds and improved understanding of MITs microscopic mechanisms, which in turn will revolutionize microelectronics science to provide energy-saving solutions. Our research has improved both creativity and efficiency in transforming rare-event discoveries of new functional materials to persistent innovations. In addition to open-sourcing the online MIT database and the classification model, the LVGP open source code has been downloaded more than 15,000 times within two years. More than 40 MIT compounds have been identified and many have been pursued experimentally via collaborators. The research results are published in close to 20 collaborative papers in high-impact journals, such as Chem. Mater., Appl. Phys. Rev., Sci. Rep., among others of design space.

36 MATERIALS SCIENCE↗

Evaluating Variable-Impedance Magnetically-Insulated Transmission Lines as a Risk-Mitigation Measure for Next-Generation Pulsed Power

This project has produced the first detailed characterizations of power flow resulting from applying the “variable-impedance MITL” concept to real-life systems in Sandia’s pulsed power program (Z and next-generation pulsed power (NGPP)). We present simulation results and analyses for constant-impedance versions of both Z and NGPP and survey the operational viability of several variable-impedance re-designs in the parameter space of linear tapers. Circuit modeling (SCREAMER/Bertha) was used to pinpoint promising candidate designs, and EM-PIC (Empire) simulations were used to evaluate these candidates more rigorously. This approach was particularly successful in the Z regime which resulted in the identification of several viable variable-impedance MITL designs for each level. The approach was more challenged in the operating space NGPP occupies, producing data points that speak to a more restrictive design space due to anode plasma turn-on. In the end, we were able to converge on one viable variable-impedance design for the highest inductance line (level “F”) and one for the highest current line (level “A”). Altogether, the body of simulation evidence presented in this report suggest there does exist flexibility in operating space for magnetically-insulated transmission lines (MITLs) having variable geometric impedance to be a potential enabling technology for safely increasing current delivery (and potentially lowering stack voltage) in pulsed-power drivers by manipulating electron losses; however, operating points for a particular design must be carefully screened. Circuit and EM-PIC modeling provided consistent verdicts in safe operating regimes for operational viability, but additional physics such as anode plasma turn-on which is included in Empire but not in SCREAMER/Bertha was found to be a critical factor affecting power flow that lead to different assessments between the codes. It is not always the case that the occurrence of anode plasma caused a design to fail (some designs turned on anode plasma yet still delivered load currents meeting design targets); the details matter such as how early in the pulse anode surfaces break down (and how large a region). However, in every case that it did fail it was found that the feedback from anode plasma was the cause (i.e., turning off the anode plasma model in Empire restored agreement with the circuit model prediction). As circuit simulations represent an efficient and practical means of surveying design space compared to more computationally-expensive approaches such as EM-PIC, it could be prudent to invest in the research and development of models to include the effects of anode plasma such as ion emission in circuit codes. The variable-impedance MITL design is a new concept that enables controlled manipulation of the initial electron losses in the outer MITL and can be tested on Z today. We encourage follow-on work to explore further optimization (including alternative variable-impedance profiles, e.g., having constant dZ/dR), and to confirm the major findings presented in this report by fielding test hardware on actual Z shots.

24 POWER TRANSMISSION AND DISTRIBUTION↗

New Blue and Red Variable Stars in NGC 247

Images recorded with the Gemini Multi-Object Spectrograph (GMOS) on Gemini South are combined with archival images from other facilities to search for variable stars in the southern half of the nearby disk galaxy NGC 247. Fifteen new periodic and nonperiodic variables are identified. These include three Cepheids with periods <25 days, four semiregular variables, one of which shows light variations similar to those of R CrB stars, five variables with intrinsic visible/red brightnesses and colors that are similar to those of luminous blue variables (LBVs), and three fainter blue variables, one of which may be a noneclipsing close binary system. The S Doradus instability strip defines the upper envelope of a distinct sequence of objects on the (i, g-i) color–magnitude diagram (CMD) of NGC 247. The frequency of variability of an amplitude ≥0.1 magnitude in the part of the CMD that contains LBVs over the seven-month period when the GMOS images were recorded is ~0.2. The light curve of the B[e] supergiant J004702.18–204739.9, which is among the brightest stars in NGC 247, is also examined. Low-amplitude variations on day-to-day timescales are found, coupled with a systematic trend in mean brightness over a six-month time interval.

47 OTHER INSTRUMENTATION↗

The Variability of the Black Hole Image in M87 at the Dynamical Timescale

The black hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5–61 days) is comparable to the 6 day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure-phase measurements on all six linearly independent nontrivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of ∼3°–5°. The only triangles that exhibit substantially higher variability (∼90°–180°) are the ones with baselines that cross the visibility amplitude minima on the u–v plane, as expected from theoretical modeling. We used two sets of general relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗

Variability-selected Intermediate-mass Black Hole Candidates in Dwarf Galaxies from ZTF and WISE

While it is difficult to observe the first black hole seeds in the early universe, we can study intermediate-mass black holes (IMBHs) in local dwarf galaxies for clues about their origins. In this paper we present a sample of variability-selected active galactic nuclei (AGN) in dwarf galaxies using optical photometry from the Zwicky Transient Facility (ZTF) and forward-modeled mid-IR photometry of time-resolved Wide-field Infrared Survey Explorer (WISE) co-added images. We found that 44 out of 25,714 dwarf galaxies had optically variable AGN candidates and 148 out of 79,879 dwarf galaxies had mid-IR variable AGN candidates, corresponding to active fractions of 0.17% ± 0.03% and 0.19% ± 0.02%, respectively. We found that spectroscopic approaches to AGN identification would have missed 81% of our ZTF IMBH candidates and 69% of our WISE IMBH candidates. Only nine candidates have been detected previously in radio, X-ray, and variability searches for dwarf galaxy AGN. The ZTF and WISE dwarf galaxy AGN with broad Balmer lines have virial masses of 10 5 M ⊙ < M BH < 10 7 M ⊙ , but for the rest of the sample, BH masses predicted from host galaxy mass range between 10 5.2 M ⊙ < M BH < 10 7.25 M ⊙ . We found that only 5 of 152 previously reported variability-selected AGN candidates from the Palomar Transient Factory in common with our parent sample were variable in ZTF. We also determined a nuclear supernova fraction of 0.05% ± 0.01% yr -1 for dwarf galaxies in ZTF. Our ZTF and WISE IMBH candidates show the promise of variability searches for the discovery of otherwise hidden low-mass AGN.

79 ASTRONOMY AND ASTROPHYSICS↗

Characterizing and Mitigating Intraday Variability: Reconstructing Source Structure in Accreting Black Holes with mm-VLBI

The extraordinary physical resolution afforded by the Event Horizon Telescope has opened a window onto the astrophysical phenomena unfolding on horizon scales in two known black holes, M87* and Sgr A*. However, with this leap in resolution has come a new set of practical complications. Sgr A* exhibits intraday variability that violates the assumptions underlying Earth aperture synthesis, limiting traditional image reconstruction methods to short timescales and data sets with very sparse (u, v) coverage. We present a new set of tools to detect and mitigate this variability. We develop a data-driven, model-agnostic procedure to detect and characterize the spatial structure of intraday variability. This method is calibrated against a large set of mock data sets, producing an empirical estimator of the spatial power spectrum of the brightness fluctuations. We present a novel Bayesian noise modeling algorithm that simultaneously reconstructs an average image and statistical measure of the fluctuations about it using a parameterized form for the excess variance in the complex visibilities not otherwise explained by the statistical errors. These methods are validated using a variety of simulated data, including general relativistic magnetohydrodynamic simulations appropriate for Sgr A* and M87*. We find that the reconstructed source structure and variability are robust to changes in the underlying image model. We apply these methods to the 2017 EHT observations of M87*, finding evidence for variability across the EHT observing campaign. The variability mitigation strategies presented are widely applicable to very long baseline interferometry observations of variable sources generally, for which they provide a data-informed averaging procedure and natural characterization of inter-epoch image consistency.

79 ASTRONOMY AND ASTROPHYSICS↗

Variability Analysis of FMVSS-121 Air Brake Systems: 60-mi/hr Service Brake System Performance Data for Truck Tractors

In support of the Federal Motor Carrier Safety Administration’s (FMCSA’s) ongoing interest in connected and automated commercial vehicles, this report summarizes analyses conducted to quantify variability in stopping distance tests conducted on commercial truck tractors. The data used were retrieved from tests performed under the controlled conditions specified for FMVSS-121 air brake system compliance testing. The report explores factors affecting the variability of the service brake stopping distance as defined by 49 CFR 571.121, S5.3.1 Stopping Distance—trucks and buses stopping distance. Variables examined in this analysis include brake type, weight, wheelbase, and tractor antilock braking system (ABS). This analysis uses existing test data collected between 2010 and 2019. Several of the examined parameters affected both tractor stopping distance and stopping distance variability. First, the average stopping distance for disc/disc brakes was shorter than either drum/drum or disc/drum brakes, but brake type did not have a statistically significant effect on stopping distance variability. Second, tractors with a gross vehicle weight rating (GVWR) of 45,000-50,000 lb had shorter stopping distances than any other examined weight category. Weight was found to have a significant effect on stopping distance variability; the 50,000-55,000-lb GVWR range had more variability than both the next lower (45,000-50,000 lb) and next higher (55,000-60,000 lb) range.

Lascurain, Mary Beth↗

Simulation of PV Variability as a Function of PV Generation and Plant Size: Preprint

The deployment of photovoltaic (PV) systems continues to show significant expansion; however, this growth has brought added attention to issues around the variability of the solar resource. Both spatial and temporal variability exist. Temporal scales can range from the sub-second to multiyear, whereas spatial scales can range from a few meters to tens of kilometers. There are multiple methods described in the literature to quantify PV variability at various spatial and temporal scales. This study focuses on short-term temporal variability and uses similar approaches with the addition of PV plant size a parameter to quantify variability. The method employed here incorporates the normalization of clear- and cloudy-sky conditions and PV plant size to quantify nominal variability metrics. The distribution and fluctuations of these metrics provide relevant information that is useful for system operations. The National Solar Radiation Database (NSRDB) is used to simulate PV variability as a function of PV generation and plant size. Hypothetical but realistic system information at 33 locations is used to model PV generation by feeding NSRDB solar irradiance data to the National Renewable Energy Laboratory’s System Advisor Model (SAM). Over the selected region, it is found that the aggregated ramp rates for the 1-minute data are associated with standard deviations ranging from 0.002–0.055 on a daily basis; however, hourly intervals induce higher aggregated ramp rates than the other timescales. Even though minute-to-minute variations are significant for the 1-minute timescale, the standard deviation aggregated into a daily metric is smaller because of the cancellation of values.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Coastal Sea Level Observations Record the Twentieth-Century Enhancement of Decadal Climate Variability

Changes in the amplitude of decadal climate variability over the twentieth century have been noted, with most evidence derived from tropical Pacific sea surface temperature records. However, the length, spatial coverage, and stability of most instrumental records are insufficient to robustly identify such nonstationarity, or resolve its global spatial structure. Here, it is found that the long-term, stable, observing platform provided by tide gauges reveals a dramatic increase in the amplitude and spatial coherence of decadal (11–14-yr period) coastal sea level (ζ) variability between 1960 and 2000. During this epoch, western North American ζ was approximately out of phase with ζ in Sydney, Australia, and led northeastern U.S. ζ by approximately 1–2 years. The amplitude and timing of changes in decadal ζ variability in these regions are consistent with changes in atmospheric variability. Specifically, central equatorial Pacific wind stress and Labrador Sea heat flux are highly coherent and exhibit contemporaneous, order-of-magnitude increases in decadal power. These statistical relationships have a mechanistic underpinning: Along the western North American coastline, equatorial winds are known to drive rapidly propagating ζ signals along equatorial and coastal waveguides, while a 1–2-yr lag between Labrador Sea heat fluxes and northeastern United States ζ is consistent with a remotely forced, buoyancy-driven, mechanism. Tide gauges thus provide strong independent support for an increase in interbasin coherence on decadal time scales over the second half of the twentieth century, with implications for both the interpretation and prediction of climate and sea level variability.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the influence of natural climate variability on in situ measurements of seasonal total and extreme daily precipitation

Abstract While various studies explore the relationship between individual sources of climate variability and extreme precipitation, there is a need for improved understanding of how these physical phenomena simultaneously influence precipitation in the observational record across the contiguous United States. In this work, we introduce a single framework for characterizing the historical signal (anthropogenic forcing) and noise (natural variability) in seasonal mean and extreme precipitation. An important aspect of our analysis is that we simultaneously isolate the individual effects of seven modes of variability while explicitly controlling for joint inter-mode relationships. Our method utilizes a spatial statistical component that uses in situ measurements to resolve relationships to their native scales; furthermore, we use a data-driven procedure to robustly determine statistical significance. In Part I of this work we focus on natural climate variability: detection is mostly limited to DJF and SON for the modes of variability considered, with the El Niño/Southern Oscillation, the Pacific–North American pattern, and the North Atlantic Oscillation exhibiting the largest influence. Across all climate indices considered, the signals are larger and can be detected more clearly for seasonal total versus extreme precipitation. We are able to detect at least some significant relationships in all seasons in spite of extremely large (> 95%) background variability in both mean and extreme precipitation. Furthermore, we specifically quantify how the spatial aspect of our analysis reduces uncertainty and increases detection of statistical significance while also discovering results that quantify the complex interconnected relationships between climate drivers and seasonal precipitation.

54 ENVIRONMENTAL SCIENCES↗

Strength and ductility of additively manufactured 316L stainless steel: Impact of neutron irradiation and data variability

Here, this article presents the mechanical properties of additively manufactured (AM) 316L stainless steel processed via the laser powder bed fusion (LPBF) method, focusing on the effects of neutron irradiation on mechanical properties and the variability in strength and ductility data. The rapid melting-solidification process and multiple heating-cooling cycles inherent in LPBF typically result in a fine, metastable microstructure with significant local variability. AM 316L builds of varying thicknesses were fabricated, and SS-J3 miniature tensile specimens were machined from six different locations. These specimens were irradiated with fast neutrons to doses of 2 and 10 dpa at target temperatures of 300 °C and 600 °C. Post-irradiation tensile tests were conducted at room temperature, 300 °C, and 600 °C. Compared to conventional 316L stainless steel, AM 316L exhibited higher initial strength but lower ductility. Irradiation at 300 °C caused significant hardening and prompt necking at yield, with limited uniform ductility, although embrittlement was not observed up to 10 dpa. While neutron irradiation, particularly at 600 °C, increased the variability in strength and ductility data, no clear dependence of mechanical properties on build thickness or sampling location was found—contrary to the conventional perception that AM materials may exhibit high property variability. Furthermore, we observed that the variability in property data for LPBF-processed 316L was relatively low compared to that of wrought 316L stainless steel. This reduced variability in AM 316L steel may be attributed to its highly metastable, stress-containing microstructure, which is discussed in the context of general tensile property variations.

Additively manufactured 316L stainless steel↗

Detecting multi‐scale riverine topographic variability and its influence on Chinook salmon habitat selection

Abstract Quantifying geomorphic conditions that impact riverine ecosystems is critical in river management due to degraded riverine habitat, changing flow and thermal conditions, and increasing anthropogenic pressure. Geomorphic complexity at different scales directly impacts habitat heterogeneity and affects aquatic biodiversity resilience. Here we showed that the combination of continuous spatial survey at high resolution, topobathymetric light detection and ranging (LiDAR), and continuous wavelet analysis can help identify and characterize that complexity. We used a continuous wavelet analysis on 1‐m resolution topobathymetry in three rivers in the Salmon River Basin, Idaho (USA), to identify different scales of topographic variability and the potential effects of this variability on salmonid redd site selection. On each river, wavelet scales characterized the topographic variability by portraying repeating patterns in the longitudinal profile. We found three major representative spatial wavelet scales of topographic variability in each river: a small wavelet scale associated with local morphology such as pools and riffles, a mid‐wavelet scale that identified larger channel unit features, and a large wavelet scale related to valley‐scale controls. The small wavelet scale was used to identify pools and riffles along the entire lengths of each river as well as areas with differing riffle‐pool development. Areas along the rivers with high local topographic variability (high wavelet power) at all wavelet scales contained the largest features (i.e., deepest or longest pools) in the systems. By comparing the wavelet power for each wavelet scale to Chinook salmon redd locations, we found that higher small‐scale wavelet power, which is related to pool‐riffle topography, is important for redd site selection. The continuous wavelet methodology objectively identified scales of topographic variability present in these rivers, performed efficient channel‐unit identification, and provided geomorphic assessment without laborious field surveys.

Duffin, Jenna↗

Assessment of rainfall variability and future change in Brazil across multiple timescales

Abstract Rainfall variability change under global warming is a crucial issue that may have a substantial impact on society and the environment, as it can directly impact biodiversity, agriculture, and water resources. Observed precipitation trends and climate change projections over Brazil indicate that many sectors of society are potentially highly vulnerable to the impacts of climate change. The purpose of this study is to assess model projections of the change in rainfall variability at various temporal scales over sub‐regions of Brazil. For this, daily data from 30 CMIP5 models for historical (1900–2005) and future (2050–2100) experiments under a high‐emission scenario are used. We assess the change in precipitation variability, applying a band‐pass filter to isolate variability on daily, weekly, monthly, intra‐seasonal, and El Nino Southern Oscillation (ENSO) time scales. For historical climate, simulated precipitation is evaluated against observations to establish model reliability. The results show that models largely agree on increases in variability on all timescales in all sub‐regions, except on ENSO timescales where models do not agree on the sign of future change. Brazil will experience more rainfall variability in the future that is, drier or more frequent dry periods and wetter wet periods on daily, weekly, monthly, and intra‐seasonal timescales, even in sub‐regions where future changes in mean rainfall are currently uncertain. This may provide useful information for climate change adaptation across, for example, the agriculture and water resource sectors in Brazil.

Alves, Lincoln M.↗

High resolution variability in wet deposition in the southeastern United States

Rainwater chemistry is determined by atmospheric pollutants and particles which vary spatially and temporally. Industrial and agricultural activities and meteorological events (e.g. sea breezes, severe weather, blowing dust) alter atmospheric particle and trace gas compositions. These gases and particles are scavenged by cloud and rain droplets that drive wet deposition. During an Intensive Operation Period (IOP) from April to October 2021, rainwater was collected at higher frequency intervals, usually daily, after precipitation events at three locations on the Savannah River Site (SRS). The farthest locations were separated by approximately 20 km. The mean concentration (μeq/L) of seven ions followed the Cl⁻ > SO 4 2− > Na⁺ > NO 3 ⁻ > K⁺ > Mg 2+ > Ca 2+ downward trend. Ion concentrations were compared to background ion concentrations from the National Atmospheric Deposition Program (NADP). The high frequency monthly averaged SRS data compared well with the monthly averaged NADP background but demonstrated extensive variability. In some months in 2021, the high frequency data compared better to the NADP site near the coast while in other months inland sites compared better. Strong spatial variability for ion concentrations was observed across SRS which was attributed to localized impacts in rainfall spatial variability. High frequency measurements allowed for comparison to regional weather patterns indicating influences from the Atlantic Ocean, Gulf of Mexico, and cities. This can account for spatial variability in the wet deposition flux. Sea breezes, Saharan dust, and anthropogenic sources were shown to impact wet deposition flux variability. Higher frequency precipitation chemistry sampling at numerous locations better captures ion concentration variability and improves measurement representativeness.

54 ENVIRONMENTAL SCIENCES↗

An embedded variable step IMEX scheme for the incompressible Navier–Stokes equations

This report presents a series of implicit–explicit (IMEX) variable stepsize algorithms for the incompressible Navier–Stokes equations (NSE). Here, IMEX means the nonlinear term is treated fully explicitly while the remaining terms are treated implicitly. With the advent of new computer architectures there has been growing demand for low memory solvers of this type. The addition of time adaptivity improves the accuracy and greatly enhances the efficiency of the algorithm. We prove energy stability of an embedded first–second order IMEX pair. For the first order member of the pair, we prove stability for variable stepsizes, and analyze convergence. We believe this to be the first proof of this type for a variable stepsize IMEX scheme for the incompressible NSE. We then define and test a variable stepsize, variable order IMEX scheme using these methods. Our work contributes several firsts for IMEX NSE schemes, including an energy argument and error analysis of a two-step, variable stepsize method, and embedded error estimation for an IMEX multi-step method.

42 ENGINEERING↗

A machine learning approach for identifying variables associated with risk of developing neutralizing antidrug antibodies to factor VIII

A key unmet need in the management of hemophilia A (HA) is the lack of clinically validated markers that are associated with the development of neutralizing antibodies to Factor VIII (FVIII) (commonly referred to as inhibitors). This study aimed to identify relevant biomarkers for FVIII inhibition using Machine Learning (ML) and Explainable AI (XAI) using the My Life Our Future (MLOF) research repository. The dataset includes biologically relevant variables such as age, race, sex, ethnicity, and the variants in the F8 gene. In addition, we previously carried out Human Leukocyte Antigen Class II (HLA-II) typing on samples obtained from the MLOF repository. Using this information, we derived other patient-specific biologically and genetically important variables. These included identifying the number of foreign FVIII derived peptides, based on the alignment of the endogenous FVIII and infused drug sequences, and the foreign-peptide HLA-II molecule binding affinity calculated using NetMHCIIpan. The data were processed and trained with multiple ML classification models to identify the top performing models. The top performing model was then chosen to apply XAI via SHAP, (SHapley Additive exPlanations) to identify the variables critical for the prediction of FVIII inhibitor development in a hemophilia A patient. Using XAI we provide a robust and ranked identification of variables that could be predictive for developing inhibitors to FVIII drugs in hemophilia A patients. These variables could be validated as biomarkers and used in making clinical decisions and during drug development. The top five variables for predicting inhibitor development based on SHAP values are: (i) the baseline activity of the FVIII protein, (ii) mean affinity of all foreign peptides for HLA DRB 3, 4, & 5 alleles, (iii) mean affinity of all foreign peptides for HLA DRB1 alleles), (iv) the minimum affinity among all foreign peptides for HLA DRB1 alleles, and (v) F8 mutation type.

60 APPLIED LIFE SCIENCES↗

A novel conditional generative model for efficient ensemble forecasts of state variables in large-scale geological carbon storage

Integrating monitoring data to efficiently update reservoir pressure and CO 2 plume distribution forecasts presents a significant challenge in geological carbon storage (GCS) applications. Inverse modeling techniques are commonly used to fuse observational data and refine reservoir model parameters, thereby improving state variable forecasts. However, these techniques often rely on linear or Gaussian assumptions, which can limit their effectiveness in accurately predicting state variables. Moreover, simulating large-scale three-dimensional (3D) GCS problems is computationally expensive, making iterative runs in inverse problems prohibitive. To address these challenges, we propose a conditional generative model utilizing the score-based diffusion method for real-time 3D pressure and saturation field distribution predictions. Our approach involves solving the score function with a mini-batch-based Monte Carlo estimator to generate labeled data. This data is subsequently employed to train a fully connected neural network, enabling it to learn the conditional sample generator within a supervised learning framework. This method enables the rapid generation of a large ensemble of predictions, facilitating comprehensive uncertainty quantification of state variables. Here we applied our method to forecast the dynamic 3D distributions of pressure and saturation fields over a 30-year injection period. The statistical assessment with low root mean square error (RMSE) values demonstrates that our method can accurately predict the spatiotemporal distributions of both pressure and saturation fields. Moreover, the developed conditional generative model shows high computational efficiency by generating 100 ensemble forecasts of 3D state variables in less than 10 min. The consistency between ensemble averages and ground truth values further illustrates the model’s capability to capture state variable dynamics during the CO 2 plume injection process. Notably, the ground truth values fall within the ensemble forecasts, indicating that our uncertainty quantification effectively captures variability and potential noise in the observations. Thus, the developed conditional generative model proves to be a more efficient, accurate, and practical tool for GCS applications, facilitating timely risk analysis and informed decision-making.

58 GEOSCIENCES↗

A data-driven global soil heterotrophic respiration dataset and the drivers of its inter-annual variability

Soil heterotrophic respiration (SHR), one of the primary carbon fluxes from terrestrial ecosystems to the atmosphere, is important for carbon-climate feedbacks because of its sensitivity to available litter and soil carbon, climatic conditions, and nutrient availability. However, until recently limited SHR data were available, and most published global SHR estimates have either a short time span, coarse spatial resolution, or reply on overly-simple model formulations. To better understand and quantify the global distribution of SHR and its sensitivity to climate variability, we produced a new global SHR dataset using Random Forest algorithms, up-scaling 455 point data from the Global Soil Respiration Database (SRDB 4.0) with gridded fields of climatic, edaphic and productivity as explanatory variables. We estimated a global total SHR of 46.8 Pg C yr-1 over 1985-2013 (95% confidence interval: 38.6-56.3 Pg C yr-1), with a significant increasing trend of 0.03 Pg C yr-2 during this period. We found that the choice of soil moisture datasets contributes more to the difference among these data-driven SHR members rather than that of productivity, temperature and precipitation data sources. We also analyzed the influence of climatic variables on the inter-annual variability (IAV) of our SHR product. Water availability was the dominant driver of IAV at global scales, although the inferred sensitivity depends on the choice of the soil moisture gridded dataset. At the ecosystem scale, temperature strongly controls the IAV of SHR in tropical forests, while water availability dominates in extra-tropical forest and semi-arid regions. Our machine-learning gridded SHR dataset and outputs from process-based land surface models (TRENDYv6) show agreement for a strong association between water variability and SHR IAV at the global scale, but the two approaches lead to different temporal trend globally and different controlling variables for IAV at the ecosystem scale. Our study provides evidence for the pervasive and important role of water availability in driving SHR, indicating both a direct effect limiting decomposition rates and an indirect effect through the amount of fresh organic matter made available to SHR from productivity. In consideration of potential limitations and uncertainties remaining in our data-driven SHR datasets, we call for a more scientifically designed observation network for SHR, more observation data compilation, and increased use of deep learning methods making maximum use of observation data in hand. This will benefit process-based models, and improve our understanding of SHR response to future anomalous environmental conditions.

Yao, Yitong↗