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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Prognostic Power of Extreme Rainfall Scaling Formulas Across Space and Time Scales

Some studies documenting changes in extreme precipitation use scaling formulas to approximate the large percentiles of the rainfall distribution from average dynamical and thermodynamical variables, called predictors. Here we instead assess the performance of these formulas as approximations to the rain rates in individual events. We evaluate the accuracies of the scaling relationships as functions of spatial and temporal scales by analyzing tropical rainfall in a superparameterized model. Relationships using full vertical profiles of the predictors are more accurate than those using their values at specific vertical levels because they better characterize the specific dynamics of each event. Both types of scaling relationships perform well over a range of length scales from 200 to 2,000 km and time scales from an hour to a week, and their precision is higher in the case of simulations with superparameterization than with parameterized convection. Uncertainties emerging from the local use of the scaling relationships suggest that they may only characterize the intensification of individual extremes for a warming of 4–5 K or larger. Finally, we argue that these formulas can be used to reconstruct the tail of the rainfall distribution directly from its predictors without prior information on P. While scalings have been used as diagnostic equations conditioned on the occurrence of extreme rainfall, they are actually able to mimic the prognostic behavior of climate model parameterizations on a variety of scales when estimating the intensity, frequency, and spatial patterns of extremes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Effects of downscaled high-resolution meteorological data on the PSCF identification of emission sources

The Potential Source Contribution Function (PSCF) model has been successfully used for identifying regions of emission source at a long distance in this study, the PSCF model relies on backward trajectories calculated by the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model. In this study, we investigated the impacts of grid resolution and Planetary Boundary Layer (PBL) parameterization (e.g., turbulent transport of pollutants) on the PSCF analysis. The Mellor-Yamada-Janjic (MYJ) and Yonsei University (YUS) parameterization schemes were selected to model the turbulent transport in the PBL within the Weather Research and Forecasting (WRF version 3.6) model. Two separate domain grid sizes (83 and 27 km) were chosen in the WRF downscaling in generating the wind data for driving the HYSPLIT calculation. The effects of grid size and PBL parameterization are important in incorporating the influ- ence of regional and local meteorological processes such as jet streaks, blocking patterns, Rossby waves, and terrain-induced convection on the transport of pollutants by a wind trajectory. We found high resolution PSCF did discover and locate source areas more precisely than that with lower resolution meteorological inputs. The lack of anticipated improvement could also be because a PBL scheme chosen to produce the WRF data was only a local parameterization and unable to faithfully duplicate the real atmosphere on a global scale. The MYJ scheme was able to replicate PSCF source identification by those using the Reanalysis and discover additional source areas that was not identified by the Reanalysis data. In conclusion, a potential benefit for using high-resolution wind data in the PSCF modeling is that it could discover new source location in addition to those identified by using the Reanalysis data input.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Predicting biomass of hyperdiverse and structurally complex central Amazonian forests - a virtual approach using extensive field data

Notice on corrigendum: This paper has a corresponding corrigendum published. Please read the corrigendum first. Old-growth forests are subject to substantial changes in structure and species composition due to the intensification of human activities, gradual climate change and extreme weather events. Trees store ca. 90 % of the total aboveground biomass (AGB) in tropical forests and precise tree biomass estimation models are crucial for management and conservation. In the central Amazon, predicting AGB at large spatial scales is a challenging task due to the heterogeneity of successional stages, high tree species diversity and inherent variations in tree allometry and architecture. We parameterized generic AGB estimation models applicable across species and a wide range of structural and compositional variation related to species sorting into height layers as well as frequent natural disturbances. We used 727 trees (diameter at breast height ≥ 5 cm) from 101 genera and at least 135 species harvested in a contiguous forest near Manaus, Brazil. Sampling from this data set we assembled six scenarios designed to span existing gradients in floristic composition and size distribution in order to select models that best predict AGB at the landscape level across successional gradients. We found that good individual tree model fits do not necessarily translate into reliable predictions of AGB at the landscape level. When predicting AGB (dry mass) over scenarios using our different models and an available pantropical model, we observed systematic biases ranging from -31 % (pantropical) to +39 %, with root-mean-square error (RMSE) values of up to 130 Mg ha -1 (pantropical). Our first and second best models had both low mean biases (0.8 and 3.9 %, respectively) and RMSE (9.4 and 18.6 Mg ha -1 ) when applied over scenarios. Predicting biomass correctly at the landscape level in hyperdiverse and structurally complex tropical forests, especially allowing good performance at the margins of data availability for model construction/calibration, requires the inclusion of predictors that express inherent variations in species architecture. The model of interest should comprise the floristic composition and size-distribution variability of the target forest, implying that even generic global or pantropical biomass estimation models can lead to strong biases. Reliable biomass assessments for the Amazon basin (i.e., secondary forests) still depend on the collection of allometric data at the local/regional scale and forest inventories including species-specific attributes, which are often unavailable or estimated imprecisely in most regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Corrigendum to "Predicting biomass of hyperdiverse and structurally complex central Amazonian forests — a virtual approach using extensive field data" Published in Biogeosciences, 13, 1553-1570, 2016

Old-growth forests are subject to substantial changes in structure and species composition due to the intensification of human activities, gradual climate change and extreme weather events. Trees store ca. 90% of the total aboveground biomass (AGB) in tropical forests and precise tree biomass estimation models are crucial for management and conservation. In the central Amazon, predicting AGB at large spatial scales is a challenging task due to the heterogeneity of successional stages, high tree species diversity and inherent variations in tree allometry and architecture. We parameterized generic AGB estimation models applicable across species and a wide range of structural and compositional variation related to species sorting into height layers as well as frequent natural disturbances. We used 727 trees (diameter at breast height ≥5 cm) from 101 genera and at least 135 species harvested in a contiguous forest near Manaus, Brazil. Sampling from this data set we assembled six scenarios designed to span existing gradients in floristic composition and size distribution in order to select models that best predict AGB at the landscape level across successional gradients. We found that good individual tree model fits do not necessarily translate into reliable predictions of AGB at the landscape level. We observed systematic biases ranging from -31% (pantropical) to +39 %, with root-mean-square error (RMSE) values of up to 130-Mg ha -1 (pantropical), when predicting AGB (dry mass) over scenarios using our different models and an available pantropical model. Our first and second best models had both low mean biases (0.8 and 3.9 %, respectively) and RMSE (9.4 and 18.6 Mg ha -1 ) when applied over scenarios. Predicting biomass correctly at the landscape level in hyperdiverse and structurally complex tropical forests, especially allowing good performance at the margins of data availability for model construction/calibration, requires the inclusion of predictors that express inherent variations in species architecture. Furthermore, the model of interest should comprise the floristic composition and size-distribution variability of the target forest, implying that even generic global or pantropical biomass estimation models can lead to strong biases. Reliable biomass assessments for the Amazon basin (i.e., secondary forests) still depend on the collection of allometric data at the local/regional scale and forest inventories including species-specific attributes, which are often unavailable or estimated imprecisely in most regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Assessment of Precipitating Marine Stratocumulus Clouds in the E3SMv1 Atmosphere Model: A Case Study from the ARM MAGIC Field Campaign

This paper presents a process-oriented evaluation of precipitating stratocumulus and its transition to cumulus in version 1 of the Energy Exascale Earth System Model (E3SMv1) using comprehensive case-study observations from a field campaign of the Atmospheric Radiation Measurement program (ARM). The E3SMv1 single-column model (SCM) of the marine boundary layer and its low clouds and precipitation are compared to observations including subcloud drizzle retrievals from a combination of Doppler radar and lidar backscatter measurements. Additionally, the SCM is also compared to a large-eddy simulation (LES) of the same case. The combination of advanced remote sensing observations and LES is a powerful framework to evaluate the physical parameterizations of large-scale models. Given the observed large-scale environment, the E3SMv1 SCM realistically represents the evolution of clouds and boundary layer structure during the stratocumulus-to-cumulus transition. The model well simulates the liquid water path and its diurnal cycle in the stratocumulus period as well as the two-layer vertical thermodynamic structure and lower cloud fraction in the transition period. E3SMv1’s success in simulating the cloud in the stratocumulus period permitted examination of its precipitation processes. Here problems were identified with E3SMv1 producing an unrealistically small subcloud precipitation fraction, an unrealistic double peak in the vertical profiles of precipitation mass, and drizzle that evaporates too close to the surface. Further model diagnostics determined that these unrealistic characteristics resulted from an overly long microphysics time step and an unrealistic parameterization of the precipitation fraction. These results imply that careful consideration of these issues is needed in order to better simulate precipitation processes in marine stratocumulus.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Building a predictive model for polycyclic aromatic hydrocarbon dosimetry in organotypically cultured human bronchial epithelial cells using benzo[ a ]pyrene

The airway epithelium is a primary route of exposure for inhaled toxicants, and organotypic culture models represent an important advancement for toxicity testing compared to simple in vitro models that may lack metabolic capability and multicellular structure/communication associated with the bronchial epithelium in vivo. A quantitative understanding of chemical dosimetry is key for interpreting and extrapolating study results; however, dosimetry is understudied in organotypic models limiting ability to predict toxicity. We developed a dosimetry model for primary human bronchial epithelial cells (HBECs) cultured at the air-liquid interface (ALI) using benzo[a]pyrene (BAP), a representative polycyclic aromatic hydrocarbon. Dose and time course evaluation of metabolite formation and enzyme activity and expression were utilized to parameterize a cellular dosimetry model to improve the utility of ALI-HBECs for assessing chemical risk. Dosimetry analysis demonstrated absorption of BAP into cells and an increase in Phase 1 and 2 metabolites over time that correlated with regulation of metabolizing enzymes. BAP was cleared from cells by 48 h after exposure, and the primary metabolites generated in ALI-HBECs were BAP-3-phenol, BAP-4,5-dihydrodiol, BAP-7,8-dihydrodiol, BAP-9,10-dihydrodiol, BAP-7,8,9,10-tetrol, BAP-3-phenol-glucuronide, BAP-4,5-dihydrodiol-glucuronide, and BAP-9,10-dihydrodiol-glucuronide. The resulting dosimetry model described BAP and 7,8-dihydrodiol toxicokinetics in ALI-HBECs and suggested active excretion of 7,8-dihydrodiol. Overall, this study demonstrates metabolic competency of ALI-HBECs for BAP metabolism, demonstrates the usefulness of complex in vitro systems for human-relevant toxicity data, and exhibits how in silico models can be utilized for understanding the dosimetry of test compounds to aid in in vitro to human extrapolation of toxicity data for risk assessments.

Benzo[a]pyrene↗
DOE OSTI2016

Final report of “A Detailed Study of the Physical Mechanisms Controlling CO2-Brine Capillary Trapping in the Subsurface” (University of Arizona, DE-SC0006696)

Carbon capture and storage (CCS) of carbon dioxide emissions generated by production or combustion of fossil fuels is a technologically viable means to reduce the build-up of CO2 in the atmosphere and oceans. Using advantages of scale and location, CCS is particularly suitable for large point sources near ubiquitous deep saline aquifers, depleted gas reservoirs, or at production reservoirs for enhanced oil recovery (EOR). In the BES-funded research project, Oregon State University (OSU) carried out capillary trapping experiments with proxy fluids that mimic the properties of the scCO2/brine system under ambient temperatures and pressures, and successfully developed a unique and novel x-ray compatible, high-pressure, elevated temperature setup to study the scCO2/brine system under challenging reservoir conditions. Both methodologies were applied to a variety of porous media, including synthetic (glass bead) and geologic (Bentheimer sandstone) materials. The University of Arizona (UA) developed pore-scale lattice Boltzmann (LB) models which are able to handle the experimental conditions for proxy fluids, as well as the scCO2/brine system, that are capable of simulating permeability in volumes of tens of millions of fluid elements. We reached the following summary findings (main institute indicated): 1. (OSU/UA) To understand capillary trapping in a multiphase fluid-porous medium system, the system must be analyzed from a pore-scale force balance perspective; trapping can be enhanced by manipulating wetting and nonwetting phase fluid properties. 2. (OSU) Pore-scale fluid connectivity and topology has a clear and direct effect on nonwetting phase capillary trapping efficiency. 3. (OSU) Rock type and flow regime also have a pronounced effects on capillary trapping. 4. (OSU/UA) There is a predictable relationship between NWP connectivity and NWP saturation, which allows for development of injection strategies that optimize trapping. The commonly used Land model (Land, 1968) does not predict amount of trapped NWP accurately. 5. (UA) There are ambiguities regarding the segmentation of large-volume gray-scale CT data into pore-volumes suitable for pore-scale modeling. Simulated permeabilities vary by three orders of magnitude and do not resemble observed values very well. Small-volume synchrotron-based CT data (such as produced by OSU) does not suffer significantly from segmentation ambiguities. 6. (UA) A standard properly parameterized Shan-Chen model LB model is useful for simulating porous media with proxy fluids as well as the scCO2/brine system and produces results that are consistent with tomographic observations. 7. (UA) A LB model with fluid-interactions defined by a (modified) Peng-Robinson Equation of State is able to handle the scCO2/brine system with variable solid phase wettability. This model is numerically stable at temperatures between 0 and 250 °C and pressures between 3 and 50 MPa, and produces appropriate densities above the critical point of CO2 and exhibits three-phase separation below. Based on above findings OSU and UA have proposed continued experimentation and pore-scale modeling of the scCO2/brine system. The reported research has extensively covered capillary trapping using proxy fluids, but due to limited beam-time availability we were unable to apply our high-pressure CO2 setup to sufficient variation in fluid properties, and initial scCO2 connectivity. New data will also allow us to test, calibrate and apply our LB models to reservoir conditions beyond those that are currently feasible experimentally. Such experiments and simulations will also allow us to provide information how suitable proxy fluids are for the scCO2/brine system. We believe it would be worthwhile to pursue the following new research questions: 1. What are the fundamental differences in the physics underlying capillary trapping at ambient vs. supercritical conditions? 2. Do newly developed pore-scale trapping interactions and relationships translate to continuum scales? A motivation for these questions was elaborated in “Capillary Trapping of Super-Critical CO2: Linking Pore and Continuum Scales to Verify new Relationships” that was submitted to DOE-BES in 2015.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2017

Battery state of health monitoring by estimation of the number of cyclable Li-ions

Herein this paper introduces a method to monitor battery state of health (SOH) by estimating the number of cyclable Li-ions, a health-relevant electrochemical variable. SOH monitoring is critical to battery management in balancing the trade-off between maximizing system performance and minimizing battery degradation. The decrease of cyclable Li-ions indicates the effect on the SOH of degradation mechanisms that consume cyclable Li-ions. The unavailability of the number of cyclable Li-ions through non-invasive measurements makes its estimation necessary for in-situ SOH monitoring. In this paper, the extended Kalman filter (EKF) is used to estimate the number of cyclable Li-ions as an unknown battery parameter. The single particle model (SPM), a simplified battery electrochemical model, is used as the model in the EKF to achieve a computational complexity suitable for on-line estimation. Simulations are performed under typical electric vehicle current trajectories using an example parameter set for a hybrid-electric-vehicle battery. In the simulations, the battery is represented by the Doyle–Fuller–Newman (DFN) model, an electrochemical model with higher fidelity than the SPM. To comply with the practice, instead of using the same parameters as the DFN model in the SPM, parameterization of the SPM is performed before estimation of the number of cyclable Li-ions. The simulations show high estimation accuracy of the number of cyclable Li-ions using the EKF, even with the structural and parametric differences between the DFN model and the SPM, under both the ideal conditions and various non-ideal conditions (i.e., SOC estimation error, additional modeling error, and measurement noise).

25 ENERGY STORAGE↗
DOE OSTI2020

Molecular models for creep in oriented polyethylene fibers

Highly oriented and crystalline polyetheylene (PE) fibers have a large failure stress under rapid tensile loading but exhibit significant creep at much smaller stresses that limits applications. A possible mechanism is slip of chains due to stress-enhanced, thermally activated nucleation of dislocations at chain ends in crystalline regions. Molecular dynamics simulations are used to parameterize a Frenkel–Kontorova model that provides analytic expressions for the limiting stress and activation energy for dislocation nucleation as a function of stress. Results from four commonly used hydrocarbon potentials are compared to show that the qualitative behavior is robust and estimate quantitative uncertainties. In all cases, the results can be described by an Eyring model with values of the zero-stress activation energy Ea0≈1.5 eV and activation volume V* ≈ 45 Å3 that are consistent with the experimental results for increasingly crystalline materials. The limiting yield stress is ∼8 GPa. These results suggest that activated dislocation nucleation at chain ends is an important mechanism for creep in highly oriented PE fibers.

O’Connor, Thomas C. (ORCID:0000000293930295)↗
DOE OSTI2019

Improved Spatiotemporal Representativeness and Bias Reduction of Satellite-Based Evapotranspiration Retrievals via Use of In Situ Meteorology and Constrained Canopy Surface Resistance

Evapotranspiration (ET) is a key component of the atmospheric and terrestrial water and energy budgets. Satellite-based vegetation index approaches have used remotely-sensed vegetation and reanalysis meteorological properties with surface energy balance models to estimate global ET (MOD16 ET). We reconstructed satellite retrievals using in situ meteorology (Argonne-ET) and evaluated them using a dense network of surface turbulent heat flux measurements. Argonne-ET resolves spatial heterogeneity of ET across the U.S. Southern Great Plains that is not well characterized by MOD16; MOD16 ET exhibits spatial autocorrelation (Moran's I (MI) significantly > 0 in May – Aug.), whereas in situ and Argonne-ET exhibit characteristics of a random spatial process (MI not significantly different from 0). The skill in resolving ET temporal variability is not significantly different between MOD16 and Argonne-ET (correlation coefficient = 0.75 and 0.72, respectively). However, the root mean square errors were significantly lower for Argonne-ET (36 W m -2 ) than MOD16 (43 W m -2 ), and MOD16 exhibits substantial bias in annual ET relative to in situ measurements (-38%). This is attributed to overestimation in the dry canopy surface resistance (rs) parameterization. Using rs constrained to the range of typical measured values, Argonne-ET substantially reduces the bias in the annual ET (+1%). The improved ET estimates are critical for regional water budget analyses. The methodology presented herein also demonstrates the ability to retrieve high temporally resolved (30-min; cf. 8-day MOD16) ET that can be used for development of processed-based diagnostics of model biases, and to elucidate avenues to improve ET model parameterizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Improving the particle dry deposition scheme in the CMAQ photochemical modeling system

Dry deposition of atmospheric aerosols in large-scale models is a critical, but highly uncertain, sink process with a strong dependence on particle size, meteorological conditions, and land surface properties. This study investigates the particle dry deposition scheme implemented in the standard Community Multiscale Air Quality (CMAQ) model v5.2.1, characterizes its underlying parameterized components with comparison to a similar scheme in a contemporary regional-scale model, and proposes two updated schemes that are then evaluated with available ambient particle deposition velocity (V d ) measurements. Both updated schemes reduce the surprisingly strong dependence of deposition velocity on the aerosol mode width, with one scheme further introducing a dependence on vegetation coverage that is broadly consistent with variability in observations between vegetated and non-vegetated surfaces. Compared to the base scheme, the updated scheme with vegetation dependence increases V d for submicron particles and decreases it for larger particles by an average of 37% and –66%, respectively. This scheme performs statistically better than the base scheme, reducing fractional biases by 56%–97% for vegetated land-use types and has roughly equivalent performance over water. The base and updated schemes are tested with three annual CMAQ (v5.2.1) simulations for the year 2011; predicted ambient aerosol concentrations are evaluated with routine monitoring network observations and predicted dry deposition fluxes are evaluated with data from the Clean Air Status and Trends Network (CASTNET). The updated scheme with vegetation dependence reduces negative fractional biases for PM 10 by 41% and positive fractional biases for PM 2.5 organic carbon by 15%. This scheme has been incorporated into the most recent publicly accessible versions of CMAQ (v5.3 and beyond) to replace the scheme used in previous versions of CMAQ (v4.5 through v5.2.1).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Using ARM Observations to Evaluate Process-Interactions in MCS Simulations Across Scales (Final Progress Report)

This project, funded by DOE Atmospheric System Research (DE-SC0020050), focused on improving the representation of mesoscale convective systems (MCSs) in numerical weather and climate models by leveraging high-resolution observations from the DOE Atmospheric Radiation Measurement (ARM) program. The research aimed to evaluate model sensitivities to grid spacing, microphysics, and planetary boundary layer (PBL) schemes, with a particular emphasis on improving convection parameterization for high-resolution modeling. Findings from this work highlight several key advancements. Model validation against ARM radar wind profiler data from the Southern Great Plains (SGP) and Manaus (MAO) sites revealed systematic biases in simulated convective mass flux profiles, leading to the development of an observationally constrained evaluation framework for diagnosing and improving model performance. Sensitivity analyses demonstrated that the representation of Amazonian MCSs was highly dependent on PBL scheme selection, while mid-latitude MCSs were more strongly influenced by microphysics parameterizations. A series of high-resolution WRF simulations, ranging from 4 km to 125 m grid spacing, provided insight into the behavior of convective drafts across scales. While updraft properties converged at sub-kilometer resolutions, biases in downdraft intensity persisted even at the finest resolution tested, emphasizing the need for further refinements in model physics. Additionally, comparisons of MCS vertical structures between mid-latitude and tropical environments revealed stronger updrafts and larger mass flux in mid-latitude MCSs, providing critical insights for improving climate model representations of storm-scale dynamics. The project’s findings have already contributed to advancing numerical modeling capabilities, particularly in WRF, MPAS, ICON, and DOE’s SCREAM model, by refining how convective processes are represented in high-resolution climate simulations. Results were disseminated through peer-reviewed publications, conference presentations, and ARM/ASR Research Highlights, engaging the broader scientific community. The project also provided valuable training opportunities for two postdoctoral researchers, who played central roles in model development, analysis, and dissemination of results. Their work contributed to several publications and conference presentations, helping prepare them for careers in atmospheric modeling. By improving the simulation of MCSs, this research directly supports the development of next-generation climate models capable of more accurately representing extreme precipitation and convective processes. The insights gained will inform future improvements in convective parameterization and guide the design of high-resolution weather and climate simulations, ultimately enhancing the reliability of climate projections and weather forecasts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Improved multifidelity Monte Carlo estimators based on normalizing flows and dimensionality reduction techniques

Here, we study the problem of multifidelity uncertainty propagation for computationally expensive models. In particular, we consider the general setting where the high-fidelity and low-fidelity models have a dissimilar parameterization both in terms of number of random inputs and their probability distributions, which can be either known in closed form or provided through samples. We derive novel multifidelity Monte Carlo estimators which rely on a shared subspace between the high-fidelity and low-fidelity models where the parameters follow the same probability distribution, i.e., a standard Gaussian. We build the shared space employing normalizing flows to map different probability distributions into a common one, together with linear and nonlinear dimensionality reduction techniques, active subspaces and autoencoders, respectively, which capture the subspaces where the models vary the most. We then compose the existing low-fidelity model with these transformations and construct modified models with an increased correlation with the high-fidelity model, which therefore yield multifidelity estimators with reduced variance. A series of numerical experiments illustrate the properties and advantages of our approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Sensitivity of Arctic Clouds to Ice Microphysical Processes in the NorESM2 Climate Model

Abstract Ice formation remains one of the most poorly represented microphysical processes in climate models. While primary ice production (PIP) parameterizations are known to have a large influence on the modeled cloud properties, the representation of secondary ice production (SIP) is incomplete and its corresponding impact is therefore largely unquantified. Furthermore, ice aggregation is another important process for the total cloud ice budget, which also remains largely unconstrained. In this study, we examine the impact of PIP, SIP, and ice aggregation on Arctic clouds, using the Norwegian Earth System Model, version 2 (NorESM2). Simulations with both prognostic and diagnostic PIP show that heterogeneous freezing alone cannot reproduce the observed cloud ice content. The implementation of missing SIP mechanisms (collisional breakup, drop shattering, and sublimation breakup) in NorESM2 improves the modeled ice properties, while improvements in liquid content occur only in simulations with prognostic PIP. However, results are sensitive to the description of collisional breakup. This mechanism, which dominates SIP in the examined conditions, is very sensitive to the treatment of the sublimation correction factor, a poorly constrained parameter that is included in the utilized parameterization. Finally, variations in ice aggregation treatment can also significantly impact cloud properties, mainly through their impact on collisional breakup efficiency. Overall, enhancement in ice production through the addition of SIP mechanisms and the reduction in ice aggregation (in line with radar observations of shallow Arctic clouds) result in enhanced cloud cover and decreased TOA radiation biases, compared to satellite measurements, especially during the cold months. Significance Statement Arctic clouds remain a large source of uncertainty in projections of the future climate due to the poor representation of the microphysical processes that govern their life cycle. Ice formation is among the least understood processes. While it is widely recognized that better constraints on primary ice production (PIP) are needed to improve existing parameterizations, we show that secondary ice production (SIP) and ice aggregation can have also a significant impact on ice number concentrations. Constraining ice formation through the addition of missing SIP mechanisms and reducing ice aggregation can improve the representation of the cloud macrophysical properties and enhance total cloud cover in the Arctic region, which in turn contributes to decreased TOA radiation biases in the cold months.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

A Finite-Time Ensemble Method for Mixed Layer Model Comparison

Here, this work evaluates the fidelity of various upper-ocean turbulence parameterizations subject to realistic monsoon forcing and presents a finite-time ensemble vector (EV) method to better manage the design and numerical principles of these parameterizations. The EV method emphasizes the dynamics of a turbulence closure multimodel ensemble and is applied to evaluate 10 different ocean surface boundary layer (OSBL) parameterizations within a single-column (SC) model against two boundary layer large-eddy simulations (LES). Both LES include realistic surface forcing, but one includes wind-driven shear turbulence only, while the other includes additional Stokes forcing through the wave-average equations that generate Langmuir turbulence. The finite-time EV framework focuses on what constitutes the local behavior of the mixed layer dynamical system and isolates the forcing and ocean state conditions where turbulence parameterizations most disagree. Identifying disagreement provides the potential to evaluate SC models comparatively against the LES. Observations collected during the 2018 monsoon onset in the Bay of Bengal provide a case study to evaluate models under realistic and variable forcing conditions. The case study results highlight two regimes where models disagree 1) during wind-driven deepening of the mixed layer and 2) under strong diurnal forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

A two-parameter criterion for classifying the explodability of massive stars by the neutrino-driven mechanism

Thus far, judging the fate of a massive star (either a neutron star [NS] or a black hole) solely by its structure prior to core collapse has been ambiguous. Our work and previous attempts find a nonmonotonic variation of successful and failed supernovae with zero-age main-sequence mass, for which no single structural parameter can serve as a good predictive measure. However, we identify two parameters computed from the pre-collapse structure of the progenitor, which in combination allow for a clear separation of exploding and nonexploding cases with only a few exceptions (∼1%–2.5%) in our set of 621 investigated stellar models. One parameter is M{sub 4}, defining the normalized enclosed mass for a dimensionless entropy per nucleon of s = 4, and the other is μ{sub 4}≡(dm/M{sub ⊙})/(dr/1000 km)|{sub s=4}, being the normalized mass derivative at this location. The two parameters μ{sub 4} and M{sub 4}μ{sub 4} can be directly linked to the mass-infall rate, M-dot , of the collapsing star and the electron-type neutrino luminosity of the accreting proto-NS, L{sub ν{sub e}}∝M{sub ns} M-dot , which play a crucial role in the “critical luminosity” concept for the theoretical description of neutrino-driven explosions as runaway phenomena of the stalled accretion shock. All models were evolved employing the approach of Ugliano et al. for simulating neutrino-driven explosions in spherical symmetry. The neutrino emission of the accretion layer is approximated by a gray transport solver, while the uncertain neutrino emission of the 1.1 M{sub ⊙} proto-NS core is parameterized by an analytic model. The free parameters connected to the core-boundary prescription are calibrated to reproduce the observables of SN 1987A for five different progenitor models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Improving Representation of Convective Transport for Scale-Aware Parameterization – Part I: Convection and Cloud Properties Simulated with Spectral Bin and Bulk Microphysics

The ultimate goal of this study is to improve representation of convective transport by cumulus parameterization for meso-scale and climate models. As Part I of the study, we perform extensive evaluations of cloud-resolving simulations of a squall line and mesoscale convective complexes in mid-latitude continent and tropical regions using the Weather Research and Forecasting (WRF) model with spectral-bin microphysics (SBM) and with two double-moment bulk microphysics schemes: a modified Morrison (MOR) and Milbrandt and Yau (MY2). Compared to observations, in general, SBM gives better simulations of precipitation, vertical velocity of convective cores, and the vertically decreasing trend of radar reflectivity than MOR and MY2, and therefore will be used for analysis of scale-dependence of eddy transport in Part II. The common features of the simulations for all convective systems are (1) the model tends to overestimate convection intensity in the middle and upper troposphere, but SBM can alleviate much of the overestimation and reproduce the observed convection intensity well; (2) the model greatly overestimates radar reflectivity in convective cores (SBM predicts smaller radar reflectivity but does not remove the large overestimation); and (3) the model performs better for mid-latitude convective systems than tropical system. The modeled mass fluxes of the mid latitude systems are not sensitive to microphysics schemes, but are very sensitive for the tropical case indicating strong microphysics modification to convection. Cloud microphysical measurements of rain, snow and graupel in convective cores will be critically important to further elucidate issues within cloud microphysics schemes.

Fan, Jiwen↗
DOE OSTI2025

A Decadal Hybrid GCM Simulation Using Deep‐Learning‐Based Cloud and Convection Parameterization Generalized to a Warm Climate

A critical challenge for machine‐learning (ML) parameterization in global climate models (GCMs) is to achieve stable, accurate simulations under climates not seen during training. Previous studies have demonstrated promising offline performance and year‐long online stability in aquaplanet simulations but have encountered difficulties in real geography and under climate warming. Here we report that a GCM with real geography configuration using neural‐network‐based cloud and convection parameterization, trained exclusively with present‐day climate data, successfully performs a stable, decade‐long simulation of a warm climate with +4 K sea surface temperature (SST). The neural network (NN) is based on Han et al. (2023, https://doi.org/10.1029/2022ms003508 ) with additional inputs. The simulation captures the global precipitation distribution, surface temperatures, vertical atmospheric structures, and extreme precipitation very well, closely matching simulations from both the superparameterized CAM (SPCAM) and the conventional CAM5 in the warm climate without accuracy degradation compared to those in the baseline climate. Moreover, it produces a climate response to +4 K SST in atmospheric thermodynamic states and circulations similar to those from SPCAM and CAM5. Prognostic ablation tests on NN input variables show that the NN without convective memory as input suffers from numerical instability, and the NN without considering radiative variables and land fraction as input, or with reduced training samples produce less accurate results. To our knowledge, this is the first time an ML parameterization successfully achieves online extrapolation to a warm climate without using additional warm‐climate data for training. It demonstrates the potential of ML‐driven parameterizations for credible long‐term climate projections.

Atmosphere model↗