Silvers, Levi G., and Thomas E Robinson, February 2021: Clouds and radiation in a mock-Walker circulation. Journal of Advances in Modeling Earth Systems, 13(2), DOI:10.1029/2020MS002196. Abstract
The Walker circulation connects the regions with deep atmospheric convection in the western tropical Pacific to the shallow-convection, tropospheric subsidence, and stratocumulus cloud decks of the eastern Pacific. The purpose of this study is to better understand the multi-scale interactions between the Walker circulation, cloud systems, and interactive radiation. To do this we simulate a mock-Walker Circulation with a full-physics general circulation model using idealized boundary conditions. Our experiments use a doubly-periodic domain with grid-spacing of 1, 2, 25, and 100 km. We thus span the range from General Circulation Models (GCMs) to Cloud-system Resolving Models (CRMs). Our model is derived from the Geophysical Fluid Dynamics Laboratory atmospheric GCM (AM4.0). We find substantial differences in the mock-Walker circulation simulated by our GCM-like and CRM-like experiments. The CRM-like experiments have more upper level clouds, stronger overturning circulations, and less precipitation. The GCM-like experiments have a low-level cloud fraction that is up to 20% larger. These differences leads to opposite atmospheric responses to changes in the longwave cloud radiative effect (LWCRE). Active LWCRE leads to increased precipitation for our GCMs, but decreased precipitation for our CRMs. The LWCRE leads to a narrower rising branch of the circulation and substantially increases the fraction of precipitation from the large-scale cloud parameterization. This work demonstrates that a mock-Walker circulation is a useful generalization of radiative convective equilibrium that includes a large-scale circulation.
Zhang, Gan, Levi G Silvers, Ming Zhao, and Thomas R Knutson, March 2021: Idealized aquaplanet simulations of tropical cyclone activity: Significance of temperature gradients, Hadley circulation, and zonal asymmetry. Journal of the Atmospheric Sciences, 78(3), DOI:10.1175/JAS-D-20-0079.1877-902. Abstract
Earlier studies have proposed many semiempirical relations between climate and tropical cyclone (TC) activity. To explore these relations, this study conducts idealized aquaplanet experiments using both symmetric and asymmetric sea surface temperature (SST) forcings. With zonally symmetric SST forcings that have a maximum at 10°N, reducing meridional SST gradients around an Earth-like reference state leads to a weakening and southward displacement of the intertropical convergence zone. With nearly flat meridional gradients, warm-hemisphere TC numbers increase by nearly 100 times due particularly to elevated high-latitude TC activity. Reduced meridional SST gradients contribute to a poleward expansion of the tropics, which is associated with a poleward migration of the latitudes where TCs form or reach their lifetime maximum intensity. However, these changes cannot be simply attributed to the poleward expansion of Hadley circulation. Introducing zonally asymmetric SST forcings tends to decrease the global TC number. Regional SST warming—prescribed with or without SST cooling at other longitudes—affects local TC activity but does not necessarily increase TC genesis. While regional warming generally suppresses TC activity in remote regions with relatively cold SSTs, one experiment shows a surprisingly large increase of TC genesis. This increase of TC genesis over relatively cold SSTs is related to local tropospheric cooling that reduces static stability near 15°N and vertical wind shear around 25°N. Modeling results are discussed with scaling analyses and have implications for the application of the “convective quasi-equilibrium and weak temperature gradient” framework.
Wing, Allison A., Catherine L Stauffer, Tobias Becker, Kevin A Reed, Min-Seop Ahn, Nathan P Arnold, Sandrine Bony, Mark Branson, George H Bryan, Jean-Pierre Chaboureau, Stephan R de Roode, Kulkarni Gayatri, Cathy Hohenegger, I-Kuan Hu, Fredrik Jansson, Todd R Jones, Marat Khairoutdinov, Daehyun Kim, Zane K Martin, Shuhei Matsugishi, Brian Medeiros, Hiroaki Miura, Yumin Moon, Sebastian K Muller, Tomoki Ohno, Max Popp, Thara Prabhakaran, David A Randall, Rosimar Rios-Berrios, Nicolas Rochetin, Romain Roehrig, David M Romps, James H Ruppert Jr, Masaki Satoh, Levi G Silvers, Martin S Singh, Bjorn Stevens, Lorenzo Tomassini, Chiel C van Heerwaarden, Shuguang Wang, and Ming Zhao, September 2020: Clouds and Convective Self-Aggregation in a Multi-Model Ensemble of Radiative-Convective Equilibrium Simulations. Journal of Advances in Modeling Earth Systems, 12(9), DOI:10.1029/2020MS002138. Abstract
The Radiative‐Convective Equilibrium Model Intercomparison Project (RCEMIP) is an intercomparison of multiple types of numerical models configured in radiative‐convective equilibrium (RCE). RCE is an idealization of the tropical atmosphere that has long been used to study basic questions in climate science. Here, we employ RCE to investigate the role that clouds and convective activity play in determining cloud feedbacks, climate sensitivity, the state of convective aggregation, and the equilibrium climate. RCEMIP is unique among intercomparisons in its inclusion of a wide range of model types, including atmospheric general circulation models (GCMs), single column models (SCMs), cloud‐resolving models (CRMs), large eddy simulations (LES), and global cloud‐resolving models (GCRMs). The first results are presented from the RCEMIP ensemble of more than 30 models. While there are large differences across the RCEMIP ensemble in the representation of mean profiles of temperature, humidity, and cloudiness, in a majority of models anvil clouds rise, warm, and decrease in area coverage in response to an increase in sea surface temperature (SST). Nearly all models exhibit self‐aggregation in large domains and agree that self‐aggregation acts to dry and warm the troposphere, reduce high cloudiness, and increase cooling to space. The degree of self‐aggregation exhibits no clear tendency with warming. There is a wide range of climate sensitivities, but models with parameterized convection tend to have lower climate sensitivities than models with explicit convection. In models with parameterized convection, aggregated simulations have lower climate sensitivities than unaggregated simulations.
GFDL's new CM4.0 climate model has high transient and equilibrium climate sensitivities near the middle of the upper half of CMIP5 models. The CMIP5 models have been criticized for excessive sensitivity based on observations of present‐day warming and heat uptake and estimates of radiative forcing. An ensemble of historical simulations with CM4.0 produces warming and heat uptake that are consistent with these observations under forcing that is at the middle of the assessed distribution. Energy budget‐based methods for estimating sensitivities based on these quantities underestimate CM4.0's sensitivities when applied to its historical simulations. However, we argue using a simple attribution procedure that CM4.0's warming evolution indicates excessive transient sensitivity to greenhouse gases. This excessive sensitivity is offset prior to recent decades by excessive response to aerosol and land use changes.
Mixed-phase clouds are frequently observed in the atmosphere. Here we present a parameterization for ice crystal concentration and ice nucleation rate based on parcel model simulations for mixed-phase stratocumulus clouds, in complement to a previous parameterization for stratus clouds. The parcel model uses a singular (time-independent) description for deposition nucleation and a time-dependent description for condensation nucleation and immersion freezing on mineral dust particles. The mineral dust and temperature-dependent parameterizations have been implemented in the Geophysical Fluid Dynamics Laboratory atmosphere model AM4.0 (new), while the standard AM4.0 (original) uses a temperature-dependent parameterization. Model simulations with the new and original AM4.0 show significant changes in cloud properties and radiative effects. In comparison to measurements, cloud-phase (i.e., liquid and ice partitioning) simulation appears to be improved in the new AM4.0 model. More supercooled liquid cloud is predicted in the new model, it is sustained even at temperatures lower than -25 °C unlike in the original model. A more accurate accounting of ice nucleating particles and ice crystals is essential for improved cloud phase simulation in the global atmosphere.
We describe GFDL's CM4.0 physical climate model, with emphasis on those aspects that may be of particular importance to users of this model and its simulations. The model is built with the AM4.0/LM4.0 atmosphere/land model and OM4.0 ocean model. Topics include the rationale for key choices made in the model formulation, the stability as well as drift of the pre‐industrial control simulation, and comparison of key aspects of the historical simulations with observations from recent decades. Notable achievements include the relatively small biases in seasonal spatial patterns of top‐of‐atmosphere fluxes, surface temperature, and precipitation; reduced double Intertropical Convergence Zone bias; dramatically improved representation of ocean boundary currents; a high quality simulation of climatological Arctic sea ice extent and its recent decline; and excellent simulation of the El Niño‐Southern Oscillation spectrum and structure. Areas of concern include inadequate deep convection in the Nordic Seas; an inaccurate Antarctic sea ice simulation; precipitation and wind composites still affected by the equatorial cold tongue bias; muted variability in the Atlantic Meridional Overturning Circulation; strong 100 year quasi‐periodicity in Southern Ocean ventilation; and a lack of historical warming before 1990 and too rapid warming thereafter due to high climate sensitivity and strong aerosol forcing, in contrast to the observational record. Overall, CM4.0 scores very well in its fidelity against observations compared to the Coupled Model Intercomparison Project Phase 5 generation in terms of both mean state and modes of variability and should prove a valuable new addition for analysis across a broad array of applications.
Andrews, Timothy, Jonathan M Gregory, David J Paynter, and Levi G Silvers, et al., August 2018: Accounting for Changing Temperature Patterns Increases Historical Estimates of Climate Sensitivity. Geophysical Research Letters, 45(16), DOI:10.1029/2018GL078887. Abstract
Eight atmospheric general circulation models (AGCMs) are forced with observed historical (1871–2010) monthly sea surface temperature and sea ice variations using the Atmospheric Model Intercomparison Project II data set. The AGCMs therefore have a similar temperature pattern and trend to that of observed historical climate change. The AGCMs simulate a spread in climate feedback similar to that seen in coupled simulations of the response to CO2 quadrupling. However, the feedbacks are robustly more stabilizing and the effective climate sensitivity (EffCS) smaller. This is due to a pattern effect, whereby the pattern of observed historical sea surface temperature change gives rise to more negative cloud and longwave clear‐sky feedbacks. Assuming the patterns of long‐term temperature change simulated by models, and the radiative response to them, are credible; this implies that existing constraints on EffCS from historical energy budget variations give values that are too low and overly constrained, particularly at the upper end. For example, the pattern effect increases the long‐term Otto et al. (2013, https://doi.org/10.1038/ngeo1836) EffCS median and 5–95% confidence interval from 1.9 K (0.9–5.0 K) to 3.2 K (1.5–8.1 K).
Giorgetta, M A., R Brokopf, T Crueger, M Esch, Stephanie Fiedler, J Helmert, Cathy Hohenegger, Luis Kornblueh, M Köhler, E Manzini, T Mauritsen, C Nam, T Raddatz, S Rast, D Reinert, M Sakradzija, H Schmidt, R Schneck, R Schnur, and Levi G Silvers, et al., July 2018: ICON‐A, the Atmosphere Component of the ICON Earth System Model: I. Model Description. Journal of Advances in Modeling Earth Systems, 10(7), DOI:10.1029/2017MS001242. Abstract
ICON‐A is the new icosahedral nonhydrostatic (ICON) atmospheric general circulation model in a configuration using the Max Planck Institute physics package, which originates from the ECHAM6 general circulation model, and has been adapted to account for the changed dynamical core framework. The coupling scheme between dynamics and physics employs a sequential updating by dynamics and physics, and a fixed sequence of the physical processes similar to ECHAM6. To allow a meaningful initial comparison between ICON‐A and the established ECHAM6‐LR model, a setup with similar, low resolution in terms of number of grid points and levels is chosen. The ICON‐A model is tuned on the base of the Atmospheric Model Intercomparison Project (AMIP) experiment aiming primarily at a well balanced top‐of atmosphere energy budget to make the model suitable for coupled climate and Earth system modeling. The tuning addresses first the moisture and cloud distribution to achieve the top‐of‐atmosphere energy balance, followed by the tuning of the parameterized dynamic drag aiming at reduced wind errors in the troposphere. The resulting version of ICON‐A has overall biases, which are comparable to those of ECHAM6. Problematic specific biases remain in the vertical distribution of clouds and in the stratospheric circulation, where the winter vortices are too weak. Biases in precipitable water and tropospheric temperature are, however, reduced compared to the ECHAM6. ICON‐A will serve as the basis of further development and as the atmosphere component to the coupled model, ICON‐Earth system model (ESM).
Paynter, David J., Thomas L Frölicher, Larry W Horowitz, and Levi G Silvers, February 2018: Equilibrium Climate Sensitivity Obtained from Multi-Millennial Runs of Two GFDL Climate Models. Journal of Geophysical Research: Atmospheres, 123(4), DOI:10.1002/2017JD027885. Abstract
Equilibrium climate sensitivity (ECS), defined as the long-term change in global mean surface air temperature in response to doubling atmospheric CO2, is usually computed from short atmospheric simulations over a mixed layer ocean, or inferred using a linear regression over a short-time period of adjustment. We report the actual ECS from multi-millenial simulations of two GFDL general circulation models (GCMs), ESM2M and CM3 of 3.3 K and 4.8 K, respectively. Both values are ~1 K higher than estimates for the same models reported in the Fifth Assessment Report of the Intergovernmental Panel on Climate Change obtained by regressing the Earth's energy imbalance against temperature. This underestimate is mainly due to changes in the climate feedback parameter (−α) within the first century after atmospheric CO2 has stabilized. For both GCMs it is possible to estimate ECS with linear regression to within 0.3 K by increasing CO2 at 1% per year to doubling and using years 51-350 after CO2 is constant. We show that changes in −α differ between the two GCMs and are strongly tied to the changes in both vertical velocity at 500 hPa (ω500) and estimated inversion strength (EIS) that the GCMs experience during the progression towards the equilibrium. This suggests that while cloud physics parametrizations are important for determining the strength of −α, the substantially different atmospheric state resulting from a changed SST pattern may be of equal importance.
Silvers, Levi G., David J Paynter, and Ming Zhao, January 2018: The Diversity of Cloud Responses to Twentieth Century Sea Surface Temperatures. Geophysical Research Letters, 45(1), DOI:10.1002/2017GL075583. Abstract
Low-level clouds are shown to be the conduit between the observed sea surface temperatures (SST) and large decadal fluctuations of the top of the atmosphere (TOA) radiative imbalance. The influence of low-level clouds on the climate feedback is shown for global mean time series as well as particular geographic regions. The changes of clouds are found to be important for a mid-century period of high sensitivity and a late century period of low sensitivity. These conclusions are drawn from analysis of amip-piForcing simulations using three atmospheric general circulation models (AM2.1, AM3, and AM4.0). All three models confirm the importance of the relationship between the global climate sensitivity and the eastern Pacific trends of SST and low-level clouds. However, this work argues that the variability of the climate feedback parameter is not driven by stratocumulus dominated regions in the eastern ocean basins, but rather by the cloudy response in the rest of the tropics.
In this two-part paper, a description is provided of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). This version, with roughly 100km horizontal resolution and 33 levels in the vertical, contains an aerosol model that generates aerosol fields from emissions and a “light” chemistry mechanism designed to support the aerosol model but with prescribed ozone. In Part I, the quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode – with prescribed sea surface temperatures (SSTs) and sea ice distribution – is described and compared with previous GFDL models and with the CMIP5 archive of AMIP simulations. The model's Cess sensitivity (response in the top-of-atmosphere radiative flux to uniform warming of SSTs) and effective radiative forcing are also presented. In Part II, the model formulation is described more fully and key sensitivities to aspects of the model formulation are discussed, along with the approach to model tuning.
In Part II of this two-part paper, documentation is provided of key aspects of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). The quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode has been provided in Part I. Part II provides documentation of key components and some sensitivities to choices of model formulation and values of parameters, highlighting the convection parameterization and orographic gravity wave drag. The approach taken to tune the model's clouds to observations is a particular focal point. Care is taken to describe the extent to which aerosol effective forcing and Cess sensitivity have been tuned through the model development process, both of which are relevant to the ability of the model to simulate the evolution of temperatures over the last century when coupled to an ocean model.
Popp, Max, and Levi G Silvers, November 2017: Double and Single ITCZs with and without Clouds. Journal of Climate, 30(22), DOI:10.1175/JCLI-D-17-0062.1. Abstract
A major bias in tropical precipitation over the Pacific in climate simulations stems from the models’ tendency to produce two strong distinct intertropical convergence zones (ITCZs) too often. Several mechanisms have been proposed that may contribute to the emergence of two ITCZs, but current theories cannot fully explain the bias. This problem is tackled by investigating how the interaction between atmospheric cloud-radiative effects (ACREs) and the large-scale circulation influences the ITCZ position in an atmospheric general circulation model. Simulations are performed in an idealized aquaplanet setup and the longwave and shortwave ACREs are turned off individually or jointly. The low-level moist static energy (MSE) is shown to be a good predictor of the ITCZ position. Therefore, a mechanism is proposed that explains the changes in MSE and thus ITCZ position due to ACREs consistently across simulations. The mechanism implies that the ITCZ moves equatorward if the Hadley circulation strengthens because of the increased upgradient advection of low-level MSE off the equator. The longwave ACRE increases the meridional heating gradient in the tropics and as a response the Hadley circulation strengthens and the ITCZ moves equatorward. The shortwave ACRE has the opposite effect. The total ACRE pulls the ITCZ equatorward. This mechanism is discussed in other frameworks involving convective available potential energy, gross moist stability, and the energy flux equator. It is thus shown that the response of the large-scale circulation to the shortwave and longwave ACREs is a fundamental driver of changes in the ITCZ position.