Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surrogates”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

The influence of TRP53 in the dose response of radiation-induced apoptosis, DNA repair and genomic stability in murine haematopoietic cells

Apoptotic and DNA damage endpoints are frequently used as surrogate markers of cancer risk, and have been well-studied in the Trp53+/- mouse model. We report the effect of differing Trp53 gene status on the dose response of ionizing radiation exposures (0.01-2 Gy), with the unique perspective of determining if effects of gene status remain at extended time points. Here we report no difference in the dose response for radiation-induced DNA double-strand breaks in bone marrow and genomic instability (MN-RET levels) in peripheral blood, between wild-type ( Trp53+/+ ) and heterozygous ( Trp53+/- ) mice. The dose response for Trp53+/+ mice showed higher initial levels of radiation-induced lymphocyte apoptosis relative to Trp53+/- between 0 and 1 Gy. Although this trend was observed up to 12 hours post-irradiation, both genotypes ultimately reached the same level of apoptosis at 14 hours, suggesting the importance of late-onset p53-independent apoptotic responses in this mouse model. Expected radiation-induced G1 cell cycle delay was observed in Trp53+/+ but not Trp53+/- . Although p53 has an important role in cancer risk, we have shown its influence on radiation dose response can be temporally variable. This research highlights the importance of caution when using haematopoietic endpoints as surrogates to extrapolate radiation-induced cancer risk estimation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

Recapturing Graphite-Based Fuel Element Technology for Nuclear Thermal Propulsion

ORNL is currently recapturing graphite based fuel forms for Nuclear Thermal Propulsion (NTP). This effort involves research and development on materials selection, extrusion, and coating processes to produce fuel elements representative of historical ROVER and NERVA fuel. Initially, lab scale specimens were fabricated using surrogate oxides to develop processing parameters that could be applied to full length NTP fuel elements. Progress toward understanding the effect of these processing parameters on surrogate fuel microstructure is presented.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2023

Evaluation of 134 Ce/ 134 La as a PET Imaging Theranostic Pair for 225 Ac α-Radiotherapeutics

225 Ac-targeted α-radiotherapy is a promising approach to treating malignancies, including prostate cancer. However, α-emitting isotopes are difficult to image because of low administered activities and a low fraction of suitable γ-emissions. The in vivo generator 134 Ce/ 134 La has been proposed as a potential PET imaging surrogate for the therapeutic nuclides 225 Ac and 227 Th. In this report, we detail efficient radiolabeling methods using the 225 Ac-chelators DOTA and MACROPA. These methods were applied to radiolabeling of prostate cancer imaging agents, including PSMA-617 and MACROPA-PEG 4 -YS5, for evaluation of their in vivo pharmacokinetic characteristics and comparison to the corresponding 225 Ac analogs. Methods: Radiolabeling was performed by mixing DOTA/MACROPA chelates with 134 Ce/ 134 La in NH 4 OAc, pH 8.0, at room temperature, and radiochemical yields were monitored by radio–thin-layer chromatography. In vivo biodistributions of 134 Ce-DOTA/MACROPA.NH 2 complexes were assayed through dynamic small-animal PET/CT imaging and ex vivo biodistribution studies over 1 h in healthy C57BL/6 mice, compared with free 134 CeCl 3 . In vivo, preclinical imaging of 134 Ce-PSMA-617 and 134 Ce-MACROPA-PEG 4 -YS5 was performed on 22Rv1 tumor–bearing male nu/nu-mice. Ex vivo biodistribution was performed for 134 Ce/ 225 Ac-MACROPA-PEG 4 -YS5 conjugates. Results: 134 Ce-MACROPA.NH 2 demonstrated near-quantitative labeling with 1:1 ligand-to-metal ratios at room temperature, whereas a 10:1 ligand-to-metal ratio and elevated temperatures were required for DOTA. Rapid urinary excretion and low liver and bone uptake were seen for 134 Ce/ 225 Ac-DOTA/MACROPA. NH 2 conjugates in comparison to free 134 CeCl 3 confirmed high in vivo stability. An interesting observation during the radiolabeling of tumor-targeting vectors PSMA-617 and MACROPA-PEG 4 -YS5—that the daughter 134 La was expelled from the chelate after the decay of parent 134 Ce—was confirmed through radio–thin-layer chromatography and reverse-phase high-performance liquid chromatography. Both conjugates, 134 Ce-PSMA-617 and 134 Ce-MACROPA-PEG 4 -YS5, displayed tumor uptake in 22Rv1 tumor–bearing mice. The ex vivo biodistribution of 134 Ce-MACROPA.NH 2 , 134 Ce-DOTA and 134 Ce-MACROPA-PEG 4 -YS5 corroborated well with the respective 225 Ac-conjugates. Conclusion: These results demonstrate the PET imaging potential for 134 Ce/ 134 La-labeled small-molecule and antibody agents. The similar 225 Ac and 134 Ce/ 134 La-chemical and pharmacokinetic characteristics suggest that the 134 Ce/ 134 La pair may act as a PET imaging surrogate for 225 Ac-based radioligand therapies.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Towards verifiable cancer digital twins: tissue level modeling protocol for precision medicine

Cancer exhibits substantial heterogeneity, manifesting as distinct morphological and molecular variations across tumors, which frequently undermines the efficacy of conventional oncological treatments. Developments in multiomics and sequencing technologies have paved the way for unraveling this heterogeneity. Nevertheless, the complexity of the data gathered from these methods cannot be fully interpreted through multimodal data analysis alone. Mathematical modeling plays a crucial role in delineating the underlying mechanisms to explain sources of heterogeneity using patient-specific data. Intra-tumoral diversity necessitates the development of precision oncology therapies utilizing multiphysics, multiscale mathematical models for cancer. This review discusses recent advancements in computational methodologies for precision oncology, highlighting the potential of cancer digital twins to enhance patient-specific decision-making in clinical settings. We review computational efforts in building patient-informed cellular and tissue-level models for cancer and propose a computational framework that utilizes agent-based modeling as an effective conduit to integrate cancer systems models that encode signaling at the cellular scale with digital twin models that predict tissue-level response in a tumor microenvironment customized to patient information. Furthermore, we discuss machine learning approaches to building surrogates for these complex mathematical models. These surrogates can potentially be used to conduct sensitivity analysis, verification, validation, and uncertainty quantification, which is especially important for tumor studies due to their dynamic nature.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Experimental and Computational Characterization of a Modified Sioutas Cascade Impactor for Respirable Radioactive Aerosols

Oak Ridge National Laboratory is collecting and characterizing aerosols released when spent nuclear fuel (SNF) rods are fractured in bending. An aerosol collection system was designed and tested to collect respirable sized (<10 μm aerodynamic diameter [AED]) particulates inside a hot cell facility. The setup is a modified version of the commercially available Sioutas cascade impactor, to which additional stages were added to expand the aerosol collection range from 2.5 to ~15 μm AED. To accommodate the additional stages and specific test conditions, the operating flow rate for aerosol collection was reduced, and testing was conducted by using pressure drop measurements, surrogate dust collection, and particle size characterization. The fluid flow distribution within the cascade and its stages was simulated in STAR-CCM+, and the stage-wise pressure drops obtained using the computational fluid dynamics model were then compared to experimental data. Lagrangian particle simulations were also performed, and stage-wise collection statistics were obtained from the simulation for comparison with the experimental data obtained using SNF-surrogate dust particles. The results provide valuable insights into the stage-wise particle collection characteristics of the modified cascade impactor and can also be used to improve the prediction accuracy of the manufacturer-determined analytical correlations.

aerosol modeling↗
DOE OSTI2018

Does Marine Surface Tension Have Global Biogeography? Addition for the OCEANFILMS Package

We apply principles of Gibbs phase plane chemistry across the entire ocean-atmosphere interface to investigate aerosol generation and geophysical transfer issues. Marine surface tension differences comprise a tangential pressure field controlling trace gas fluxes, primary organic inputs, and sea spray salt injections, in addition to heat and momentum fluxes. Mapping follows from the organic microlayer composition, now represented in ocean system models. Organic functional variations drive the microforcing, leading to (1) reduced turbulence and (by extension) laminar gas-energy diffusion; plus (2) altered bubble film mass emission into the boundary layer. Interfacial chemical behaviors are, therefore, closely reviewed as the background. We focus on phase transitions among two dimensional “solid, liquid, and gaseous” states serving as elasticity indicators. From the pool of dissolved organic carbon (DOC) only proteins and lipids appear to occupy significant atmospheric interfacial areas. The literature suggests albumin and stearic acid as the best proxies, and we distribute them through ecodynamic simulation. Consensus bulk distributions are obtained to control their adsorptive equilibria. We devise parameterizations for both the planar free energy and equation of state, relating excess coverage to the surface pressure and its modulus. Constant settings for the molecular surrogates are drawn from laboratory study and successfully reproduce surfactant solid-to-gas occurrence in compression experiments. Since DOC functionality measurements are rare, we group them into super-ecological province tables to verify aqueous concentration estimates. Outputs are then fed into a coverage, tension, elasticity code. The resulting two dimensional pressure contours cross a critical range for the regulation of precursor piston velocity, bubble breakage, and primary aerosol sources plus ripple damping. Concepts extend the water-air adsorption theory currently embodied in our OCEANFILMS aerosol emissions package, and the two approaches could be inserted into Earth System Models together. Uncertainties in the logic include kinetic and thermochemical factors operating at multiple scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Inverse Modeling of Hydrologic Parameters in CLM4 via Generalized Polynomial Chaos in the Bayesian Framework

In this work, generalized polynomial chaos (gPC) expansion for land surface model parameter estimation is evaluated. We perform inverse modeling and compute the posterior distribution of the critical hydrological parameters that are subject to great uncertainty in the Community Land Model (CLM) for a given value of the output LH. The unknown parameters include those that have been identified as the most influential factors on the simulations of surface and subsurface runoff, latent and sensible heat fluxes, and soil moisture in CLM4.0. We set up the inversion problem in the Bayesian framework in two steps: (i) building a surrogate model expressing the input–output mapping, and (ii) performing inverse modeling and computing the posterior distributions of the input parameters using observation data for a given value of the output LH. The development of the surrogate model is carried out with a Bayesian procedure based on the variable selection methods that use gPC expansions. Our approach accounts for bases selection uncertainty and quantifies the importance of the gPC terms, and, hence, all of the input parameters, via the associated posterior probabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Adaptive Machine Learning for Robust Diagnostics and Control of Time-Varying Particle Accelerator Components and Beams

Machine learning (ML) is growing in popularity for various particle accelerator applications including anomaly detection such as faulty beam position monitor or RF fault identification, for non-invasive diagnostics, and for creating surrogate models. ML methods such as neural networks (NN) are useful because they can learn input-output relationships in large complex systems based on large data sets. Once they are trained, methods such as NNs give instant predictions of complex phenomenon, which makes their use as surrogate models especially appealing for speeding up large parameter space searches which otherwise require computationally expensive simulations. However, quickly time varying systems are challenging for ML-based approaches because the actual system dynamics quickly drifts away from the description provided by any fixed data set, degrading the predictive power of any ML method, and limits their applicability for real time feedback control of quickly time-varying accelerator components and beams. In contrast to ML methods, adaptive model-independent feedback algorithms are by design robust to un-modeled changes and disturbances in dynamic systems, but are usually local in nature and susceptible to local extrema. In this work, we propose that the combination of adaptive feedback and machine learning, adaptive machine learning (AML), is a way to combine the global feature learning power of ML methods such as deep neural networks with the robustness of model-independent control. We present an overview of several ML and adaptive control methods, their strengths and limitations, and an overview of AML approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Validating a Tethered Balloon System and Optical Technologies for Marine Wildlife Detection and Tracking

The interactions between marine wildlife and marine energy devices are not well understood, leading to regulatory delays for device deployments and testing. Technologies that enable marine wildlife observations can help to fill data gaps and reduce uncertainties about animal–device interactions. A validation test conducted in Galveston Bay near La Porte, Texas, in December 2022 used a technology package consisting of a tethered balloon system and three independent sensor systems, including three-band visible, eight-band multispectral, and single-band thermal to detect three marine-mammal-shaped surrogates. The field campaign aimed to provide an initial step to evaluating the use of the TBS and the effectiveness of the sensor suite for marine wildlife observations and detection. From 2 December to 7 December 2022, 6 flights were conducted under varying altitudes and environmental conditions resulting in the collection of 5454 images. A subset of the images was classified and analyzed with two collection criteria including Beaufort wind force scale and TBS altitude to assess a range of observations of a surrogate from near-shore to offshore based on pixel count. The results of this validation test demonstrate the potential for using TBSs and imaging sensors for marine wildlife observations and offer valuable information for further development and application of this technology for marine energy and other blue economy sectors.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

THE MIRA–TITAN UNIVERSE: PRECISION PREDICTIONS FOR DARK ENERGY SURVEYS

Large-scale simulations of cosmic structure formation play an important role in interpreting cosmological observations at high precision. The simulations must cover a parameter range beyond the standard six cosmological parameters and need to be run at high mass and force resolution. A key simulation-based task is the generation of accurate theoretical predictions for observables using a finite number of simulation runs, via the method of emulation. Using a new sampling technique, we explore an eight-dimensional parameter space including massive neutrinos and a variable equation of state of dark energy. We construct trial emulators using two surrogate models (the linear power spectrum and an approximate halo mass function). The new sampling method allows us to build precision emulators from just 26 cosmological models and to systematically increase the emulator accuracy by adding new sets of simulations in a prescribed way. Emulator fidelity can now be continuously improved as new observational data sets become available and higher accuracy is required. Finally, using one ΛCDM cosmology as an example, we study the demands imposed on a simulation campaign to achieve the required statistics and accuracy when building emulators for investigations of dark energy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

StOKeDMD: Streaming Occupation kernel dynamic mode decomposition

Dynamic mode decomposition (DMD) has become a common technique for constructing surrogate models for dynamical systems from observed system states. The Occupation Kernel DMD (OKDMD) method proposed in (Rosenfeld et al., 2022) and (Rosenfeld et al., 2024) is a Liouville operator based method that builds surrogate models from system state trajectories. Here, this paper proposes an extension of OKDMD to the case when the system states are observed in a streaming fashion, i.e., only a small fraction of the state trajectory is available at a given time. The developed method, Streaming Occupation Kernel DMD (StOKeDMD), accommodates the streaming data input by leveraging properties of specific choices of kernel functions and occupation kernels. We apply the StoKeDMD method as a compression method for streaming data, analyze the memory complexity, and demonstrate the performance of StoKeDMD in the compression of streaming data generated from a Lorenz system and a fluid flow simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Optimal SCR Control Using Data-Driven Models

We present an optimal control solution for the urea injection for a heavy-duty diesel (HDD) selective catalytic reduction (SCR). The approach taken here is useful beyond SCR and could be applied to any system where a control strategy is desired and input-output data is available. For example, the strategy could also be used for the diesel oxidation catalyst (DOC) system. In this paper, we identify and validate a one-step ahead Kalman state-space estimator for downstream NOx using the bench reactor data of an SCR core sample. The test data was acquired using a 2010 Cummins 6.7L ISB production engine with a 2010 Cummins production aftertreatment system. We used a surrogate HDD federal test procedure (FTP), developed at Michigan Technological University (MTU), which simulates the representative transients of the standard FTP cycle, but has less engine speed/load points. The identified state-space model is then used to develop a tunable cost function that simultaneously minimizes NOx emissions and urea usage. The cost function is quadratic and univariate, thus the minimum can be computed analytically. We show the performance of the closed-loop controller in using a reduced-order discrete SCR simulator developed at MTU. Our experiments with the surrogate HDD-FTP data show that the strategy developed in this paper can be used to identify performance bounds for urea dose controllers.

diesel, emmissions, control, environment, modeling↗
DOE OSTI2018

Combustion Characteristics of PRF and TSF Ethanol Blends with RON 98 in an Instrumented CFR Engine

The CFR F1 engine is the standard testing apparatus used for rating the research octane number (RON) of gasoline fuels. Unlike the motor octane number (MON) method, where the intake port temperature after the carburetor is controlled by an electric heater, the mixture temperature can vary during the RON test due to the heat of vaporization (HoV) of the fuel. Ethanol is receiving increasing attention as a high octane and high HoV fuel component. This work presents an analysis of the combustion characteristics during the RON rating of ethanol fuel blends according to the standard ASTM D2699 method, highlighting the effects of ethanol concentration and base fuel composition. All fuels were blended to a constant RON of 98. Ethanol levels varied from 0 to 50 vol% and the base fuels were surrogate blends composed of primary reference fuels (PRF), toluene standardization fuels (TSF), and a four component gasoline surrogate. These were compared against two full boiling range gasolines, also having a RON of 98. Through the use of detailed cylinder pressure analyses, this paper provides insights into the combustion behavior of various RON 98 fuels which are commonly not captured by the standard RON knock rating method. Ethanol was found to significantly reduce the pressure transducer based knock intensity of PRF fuel blends up to a level of 30 vol%, despite all fuels having the same RON. Meanwhile, the knock intensity of the equivalent TSF fuel blends was found to be relatively insensitive to ethanol concentration. For ethanol concentrations beyond 30 vol%, the pressure transducer based knock intensity behavior did not appear to be affected by the base fuel composition.

CFR↗
DOE OSTI2020

Real Fuel Modeling for Gasoline Compression Ignition Engine

Increasing regulatory demand for efficiency has led to development of novel combustion modes such as HCCI, GCI and RCCI for gasoline light duty engines. In order to realize HCCI as a compression ignition combustion mode system, in-cylinder compression temperatures must be elevated to reach the autoignition point of the premixed fuel/air mixture. This should be co-optimized with appropriate fuel formulations that can autoignite at such temperatures. CFD combustion modeling is used to model the auto ignition of gasoline fuel under compression ignition conditions. Using the fully detailed fuel mechanism consisting of thousands of components in the CFD simulations is computationally expensive. To overcome this challenge, the real fuel is represented by few major components of create a surrogate fuel mechanism. In this study, 9 variations of gasoline fuel sets were chosen as candidates to run in HCCI combustion mode. A study detailing the development of the gasoline real fuel model was performed and various surrogates for gasoline fuel were investigated. The gasoline real fuel model will be used in subsequent CFD modelling activities for the development of an advanced mixed mode combustion system as part of the Department of Energy funded project DE-EE0008478.

Gasoline Compression Ignition, real fuel modeling,↗
DOE OSTI2020

Knock Mitigation Effectiveness of EGR across the Pressure-Temperature Domain

Exhaust gas recirculation (EGR) has been shown to enable efficiency improvements in SI engines through multiple different mechanisms, including decreasing the knock propensity at high load, which allows higher compression ratio. While many of the benefits of EGR are applicable to both low and high power density engines, including reductions in pumping work and improved specific heat ratio, the knock benefits and corresponding compression ratio increases have been limited to low power density naturally aspirated engines primarily intended for hybrid vehicle architectures. An earlier study [1] indicated that there may be a kinetic limitation for the ability of EGR to mitigate knock under these conditions, but that study only considered a small number of conditions. In this investigation, we expand on that study while also providing data for model validation for the new light-duty combustion consortium from the U.S. Department of Energy: Partnership for Advancing Combustion Engines (PACE). In this investigation, the effectiveness of EGR to mitigate knock is studied with regards to the effect of engine speed (1,500 and 3,000 rpm), changing trajectory in the pressure-temperature domain by varying the intake manifold temperature (35, 60, and 90 deg C), and by considering the effect of minor species by studying the effect of untreated EGR vs. EGR that has been treated by an automotive three-way catalyst. Additionally, to increase the relevance of these data for future modeling studies, the performance of the full boiling range gasoline was compared relative to a surrogate formulation. The study found that the fuel surrogate performs well, confirmed the kinetic limitations of EGR to mitigate knock under boost, and showed improvements in EGR performance with catalyzed EGR.

42 ENGINEERING↗
DOE OSTI2023

Multifidelity Monte Carlo estimation for efficient uncertainty quantification in climate-related modeling

Abstract. Uncertainties in an output of interest that depends on the solution of a complex system (e.g., of partial differential equations with random inputs) are often, if not nearly ubiquitously, determined in practice using Monte Carlo (MC) estimation. While simple to implement, MC estimation fails to provide reliable information about statistical quantities (such as the expected value of the output of interest) in application settings such as climate modeling, for which obtaining a single realization of the output of interest is a costly endeavor. Specifically, the dilemma encountered is that many samples of the output of interest have to be collected in order to obtain an MC estimator that has sufficient accuracy – so many, in fact, that the available computational budget is not large enough to effect the number of samples needed. To circumvent this dilemma, we consider using multifidelity Monte Carlo (MFMC) estimation which leverages the use of less costly and less accurate surrogate models (such as coarser grids, reduced-order models, simplified physics, and/or interpolants) to achieve, for the same computational budget, higher accuracy compared to that obtained by an MC estimator – or, looking at it another way, an MFMC estimator obtains the same accuracy as the MC estimator at lower computational cost. The key to the efficacy of MFMC estimation is the fact that most of the required computational budget is loaded onto the less costly surrogate models so that very few samples are taken of the more expensive model of interest. We first provide a more detailed discussion about the need to consider an alternative to MC estimation for uncertainty quantification. Subsequently, we present a review, in an abstract setting, of the MFMC approach along with its application to three climate-related benchmark problems as a proof-of-concept exercise.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A machine learning approach targeting parameter estimation for plant functional type coexistence modeling using ELM-FATES (v2.0)

Tropical forest dynamics play a crucial role in the global carbon, water, and energy cycles. However, realistically simulating the dynamics of competition and coexistence between different plant functional types (PFTs) in tropical forests remains a significant challenge. This study aims to improve the modeling of PFT coexistence in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), a vegetation demography model implemented in the Energy Exascale Earth System Model (E3SM) land model (ELM), ELM-FATES. Specifically, we explore (1) whether plant trait relationships established from field measurements can constrain ELM-FATES simulations and (2) whether machine learning (ML)-based surrogate models can emulate the complex ELM-FATES model and optimize parameter selections to improve PFT coexistence modeling. We conducted three ensembles of ELM-FATES experiments at a tropical forest site near Manaus, Brazil. By comparing the ensemble experiments without (Exp-CTR) and with (Exp-OBS) consideration of observed trait relationships, we found that accounting for these relationships slightly improves the simulations of water, energy, and carbon variables when compared to observations but degrades the simulation of PFT coexistence. Using ML-based surrogate models trained on Exp-CTR, we optimized the trait parameters in ELM-FATES and conducted another ensemble of experiments (Exp-ML) with these optimized parameters. The proportion of PFT coexistence experiments significantly increased from 21 % in Exp-CTR to 73 % in Exp-ML. After filtering the experiments that allow for PFT coexistence to agree with observations (within 15 % tolerance), 33 % of the Exp-ML experiments were retained, which is a significant improvement compared to the 1.4 % in Exp-CTR. Exp-ML also accurately reproduces the annual means and seasonal variations in water, energy, and carbon fluxes and the field inventory of aboveground biomass. This study represents a reproducible method that utilizes machine learning to identify parameter values that improve model fidelity against observations and PFT coexistence in vegetation demography models for diverse ecosystems. Our study also suggests the need for new mechanisms to enhance the robust simulation of coexisting plants in ELM-FATES and has significant implications for modeling the response and feedbacks of ecosystem dynamics to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Fortran–Python interface for integrating machine learning parameterization into earth system models

Abstract. Parameterizations in earth system models (ESMs) are subject to biases and uncertainties arising from subjective empirical assumptions and incomplete understanding of the underlying physical processes. Recently, the growing representational capability of machine learning (ML) in solving complex problems has spawned immense interests in climate science applications. Specifically, ML-based parameterizations have been developed to represent convection, radiation, and microphysics processes in ESMs by learning from observations or high-resolution simulations, which have the potential to improve the accuracies and alleviate the uncertainties. Previous works have developed some surrogate models for these processes using ML. These surrogate models need to be coupled with the dynamical core of ESMs to investigate the effectiveness and their performance in a coupled system. In this study, we present a novel Fortran–Python interface designed to seamlessly integrate ML parameterizations into ESMs. This interface showcases high versatility by supporting popular ML frameworks like PyTorch, TensorFlow, and scikit-learn. We demonstrate the interface's modularity and reusability through two cases: an ML trigger function for convection parameterization and an ML wildfire model. We conduct a comprehensive evaluation of memory usage and computational overhead resulting from the integration of Python codes into the Fortran ESMs. By leveraging this flexible interface, ML parameterizations can be effectively developed, tested, and integrated into ESMs.

54 ENVIRONMENTAL SCIENCES↗