SEARCH · Search NASA

Results for “ensemble simulation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Estimating the Likelihood of GHG Concentration Scenarios From Probabilistic Integrated Assessment Model Simulations

The climate scenarios that form the basis for current climate risk assessments have no assigned probabilities, and this impedes the analysis of future climate risks. This paper proposes an approach to estimate the probability of carbon dioxide (CO 2 ) concentration scenarios used in key climate change modeling experiments. It computes the CO 2 emissions compatible with the concentrations prescribed by Coupled Model Intercomparison Project Phase 5 (CMIP5) and CMIP6 experiments. The distribution of these compatible cumulative emissions is interpreted as the likelihood of future emissions given a concentration pathway. Using Bayesian analysis, the probability of each pathway can be estimated from a probabilistic sample of future emissions. The approach is demonstrated with five probabilistic CO 2 emission simulation ensembles from four Integrated Assessment Models (IAM), leading to independent estimates of the likelihood of the CO 2 concentration of Representative Concentration Pathways (RCP) and Shared Socioeconomic Pathways (SSP). Results suggest that SSP5-8.5 is unlikely for the second half of the 21st century, but offer no clear consensus on which of the remaining scenarios is most likely. Estimates of likelihoods of CO 2 concentrations associated with RCP and SSP scenarios are affected by sampling errors, differences in emission sources simulated by the IAMs, and a lack of a common experimental framework for IAM simulations. These shortcomings, along with a small IAM ensemble size, limit the applicability of the results presented here. Novel joint IAM and the Earth System Model experiments are needed to deliver actionable probabilistic climate risk assessments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

GrainNN: A neighbor-aware long short-term memory network for predicting microstructure evolution during polycrystalline grain formation

High fidelity simulations of grain formation in alloys are an indispensable tool for process-to-mechanical-properties characterization. Such simulations, however, can be computationally expensive as they require fine spatial and temporal discretizations. Their cost becomes an obstacle to parametric studies and ensemble runs and ultimately makes downstream tasks like optimal control and uncertainty quantification challenging. To enable such downstream tasks, we introduce GrainNN, an efficient and accurate reduced-order model for epitaxial grain growth in additive manufacturing conditions. GrainNN is a sequence-to-sequence long-short-term-memory (LSTM) deep neural network that evolves the dynamics of manually crafted features. Its innovations are (1) an attention mechanism with grain-microstructure-specific transformer architecture; and (2) an overlapping combination of several clones of the network to generalize to grain configurations that are different from those used for training. This design enables GrainNN to predict grain formation for unseen physical parameters, grain number, domain size and geometry. Furthermore, GrainNN not only reconstructs the quantities of interest but also can be pointwise accurate. In our numerical experiments, we use a polycrystalline phase field method to both generate the training data and assess GrainNN. For multiparametric, ensemble simulations with many grains, GrainNN can be orders of magnitude faster than phase field simulations, while delivering 5%–15% pointwise error. Additionally, this speedup includes the cost of the phase field simulations for generating training data.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Effect of Magnetic Field on the Dynamics of Movement of Charged Particles in Cluster Systems

The results of numerical and analytical studies of the influence of a constant magnetic field on the dynamics of the thermal motion of charged particles in an isotropic electrostatic trap are presented. Simulations were performed for cluster systems of particles with Coulomb interaction in a wide range of their parameters. The analysis of the spectral density of displacement of particles and mass transfer processes in simulated ensembles is performed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Estimating Changes in Temperature Distributions in a Large Ensemble of Climate Simulations Using Quantile Regression

Understanding future changes in extreme temperature events in a transient climate is inherently challenging. A single model simulation is generally insufficient to characterize the statistical properties of the evolving climate, but ensembles of repeated simulations with different initial conditions greatly expand the amount of data available. In this paper, we present a new approach for using ensembles to characterize changes in temperature distributions based on quantile regression that more flexibly characterizes seasonal changes. Specifically, our approach uses a continuous representation of seasonality rather than breaking the dataset into seasonal blocks; that is, we assume that temperature distributions evolve smoothly both day to day over an annual cycle and year to year over longer secular trends. To demonstrate our method’s utility, we analyze an ensemble of 50 simulations of the Community Earth System Model (CESM) under a scenario of increasing radiative forcing to 2100, focusing on North America. As previous studies have found, we see that daily temperature bulk variability generally decreases in wintertime in the continental mid- and high latitudes (>40°). A more subtle result that our approach uncovers is that differences in two low quantiles of wintertime temperatures do not shrink as much as the rest of the temperature distribution, producing a more negative skew in the overall distribution. Although the examples above concern temperature only, the technique is sufficiently general that it can be used to generate precise estimates of distributional changes in a broad range of climate variables by exploiting the power of ensembles.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Accelerated Probabilistic Marching Cubes by Deep Learning for Time-Varying Scalar Ensembles

Visualizing the uncertainty of ensemble simulations is challenging due to the large size and multivariate and temporal features of en-semble data sets. One popular approach to studying the uncertainty of ensembles is analyzing the positional uncertainty of the level sets. Probabilistic marching cubes is a technique that performs Monte Carlo sampling of multivariate Gaussian noise distributions for positional uncertainty visualization of level sets. However, the technique suffers from high computational time, making interactive visualization and analysis impossible to achieve. This paper introduces a deep-learning-based approach to learning the level-set uncertainty for two-dimensional ensemble data with a multivariate Gaussian noise assumption. We train the model using the first few time steps from time-varying ensemble data in our workflow. We demonstrate that our trained model accurately infers uncertainty in level sets for new time steps and is up to 170X faster than that of the original probabilistic model with serial computation and 10X faster than that of the original parallel computation.

Han, Mengjiao↗

NASA NTRSDate not supplied

The Role of Internal Variability and Feedbacks Controlling AMOC Stability

A bi-stable mode of the Atlantic Meridional Overturning Circulation (AMOC) is found in a 10-member ensemble simulation of the SSP2-4.5 scenario using the NASA GISS-E2-1-G climate model. Local feedbacks in the subpolar North Atlantic region in conjunction with internal variability in sea-ice transport and melt play a critical role in causing the divergent behavior of the AMOC in the ensemble members. While other fully coupled models have demonstrated the important role of surface freshening in leading to AMOC shutdown, either through hosing experiments or increased precipitation and greenhouse gas warming at high latitudes, in the GISS simulations, there are no external freshwater perturbations. This is the first time that a CMIP-class model has shown such a bifurcation across an initial condition ensemble.

Atlantic Meridional Overturning Circulation↗

DOE OSTI2023

Ice Supersaturation Variability in Cirrus Clouds: Role of Vertical Wind Speeds and Deposition Coefficients

Abstract Aircraft measurements reveal ice supersaturation statistics in cirrus (ISSs) with broad maxima around ice saturation and pronounced variance. In this study, processes shaping ISSs in midlatitude and tropical upper tropospheric conditions are systematically investigated. Water vapor deposition and sublimation of size‐resolved ice crystal populations are simulated in an air parcel framework. Mesoscale temperature fluctuations (MTFs) due to gravity waves force the temporal evolution of supersaturation. Various levels of background wave forcing and cirrus thickness are distinguished in stochastic ensemble simulations. Kinetic limitations to ice mass growth are brought about by supersaturation‐dependent deposition coefficients that represent efficient and inefficient growth modes as a function of ice crystal size. The simulations identify a wide range of deposition coefficients in cirrus, but most values stay above 0.01 such that kinetic limitations to water uptake remain moderate. Supersaturation quenching times are long, typically 0.5–2 hr. The wave forcing thus causes a remarkably large variability in ISSs and cirrus microphysical properties except in the thickest cirrus, producing ensemble‐mean ISSs in line with in‐situ measurements. ISS variance is controlled by MTFs and increases with decreasing cirrus integral radii. In comparison, the impact of ice crystal growth rates on ISSs is small. These results contribute to efforts directed at identifying and solving issues associated with ice‐supersaturated areas and non‐equilibrium cirrus physics in global models.

Meteorology & Atmospheric Sciences↗

DOE OSTI2024

Probabilistic projections of the Amery Ice Shelf catchment, Antarctica, under conditions of high ice-shelf basal melt

Abstract. Antarctica's Lambert Glacier drains about one-sixth of the ice from the East Antarctic Ice Sheet and is considered stable due to the strong buttressing provided by the Amery Ice Shelf. While previous projections of the sea-level contribution from this sector of the ice sheet have predicted significant mass loss only with near-complete removal of the ice shelf, the ocean warming necessary for this was deemed unlikely. Recent climate projections through 2300 indicate that sufficient ocean warming is a distinct possibility after 2100. This work explores the impact of parametric uncertainty on projections of the response of the Lambert–Amery system (hereafter “the Amery sector”) to abrupt ocean warming through Bayesian calibration of a perturbed-parameter ice-sheet model ensemble. We address the computational cost of uncertainty quantification for ice-sheet model projections via statistical emulation, which employs surrogate models for fast and inexpensive parameter space exploration while retaining critical features of the high-fidelity simulations. To this end, we build Gaussian process (GP) emulators from simulations of the Amery sector at a medium resolution (4–20 km mesh) using the Model for Prediction Across Scales (MPAS)-Albany Land Ice (MALI) model. We consider six input parameters that control basal friction, ice stiffness, calving, and ice-shelf basal melting. From these, we generate 200 perturbed input parameter initializations using space filling Sobol sampling. For our end-to-end probabilistic modeling workflow, we first train emulators on the simulation ensemble and then calibrate the input parameters using observations of the mass balance, grounding line movement, and calving front movement with priors assigned via expert knowledge. Next, we use MALI to project a subset of simulations to 2300 using ocean and atmosphere forcings from a climate model for both low- and high-greenhouse-gas-emission scenarios. From these simulation outputs, we build multivariate emulators by combining GP regression with principal component dimension reduction to emulate multivariate sea-level contribution time series data from the MALI simulations. We then use these emulators to propagate uncertainty from model input parameters to predictions of glacier mass loss through 2300, demonstrating that the calibrated posterior distributions have both greater mass loss and reduced variance compared to the uncalibrated prior distributions. Parametric uncertainty is large enough through about 2130 that the two projections under different emission scenarios are indistinguishable from one another. However, after rapid ocean warming in the first half of the 22nd century, the projections become statistically distinct within decades. Overall, this study demonstrates an efficient Bayesian calibration and uncertainty propagation workflow for ice-sheet model projections and identifies the potential for large sea-level rise contributions from the Amery sector of the Antarctic Ice Sheet after 2100 under high-greenhouse-gas-emission scenarios.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

On the additivity of climate responses to the volcanic and solar forcing in the early 19th century

The early 19th century was the coldest period over the past 500 years, when strong tropical volcanic events and a solar minimum coincided. The 1809 unidentified eruption and the 1815 Tambora eruption happened consecutively during the Dalton minimum of solar irradiance; however, the relative role of the two forcing (volcano and solar) agents is still unclear. In this study, we examine the responses from a set of early 19th century simulations with combined and separated volcanic and solar forcing agents, as suggested in the protocol for the past1000 experiment of the Paleoclimate Modelling Intercomparison Project – Phase 4 (PMIP4). From 20-member ensemble simulations with the Max Planck Institute Earth system model (MPI-ESM1.2-LR), we find that the volcano- and solar-induced surface cooling is additive in the global mean/large scale, regardless of combining or separating the forcing agents. The two solar reconstructions (SATIRE (Spectral and Total Irradiance REconstruction-Millennia model) and PMOD (Physikalisch-Meteorologisches Observatorium Davos)) contribute to a cooling before and after 1815 of ~0.05 and ~0.15 K monthly average near-surface air cooling, respectively, indicating a limited solar contribution to the early 19th century cold period. The volcanic events provide the main cooling contributions, inducing a surface cooling that peaks at ~0.82 K for the 1809 event and ~1.35 K for Tambora. After the Tambora eruption, the temperature in most regions increases toward climatology largely within 5 years, along with the reduction of volcanic forcing. In the northern extratropical oceans, the temperature increases slowly at a constant rate until 1830, which is related to the reduction of seasonality and the concurrent changes in Arctic sea-ice extent. The albedo feedback of Arctic sea ice is found to be the main contributor to the Arctic amplification of the cooling signal. Several non-additive responses to solar and volcanic forcing happen on regional scales. In the atmosphere, the stratospheric polar vortex tends to strengthen when combining both volcano and solar forcing, even though the two forcing agents separately induce opposite-sign changes in stratospheric temperatures and zonal winds. In the ocean, when combining the two forcings, additional surface cold water propagates to the northern extratropics from the additional solar cooling in the tropics, which results in regional cooling along the propagation. Overall, this study not only quantifies the surface responses from combinations of the volcano and solar forcing, but also highlights the components that cannot be simply added from the responses of the individual forcing agents, indicating that a relatively small forcing agent (such as solar in early 19th century) can impact the response from the large forcing (such as the 1815 Tambora eruption) when considering regional climates.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Transfer learning of neural surrogates on multifidelity groundwater simulations

Multifidelity data used in the paper published in Advances in Water Resources 206 (2025) 105140, https://doi.org/10.1016/j.advwatres.2025.105140 The code used to process the data is openly available on GitHub at https://github.com/Model-Reduction-and-UQ-Group/Transfer_Learning_K_reconstruction Computationally inexpensive surrogates of process-based models, such as deep neural networks, enable ensemble-based computations used in risk assessment, data assimilation, etc. However, generation of large datasets required to train a neural network can be as expensive as the ensemble simulations themselves. We ameliorate this challenge by using data from multifidelity (MF) groundwater simulations and transfer learning (TL) to reduce data generation costs while maintaining model accuracy. As a computational example, we train a deep convolutional neural network (CNN) to reconstruct permeability fields from saturation maps derived from a multiphase flow model. Starting with very low- and low-fidelity data generated on increasingly coarse meshes, we pretrain the CNN, followed by output-layer training and fine-tuning using only a limited number of high-fidelity samples. We demonstrate the surrogate’s robustness when interpreting low-quality inputs—such as interpolated maps or data affected by noise—which has strong implications for the applicability in practical hydrogeological scenarios. This multilevel MF-TL strategy achieves a favorable trade-off between computational efficiency and predictive accuracy, significantly outperforming high-fidelity-only approaches under the same computational budget.

Chiofalo, Alessia [University of Bologna] (ORCID:0↗

DOE OSTI2021

Supersaturation, buoyancy, and deep convection dynamics

Motivated by recent discussions concerning differences of convective dynamics in polluted and pristine environments, the so-called convective invigoration in particular, this paper provides an analysis of factors affecting convective updraft buoyancy, such as the in-cloud supersaturation, condensate and precipitation loading, and entrainment. We use the deep convective period from simulations of daytime convection development over land discussed in our previous publications. An entraining parcel framework is used in the theoretical analysis. We show that for the specific case considered here, finite (positive) supersaturation noticeably reduces pseudo-adiabatic parcel buoyancy and cumulative convective available potential energy (cCAPE) in the lower troposphere. This comes from keeping a small fraction of the water vapor in a supersaturated state and thus reducing the latent heating. Such a lower-tropospheric impact is comparable to the effects of condensate loading and entrainment in the idealized parcel framework. For the entire tropospheric depth, loading and entrainment have a much more significant impact on the total CAPE. For the cloud model results, we compare ensemble simulations applying either a bulk microphysics scheme with saturation adjustment or a more comprehensive double-moment scheme with supersaturation prediction. We compare deep convective updraft velocities, buoyancies, and supersaturations from all ensembles. In agreement with the parcel analysis, the saturation-adjustment scheme provides noticeably stronger updrafts in the lower troposphere. For the simulations predicting supersaturation, there are small differences between pristine and polluted conditions below the freezing level that are difficult to explain by standard analysis of the in-cloud buoyancy components. By applying the piggybacking technique, we show that the lower-tropospheric buoyancy differences between pristine and polluted simulations come from a combination of temperature (i.e., latent heating) and condensate loading differences that work together to make polluted buoyancies and updraft velocities slightly larger when compared to their pristine analogues. Overall, the effects are rather small and contradict previous claims of a significant invigoration of deep convection in polluted environments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Convection-Permitting Ensembles of an Isolated Mountain Thunderstorm during RELAMPAGO/CACTI

Abstract The north–south-oriented Sierras de Córdoba (SDC) ridge in central Argentina is noted for initiating thunderstorms that may grow into intense mesoscale convective systems (MCSs). It also initiates more isolated, shorter-lived cells under weaker synoptic forcing. These cells are less impactful than MCSs but may be difficult to predict in convective-scale numerical weather prediction (NWP) due to their strong sensitivities to subgrid and partially resolved processes. To study the mechanisms and predictability of such cells, convection-permitting ensemble simulations were conducted of an isolated, diurnally forced SDC thunderstorm during Cloud, Aerosol, and Complex Terrain Interactions (CACTI)/Remote Sensing of Electrification, Lightning, and Mesoscale/Microscale Processes with Adaptive Ground Observations (RELAMPAGO). The rich observational data facilitated detailed ensemble verification, where dry biases in the surface energy balance and soil moisture were identified. These biases promoted rapid removal of convective inhibition and an early onset of precipitating cells over the SDC that were shallower and weaker than the observed cell. Correction, and then overcorrection, of the soil moisture bias in two successive ensembles was required to rectify the surface energy balance and improve the representation of the SDC cell. Nevertheless, substantial ensemble variability in convective precipitation was found, with some members producing more widespread convection than observed and others producing no deep convection at all. This variability was largely explained by a combination of thermodynamic and dynamic mechanisms, dominated by a positive sensitivity of convective precipitation to preconvective moist instability over the ridge. Secondary sensitivities were found to low-level upward mass flux and midlevel cross-barrier winds, the latter of which caused gravity waves with elevated downdrafts that tended to suppress incipient clouds.

Lopez, Andres [Department of Atmospheric and Ocean↗

NASA NTRSDate not supplied

AMOC Stability in GISS-E2-1-G: A Case of Stochastic Bifurcation

A 10-member ensemble simulation using the NASA GISS-E2-1-G climate model shows a clear bifurcation in the Atlantic Meridional Overturning Circulation (AMOC) strength under the SSP2-4.5 future scenario, leading to 8 stable AMOC-on and 2 stable AMOC-off climate states after 400 years of integration, despite identical forcing of each ensemble member. Ocean models have previously shown a tipping point where sufficient forcing (e.g. by freshwater inputs) can shut down the AMOC. A variety of fully coupled models have demonstrated this, either through hosing experiments or increased precipitation and warming at high latitudes due to increased concentrations of greenhouse gases. In the GISS simulations, there are no external freshwater perturbations. The bifurcation arises freely in the coupled system and is the result of stochastic variability associated with sea-ice transport and melting in the Irminger Sea following a slowing of the GHG concentration increase at the end of the twenty-first century. We believe this is the first time that a CMIP-class model has shown such a bifurcation across an initial condition ensemble.

GISS ModelE2.1↗

NASA NTRSDate not supplied

Global Climate Impacts of Greenland and Antarctic Glacial Melt: A Comparative Study

GISS ModelE2.1↗

DOE OSTI2022

Data and scripts associated with the manuscript "Using Mutual Information for Global Sensitivity Analysis on Watershed Modeling"

This package contains the data and scripts used in "Using Mutual Information for Global Sensitivity Analysis on Watershed Modeling" (Jiang et al., 2022). The ARW_SWAT.zip file contains the SWAT simulation and the sensitivity analysis result in the American River Watershed. The Portage_SWAT.zip file contains the SWAT simulation and the sensitivity analysis result in the Portage River Watershed. The src folder contains the source code of sensitivity analysis and plotting functions. The notebooks folder contains the Jupyter notebooks for performing sensitivity analysis and other post analyses.Global Sensitivity Analysis (GSA) often is applied to assess the sensitivity of model outputs to their inputs using ensemble simulations. However, increasing model complexity and the associated computational cost have limited the use of most GSA approaches for process-based watershed models. We propose to use Mutual Information (MI) as a computationally efficient GSA method for watershed modeling. Such MI computed from several hundred realizations usually can capture nonlinear relationships between inputs and outputs of interest. We perform MI-based watershed sensitivity analyses in studies of the Portage River Watershed in Ohio and the American River Watershed in Washington. In these studies, MI is used to evaluate the sensitivity of river discharges simulated by the Soil and Water Assessment Tool to no less than 20 parameters for each watershed. Our MI-based sensitivity analyses achieved convergence with about 300~500 realizations, a small fraction of the ensemble size (i.e., several thousand) required by the Sobol method. Nevertheless, MI yields similar sensitivity ranking compared to the Sobol method, especially for sensitive parameters. Our study thus sheds new light on the use of MI as an affordable GSA method for computationally intensive models such as the hyper-resolution, watershed hydrobiogeochemical models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

Mechanisms Governing the Development of the North Atlantic Warming Hole in the CESM-LE Future Climate Simulations

A warming deficit in North Atlantic sea surface temperatures is a striking feature in global climate model future projections. This North Atlantic warming hole has been related to a slowing of the Atlantic meridional overturning circulation (AMOC); however, the detailed mechanisms involved in its generation remain an open question. An analysis of the Community Earth System Model Large Ensemble simulations is conducted to obtain further insight into the development of the warming hole and its relationship to the AMOC. It is shown that increasing freshwater fluxes through the Arctic gates lead to surface freshening and reduced Labrador Sea deep convection, which in turn act to cool Labrador Sea sea surface temperatures. Furthermore, the resulting changes in surface ocean circulation lead to enhanced transport of cooled Labrador Sea surface waters into the interior of the subpolar gyre and a more zonal orientation of the North Atlantic Current. As a result, there is an increase in ocean advective heat flux divergence within the center of the subpolar gyre, causing this warming deficit in North Atlantic sea surface temperatures. These local changes to the ocean circulation affect the AMOC and lead to its slowdown.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

Community climate simulations to assess avoided impacts in 1.5 and 2 °C futures

The Paris Agreement of December 2015 stated a goal to pursue efforts to keep global temperatures below 1.5 °C above preindustrial levels and well below 2 °C. The IPCC was charged with assessing climate impacts at these temperature levels, but fully coupled equilibrium climate simulations do not currently exist to inform such assessments. Here, we produce a set of scenarios using a simple model designed to achieve long-term 1.5 and 2 °C temperatures in a stable climate. These scenarios are then used to produce century-scale ensemble simulations using the Community Earth System Model, providing impact-relevant long-term climate data for stabilization pathways at 1.5 and 2 °C levels and an overshoot 1.5 °C case, which are realized (for the 21st century) in the coupled model and are freely available to the community. We also describe the design of the simulations and a brief overview of their impact-relevant climate response. Exceedance of historical record temperature occurs with 60 % greater frequency in the 2 °C climate than in a 1.5 °C climate aggregated globally, and with twice the frequency in equatorial and arid regions. Extreme precipitation intensity is statistically significantly higher in a 2.0 °C climate than a 1.5 °C climate in some specific regions (but not all). The model exhibits large differences in the Arctic, which is ice-free with a frequency of 1 in 3 years in the 2.0 °C scenario, and 1 in 40 years in the 1.5 °C scenario. Significance of impact differences with respect to multi-model variability is not assessed.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Aerosol Direct Radiative and Cloud Adjustment Effects on Surface Climate over Eastern China: Analyses of WRF Model Simulations

The WRF-simulated changes in clouds and climate due to the increased anthropogenic aerosols for the summers of 2002–08 (vs the 1970s) over eastern China were used to offline calculate the radiative forcings associated with aerosol–radiation (AR) and aerosol–cloud–radiation (ACR) interactions, which subsequently facilitated the interpretation of surface temperature changes. During this period, the increases of aerosol optical depth (ΔAOD) averaged over eastern China range from 0.18 in 2004 to 0.26 in 2007 as compared to corresponding cases in the 1970s, and the multiyear means (standard deviations) of AR and ACR forcings at the surface are −6.7 (0.58) and −3.5 (0.63) W m −2 , respectively, indicating the importance of cloud changes in affecting both the aerosol climate forcing and its interannual variation. The simulated mean surface cooling is 0.35°C, dominated by AR and ACR with a positive (cooling) feedback associated with changes in meteorology (~10%), and two negative (warming) feedbacks associated with decreases in latent (~70%) and sensible (~20%) heat fluxes. More detailed spatial characteristics were analyzed using ensemble simulations for the year 2008. Three regions—Jing-Jin-Ji (ΔAOD ~ 0.63), Sichuan basin (ΔAOD ~ 0.31), and middle Yangtze River valley (ΔAOD ~ 0.26)—at different climate regimes were selected to investigate the relative roles of AR and ACR. While the AR forcing is closely related to ΔAOD values, the ACR forcing presents different regional characteristics owing to cloud changes. In addition, the surface heat flux feedbacks are also different between regions. The study thus illustrates that ACR forcing is useful as a diagnostic parameter to unravel the complexity of climate change to aerosol forcing over eastern China.

54 ENVIRONMENTAL SCIENCES↗