Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Implementation of a comprehensive ice crystal formation parameterization for cirrus and mixed-phase clouds in the EMAC model (based on MESSy 2.53)

A comprehensive ice nucleation parameterization has been implemented in the global chemistry-climate model EMAC to improve the representation of ice crystal number concentrations (ICNCs). The parameterization of Barahona and Nenes (2009, hereafter BN09) allows for the treatment of ice nucleation taking into account the competition for water vapour between homogeneous and heterogeneous nucleation in cirrus clouds. Furthermore, the influence of chemically heterogeneous, polydisperse aerosols is considered by applying one of the multiple ice nucleating particle parameterizations which are included in BN09 to compute the heterogeneously formed ice crystals. BN09 has been modified in order to consider the pre-existing ice crystal effect and implemented to operate both in the cirrus and in the mixed-phase regimes. Compared to the standard EMAC parameterizations, BN09 produces fewer ice crystals in the upper troposphere but higher ICNCs in the middle troposphere, especially in the Northern Hemisphere where ice nucleating mineral dust particles are relatively abundant. Overall, ICNCs agree well with the observations, especially in cold cirrus clouds (at temperatures below 205 K), although they are underestimated between 200 and 220 K. As BN09 takes into account processes which were previously neglected by the standard version of the model, it is recommended for future EMAC simulations.

58 GEOSCIENCES↗
DOE OSTI2026

Prognostic simulations of mixed-phase clouds with model AC-1D v1.0: the impact of aerosol types and freezing parameterizations on ice crystal budgets

Mixed-phase clouds at high latitudes contribute to the uncertainty in predicting cloud feedbacks and climate sensitivity, mainly due to the complexity of microphysical processes that influence the partitioning between the supercooled liquid and ice phases, and hence, cloud radiative effects on regional scales. Particularly in Arctic mixed-phase clouds, the activation of ice-nucleating particles (INPs) from various aerosol populations remains a leading source of uncertainty. We developed an aerosol-cloud one-dimensional (AC-1D) model, which provides a novel framework to prognostically treat INP and ice crystal budgets while explicitly accounting for polydisperse and multicomponent aerosol that activate INPs following different freezing parameterizations. The AC-1D model is informed by large-eddy simulations to probe the impact of INP representation on predicted ice crystal number concentrations (N i ) and ice crystal budgets in mixed-phase Arctic stratus. We apply three immersion freezing (IMF) parameterizations, two time-independent (singular) and one time-dependent (classical nucleation theory), to predict the evolution of the INP reservoir and resulting ice crystal budget from polydisperse mineral dust, organic (humic-like substances), and sea spray aerosol particle size distributions. Our analysis focuses on how variations in aerosol number concentration and cloud system parameters such as cloud cooling rate, cloud-top entrainment rate, and ice crystal fall speed influence the INP reservoir and ice crystal budgets. Furthermore, this study investigates the competitive ice nucleation dynamics in mixed aerosol environments and provides a process-level quantification of the INP budget terms, which directly controls ice crystal budgets. For all studied case scenarios, the aerosol types and associated particle size distributions significantly impact INP and N i , and the choice between a time-dependent and a singular freezing description yields orders-of-magnitude differences in the predicted INP and N i over the 10 h simulation time, reflecting typical cloud lifetimes. Our results show that the influence of cloud cooling, INP entrainment, and sedimentation varies significantly depending on the chosen freezing parameterization. These findings underscore the critical need for robust IMF parameterizations and precise cloud system observations to enhance the accuracy of models in predicting mixed-phase cloud structure and evolution.

Arctic clouds↗
DOE OSTI2016

Global climate impacts of stochastic deep convection parameterization in the NCAR CAM5

In this paper, the stochastic deep convection parameterization of Plant and Craig (PC) is implemented in the Community Atmospheric Model version 5 (CAM5) to incorporate the stochastic processes of convection into the Zhang-McFarlane (ZM) deterministic deep convective scheme. Its impacts on deep convection, shallow convection, large-scale precipitation and associated dynamic and thermodynamic fields are investigated. Results show that with the introduction of the PC stochastic parameterization, deep convection is decreased while shallow convection is enhanced. The decrease in deep convection is mainly caused by the stochastic process and the spatial averaging of input quantities for the PC scheme. More detrained liquid water associated with more shallow convection leads to significant increase in liquid water and ice water paths, which increases large-scale precipitation in tropical regions. Specific humidity, relative humidity, zonal wind in the tropics, and precipitable water are all improved. The simulation of shortwave cloud forcing (SWCF) is also improved. The PC stochastic parameterization decreases the global mean SWCF from -52.25 W/m 2 in the standard CAM5 to -48.86 W/m 2 , close to -47.16 W/m 2 in observations. The improvement in SWCF over the tropics is due to decreased low cloud fraction simulated by the stochastic scheme. Sensitivity tests of tuning parameters are also performed to investigate the sensitivity of simulated climatology to uncertain parameters in the stochastic deep convection scheme.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Systematic Parameterization of Lignin for the CHARMM Force Field

Plant cell walls have three primary components, cellulose, hemicellulose, and lignin, the latter of which is a recalcitrant, aromatic heteropolymer that provides structure to plants, water and nutrient transport through plant tissues, and a highly effective defense against pathogens. Overcoming the recalcitrance of lignin is key to effective biomass deconstruction, which would in turn enable the use of biomass as a feedstock for industrial processes. Our understanding of lignin structure in the plant cell wall is hampered by the limitations of the available lignin forcefields, which currently only account for a single linkage between lignins and lack explicit parameterization for emerging lignin structures both from natural variants and engineered lignin structures. Since polymerization of lignin occurs via radical intermediates, multiple C-O and C-C linkages have been isolated , and the current force field only represents a small subset of lignin the diverse lignin structures found in plants. In order to take into account the wide range of lignin polymerization chemistries, monomers and dimer combinations of C-, H-, G-, and S-lignins as well as with hydroxycinnamic acid linkages were subjected to extensive quantum mechanical calculations to establish target data from which to build a complete molecular mechanics force field tuned specifically for diverse lignins. This was carried out in a GPU-accelerated global optimization process, whereby all molecules were parameterized simultaneously using the same internal parameter set. By parameterizing lignin specifically, we are able to more accurately represent the interactions and conformations of lignin monomers and dimers relative to a general force field. This new force field will enables computational researchers to study the effects of different linkages on the structure of lignin, as well as construct more accurate plant cell wall models based on observed statistical distributions of lignin that differ between disparate feedstocks, and guide further lignin engineering efforts.

09 BIOMASS FUELS↗
DOE OSTI2020

The Impacts of Horizontal Grid Spacing and Cumulus Parameterization on Subseasonal Prediction in a Global Convection-Permitting Model

Monthlong simulations targeting four Madden–Julian oscillation events made with several global model configurations are verified against observations to assess the roles of grid spacing and convective parameterization on the representation of tropical convection and midlatitude forecast skill. Specifically, the performance of a global convection-permitting model (CPM) configuration with a uniform 3-km mesh is compared to that of a global 15-km mesh with and without convective parameterization, and of a variable-resolution “channel” simulation using 3-km grid spacing only in the tropics with a scale-aware convection scheme. It is shown that global 3-km simulations produce realistic tropical precipitation statistics, except for an overall wet bias and delayed diurnal cycle. The channel simulation performs similarly, although with an unrealistically higher frequency of heavy rain. The 15-km simulations with and without cumulus schemes produce too much light and heavy tropical precipitation, respectively. Without convection parameterization, the 15-km global model produces unrealistically abundant, short-lived, and intense convection throughout the tropics. Only the global CPM configuration is able to capture eastward-propagating Madden–Julian oscillation events, and the 15-km runs favor stationary or westward-propagating convection organized at the planetary scale. The global 3-km CPM exhibits the highest extratropical forecast skill aloft and at the surface, particularly during week 3 of each hindcast. Although more cases are needed to confirm these results, this study highlights many potential benefits of using global CPMs for subseasonal forecasting. Furthermore, results show that alternatives to global convection-permitting resolution—using coarser or spatially variable resolution—feature compromises that may reduce their predictive performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

“Scrutinizing entrainment and mass flux closures in shallow cumulus parameterizations using cloud-radar observations and large-eddy simulation” (Final report for DOE grant)

Warm, shallow clouds are responsible for a substantial portion of Earth System Model (ESM) uncertainty and divergence among model climate projections. A substantial amount of this uncertainty is tied to the representation of convection in ESMs, in particular assumptions underlying the convective parameterizations. This study focuses on entrainment and entrainment closures in ESM parameterizations. The completed project was a joint observational and modeling study to observationally derive estimates of entrainment and to use high-resolution modeling to scrutinize entrainment and mass flux closures in shallow convection parameterizations.

58 GEOSCIENCES↗
DOE OSTI2018

Improvement in the Modeled Representation of North American Monsoon Precipitation Using a Modified Kain–Fritsch Convective Parameterization Scheme

A commonly noted problem in the simulation of warm season convection in the North American monsoon region has been the inability of atmospheric models at the meso-β scales (10 s to 100 s of kilometers) to simulate organized convection, principally mesoscale convective systems. With the use of convective parameterization, high precipitation biases in model simulations are typically observed over the peaks of mountain ranges. To address this issue, the Kain–Fritsch (KF) cumulus parameterization scheme has been modified with new diagnostic equations to compute the updraft velocity, the convective available potential energy closure assumption, and the convective trigger function. The scheme has been adapted for use in the Weather Research and Forecasting (WRF). A numerical weather prediction-type simulation is conducted for the North American Monsoon Experiment Intensive Observing Period 2 and a regional climate simulation is performed, by dynamically downscaling. In both of these applications, there are notable improvements in the WRF model-simulated precipitation due to the better representation of organized, propagating convection. The use of the modified KF scheme for atmospheric model simulations may provide a more computationally economical alternative to improve the representation of organized convection, as compared to convective-permitting simulations at the kilometer scale or a super-parameterization approach.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Sensitivity of U.S. summer precipitation to model resolution and convective parameterizations across gray zone resolutions: Sensitivity Across Gray Zone Resolutions

Simulating summer precipitation is a significant challenge for climate models that rely on cumulus parameterizations to represent moist convection processes. Motivated by recent advances in computing that support very high-resolution modeling, this study aims to systematically evaluate the effects of model resolution and convective parameterizations across the gray zone resolutions. Simulations usingthe Weather Research and Forecasting model were conducted at grid spacings of 36 km, 12 km, and 4 km fortwo summers over the conterminous U.S. The convection-permitting simulations at 4 km grid spacing are most skillful in reproducing the observed precipitation spatial distributions and diurnal variability. Notable differences are found between simulations with the traditional Kain-Fritsch (KF) and the scale-aware Grell-Freitas (GF) convection schemes, with the latter more skillful in capturing the nocturnal timing in the Great Plains and North American monsoon regions. The GF scheme also simulates a smoother transition from convective to large-scale precipitation as resolution increases, resulting in reduced sensitivity to model resolution compared to the KF scheme. Nonhydrostatic dynamics has a positive impact on precipitation over complex terrain even at 12 km and 36 km grid spacings. With nudging of the winds toward observations, we show that the conspicuous warm biases in the Southern Great Plains are related to precipitation biases induced by large-scale circulation biases, which are insensitive to model resolution. Overall, notable improvements in simulating summer rainfall and its diurnal variability through convection-permitting modeling and scale-aware parameterizations suggest promising venues for improving climate simulations of water cycle processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Mode connectivity in the loss landscape of parameterized quantum circuits

Variational training of parameterized quantum circuits (PQCs) underpins many workflows employed on near-term noisy intermediate scale quantum (NISQ) devices. It is a hybrid quantum-classical approach that minimizes an associated cost function in order to train a parameterized ansatz. In this work we adapt the qualitative loss landscape characterization for neural networks introduced in Goodfellow et al. (2014); Li et al. (2017) and tests for connectivity used in Draxler et al. (2018) to study the loss landscape features in PQC training. We present results for PQCs trained on a simple regression task, using the bilayer circuit ansatz, which consists of alternating layers of parameterized rotation gates and entangling gates. Multiple circuits are trained with 3 different batch gradient optimizers: stochastic gradient descent, the quantum natural gradient, and Adam. We identify large features in the landscape that can lead to faster convergence in training workflows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

EZFF: Python library for multi-objective parameterization and uncertainty quantification of interatomic forcefields for molecular dynamics

Parameterization of interatomic forcefields is a necessary first step in performing molecular dynamics simulations. This is a non-trivial global optimization problem involving quantification of multiple empirical variables against one or more properties. We present EZFF, a lightweight Python library for parameterization of several types of interatomic forcefields implemented in several molecular dynamics engines against multiple objectives using genetic-algorithm-based global optimization methods. The EZFF scheme provides unique functionality such as the parameterization of hybrid forcefields composed of multiple forcefield interactions as well as built-in quantification of uncertainty in forcefield parameters and can be easily extended to other forcefield functional forms as well as MD engines.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Prognostic Validation of a Neural Network Unified Physics Parameterization

Weather and climate models approximate diabatic and sub-grid-scale processes in terms of grid-scale variables using parameterizations. Current parameterizations are designed by humans based on physical understanding, observations, and process modeling. As a result, they are numerically efficient and interpretable, but potentially oversimplified. However, the advent of global high-resolution simulations and observations enables a more robust approach based on machine learning. In this letter, a neural network-based parameterization is trained using a near-global aqua-planet simulation with a 4-km resolution (NG-Aqua). The neural network predicts the apparent sources of heat and moisture averaged onto (160 km)2 grid boxes. A numerically stable scheme is obtained by minimizing the prediction error over multiple time steps rather than single one. In prognostic single-column model tests, this scheme matches both the fluctuations and equilibrium of NG-Aqua simulation better than the Community Atmosphere Model does.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Observational Relationship between Entrainment Rate and Environmental Relative Humidity and Implications for Convection Parameterization

Entrainment rate is a critical but highly uncertain quantity in convective parameterizations; especially, the effects of environmental relative humidity on entrainment rate are controversial, or even opposite, in different studies. Analysis of aircraft observations of cumuli from the RACORO and RICO field campaigns shows that entrainment rate is positively correlated with relative humidity. Physical analysis shows that higher relative humidity promotes entrainment by reducing buoyancy in the cloud cores and by weakening downdrafts near the cloud cores. The reduced buoyancy in the cloud cores and weakened downdrafts surrounding the cores further reduce updrafts in the cloud cores; the cloud cores with smaller updrafts are more significantly affected by their environment, resulting in larger entrainment rate. The relationship between entrainment rate and relative humidity is consistent with the buoyancy sorting concept widely used in convection parameterizations. The results provide reliable in situ observations to improve parameterizations of entrainment rate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Improving the Physical Basis for Updraft Dynamics in Deep Convection Parameterizations

This article presents a new deep convective parameterization that determines cloud characteristics based on a specified cloud size distribution. The vertical profiles of cloud properties are determined by analytic equations, which formulate entrainment with an inverse relationship to cloud width. In line with recent studies of large eddy simulations (LES), cloud widths are assumed to be constant with height and vertical mass flux (M) characteristics of the clouds are therefore regulated by the vertical velocity profile. The parameterization is configured to work with existing cloud base M closure formulations, with the closure predicting the total cloud area rather than the cloud base M directly. Analytic formula are also used to connect the vertical wind shear magnitude to the cloud size distribution, wherein larger shear magnitudes result in more numerous large updrafts than weaker shear magnitudes, which is in line with recent research results. The parameterization is compared against 10 deep convective LES with varying thermodynamic and vertical wind shear profiles. Results show dramatic improvements in the prediction of normalized M, detrainment, and the properties of detrained air over the existing Zhang and McFarlane (1995) scheme. In particular, the new model is able to correctly portray the transition from a bottom-heavy M profile in weakly sheared environments, to a top-heavy M profile in strongly sheared environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Development of Heterogeneous Ice Nucleation Rate Coefficient Parameterizations From Ambient Measurements

Ice nucleating particles (INPs) are a rare subset of particles that can have an outsized impact relative to their prevalence. To simulate INP variability, models require parameterizations for the most important sources of INPs. Most parameterizations in the literature were developed from laboratory experiments that used proxies for ambient particles. Whether these laboratory models effectively represent ambient aerosol is still not well understood and parameterizations developed from field measurements may be more atmospherically relevant. Expanding on previous work in which we used an ice chamber and a single-particle mass spectrometer (SPMS) to characterize INP composition, we calculate heterogeneous ice nucleation rate coefficients for ambient particles. We find good agreement between our dust and previous measurements of airborne dust. We also use a Monte Carlo approach to assess the relative contribution of each measurement to the uncertainty and find that the biggest source is the sampling efficiency of the SPMS.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Coupling Warm Rain With an Eddy Diffusivity/Mass Flux Parameterization: 2. Sensitivities and Comparison to Observations

Abstract This is the second part of the description and validation of the coupling of a warm‐rain microphysics scheme to a stochastic multi‐plume eddy‐diffusivity/mass‐flux parameterization. Part 1 provides the model description and initial validation of the parameterization in a single column model (SCM) against large eddy simulations. This manuscript provides additional comparisons against satellite observations of 500 unique simulation cases, an analysis of model parameter‐based sensitivities in those 500 cases, and an evaluation of the idealized LES case study paradigm for parameterization development. The SCM is forced by a weather reanalysis in the subtropical northeast Pacific region characterized by the stratocumulus‐to‐cumulus cloud transition. 18 scalar outputs are examined to identify the sensitivities to 22 of the model's parameters. The model is highly sensitive to parameters that determine the behavior of convective plumes, which collectively lead to the most sensitivity in 13 of the 18 outputs. Only a few parameters contribute sensitivities that are not highly correlated with those of other parameters. While forcing the model with two canonical field campaign cases results in similar picture of which parameters lead to high sensitivity in model outputs, the sensitivity magnitudes can be overestimated compared to when the model is forced by weather reanalyzes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Parameterizing Convective Organization Effects With a Moisture-PDF Approach in Climate Models: Concept and a Regional Case Simulation

We propose a parameterization scheme of convective organization effects based on a moisture-distribution approach, which can reflect aggregation of convective cells within a model grid as well as the interaction between convection and spatial heterogeneity in free-troposphere moisture. With this concept, convective cells in an aggregated state are surrounded by air that is moister than the grid-mean condition, which provides a “shielding” effect favorable for deeper convective updrafts. Such effects are represented quantitatively via utilizing a schematic diagram about the geometric interpretation of convective cluster, dry area, and their contact area, in which both the convective cluster size and dry area fraction increase during the development of convection. Our sensitivity analyses indicate that the new scheme performs well in capturing the mean precipitation features. Excluding the convective organization effects in the model leads to a considerable reduction in the simulated precipitation magnitude. Notably, there exists an inverse relationship between background mean moisture condition and sub-grid moisture variability for a given rain rate, suggesting that with a high spatial heterogeneity in free-troposphere moisture, convection can still develop and maintain its strength under relatively dry background conditions, consistent with the convection-permitting model simulation. Overall, our new parameterization of convective organization effects can successfully reproduce the relationship between precipitation and sub-grid moisture variability, which is a missing element in traditional convection parameterization schemes but important for the simulations of precipitation variability at various scales in climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Acceleration of the Parameterization of Unified Microphysics Across Scales (PUMAS) on the Graphics Processing Unit (GPU) With Directive-Based Methods

Cloud microphysics is one of the most time-consuming components in a climate model. In this study, we port the cloud microphysics parameterization in the Community Atmosphere Model (CAM), known as Parameterization of Unified Microphysics Across Scales (PUMAS), from CPU to GPU to seek a computational speedup. The directive-based methods (OpenACC and OpenMP target offload) are determined as the best fit specifically for our development practices, which enable a single version of source code to run either on the CPU or GPU, and yield a better portability and maintainability. Their performance is first examined in a PUMAS stand-alone kernel and the directive-based methods can outperform a CPU node as long as there is enough computational burden on the GPU. A consistent behavior is observed when we run PUMAS on the GPU in a practical CAM simulation. A 3.6× speedup of the PUMAS execution time, including data movement between CPU and GPU, is achieved at a coarse horizontal resolution (8 NVIDIA V100 GPUs against 36 Intel Skylake CPU cores). This speedup further increases up to 5.4× at a high resolution (24 NVIDIA V100 GPUs against 108 Intel Skylake CPU cores), which highlights the fact that GPU favors larger problem size. This study demonstrates that using GPU in a CAM simulation can save noticeable computational costs even with a small portion of code being GPU-enabled. Therefore, we are encouraged to port more parameterizations to GPU to take advantage of its computational benefit.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Integrating State Data Assimilation and Innovative Model Parameterization Reduces Simulated Carbon Uptake in the Arctic and Boreal Region

Model representation of carbon uptake and storage is essential for accurate projection of the response of the arctic-boreal zone to a rapidly changing climate. Land model estimates of LAI and aboveground biomass that can have a marked influence on model projections of carbon uptake and storage vary substantially in the arctic and boreal zone, making it challenging to correctly evaluate model estimates of Gross Primary Productivity (GPP). To understand and correct bias of LAI and aboveground biomass in the Community Land Model (CLM), we assimilated the 8-day Moderate Resolution Imaging Spectroradiometer (MODIS) LAI observation and a machine learning product of annual aboveground biomass into CLM using an Ensemble Adjustment Kalman Filter (EAKF) in an experimental region including Alaska and Western Canada. Assimilating LAI and aboveground biomass reduced these model estimates by 58% and 72%, respectively. The change of aboveground biomass was consistent with independent estimates of canopy top height at both regional and site levels. The International Land Model Benchmarking system assessment showed that data assimilation significantly improved CLM's performance in simulating the carbon and hydrological cycles, as well as in representing the functional relationships between LAI and other variables. Here, to further reduce the remaining bias in GPP after LAI bias correction, we re-parameterized CLM to account for low temperature suppression of photosynthesis. The LAI bias corrected model that included the new parameterization showed the best agreement with model benchmarks. Combining data assimilation with model parameterization provides a useful framework to assess photosynthetic processes in LSMs.

58 GEOSCIENCES↗