Previous modeling studies indicate that even though marine cloud brightening under a susceptibility-based strategy is effective in reducing the global average surface temperature, it triggers a La Niña-like sea-surface temperature response with cooling mostly confined within lower latitudes. Here we explore a different cloud seeding strategy involving seeding of regions with low susceptibility. Simulations with the Community Earth System Model, version 2 (CESM2) reveal that because the regional forcing is weaker and more widespread, cooling is more evenly distributed over the globe. This new strategy also does not result in the La Niña-like state seen in the other strategies.

The efficiency of marine cloud brightening in cooling Earth’s surface temperature is investigated by using a medium ensemble of simulations with the Community Earth System Model version 2 (CESM2). Various cloud seeding schemes based on susceptibility are examined to determine what area extent will be required to induce 1oC cooling under SSP2-4.5. The results indicate that cloud seeding over 5% of the ocean area is capable of achieving this goal. Under this seeding scheme, cloud seeding is mainly deployed over lower latitudes where strong surface temperature and precipitation responses are induced. The simulations also reveal that the 5% cloud seeding scheme induces an overall reduction in global precipitation, with an increase over land and a decrease over the ocean.

Mesoscale convective systems (MCSs) play an important role in modulating the global hydrological cycle, general circulation, and radiative energy budget. In this study, we evaluate MCS simulations in the second version of U.S. Department of Energy (DOE) Energy Exascale Earth System Model (E3SMv2). E3SMv2 atmosphere model (EAMv2) is run at the uniform 0.25° horizontal resolution. We track MCSs consistently in the model and observations using the PyFLEXTRKR algorithm, which defines MCS based on both cloud-top brightness temperature (Tb) and surface precipitation. Results from using Tb only to define MCS, commonly used in previous studies, are also discussed. Furthermore, sensitivity experiments are performed to examine the impact of new cloud and convection parameterizations developed for EAMv3 on simulated MCSs. Our results show that EAMv2 simulated MCS precipitation is largely underestimated in the tropics and contiguous United States. This is mainly attributed to the underestimated precipitation intensity in EAMv2. In contrast, the simulated MCS frequency becomes more comparable to observations if MCSs are defined only based on cloud-top Tb. The Tb-based MCS tracking method, however, includes many cloud systems with very weak precipitation which conflicts with the MCS definition. This result illustrates the importance of accounting for precipitation in evaluating simulated MCSs. We also find that the new physics parameterizations help increase the relative contribution of convective precipitation to total precipitation in the tropics, but the simulated MCS properties are overall not significantly improved. This suggests that simulating MCSs will remain a challenge for the next version of E3SM.

Current global actions to reduce greenhouse gas emissions are very likely to be insufficient to meet the climate targets outlined under the Paris Agreement. This motivates research on possible methods for intervening in the Earth system to minimize climate risk while decarbonization efforts continue. One such hypothetical climate intervention is stratospheric aerosol injection (SAI), where reflective particles would be released into the stratosphere to cool the planet by reducing solar insolation. The climate response to SAI is not well understood, particularly on short-term time horizons frequently used by decision makers and planning practitioners to assess climate information. This knowledge gap limits informed discussion of SAI outside the scientific community. We demonstrate two framings to explore the climate response in the decade after SAI deployment in modeling experiments with parallel SAI and no-SAI simulations. The first framing, which we call a snapshot around deployment, displays change over time within the SAI scenarios and applies to the question “What happens before and after SAI is deployed in the model?” The second framing, the intervention impact, displays the difference between the SAI and no-SAI simulations, corresponding to the question “What is the impact of a given intervention relative to climate change with no intervention?” We apply these framings to annual mean 2-meter temperature, precipitation, and a precipitation extreme in the first two experiments to use large ensembles of Earth system models that comprehensively represent both the SAI injection process and climate response, and connect these results to implications for other climate variables.

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.

By injecting SO2 into the stratosphere at four latitudes (30°, 15° N/S), it might be possible not only to reduce global mean surface temperature but also to minimize changes in the equator-to-pole and inter-hemispheric gradients of temperature, further reducing some of the impacts arising from climate change relative to equatorial injection. This can happen only if the aerosols are transported to higher latitudes by the stratospheric circulation, ensuring that a greater part of the solar radiation is reflected back to space at higher latitudes, compensating for the reduced sunlight. However, the stratospheric heating produced by these aerosols modifies the circulation and strengthens the stratospheric polar vortex which acts as a barrier to the transport of air toward the poles. We show how the heating results in a feedback where increasing injection rates lead to stronger high-latitudinal transport barriers. This implies a potential limitation in the high-latitude aerosol burden and subsequent cooling.

This work documents version two of the Department of Energy’s Energy Exascale Earth System Model (E3SM). E3SM version 2 (E3SMv2) is a significant evolution from its predecessor E3SMv1, resulting in a model that is nearly twice as fast and with a simulated climate that is improved in many metrics. We describe the physical climate model in its lower horizontal resolution configuration consisting of 110 km atmosphere, 165 km land, 0.5° river routing model, and an ocean and sea ice with mesh spacing varying between 60 km in the mid-latitudes and 30 km at the equator and poles. The model performance is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima (DECK) simulations augmented with historical simulations as well as simulations to evaluate impacts of different forcing agents. The simulated climate is generally realistic, with notable improvements in clouds and precipitation compared to E3SMv1. E3SMv1 suffered from an excessively high equilibrium climate sensitivity (ECS) of 5.3 K. In E3SMv2, ECS is reduced to 4.0 K which is now within the plausible range based on a recent World Climate Research Programme (WCRP) assessment. However, E3SMv2 significantly underestimates the global mean surface temperature in the second half of the historical record. An analysis of single-forcing simulations indicates that correcting the historical temperature bias would require a substantial reduction in the magnitude of the aerosol-related forcing.

Model dependence in simulated responses to stratospheric aerosol injection (SAI) is a major uncertainty surrounding the potential implementation of this solar climate intervention strategy. We identify large differences in the aerosol mass latitudinal distributions between two recently produced climate model SAI large ensembles, despite using similar climate targets and controller algorithms, with the goal of understanding the drivers of such differences. Using a hierarchy of recently produced simulations, we identify three main contributors including: 1) the rapid adjustment of clouds and rainfall to elevated levels of carbon dioxide, 2) the associated low-frequency dynamical responses in the Atlantic Meridional Overturning Circulation, and 3) the contrasts in future climate forcing scenarios. Each uncertainty is unlikely to be significantly narrowed over the likely timeframe of a potential SAI deployment if a 1.5C target is to be met. The results thus suggest the need for significant flexibility in climate intervention deployment in order to account for these large uncertainties in the climate system response.

Mesoscale organization of convection is typically not represented in global circulation models, and hence its influence on the global circulation is not accounted for. A parameterization aiming at representing the dynamical and physical effects of the circulation associated with organized convection, referred to as the multiscale coherent structure parameterization (MCSP), is implemented in the Energy Exascale Earth System Model version 1 (E3SMv1). Simulations are conducted to assess its impact on the simulated climate. Besides E3SMv1 simulations, we performed high-resolution (1 km) simulations using the Weather Research and Forecasting (WRF) Model to determine the temperature tendencies induced by mesoscale convective systems embedded in deep convection. We tuned the free parameters of the MCSP based on the WRF simulations. We found that the MCSP enhances Kevin wave spectra in E3SMv1, improves the representation of the Madden-Julian Oscillation, and reduces model precipitation biases over the tropical Pacific.

The week 3-6 averaged winter North Atlantic Oscillation (NAO) predictive skill in a state-of-the-art coupled climate prediction system is attributed to two principle sources: upper and lower boundary conditions linked to the stratosphere and El Niño-Southern Oscillation (ENSO), respectively. A 20-member ensemble of 45-day reforecasts over 1999-2015 is utilized, together with uninitialized simulations with the atmospheric component of the prediction system forced with observed radiative forcing and lower boundary conditions. NAO forecast skill for lead times out to six weeks is higher following extreme stratospheric polar vortex conditions (weak and strong vortex events) compared to neutral states. Enhanced week 3-6 NAO predictive skill for weak vortex events results primarily from stratospheric downward coupling to the troposphere, while enhanced skill for strong vortex events can be partly attributed to lower boundary forcing related to the ENSO phenomenon. Implications for forecast system development and improvement are discussed.

The impacts of Stratospheric Aerosol Injection strategies on the Southern Annular Mode (SAM) are analysed with the Community Earth System Model (CESM). Using a set of simulations with fixed single-point SO2 injections we demonstrate the first-order dependence of the SAM response on the latitude of injection, with the northern hemispheric and equatorial injections driving a response corresponding to a positive phase of SAM and the southern hemispheric injections driving a negative phase of SAM. We further demonstrate that the results can to first order explain the differences in the SAM responses diagnosed from the two recent large ensembles of geoengineering simulations utilising more complex injection strategies - GLENS and ARISE-SAI - as driven by the differences in the simulated sulfate aerosol distributions. Our results point to the meridional extent of aerosol-induced lower stratospheric heating as an important driver of the sensitivity of the SAM response to the injection location.

A framework to enable Earth system predictability research on the subseasonal timescale is developed with the Community Earth System Model, version 2 (CESM2) using two model configurations that differ in their atmospheric components. One configuration uses the Community Atmosphere Model, version 6 (CAM6) with its top near 40 km, referred to as CESM2(CAM6). The other employs the Whole Atmosphere Community Climate Model, version 6 (WACCM6) whose top extends to ~ 140 km in the vertical and it includes fully interactive tropospheric and stratospheric chemistry (CESM2(WACCM6)). Both configurations were used to carry out subseasonal reforecasts for the time period 1999 to 2020 following the Subseasonal Experiment’s (SubX) protocol. CESM2(CAM6) and CESM2(WACCM6) show very similar subseasonal prediction skill of 2-meter temperature, precipitation, the Madden-Julian Oscillation (MJO), and North Atlantic Oscillation (NAO) to the Community Earth System Model, version 1 with the Community Atmosphere Model, version 5 (CESM1(CAM5)) and to operational models. CESM2(CAM6) and CESM2(WACCM6) reforecast sets provide a comprehensive dataset for predictability research of multiple Earth system components, including three-dimensional output for many variables, and output specific to the mesosphere and lower-thermosphere (MLT) region. We show that MLT variability can be predicted ~ 10 days in advance of sudden stratospheric warming events. Weekly real-time forecasts with CESM2(WACCM6) contribute to the multi-model mean ensemble forecast used to issue the NOAA weeks 3-4 outlooks. As a freely available community model, both CESM2 configurations can be used to carry out additional experiments to elucidate sources of subseasonal predictability.