Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model parameterization”

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,261 records · Page 70

High Latitude Climate Modeling: ARM Takes Us Beyond Case Studies (Final Report)

The main thrust of this project was to devise a method by which the majority of North Slope of Alaska (NSA) meteorological and radiometric data, collected on a daily basis, could be used to evaluate and improve global climate model (GCM) simulations and their parameterizations, particularly for cloud microphysics. Although the standard ARM Program sensors for a less complete suite of instruments for cloud and aerosol studies than the instruments on an intensive field program such as the 2008 Indirect and Semi-Direct Aerosol Campaign (ISDAC), the advantage they offer lies in the long time base and large volume of data that covers a wide range of meteorological and climatological conditions. The challenge has been devising a method to interpret the NSA data in a practical way, so that a wide variety of meteorological conditions in all seasons can be examined with climate models. If successful, climate modelers would have a robust alternative to the usual “case study” approach (i.e., from intensive field programs only) for testing and evaluating their parameterizations’ performance. Understanding climate change on regional scales requires a broad scientific consideration of anthropogenic influences that goes beyond greenhouse gas emissions to also include aerosol-induced changes in cloud properties. For instance, it is now clear that on small scales, human-induced aerosol plumes can exert microclimatic radiative and hydrologic forcing that rivals that of greenhouse gas–forced warming. This project has made significant scientific progress by investigating what causes successive versions of climate models continue to exhibit errors in cloud amount, cloud microphysical and radiative properties, precipitation, and radiation balance, as compared with observations and, in particular, in Arctic regions. To find out what is going wrong, we have tested the models' cloud representation over the full range of meteorological conditions found in the Arctic using the ARM North Slope of Alaska (NSA) data.

54 ENVIRONMENTAL SCIENCES↗

Decadal Variability of Southeast US Rainfall in an Eddying Global Coupled Model

Abstract Ocean variability is a dominant source of remote rainfall predictability, but in many cases the physical mechanisms driving this predictability are not fully understood. This study examines how ocean mesoscales (i.e., the Gulf Stream SST front) affect decadal Southeast US (SEUS) rainfall, arguing that the local imprint of large‐scale teleconnections is sensitive to resolved mesoscale features. Based on global coupled model experiments with eddying and eddy‐parameterizing ocean, we find that a resolved Gulf Stream improves localized rainfall and remote circulation response in the SEUS. The eddying model generally improves the air‐sea interactions in the Gulf Stream and the North Atlantic Subtropical High that modulate SEUS rainfall over decadal timescales. The eddy‐parameterizing simulation fails to capture the sharp SST gradient associated with the Gulf Stream and overestimates the role of tropical Pacific SST anomalies in the SEUS rainfall.

54 ENVIRONMENTAL SCIENCES↗

One Stomatal Model to Rule Them All? Toward Improved Representation of Carbon and Water Exchange in Global Models

Abstract Stomatal conductance schemes that optimize with respect to photosynthetic and hydraulic functions have been proposed to address biases in land‐surface model (LSM) simulations during drought. However, systematic evaluations of both optimality‐based and alternative empirical formulations for coupling carbon and water fluxes are lacking. Here, we embed 12 empirical and optimization approaches within a LSM framework. We use theoretical model experiments to explore parameter identifiability and understand how model behaviors differ in response to abiotic changes. We also evaluate the models against leaf‐level observations of gas‐exchange and hydraulic variables, from xeric to wet forest/woody species spanning a mean annual precipitation range of 361–3,286 mm yr −1 . We find that models differ in how easily parameterized they are, due to: (a) poorly constrained optimality criteria (i.e., resulting in multiple solutions), (b) low influence parameters, (c) sensitivities to environmental drivers. In both the idealized experiments and compared to observations, sensitivities to variability in environmental drivers do not agree among models. Marked differences arise in sensitivities to soil moisture (soil water potential) and vapor pressure deficit. For example, stomatal closure rates at high vapor pressure deficit range between −45% and +70% of those observed. Although over half the new generation of stomatal schemes perform to a similar standard compared to observations of leaf‐gas exchange, two models do so through large biases in simulated leaf water potential (up to 11 MPa). Our results provide guidance for LSM development, by highlighting key areas in need for additional experimentation and theory, and by constraining currently viable stomatal hypotheses.

54 ENVIRONMENTAL SCIENCES↗

Uncertainty quantification and global sensitivity analysis of the Los Alamos sea ice model

Changes in the high-latitude climate system have the potential to affect global climate through feedbacks with the atmosphere and connections with midlatitudes. Sea ice and climate models used to understand these changes have uncertainties that need to be characterized and quantified. We present a quantitative way to assess uncertainty in complex computer models, which is a new approach in the analysis of sea ice models. We characterize parametric uncertainty in the Los Alamos sea ice model (CICE) in a standalone configuration and quantify the sensitivity of sea ice area, extent, and volume with respect to uncertainty in 39 individual model parameters. Unlike common sensitivity analyses conducted in previous studies where parameters are varied one at a time, this study uses a global variance-based approach in which Sobol' sequences are used to efficiently sample the full 39-dimensional parameter space. We implement a fast emulator of the sea ice model whose predictions of sea ice extent, area, and volume are used to compute the Sobol' sensitivity indices of the 39 parameters. Main effects and interactions among the most influential parameters are also estimated by a nonparametric regression technique based on generalized additive models. A ranking based on the sensitivity indices indicates that model predictions are most sensitive to snow parameters such as snow conductivity and grain size, and the drainage of melt ponds. Lastly, it is recommended that research be prioritized toward more accurately determining these most influential parameter values by observational studies or by improving parameterizations in the sea ice model.

54 ENVIRONMENTAL SCIENCES↗

Spatiotemporal 4D Whole-cell Modeling of a Minimal Autotroph Reveals Central Carbon Metabolism Regulated Locally by Protein Megacomplexes via Post-translational Modifications under Light Disturbance

Photosynthetic microorganisms rely on multiple pathways in central carbon metabolism to adapt to fluctuating light and energy availability across diel cycles. Mechanistic insight into the regulatory dynamics of this adaptation requires integrating processes spanning disparate timescales, from rapid redox-dependent post-translational modifications (PTMs) to slower changes in protein expression and metabolic pathway usage. To address this complexity beyond genome-based inference and traditional modeling, we develop a whole-cell four-dimensional (3D + time) model of the marine cyanobacterium Prochlorococcus marinus MED4 that explicitly represents the spatial organization of enzymatic and molecular processes in central carbon metabolism under light perturbation. We employ a perturbation-based research design to experimentally generate time-series, multi-omics measurements that provide molecular descriptors and cryo-ET derived 3D segmented volumes as constraints for this dynamic 4D framework. The integration of experiments and modeling across defined light regimes enables quantitative validation of system-level responses and forecasting under distinct light disturbances. We test the hypothesis that light-dependent redox PTMs regulating the structural assembly of a protein megacomplex, the “dark complex,” modulate metabolic flux at a conserved regulatory node of the Calvin–Benson cycle (CBC) in cyanobacteria. Our model shows that subcellular spatial organization buffers rapid light-induced changes in thylakoid reaction rates, which are followed by redox-PTM-mediated sequestration or release of CBC enzymes in the dark complex, ultimately impacting carbon fixation dynamics within carboxysomes. Comparison with an equivalently parameterized well-mixed stochastic model demonstrates that post-translational regulation not only buffers transcriptional noise and diffusion-driven fluctuations but also stabilizes phenotypic outcomes, underscoring the importance of spatial heterogeneity in phenotypic robustness. This ability to probe adaptive, spatiotemporally resolved mechanisms in photosynthetic machinery and central carbon metabolism addresses a critical gap in genotype-to-phenotype inference and expands modeling and design capabilities for understudied or genetically intractable autotrophs such as P. marinus MED4.

Johnson, Connah G.↗

Deep Learning Emulation of Subgrid-Scale Processes in Turbulent Shear Flows

Deep neural networks (DNNs) are developed from a data set obtained from the dynamic Smagorinsky model to emulate the subgrid-scale (SGS) viscosity (ν sgs ) and diffusivity (κ sgs ) for turbulent stratified shear flows encountered in the oceans and the atmosphere. These DNNs predict ν sgs and κ sgs from velocities, strain rates, and density gradients such that the evolution of the kinetic energy budget and density variance budget terms is similar to the corresponding values obtained from the original dynamic Smagorinsky model. These DNNs also compute ν sgs and κ sgs ~2–4 times quicker than the dynamic Smagorinsky model resulting in a ~2–2.5 times acceleration of the entire simulation. This study demonstrates the feasibility of deep learning in emulating the subgrid-scale (SGS) phenomenon in geophysical flows accurately in a cost-effective manner. In a broader perspective, deep learning-based surrogate models can present a promising alternative to the traditional parameterizations of the subgrid-scale processes in climate models.

58 GEOSCIENCES↗

Making Atmospheric Convective Parameterizations Obsolete with Machine Learning Emulation

Parameterizations of moist convection in atmospheric models are notoriously problematic, and while global cloud resolving models (GCRM) are often touted as the ultimate solution, the computational cost is a considerable hurdle to overcome. Machine learning emulation of GCRMs for predictive modelling can leverage the DOE’s computational resource investments and allow widespread use of GCRMs such that traditional parameterizations become obsolete for most applications.

54 ENVIRONMENTAL SCIENCES↗

The evolution of process-based hydrologic models: historical challenges and the collective quest for physical realism

The diversity in hydrologic models has historically led to great controversy on the correct approach to process-based hydrologic modeling, with debates centered on the adequacy of process parameterizations, data limitations and uncertainty, and computational constraints on model analysis. Here, we revisit key modeling challenges on requirements to (1) define suitable model equations, (2) define adequate model parameters, and (3) cope with limitations in computing power. We outline the historical modeling challenges, provide examples of modeling advances that address these challenges, and define outstanding research needs. We also illustrate how modeling advances have been made by groups using models of different type and complexity, and we argue for the need to more effectively use our diversity of modeling approaches in order to advance our collective quest for physically realistic hydrologic models.

54 ENVIRONMENTAL SCIENCES↗

Determining Catalytically Relevant Surfaces through Coverage-Dependent Lattice Gas Models: Carbon Adsorption on Fe(100)

Here, we have quantified the C–C lateral interactions on Fe(100) using a density functional theory (DFT)-parameterized lattice gas cluster expansion (LG CE) model trained using 265 unique configurations spanning a C coverage from 0 to 1 monolayer (ML). Our LG CE model shows high predictive accuracy with a leave-multiple-out cross-validation score of 10.2 and 16.6 meV/site for systems with and without the top two layers of Fe atoms fixed, respectively. Electronic ground-state structures identified from the lattice gas model (including the structures at 0 and 1 monolayers) were further used to generate ab initio phase diagrams under a range of temperatures and pressures. At low temperatures (<400 K), we found that the 1.0 monolayer structure is dominant, whereas at higher temperatures (>500 K), the 0.88 ML structure is most likely to form on the Fe surface. Interestingly, our model identified a c (2 × 2) ordered structure at 1/2 ML, which correlates well with previous DFT studies for carbon adsorption on iron surfaces and matches with the experimentally observed low-energy electron diffraction structure. Overall, the DFT-parameterized energies for the C/Fe system including effects of coverage and configurational space can further help in developing multiscale models for various heterogeneous reactions involving C–C and C–Fe interactions.

08 HYDROGEN↗

Non-Invasive Electrode Impedance Estimation for Optimized Charge Profile Parameterization of Lithium-Ion Batteries

This work presents a non-invasive method for parameterizing a physically motivated equivalent circuit model of lithium-ion batteries using operando electrochemical impedance spectroscopy and time-domain data. The proposed model consists exclusively of linear circuit elements, enabling computationally efficient simulation and real-time implementation on battery management system chips. By integrating frequency- and time-domain measurements, the model accurately estimates internal states such as the negative electrode potential, achieving a root mean square error of 12.3 mV during fast charging validation. Parameterization requires only rate tests with sinusoidal perturbations at three different ambient temperatures, making the approach experimentally accessible. The model reveals key insights into battery behavior, including rate-dependent overpotentials primarily governed by charge transfer kinetics at the positive electrode, and temperature-dependent impedance contributions from both charge transfer and solid-state diffusion processes. Validation using reference electrodes confirms the model’s ability to detect lithium plating onset and reproduce impedance behavior across a wide range of operating conditions. The approach enables in situ optimization of fast charging profiles and lays the foundation for future extensions incorporating aging effects and plating dynamics.

25 ENERGY STORAGE↗

Accelerating Radiation Computations for Dynamical Models With Targeted Machine Learning and Code Optimization

Abstract Atmospheric radiation is the main driver of weather and climate, yet due to a complicated absorption spectrum, the precise treatment of radiative transfer in numerical weather and climate models is computationally unfeasible. Radiation parameterizations need to maximize computational efficiency as well as accuracy, and for predicting the future climate many greenhouse gases need to be included. In this work, neural networks (NNs) were developed to replace the gas optics computations in a modern radiation scheme (RTE+RRTMGP) by using carefully constructed models and training data. The NNs, implemented in Fortran and utilizing BLAS for batched inference, are faster by a factor of 1–6, depending on the software and hardware platforms. We combined the accelerated gas optics with a refactored radiative transfer solver, resulting in clear‐sky longwave (shortwave) fluxes being 3.5 (1.8) faster to compute on an Intel platform. The accuracy, evaluated with benchmark line‐by‐line computations across a large range of atmospheric conditions, is very similar to the original scheme with errors in heating rates and top‐of‐atmosphere radiative forcings typically below 0.1 K day −1 and 0.5 W m −2 , respectively. These results show that targeted machine learning, code restructuring techniques, and the use of numerical libraries can yield material gains in efficiency while retaining accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Planck 2015 Constraints on the Non-flat XCDM Inflation Model

We examine the Planck 2015 cosmic microwave background (CMB) anisotropy data by using a physically consistent energy density inhomogeneity power spectrum generated by quantum-mechanical fluctuations during an early epoch of inflation in the non-flat XCDM model. Here dark energy is parameterized using a fluid with a negative equation of state parameter but with the speed of fluid acoustic inhomogeneities set to the speed of light. We find that the Planck 2015 data in conjunction with baryon acoustic oscillation distance measurements are reasonably well fit by a closed-XCDM model in which spatial curvature contributes a percent of the current cosmological energy density budget. In this model, the measured non-relativistic matter density parameter and Hubble constant are in good agreement with values determined using most other data. Depending on cosmological parameter values, the closed-XCDM model has reduced power, relative to the tilted, spatially flat ΛCDM case, and can partially alleviate the low multipole CMB temperature anisotropy deficit and can help partially reconcile the CMB anisotropy and weak lensing σ 8 constraints, at the expense of somewhat worsening the fit to higher multipole CMB temperature anisotropy data. Furthemore, the closed-XCDM inflation model does not seem to improve the agreement much, if at all, compared to the closed ΛCDM inflation case, even though it has one additional free parameter.

79 ASTRONOMY AND ASTROPHYSICS↗

Physics-Informed CoKriging: A Gaussian-Process-Regression-Based Multifidelity Method for Data-Model Convergence

In this work, we propose a new Gaussian process regression (GPR)-based multifidelity method: physics-informed CoKriging (CoPhIK). In CoKriging-based multifidelity methods, the quantities of interest are modeled as linear combinations of multiple parameterized stationary Gaussian processes (GPs), and the hyperparameters of these GPs are estimated from data via optimization. In CoPhIK, we construct a GP representing low-fidelity data using physics-informed Kriging (PhIK), and model the discrepancy between low- and highfidelity data using a parameterized GP with hyperparameters identified via optimization. Our approach reduces the cost of optimization for inferring hyperparameters by incorporating partial physical knowledge. We prove that the physical constraints in the form of deterministic linear operators are satisfied up to an error bound. Furthermore, we combine CoPhIK with a greedy active learning algorithm for guiding the selection of additional observation locations. The eciency and accuracy of CoPhIK are demonstrated for reconstructing the partially observed modified Branin function, reconstructing the sparsely observed state of a steady state heat transport problem, and learning a conservative tracer distribution from sparse tracer concentration measurements.

Yang, Xiu↗

Building Intelligent Cyberinfrastructure to Learn Iteratively from both Observations and Models for Understanding Watershed Dynamics

Focal Area(s): Predictive modeling through the use of AI techniques and AI-derived model components; the use of AI and other tools to design a prediction system comprising of a hierarchy of models (e.g., AI-driven model/component/parameterization selection). Science Challenge: Watershed processes, such as the fate and transport of sediment, carbon and nutrients across landscapes and their fluxes to water bodies (e.g., streams, rivers and lakes), have important implications for global and regional carbon and nutrient dynamics, biogeochemical functioning of terrestrial ecosystems, and soil functions. The magnitude of lateral surface/subsurface transport and fluxes of sediment, carbon and nutrients are key factors controlling the vulnerability of watersheds to climate extremes such as droughts, wildfires, and floods. Recent field observations and other scientific evidence suggest that the magnitudes of lateral transport and fluxes of sediment, carbon and nutrients are governed primarily by the spatial and vertical heterogeneity of landscape and soil properties and by pedogenic processes. However, the current generation of land surface and watershed models do not mechanistically couple the terrestrial and hydrologic systems, nor do they represent sufficiently the spatial and vertical heterogeneity of land surface and subsurface properties. On the other hand, increasing complexity of coupled watershed and land surface models requires more data to parameterize, calibrate and validate. Remote sensing (RS) provides a means to acquire spatial data and characterize their heterogeneity at the watershed scale, overcoming a major limitation associated with conventional point measurements. To improve the representation of land-surface and surface/subsurface process coupling and sub-grid heterogeneity in watershed models, it is essential to build our predictive understanding by learning from both the multi-scale multi-process modeling and diverse multi-scale data while leveraging powerful artificial intelligence (AI) techniques.

54 ENVIRONMENTAL SCIENCES↗

Improving the representation of shallow cumulus convection with the simplified-higher-order-closure–mass-flux (SHOC+MF v1.0) approach

Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.

54 ENVIRONMENTAL SCIENCES↗

The Impact of Model Resolution on the Southern Hemisphere in Community Climate System Model Version 4 Idealized Climate Simulations

Model resolution plays a large role in accurately simulating the Southern Hemisphere circulation in both the ocean and atmosphere. Resolving the mesoscale field is important as it has been shown to have a significant impact on the large-scale climate in eddy-rich regions, which are regions of large CO2 absorption. The presence of ocean and atmospheric mesoscale features can affect sea surface temperatures, the strength and location of storm tracks, and many other air-sea processes. Additionally, with an improvement in resolution, the eddy kinetic energy in the ocean can be expected to change considerably. The significance model resolution has on the Southern Hemisphere is examined using the Community Climate System Model, Version 4, eddy-parameterizing and eddy-resolving simulations. The CO2 concentrations and ozone levels are specified independently to better understand how the mesoscale field responds to extreme changes in external forcing and the resulting climate impacts. Overall, in the eddy-parameterizing simulations, the ozone forcing is found to be more important than the changes in CO2 concentrations. However, in the case of the eddy-resolving simulations, the CO2 concentrations are found to be more dominant, especially in eddy-rich regions. These results demonstrate the need for an increase in model resolution for climate prediction.

Daher, Houraa (ORCID:0000000200177346)↗

The Ice Particle and Aggregate Simulator (IPAS). Part III: Verification and Analysis of Ice–Aggregate and Aggregate–Aggregate Collection for Microphysical Parameterization

Abstract The Ice Particle and Aggregate Simulator (IPAS) is used to theoretically represent the aggregation process of ice crystals. Aggregates have a variety of formations based on initial ice particle size, shape, and falling orientation, all of which influence water phase partitioning. Aggregate dimensional properties and density changes are calculated for monomer–monomer (MON–MON), monomer–aggregate (MON–AGG), and aggregate–aggregate (AGG–AGG) collection to be used by ice-microphysical models for improvement in aggregation parameterizations. Aggregates are chosen from a database of 9 744 000 preformed combinations to be further collected (see Part II). AGG–AGG collection results in more extreme and a smaller range of aggregate aspect ratios than MON–AGG collection. A majority of aggregates are closer to prolate than oblate spheroids, regardless of collection type, except for quasi-horizontally oriented particles that have extreme aspect ratios to begin with. MON–AGG collection frequently results in an increase in density upon collection, whereas MON–MON and AGG–AGG collection almost always result in particle density decreases, with extreme reductions near 99% for MON–MON collection. MON–MON collection results in the greatest decreases in density but then quickly becomes unaffected by the addition of more monomers due to inherent size differences between monomers and aggregates. Finally, a holistic analysis to in situ observations of cloud particle images is presented. IPAS 2D aspect ratios surround a median value of 0.6 and closely follow that of previous studies while varying by no more than ≈12% on average from observed aggregates.

54 ENVIRONMENTAL SCIENCES↗