Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “perturbed parameter ensemble”

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 37 records · Page 2

Uncertainty quantification and optimization of precipitating hydrometeor parameters for winter precipitation in a cloud microphysics scheme

The precipitating hydrometeor parameters used in cloud microphysics schemes carry inherent uncertainties. The quantification of these uncertainties, together with parameter optimization, can significantly improve precipitation forecasts. This study investigates the effects of 13 parameters in the Weather Research and Forecasting (WRF) Double-Moment 6-class (WDM6) microphysics scheme, which define the hydrometeor characteristics such as fall velocity–diameter and mass–diameter relationships, as well as the shape parameter of the drop size distribution for precipitating particles such as rain, snow, and graupel on simulated winter precipitation. A comparison between the model's pre-defined parameters and observations from the International Collaborative Experiments for the PyeongChang 2018 Olympic and Paralympic winter games (ICE-POP 2018) field campaign reveals that the fall velocity–diameter relationship for rain, the mass–diameter relationships for snow and graupel, and the shape parameters for all precipitating particles in the WDM6 scheme deviate from the median values observed by the two-dimensional video disdrometer (2DVD). To quantify parameter sensitivities, a perturbed parameter ensemble (PPE) of 256 simulations was conducted within parameter ranges constrained by 2DVD observations for three winter precipitation cases. Bayesian optimization was then applied to identify parameter sets that minimized the root mean square error (RMSE) for each case, achieving reductions of up to 30.2 %. These results demonstrate that ensemble-based uncertainty quantification and parameter optimization can help identify key parameters and provide a pathway to improving precipitation simulation performance. In addition, measurement sites can be strategically selected based on regions that show high sensitivity to variations in hydrometeor characteristic parameters.

Bayesian optimization↗

Emergent Relationships Among Sea Ice, Longwave Radiation, and the Beaufort High Circulation Exposed Through Parameter Uncertainty Analysis

Using a perturbed parameter ensemble of a coupled climate model, emerging relationships are identified between sea ice area, net surface longwave radiation, and the atmospheric circulation over the Beaufort gyre. There is a strong positive correlation between sea ice area and the net longwave radiation over the ocean-ice surface during the melting season and a negative correlation during the freezing season. The mechanisms responsible for the longwave radiation balance at the surface are mainly driven by sea ice variations in the freezing season and by clouds in the melting season. A strong positive (negative) correlation is also found between the fall (summer) total sea ice area in the Arctic and the sea level pressure over the Beaufort High region. It is argued that as sea ice coverage is lost, static stability losses are severe in fall, resulting in enhanced evaporation, vertical motions, and weakening of the general large-scale anticyclonic circulation of the Beaufort High.

54 ENVIRONMENTAL SCIENCES↗

Time Evolution and Diurnal Variability of the Parametric Sensitivity of Turbine-Height Winds in the MYNN-EDMF Parameterization

The Mellor-Yamada-Nakanishi-Niino (MYNN) parameterization applied in the Weather Research and Forecasting (WRF) model has been augmented to include the Eddy-Diffusion Mass-Flux (EDMF) approach to better represent transport by boundary-layer eddies. This change includes the addition of new parameters associated with convective updrafts and boundary-layer clouds. These additions lead to new parametric sensitivities in the hub-height wind speed compared to simulations using the standard MYNN parameterization that were explored previously. This work builds on efforts focused on WRF’s MYNN parameterization by examining the sensitivity of hub-height wind speed to parameters in the MYNN-EDMF parameterization as a function of simulation lead time. Summer and winter periods were selected from the second Wind Forecast Improvement Project (WFIP2). Five sets of simulations were completed for each season, with lead times ranging from 2 to 6 days. The results show that the sensitivity to the new parameters associated with the EDMF scheme is generally small compared to other parameters, with the largest impact in convective boundary layers. The spread in the perturbed parameter ensembles was found to grow quickly over the first 8 to 19 hours in the summer simulations and 17 to 24 hours in the winter simulations with little change after that, regardless of the simulation length. A strong diurnal cycle in the parameter sensitivity was also found associated with changes in the atmospheric stability, as well as an increase in the sensitivity to the entrainment parameter used in the EDMF parameterization that is associated with increasing fractional areas covered by plumes.

54 ENVIRONMENTAL SCIENCES↗

Matrix Approach to Accelerate Spin-Up of CLM5

Numerical models have been developed to investigate and understand responses of biogeochemical cycle to global changes. Steady state, when a system is in dynamic equilibrium, is generally required to initialize these model simulations. However, the spin-up process that is used to achieve steady state pose a great burden to computational resources, limiting the efficiency of global modeling analysis on biogeochemical cycles. This study introduces a new Semi-Analytical Spin-Up (SASU) to tackle this grand challenge. We applied SASU to Community Land Model version 5 and examined its computational efficiency and accuracy. At the Brazil site, SASU is computationally 7 times more efficient than (or saved up to 86% computational cost in comparison with) the traditional native dynamics (ND) spin-up to reach the same steady state. Globally, SASU is computationally 8 times more efficient than the accelerated decomposition spin-up and 50 times more efficient than ND. In summary, SASU achieves the highest computational efficiency for spin-up on site and globally in comparison with other spin-up methods. It is generalizable to wide biogeochemical models and thus makes computationally costly studies (e.g., parameter perturbation ensemble analysis and data assimilation) possible for a better understanding of biogeochemical cycle under climate change.

54 ENVIRONMENTAL SCIENCES↗

Buffering of Aerosol-Cloud Adjustments by Coupling Between Radiative Susceptibility and Precipitation Efficiency

Aerosol-cloud interactions (ACI) in warm clouds are the primary source of uncertainty in effective radiative forcing (ERF) during the historical period and, by extension, inferred climate sensitivity. The ERF due to ACI (ERFaci) is composed of the radiative forcing due to changes in cloud microphysics and cloud adjustments to microphysics. Here, we examine the processes that drive ERFaci using a perturbed parameter ensemble (PPE) hosted in CAM6. Observational constraints on the PPE result in substantial constraints in the response of cloud microphysics and macrophysics to anthropogenic aerosol, but only minimal constraint on ERFaci. Examination of cloud and radiation processes in the PPE reveal buffering of ERFaci by the interaction of precipitation efficiency and radiative susceptibility.

54 ENVIRONMENTAL SCIENCES↗

The Interaction Between Climate Forcing and Feedbacks

A Perturbed Parameter Ensemble (PPE) with the Community Atmosphere Model version 6 (CAM6) is used to better understand the sensitivity of aerosol forcing and cloud feedbacks to changes in model processes. Aerosol forcing through aerosol-cloud interactions is mostly negative (a cooling) due to shortwave radiation, while feedbacks are positive or negative in different regions due to contrasting longwave and shortwave effects. Both forcing and feedbacks are related to the mean climate state. Higher magnitude cloud radiative effects generally mean larger magnitude net negative forcing and larger magnitude net positive feedback. Aerosol forcing is broadly related to the susceptibility of clouds to drop number. Feedbacks also related to susceptibility, but to a lesser extent and in different regions to aerosol forcing. Aerosol forcing and cloud feedbacks are anti-correlated in the CAM6 PPE such that stronger negative forcing is associated with stronger positive feedbacks. Even the processes governing forcing and feedback sensitivity in the PPE are similar. These include the warm rain formation process, ice loss processes and deep convective intensity.

54 ENVIRONMENTAL SCIENCES↗

Process‐Oriented Calibration of a Turbulence Scheme in the DOE's Global Storm‐Resolving Model Using Machine Learning

A process‐oriented calibration framework is developed for the Simplified Higher‐Order Closure (SHOC) turbulence scheme in DOE's Simple Cloud Resolving E3SM Atmospheric Model (SCREAM). This framework leverages machine learning surrogates and observational constraints to efficiently calibrate SHOC adjustable parameters across two convective regimes: clear‐sky dry convective boundary layer and fair‐weather shallow cumulus clouds from ARM observations. We use perturbed‐parameter ensembles of a doubly periodic version of SCREAM to train surrogates and apply Markov Chain Monte Carlo sampling guided by cost functions based on benchmarking large‐eddy simulations and observations to identify optimized parameter sets that perform well in both regimes. The calibrated SHOC parameters substantially improve boundary‐layer turbulence and cloud boundaries, and modeled cloud fraction and radiative effects align better with observations than the default. These results demonstrate that combining multiple process‐specific convective regimes with machine‐learning surrogates can reduce parametric uncertainties and yield a model more faithful to cloud–turbulence interactions.

58 GEOSCIENCES↗

Sensitivity of Anthropogenic Cloud Droplet Number Change to Preindustrial Emission Inventories and Physics Parameterizations

Understanding the preindustrial (PI) to present-day (PD) change in cloud droplet number concentration (ΔNd) is important to constrain the anthropogenic influence on clouds and aid future model projections. Perturbed parameter ensemble simulations of PI and PD conditions in the Energy Exascale Earth System Model (E3SMv3) reveal two locations with consistent negative ΔNd signal across parameter space–the Southeastern US and the United Kingdom–despite a trend of positive ΔNd globally from increased fossil fuel consumption. These negative ΔNd signals are driven by higher anthropogenic biofuel emissions in the PI. By varying the PI aerosol mass and/or number emissions across these data sets' uncertainty ranges, we can change the magnitude and sign of ΔNd and the effective radiative forcing from aerosol-cloud interactions (ERFaci) regionally and globally, with near-zero ΔNd leading to near-zero ERFaci. These results highlight the need to constrain PI emission uncertainty to better understand Earth system responses to aerosol-cloud interactions.

Brown, Hunter Y. [University of Wyoming, Laramie, ↗

Matilda v1.0: An R package for probabilistic climate projections using a reduced complexity climate model

A primary advantage to using reduced complexity climate models (RCMs) has been their ability to quickly conduct probabilistic climate projections, a key component of uncertainty quantification in many impact studies and multisector systems. Providing frameworks for such analyses has been a target of several RCMs used in studies of the future co-evolution of the human and Earth systems. In this paper, we present Matilda, an open-science R software package that facilitates probabilistic climate projection analysis, implemented here using the Hector simple climate model in a seamless and easily applied framework. The primary goal of Matilda is to provide the user with a turn-key method to build parameter sets from literature-based prior distributions, run Hector iteratively to produce perturbed parameter ensembles (PPEs), weight ensembles for realism against observed historical climate data, and compute probabilistic projections for different climate variables. This workflow gives the user the ability to explore viable parameter space and propagate uncertainty to model ensembles with just a few lines of code. The package provides significant freedom to select different scoring criteria and algorithms to weight ensemble members, as well as the flexibility to implement custom criteria. Additionally, the architecture of the package simplifies the process of building and analyzing PPEs without requiring significant programming expertise, to accommodate diverse use cases. We present a case study that provides illustrative results of a probabilistic analysis of mean global surface temperature as an example of the software application.

54 ENVIRONMENTAL SCIENCES↗

Aerosol-Cloud Interactions and the Intersection of Climate Forcing and Feedback

The goal of this project has been to understand at the process level how the interaction of aerosols and cloud microphysics may alter the cloud response to climate change from the scale of turbulence to the global scale. Atmospheric aerosols (suspended particulates) have a unique role to play in the climate system. Anthropogenic aerosols may be an air pollution hazard and degrade human health. Their presence also results in the direct scattering or absorption of sunlight, with small net impacts on climate. But aerosols are also a key part of cloud formation, as the locations on which cloud droplets and ice crystals form. Understanding these Aerosol Cloud Interactions (ACI) is a critical uncertainty for climate models, including the Energy Exascale Earth System Model (E3SM), and affects how well models like E3SM can simulate the past and future. This work aims to understand cloud microphysics as well as the interaction of aerosol forcing and cloud feedback through improving our understanding of aerosol-cloud interactions (ACI) by looking at their uncertainty in critical environments at the process level. It will also contribute to the development of E3SM. Specific goals include: first, high level analysis of Perturbed Parameter Ensemble (PPE) simulations already available from the Community Earth System Model (CESM) and E3SM. Second, detailed analysis of recent changes to temperature and the role of aerosols.

54 ENVIRONMENTAL SCIENCES↗

Final report: Insights from ARM observations into aerosol processing and transport by extratropical cyclones and aerosol effects on cyclone clouds

This project advanced understanding of aerosol–cloud interactions in extratropical cyclones using observations from the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) observatory and ACE-ENA field campaigns, along with satellite datasets and numerical modeling. We have quantified the susceptibility of frontal clouds to aerosols, showing that this susceptibility is similar to that in non-frontal clouds despite substantial (factor 2) differences in droplet concentrations between these cloud types, and also significant differences in cloud albedo. A key finding was that Aitken-mode (here, 60-100nm diameter) aerosols play a disproportionately important role in droplet activation in frontal clouds, and concentrations of these particles are frequently strongly biased in climate models. Activation of aerosols depends strongly on updraft speeds, and via collaboration with Argonne National Laboratory we contributed to quantifying these better with the ARM radar wind profiler at the ENA site. Our studies in the North Atlantic have analogs in the Southern Ocean, and to draw these out we examined data from the MARCUS and SOCRATES field campaigns. We identified new particle formation within extratropical cyclones as a potentially important source of cloud condensation nuclei in these pristine marine environments. By combining ARM observations with perturbed parameter ensemble simulations, we demonstrated that surface observations can effectively constrain aerosol–cloud adjustments and reduce uncertainty in simulated cloud liquid water path responses by approximately 15%, yielding stronger constraints on historical aerosol cooling effects. The work further established precipitation processes as a dominant control on aerosol–cloud interactions and Earth system predictability.

Gordon, Hamish [Department of Chemical Engineering↗

Aircraft in-situ measurements from SOCRATES constrain the anthropogenic perturbations of cloud droplet number

Aerosol-cloud interactions (ACI) in warm clouds alter reflected shortwave radiation by influencing cloud microphysical and macrophysical properties. The variable of state controlling ACI is the cloud droplet number concentration (N d ). Here, we examine the perturbations in N d due to anthropogenic aerosols (ΔN d, PD-PI ) using a perturbed parameter ensemble (PPE) hosted in the sixth Community Atmosphere Model (CAM6). Surrogate models are created for the CAM6 PPE outputs and are used to generate 1 million model variants of CAM6 by sampling 45 sources of parameter uncertainty. The range of uncertain physical parameters related to ACI are constrained with observations of aerosol and cloud properties from SOCRATES. The likely range of uncertain parameters and the associated range of ΔN d, PD-PI are more strongly constrained with observations of N d relative to observations of cloud condensation nuclei. We conduct sensitivity tests of how constraints on ΔN d, PD-PI are affected by systematic uncertainties in observations and our limitations in our surrogate models created for CAM6 PPE outputs. Based on this, we provide guidance on the impact of reducing systematic uncertainty in airborne microphysical observations and in surrogate models.

Song, Ci [Univ. of Wyoming, Laramie, WY (United St↗

A machine learning approach to emulation and biophysical parameter estimation with the Community Land Model, version 5

Abstract. Land models are essential tools for understanding and predicting terrestrial processes and climate–carbon feedbacks in the Earth system, but uncertainties in their future projections are poorly understood. Improvements in physical process realism and the representation of human influence arguably make models more comparable to reality but also increase the degrees of freedom in model configuration, leading to increased parametric uncertainty in projections. In this work we design and implement a machine learning approach to globally calibrate a subset of the parameters of the Community Land Model, version 5 (CLM5) to observations of carbon and water fluxes. We focus on parameters controlling biophysical features such as surface energy balance, hydrology, and carbon uptake. We first use parameter sensitivity simulations and a combination of objective metrics including ranked global mean sensitivity to multiple output variables and non-overlapping spatial pattern responses between parameters to narrow the parameter space and determine a subset of important CLM5 biophysical parameters for further analysis. Using a perturbed parameter ensemble, we then train a series of artificial feed-forward neural networks to emulate CLM5 output given parameter values as input. We use annual mean globally aggregated spatial variability in carbon and water fluxes as our emulation and calibration targets. Validation and out-of-sample tests are used to assess the predictive skill of the networks, and we utilize permutation feature importance and partial dependence methods to better interpret the results. The trained networks are then used to estimate global optimal parameter values with greater computational efficiency than achieved by hand tuning efforts and increased spatial scale relative to previous studies optimizing at a single site. By developing this methodology, our framework can help quantify the contribution of parameter uncertainty to overall uncertainty in land model projections.

54 ENVIRONMENTAL SCIENCES↗

Evaluating uncertainty in aerosol forcing of tropical precipitation shifts

Abstract. An observed southward shift in tropical rainfall over land between 1950 and 1985, followed by a weaker recovery post-1985, has been attributed to anthropogenic aerosol radiative forcing and cooling of the Northern Hemisphere relative to the Southern Hemisphere. We might therefore expect models that have a strong historic hemispheric contrast in aerosol forcing to simulate a further northward tropical rainfall shift in the near-term future when anthropogenic aerosol emission reductions will predominantly warm the Northern Hemisphere. We investigate this paradigm using a perturbed parameter ensemble (PPE) of transient coupled ocean–atmosphere climate simulations that span a range of aerosol radiative forcing comparable to multi-model studies. In the 20th century, in our single-model ensemble, we find no relationship between the magnitude of pre-industrial to 1975 inter-hemispheric anthropogenic aerosol radiative forcing and tropical precipitation shifts. Instead, tropical precipitation shifts are associated with major volcanic eruptions and are strongly affected by internal variability. However, we do find a relationship between the magnitude of pre-industrial to 2005 inter-hemispheric anthropogenic aerosol radiative forcing and future tropical precipitation shifts over 2006 to 2060 under scenario RCP8.5. Our results suggest that projections of tropical precipitation shifts will be improved by reducing aerosol radiative forcing uncertainty, but predictive gains may be offset by temporary shifts in tropical precipitation caused by future major volcanic eruptions.

54 ENVIRONMENTAL SCIENCES↗

A Multiyear Gridded Data Ensemble of Surface Biogenic Carbon Fluxes for North America: Evaluation and Analysis of Results

Accurate and fine-scale estimates of biogenic carbon fluxes are critical for measuring and monitoring the biosphere's responses and feedback to the climate system. Currently available data products from flux towers and model-intercomparison projects struggle to adequately represent spatiotemporal dynamics of surface biogenic carbon fluxes, and to quantify their uncertainties, which also are crucial to atmospheric inversion systems. To address these gaps, we introduce a new perturbed-parameter model ensemble with the CASA model to estimate surface biogenic carbon fluxes at monthly and 3-hourly scales for North America at ~500-m and 5-km resolutions. We first use the Extended Fourier Amplitude Sensitivity Testing to choose the three most sensitive parameters to be perturbed, maximum light use efficiency (E max ), optimal temperature of photosynthesis (T opt ), and temperature response of respiration (Q 10 ). The initial range for each parameter is broadly sampled for the L1 ensemble, but then we pruned E max with site-level primary productivity to derive an L2 ensemble with narrower uncertainty ranges. Ensembles are strongly correlated with site-level results at both monthly and 3-hourly scales, and the spread across L1/L2 ensemble members encompasses the range of AmeriFlux observations. Monthly variability in the L2 ensemble mean is 85% of the observed variability. The L2 ensemble outperforms diverse data products with the highest Taylor skill scores at diurnal to annual scales. Finally, the ensemble's seasonality agrees well with other models for most biome types and in high and middle latitudes, but inconsistencies are found in subtropical and tropical ecoregions and for annual totals over North America.

54 ENVIRONMENTAL SCIENCES↗

A Process Model for ITCZ Narrowing under Warming Highlights Clear-Sky Water Vapor Feedbacks and Gross Moist Stability Changes in AMIP Models

Tropical areas with mean upward motion—and as such the zonal-mean intertropical convergence zone (ITCZ)—are projected to contract under global warming. To understand this process, a simple model based on dry static energy and moisture equations is introduced for zonally symmetric overturning driven by sea surface temperature (SST). Processes governing ascent area fraction and zonal mean precipitation are examined for insight into Atmospheric Model Intercomparison Project (AMIP) simulations. Bulk parameters governing radiative feedbacks and moist static energy transport in the simple model are estimated from the AMIP ensemble. Uniform warming in the simple model produces ascent area contraction and precipitation intensification—similar to observations and climate models. Contributing effects include stronger water vapor radiative feedbacks, weaker cloud-radiative feedbacks, stronger convection-circulation feedbacks, and greater poleward moisture export. The simple model identifies parameters consequential for the inter-AMIP-model spread; an ensemble generated by perturbing parameters governing shortwave water vapor feedbacks and gross moist stability changes under warming tracks inter-AMIP-model variations with a correlation coefficient ~0.46. Here, the simple model also predicts the multimodel mean changes in tropical ascent area and precipitation with reasonable accuracy. Furthermore, the simple model reproduces relationships among ascent area precipitation, ascent strength, and ascent area fraction observed in AMIP models. A substantial portion of the inter-AMIP-model spread is traced to the spread in how moist static energy and vertical velocity profiles change under warming, which in turn impact the gross moist stability in deep convective regions—highlighting the need for observational constraints on these quantities.

54 ENVIRONMENTAL SCIENCES↗

Data for Physiological Controls on Carbon Fluxes and Biomass Production in Miscanthus: Insights From a Process- Based Agroecosystem Model

Biomass crops serve as essential feedstocks for renewable energy and bioproducts and play a critical role in achieving lower emissions in the transportation sector. However, dedicated perennial biomass crops such as Miscanthus × giganteus (Miscanthus) remain underrepresented in process- based agroecosystem models, limiting robust evaluation of their economic and environmental performance. In this study, we developed a data- constrained representation of the sterile triploid Miscanthus (IL clone) within the process- based model ecosys, integrating global sensitivity analysis, ensemble simulation, and parameter calibration. Planting, harvesting, and fertilization practices consistent with field management were incorporated, and phenology was constrained using PhenoCam- derived Green Chromatic Coordinate (GCC) data. Using the Morris global sensitivity analysis method, we identified 11 key physiological parameters governing plant carbon, water, and nutrient relations, particularly processes associated with CO2 assimilation. We then conducted ensemble simulations by perturbing these parameters and calibrated the model against eddy covariance fluxes and field- measured biomass. Building on the calibrated operating state, parameter- response analyses show that different photosynthetic processes influence productivity in different ways. Protein allocation determines whether productivity increases toward a higher level, whereas electron transport capacity controls additional gains once protein allocation approaches saturation. These findings demonstrate that parameter importance depends on physiological context and on which photosynthetic processes remain limiting. Calibration and validation against observations show that ecosys can reliably reproduce carbon and water fluxes, as well as both above- and belowground biomass, with post- calibration GPP R2 improving from 0.67 to 0.95 during the calibration period and remaining high during validation (R2 = 0.95). These results provide a mechanistic foundation for regional simulations and sustainable bioenergy assessments.

Carbon↗

Physiological Controls on Carbon Fluxes and Biomass Production in Miscanthus: Insights From a Process‐Based Agroecosystem Model

Biomass crops serve as essential feedstocks for renewable energy and bioproducts and play a critical role in achieving lower emissions in the transportation sector. However, dedicated perennial biomass crops such as Miscanthus × giganteus (Miscanthus) remain underrepresented in process-based agroecosystem models, limiting robust evaluation of their economic and environmental performance. In this study, we developed a data-constrained representation of the sterile triploid Miscanthus (IL clone) within the process-based model ecosys, integrating global sensitivity analysis, ensemble simulation, and parameter calibration. Planting, harvesting, and fertilization practices consistent with field management were incorporated, and phenology was constrained using PhenoCam-derived Green Chromatic Coordinate (GCC) data. Using the Morris global sensitivity analysis method, we identified 11 key physiological parameters governing plant carbon, water, and nutrient relations, particularly processes associated with CO 2 assimilation. We then conducted ensemble simulations by perturbing these parameters and calibrated the model against eddy covariance fluxes and field-measured biomass. Building on the calibrated operating state, parameter-response analyses show that different photosynthetic processes influence productivity in different ways. Protein allocation determines whether productivity increases toward a higher level, whereas electron transport capacity controls additional gains once protein allocation approaches saturation. These findings demonstrate that parameter importance depends on physiological context and on which photosynthetic processes remain limiting. Calibration and validation against observations show that ecosys can reliably reproduce carbon and water fluxes, as well as both above- and belowground biomass, with post-calibration GPP R 2 improving from 0.67 to 0.95 during the calibration period and remaining high during validation (R 2 = 0.95). These results provide a mechanistic foundation for regional simulations and sustainable bioenergy assessments.

ecosys↗