A Perturbed Parameter Ensemble (PPE) with the Community Atmosphere Model version 6 (CAM6) is used to better understand the sensitivity of simulated clouds to both aerosol forcing and cloud feedbacks and the interactions between them. 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 net forcing and larger net feedback. Aerosol forcing is broadly related to the susceptibility of clouds to drop number. Feedbacks are less related to susceptibility, and in different regions. 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.

Ice nucleation in mixed-phase clouds has recently been identified as a critical factor in projections of future climate. Here we explore how this process influences climate sensitivity using the Community Earth System Model 2 (CESM2). We find that ice nucleation affects simulated cloud feedbacks over most regions and levels of the troposphere, not just extratropical low clouds. Ice nucleation’s impact on climate sensitivity is found to primarily operate through this process’s role setting global-scale cloud phase. Conversely, whether ice nucleation is treated as aerosol-sensitive is of limited importance. In satellite-constrained model experiments, dissimilar ice nucleation realizations all result in a strongly positive total cloud feedback, as in the default CESM2. A microphysics update from CESM1 to CESM2 had substantially weakened ice nucleation, due partly to a model issue. Our findings suggest that this contributed to increased climate sensitivity by reducing global cloud phase bias, resulting in more realistic mixed-phase clouds.

A single column model with parameterized large-scale dynamics is used to better understand the response of steady-state tropical precipitation to relative sea surface temperature under various representations of radiation, convection, and circulation. The large-scale dynamics are parametrized via the weak temperature gradient (WTG), damped gravity wave (DGW), and spectral weak temperature gradient (Spectral WTG) method in NCAR’s Single Column Atmosphere Model (SCAM6). Radiative cooling is either specified or interactive, and the convective parameterization is run using two different values of a parameter that controls the degree of convective inhibition. Results are interpreted in the context of the Global Atmospheric System Studies (GASS) Intercomparison (Daleu et al. 2016). Using the settings given in Daleu et al. (2016), SCAM6 under the WTG and DGW methods produces erratic results, suggestive of numerical instability. However, when key parameters are changed to weaken the strength with which the circulation acts to eliminate tropospheric temperature variations, SCAM6 performs comparably to single column models in the GASS Intercomparison. The Spectral WTG method is less sensitive to changes in convection and radiation than are the other two methods, performing at least qualitatively similarly across all configurations considered. Under all three methods, circulation strength, represented in 1D by grid-scale vertical velocity, is decreased when barriers to convection are reduced. This effect is most extreme under specified radiative cooling, and is shown to come from increased static stability in the column’s reference radiative-convective equilibrium profile. This argument can be extended to interactive radiation cases as well, though perhaps less conclusively.

The simulation of ice sheet-climate interaction such as surface mass balance fluxes are sensitive to model grid resolution. Here we simulate the multicentury evolution of the Greenland Ice Sheet (GrIS) and its interaction with the climate using the Community Earth System Model version 2.2 (CESM2.2) including an interactive GrIS component (the Community Ice Sheet Model v2.1 [CISM2.1]) under an idealized warming scenario (atmospheric CO2 increases by 1% yr−1 until quadrupling the pre-industrial level and then is held fixed). A variable-resolution (VR) grid with 1/4◦ regional refinement over broader Arctic and 1◦ resolution elsewhere is applied to the atmosphere and land components, and the results are compared to conventional 1◦ lat-lon grid simulations to investigate the impact of grid refinement. An acceleration of GrIS mass loss is found at around year 110, caused by rapidly increasing surface melt as the ablation area expands with associated albedo feedback and increased turbulent fluxes. Compared to the 1◦ runs, the VR run features slower melt increase, especially over Western and Northern Greenland, which slope gently towards the peripheries. This difference pattern originates primarily from the weaker albedo feedback in the VR run, complemented by its smaller cloud longwave radiation. The steeper VR Greenland surface topography favors slower ablation zone expansion, thus leading to its weaker albedo feedback. The sea level rise contribution from the GrIS in the VR run is 53 mm by year 150 and 831 mm by year 350, approximately 40% and 20% smaller than the 1◦ runs, respectively.

Simulating whole atmosphere dynamics, chemistry, and physics is computationally expensive. It can require high vertical resolution throughout the middle and upper atmosphere, as well as a comprehensive chemistry and aerosol scheme coupled to radiation physics. An unintentional outcome of the development of one of the most sophisticated and hence computationally expensive model configurations is that it often excludes a broad community of users with limited computational resources. Here, we analyze two configurations of the Community Earth System Model Version 2, Whole Atmosphere Community Climate Model Version 6 (CESM2(WACCM6)) with simplified “middle atmosphere” chemistry at nominal 1 and 2 degree horizontal resolutions. Using observations, a reanalysis, and direct model comparisons, we find that these configurations generally reproduce the climate, variability, and climate sensitivity of the 1 degree nominal horizontal resolution configuration with comprehensive chemistry. While the background stratospheric aerosol optical depth is elevated in the middle atmosphere configurations as compared to the comprehensive chemistry configuration, it is comparable between all configurations during volcanic eruptions. For any purposes other than those needing an accurate representation of tropospheric organic chemistry and secondary organic aerosols, these simplified chemistry configurations deliver reliable simulations of the whole atmosphere that require 35% to 86% fewer computational resources at nominal 1 and 2 degree horizontal resolution, respectively.