the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
![](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-avatar-thumb150.png)
RCEMIP-II: mock-Walker simulations as phase II of the radiative–convective equilibrium model intercomparison project
Levi G. Silvers
Kevin A. Reed
The radiative–convective equilibrium (RCE) model intercomparison project (RCEMIP) leveraged the simplicity of RCE to focus attention on moist convective processes and their interactions with radiation and circulation across a wide range of model types including cloud-resolving models (CRMs), general circulation models (GCMs), single-column models, global cloud-resolving models, and large-eddy simulations. While several robust results emerged across the spectrum of models that participated in the first phase of RCEMIP (RCEMIP-I), two points that stand out are (1) the strikingly large diversity in simulated climate states and (2) the strong imprint of convective self-aggregation on the climate state. However, the lack of consensus in the structure of self-aggregation and its response to warming is a barrier to understanding. Gaining a deeper understanding of convective aggregation and tropical climate will require reducing the degrees of freedom with which convection can vary. Therefore, we propose phase II of RCEMIP (RCEMIP-II) that utilizes a prescribed sinusoidal sea surface temperature (SST) pattern to provide a constraint on the structure of convection and move one critical step up the model hierarchy. This so-called “mock-Walker” configuration generates features that resemble observed tropical circulations. The specification of the mock-Walker protocol for RCEMIP-II is described, along with example results from one CRM and one GCM. RCEMIP-II will consist of five required simulations: three simulations with the same three mean SSTs as in RCEMIP-I but with an SST gradient and two additional simulations at one of the mean SSTs with different values of the SST gradients. We also test the sensitivity to the imposed SST gradient and the domain size. Under weak SST gradients, unforced self-aggregation emerges across the entire domain, similar to what was found in RCEMIP. As the SST gradient increases, the convective region narrows and is more confined to the warmest SSTs. At warmer mean SSTs and stronger SST gradients, low-frequency variability in the convective aggregation emerges, suggesting that simulations of at least 200 d may be needed to achieve robust equilibrium statistics in this configuration. Simulations with different domain sizes generally have similar mean statistics and convective structures, depending on the value of the SST gradient. The prescribed SST boundary condition is the only difference in the set-up between RCEMIP-II and RCEMIP-I, which enables comparison between the two; however, we also welcome participation in RCEMIP-II from models that did not participate in RCEMIP-I.
Please read the corrigendum first before continuing.
-
Notice on corrigendum
The requested paper has a corresponding corrigendum published. Please read the corrigendum first before downloading the article.
-
Article
(15006 KB)
-
The requested paper has a corresponding corrigendum published. Please read the corrigendum first before downloading the article.
- Article
(15006 KB) - Full-text XML
- Corrigendum
- Companion paper
- BibTeX
- EndNote
On Earth, the tropics play an important role in climate through coupled interactions of clouds, circulation, and radiative fluxes. As such, tropical regions influence the global energy balance, and tropical variability has far-reaching global effects on weather patterns and extremes. Thus, understanding the underlying mechanisms that connect climate, clouds, and circulation remains paramount in climate science, particularly in the context of global warming (Bony et al., 2015). Global models, such as general circulation models, have evolved over recent decades to advance our understanding and simulation of the climate system, but biases and uncertainty in climate projections remain. The links between tropical circulation and clouds on a wide range of space scales and timescales are complicated by a host of scale interactions that are challenging to effectively represent in comprehensive Earth system models.
Given the complexities of these interactions, theory and idealized models of the tropical atmosphere are important tools for advancing understanding. Idealization allows for experiments that target specific features of interest while retaining the fundamental properties of the tropical climate. Furthermore, the simplicity and flexibility of idealized models provide an opportunity for developing new understanding. One line of idealized modeling has focused on radiative–convective equilibrium (RCE), which approximates the tropical atmosphere as a statistical balance between radiative cooling and convective heating. RCE represents the simplest approximation of not only the tropics but also global climate more broadly. There has been a long history exploring RCE in idealized models for advancing our understanding of the tropics (e.g., Manabe and Strickler, 1964; Held et al., 1993; Tompkins and Craig, 1998; Satoh et al., 2016; Held et al., 2007; Popke et al., 2013).
The RCE model intercomparison project (RCEMIP; Wing et al., 2018) was a coordinated international effort to standardize the boundary conditions and forcing assumptions of RCE simulations to allow for the intercomparison across various atmospheric model types. RCEMIP included cloud-resolving models (CRMs), large-eddy simulation (LES) models, global cloud-resolving models (GCRMs), and general circulation models (GCMs). The initial phase of RCEMIP (hereafter RCEMIP-I) focused on changes in clouds and convective activity with surface warming, as well as associated cloud feedbacks, the implications for climate sensitivity, and the role of aggregation in tropical climate. Results of RCEMIP-I have revealed robust behaviors across the model hierarchy. For example, in response to sea surface temperature (SST) warming, deep-convective anvil clouds rise, slightly warm, thin, and decrease in extent (Wing et al., 2020a; Stauffer and Wing, 2022, 2023), and the large-scale circulation weakens (Silvers et al., 2023). Self-aggregation occurs across nearly all large-domain RCEMIP-I simulations, and it consistently warms and dries the mean state with a reduction in the extent of high clouds (Wing et al., 2020a; Stauffer and Wing, 2022). Cloud-radiative feedbacks are consistently found to be the most important mechanism driving self-aggregation (Pope et al., 2023), and variability in the degree of self-aggregation modulates cloud feedbacks and climate sensitivity (Becker and Wing, 2020; Stauffer, 2023). The fact that these common behaviors emerged across the wide variety of model numerics and physics present in RCEMIP indicates that they are a result of fundamental physical mechanisms. However, the intercomparison also revealed substantial disagreement in the representation of mean profiles of temperature, humidity, and cloudiness; a wide range of static stability and climate sensitivities; a large variation in the degree of convective aggregation; and no consensus in the response of aggregation to warming (Wing et al., 2020a; Becker and Wing, 2020; Wing and Singh, 2024; Silvers et al., 2023).
RCEMIP-I prescribed homogeneous thermal forcings at the boundaries, which consisted of uniform SST and insolation. The simplification of the boundary conditions demonstrates how sensitive the simulation of convection and its complex interactions with radiation and circulation are to model design. The divergent behavior in RCEMIP-I reveals dependencies on representations of convection, microphysics, turbulence, and dynamical cores that may have been masked in other intercomparisons by dynamical constraints. To build on this work, we propose a second phase of RCEMIP (hereafter RCEMIP-II) in which standardized heterogeneities in the SST are prescribed to explore robust behaviors in the tropical system when strict RCE is relaxed to include more realistic, but still idealized, circulations.
A “mock-Walker” circulation (Raymond, 1994) confines deep convection to regions of low-level convergence and large-scale ascent, as determined by the SST pattern (Lindzen and Nigam, 1987; Bretherton and Sobel, 2002; Back and Bretherton, 2009). The general features resemble observed tropical circulations (Grabowski et al., 2000; Larson and Hartmann, 2003), reflecting an interplay of convection, radiation, humidity, and large-scale circulation that is fundamental to tropical deep-convective regions (Bretherton et al., 2006). While the SST pattern fixes the location of the circulation and constrains its strength, interactions between convection, surface gustiness, water vapor, and radiation modify the large-scale flow and the strength and spatial extent of convection (Tompkins, 2001; Liu and Moncrieff, 2008; Wofsy and Kuang, 2012; Silvers and Robinson, 2021). Cloud-radiative feedbacks narrow the convective region in a manner reminiscent of self-aggregation (Grabowski et al., 2000; Bretherton and Sobel, 2002; Liu and Moncrieff, 2008), similar to the zonal contraction of convection in spherical simulations with zonally uniform and meridionally varying SST (Müller and Hohenegger, 2020).
The prescribed SST gradient in the RCEMIP-II mock-Walker simulations will drive a large-scale circulation that provides a partial dynamical constraint on the structure of convection compared to the strict RCE set-up in RCEMIP-I in which there were no external constraints on the location or spatial pattern of convection. The “forced” aggregation in RCEMIP-II provides a common null expectation for the convective structure while still allowing for “un-forced” self-aggregation intrinsic to the model to emerge via radiative–convective feedbacks (e.g., Grabowski et al., 2000; Bretherton and Sobel, 2002; Liu and Moncrieff, 2008; Müller and Hohenegger, 2020). By varying the strength of the SST gradient with warming, the impact of aggregation on climate sensitivity will be able to be attributed to the forced vs. unforced aggregation. An additional motivation for the simulations described below is that mock-Walker simulations are a valuable component of the model hierarchy for tropical dynamics (Jeevanjee et al., 2017) and can be performed by both CRMs and GCM-like models in limited-area planar domains and GCMs on the global sphere. By changing just one parameter compared to the RCEMIP-I simulations (the analytic SST boundary condition), the proposed RCEMIP-II simulations represent a single, clean move up the model hierarchy of model complexity that is closer to the observed tropics than strict RCE.
This paper is organized as follows. Section 2 discusses the motivation for the experimental design of RCEMIP-II in more detail and provides an overview of the scientific objectives. The experimental design, which builds on RCEMIP-I to add an SST pattern to enable an intercomparison of idealized tropical circulations across the model hierarchy, is described in detail in Sect. 3. Section 4 presents sample results for RCEMIP-II to demonstrate the set-up, and Sect. 5 provides some next steps and brief discussion.
RCEMIP was designed to explore questions about tropical clouds, convective aggregation, and climate, with a particular focus on the following three scientific themes:
-
the robustness of the simulated mean state across the spectrum of models;
-
the response of convective clouds to warming and climate sensitivity; and
-
the dependence of convective self-aggregation on temperature.
While RCEMIP-I dictated common domain configurations, grid spacing, trace gas concentrations, insolation, and SST boundary conditions, it purposefully minimized the code changes (such as modifications to physical schemes) needed to run the simulations, both to ensure broad participation by the community and to reveal the true spread in RCE states using a model's “out-of-the-box” suite of physical schemes (Wing et al., 2018). As summarized in the prior section, several common behaviors emerged from RCEMIP-I despite the great diversity in model physics and numerics, providing strong evidence for fundamental physical mechanisms that are not dependent on the details of physical parameterizations. The model diversity also provided an opportunity to use simple theory to explain the intermodel spread (Wing and Singh, 2024). However, two other points that stand out in considering the RCEMIP-I results are (1) how strikingly large the spread in simulated RCE states is and (2) how strong of an imprint convective self-aggregation has on the climate state (Wing et al., 2020a). In particular, the wide range in the degree of self-aggregation and the lack of consensus in its temperature dependence are barriers to understanding.
Our vision for RCEMIP was always that the initial simulations would serve as a starting point, but deep understanding would require performing additional simulations to address issues such as the robustness of the results to experimental design, the sensitivity to model physics and dynamics, and the impact of other factors such as ocean–atmosphere interactions or rotation. In considering possibilities for RCEMIP-II, we sought a protocol that follows the following four principles in the spirit of the design of RCEMIP-I:
-
the ability to directly compare limited-area models with explicit convection and global climate models with parameterized convection;
-
ease of implementation, to encourage the broadest possible participation;
-
a continued investigation of the above three themes of RCEMIP while moving a step up the model hierarchy of complexity; and
-
the provision of an external constraint on convection.
The proposed mock-Walker configuration follows the above philosophy. As described in detail below, it maintains an identical set-up to that of RCEMIP-I with the exception of a simple, prescribed SST pattern. This is easy to implement and allows for direct comparison with RCEMIP-I. Care was taken to maintain consistency between the CRM and GCM domains as much as possible. These characteristics satisfy principles 1 and 2.
From its inception, RCEMIP has been motivated in part by a desire to better understand how the balance between convection and radiation interacts with large-scale circulations (Wing et al., 2018). However, the only large-scale circulations present in RCEMIP-I are those generated by self-aggregation (Silvers et al., 2023). One of our motivations for selecting mock-Walker simulations for RCEMIP-II is a desire to broaden the range of dynamical regimes and cloud types that can be simulated by moving one step up the model hierarchy from RCE (principle 3). In order to do so, we need to relax the idealization of the RCE boundary conditions. Instead of uniform SST, a prescribed sinusoidal SST pattern provides a more realistic (i.e., heterogeneous) boundary condition with a clearer tie to observations. Interactions between the convection and a large-scale circulation that is forced by SST anomalies have direct analogs on Earth to the Intertropical Convergence Zone (ITCZ), the Walker circulation, and the Hadley circulation (in a non-rotating context). The presence of subsiding circulations consistently occurring over regions of cooler SST also allows for the possibility of simulations that include stratocumulus clouds. The SST gradient, combined with an overturning circulation, also allows for the possibility of modeling the transition between shallow- and deep-convective clouds. The mock-Walker configuration thus allows investigation into the three themes of RCEMIP in a framework focused on cloud–circulation coupling that is one step up the model hierarchy from RCE, satisfying principle 3. In particular, varying the SST gradient with mean warming will permit an investigation of the “pattern effect”, supporting the second theme of RCEMIP in a new way compared to RCEMIP-I. The pattern of SST warming, that is, whether warming is enhanced in warm regions where deep convection occurs or in cold regions where there is subsidence and low clouds, is known to influence climate feedbacks (Andrews et al., 2022; Fueglistaler and Silvers, 2021). Prior work using GCMs has shown that if warming is focused on subsidence regions with cold SST, then this has little remote effect and results in a negative lapse rate feedback and large, positive low-cloud feedback. However, if the warming is instead enhanced in regions of deep convection, the warming is communicated through the free troposphere, resulting in increased stability that leads to a stronger negative lapse rate feedback and a weaker (less positive/more negative) low-cloud feedback through increases in low-level cloud cover in the cold-SST region (Andrews and Webb, 2018; Ceppi and Gregory, 2017; Dong et al., 2019; Zhou et al., 2016). An RCEMIP-II with mock-Walker simulations will allow us to investigate these cloud feedbacks on climate in a setting with both CRMs and GCMs, satisfying principles 1 and 3.
The prescribed SST gradient in the mock-Walker simulations provides forcing for low-level convergence towards the warmest SSTs (Lindzen and Nigam, 1987; Bretherton and Sobel, 2002; Back and Bretherton, 2009) and will drive a large-scale circulation that provides a dynamical constraint on the location and spatial pattern of convection relative to uniform SSTs, satisfying principle 4. The ability of SST gradients to at least partially dynamically constrain and organize convection is clear in the test simulations described below, as well as numerous previous studies (Grabowski et al., 2000; Tompkins, 2001; Bretherton et al., 2006; Lutsko and Cronin, 2018; Silvers and Robinson, 2021). The extent to which the prescribed SST gradient constrains convection and circulation in an environment of complex interactions among moist convective processes, radiation, and microphysics will be a subject of investigation across the RCEMIP-II ensemble.
Other possibilities were considered for RCEMIP-II, such as rotation, interactive SST with a slab mixed-layer ocean, and simplified physics. However, none of these options satisfy the above four principles as well as the mock-Walker configuration. It is not possible to satisfy principle 1 while adding rotation, since f-plane RCE simulations in a limited-area CRM domain are quite different than realistic rotation on the GCM sphere with uniform thermal forcing (Merlis and Held, 2019; Sobel et al., 2021). There is no obvious CRM analog to the latter other than a global CRM or perhaps a large beta plane, but either would be more computationally expensive and in opposition to principle 2. While there are many important questions that could be investigated in simulations with rotation, they are likely to be tropical-cyclone-focused rather than supporting the current themes of RCEMIP, thus opposing principle 3. Silvers et al. (2024) provide an example of rotating RCE experiments using a GCM in the context of RCEMIP-I as the control configuration. Simulations with interactive SST involve jumping much further up the model hierarchy and add significant complexity, which complicates interpretation (Coppin and Bony, 2017, 2018; Drotos et al., 2020; Hartmann and Dygert, 2022; Dygert and Hartmann, 2023). Slab mixed-layer oceans of even relatively shallow depth take many hundreds of days to equilibrate (Cronin and Emanuel, 2013), greatly increasing the computational expense and opposing principle 2. Simplified physics (i.e., microphysics and radiation) is more challenging to implement in most models, opposing principle 2, and would move down the hierarchy towards more idealization and could remove some phenomena of interest, opposing principle 3.
By utilizing a mock-Walker configuration, RCEMIP-II will satisfy all four of the above principles and will facilitate deeper understanding of convective aggregation and its role in climate – including the hydrological cycle, cloud feedbacks, and climate sensitivity.
The experimental design of the RCEMIP-II simulations follows the philosophy set out by RCEMIP-I (Wing et al., 2018) and is in line with the principles described above – that is, a small set of experiments designed to maximize the utility of the simulations in answering the above questions while minimizing the effort required of modeling groups. To facilitate comparison with the RCEMIP-I simulations, an identical configuration to that described by Wing et al. (2018) is used, except for the analytic SST boundary condition. Participation in RCEMIP-I, while beneficial for comparison, is not required to participate in RCEMIP-II. The experimental design for the RCEMIP-II mock-Walker simulations is motivated by similar model configurations in recent studies of tropical convection and circulation (e.g., Silvers and Robinson, 2021; Lutsko and Cronin, 2018; Müller and Hohenegger, 2020; Lutsko and Cronin, 2023).
3.1 Basic model set-up
A non-rotating aquaplanet model configuration is to be used (i.e., the Coriolis parameter, f, or Earth's angular velocity, Ω, are set to zero) with no sea ice and no land. Recommended geophysical constants and parameters are provided in Table 1 of Wing et al. (2018).
Models which participated in RCEMIP-I should ideally use the same model version and configuration as RCEMIP-I for their RCEMIP-II mock-Walker simulations to ensure that the SST boundary condition is the only thing that is different. If this is not possible due to model development in the intervening years, the RCE_large300
simulation from RCEMIP-I should be repeated, if possible with the new version of the model, to represent a reference point for comparison.
The RCEMIP-II simulations are to be initialized by the same sounding that was used to initialize a model's RCEMIP-I RCE_large
simulation at the corresponding mean SST. That is, they should be initialized by a horizontally averaged equilibrium sounding from the corresponding RCE_small
simulation.
We welcome participation in RCEMIP-II from models that did not participate in RCEMIP-I. Such models should at minimum complete the RCE_small
simulations described in Wing et al. (2018) to derive a sounding from which to initialize their RCEMIP-II simulations. We encourage but do not require that they complete the RCEMIP-I RCE_large
simulations to serve as a reference.
3.2 Domain configuration
The domain configuration follows the RCE_large
set-up described in Wing et al. (2018) for CRMs, GCMs, or GCRMs, which are reviewed here. Single-column models (SCMs) are not eligible to participate in RCEMIP-II. While it is likely prohibitively computationally expensive for models at LES resolutions (200 m grid spacing was used in RCEMIP-I) to perform simulations on the RCEMIP-II domain, they are welcome to do so if the necessary computational resources are available.
3.2.1 CRMs
CRMs (models with an explicit convection run on a limited-area planar domain) are to employ a three-dimensional domain with doubly periodic lateral boundary conditions. A horizontal grid spacing of 3 km is to be used with an elongated channel geometry of ∼6000 km in the long (x) direction and ∼400 km in the short (y) direction. The domain should be as close to 6000 km long as possible, given the numerical limitations of a given model. The vertical grid should be at least 74 vertical levels, have a model top no lower than 33 km, and have a sponge layer in the top model layers to damp gravity waves. Table 3 in Wing et al. (2018) provides the recommended vertical grid. The simulations are to be performed for at least 200 d, or longer if needed, to reach equilibrium.
3.2.2 GCMs
GCMs (models with parameterized convection run on the global sphere) should employ whichever dynamical core and grid are standard for each model and the horizontal resolution, vertical coordinate, and grid of their Coupled Model Intercomparison Project Phase 6 (CMIP6) configuration. We have chosen not to constrain the GCMs to a common horizontal and vertical grid because the physical parameterizations are sensitive to particular configurations. The simulations are to be performed for at least 1000 d but should be run longer if needed to reach equilibrium.
3.2.3 GCRMs
GCRMs (models with an explicit convection run on a sphere) should ideally be run with the same grid spacing as CRMs (3 km) and the same domain size as GCMs (the real Earth radius, RE). However, a reduced Earth radius may be used to reduce the computational expense. A radius of will yield the correct wavelength of the SST pattern (see below). The simulations are to be performed for at least 200 d, or longer if needed, to reach equilibrium.
3.3 Surface boundary condition
Surface enthalpy fluxes are to be calculated interactively from the resolved surface wind speed and air–sea enthalpy disequilibrium. If allowed by a model's surface layer formulation, a minimum wind speed of 1 m s−1 should be enforced. The lower boundary represents the thermodynamic state of a sea surface that is fixed in time but varies spatially according to a prescribed sinusoidal temperature pattern. We desire to keep the mean SST and the SST gradient (change in SST per unit distance) as consistent as possible between the Cartesian CRM and spherical GCM configurations.
3.3.1 CRMs
For CRMs on the Cartesian channel domain described above,
where 〈SST〉 is the mean SST, ΔSST is the difference between the maximum and minimum SST, x is the horizontal position along the long axis, and Lx is the domain length. This sets the wavelength of the SST pattern, λ, equal to Lx and places the maximum SST at to maintain periodicity at the lateral boundaries.
There are several relevant SST-related quantities that may matter for the climate of mock-Walker simulations. From the perspective of the weak temperature gradient approximation, the absolute SST contrast (ΔSST) is what ought to matter to the dynamics, and the maximum SST (where the deep convection occurs) will set the mean temperature in the free troposphere (Sobel and Bretherton, 2000; Bretherton and Sobel, 2002). Climate changes in response to SST warming are typically referenced to the domain mean SST 〈SST〉. The SST gradient (dSST/dx) might plausibly set horizontal flow speeds, based on Lindzen and Nigam (1987), and if that is the case, the SST Laplacian ought to matter for vertical motion and precipitation. Due to code and numerical optimization and limitations, the CRMs participating in RCEMIP use slightly different domain lengths, Lx. This means that it is not possible to keep all of these parameters fixed across the simulations. We have elected to fix ΔSST across models but to allow λ=Lx for each particular model after testing the alternative options. With slightly different Lx, the SST gradient will vary slightly across models. Over the range of Lx in the RCEMIP-I simulations, the resulting SST gradients differ from that with Lx=6000 km by less than 3 % (except for one model that used Lx=6480 km and which has an SST gradient that is 7.4 % smaller than that with Lx=6000 km). While not ideal, this situation is more desirable than the alternative, which is to set λ=6000 km regardless of the domain length. A fixed λ would lead to slightly different 〈SST〉 across the models, a discontinuous SST distribution at the boundaries (since λ≠Lx), and, most critically, the projection of the prescribed SST forcing onto all scales. The first two are minor issues, but while the vast majority of the power would be in the desired domain-scale wave, λ≠Lx introduces substantial noise at higher wavenumbers (contributing to the majority of the variance of ; not shown). A third option would be to adjust the value of ΔSST to maintain the same SST gradient, but this would cause problematic differences in ΔSST and the maximum SST. Therefore, we elect to use Eq. (1) to define the SST in the CRMs. Keeping ΔSST, maximum SST, and 〈SST〉 consistent is the most elegant and simplest option to implement. It maintains the quantities that plausibly matter most for dynamics, precipitation, and clouds based on weak temperature gradient arguments and avoids non-physical artifacts. We acknowledge that the SST gradients may differ slightly across models if their Lx≠6000 km, which could make it difficult to disentangle the relative contributions of Lx difference from other aspects of model physics and numerics to intermodel differences. While we cannot rule out the differences in Lx as a source of intermodel spread, we anticipate that it will be a small effect based on the test simulations we performed.
3.3.2 GCMs
For GCMs on a sphere using the observed radius of Earth,
where 〈SST〉 is the mean SST, ΔSST is the difference between the maximum and minimum SST, ϕ is latitude in degrees, and λ=54° yields a wavelength of 6004.53 km (when centered on the Equator) to approximately match the CRM configuration. Note that since these simulations are non-rotating, there is no dynamical difference between a wave in latitude and a wave in longitude. The wave in latitude proposed here maintains a consistent distance between each peak that is comparable to that in the CRM domain.
Since the simulations are non-rotating, and the Cartesian CRM domain is doubly periodic, the x and y dimensions are interchangeable, and the CRM domain should be conceptualized as being infinitely repeated in both dimensions (Fig. 1c). This means that, other than the sphericity, it is analogous to a GCM set-up with a meridional circulation forced by a zonally uniform and meridionally varying SST (Fig. 1d) which has also previously been used to study convective aggregation (e.g., Müller and Hohenegger, 2020). Setting the GCM SST pattern with zonal bands of warm and cool SSTs, following Eq. (2), ensures the closest possible match to the CRM domain and avoids any irregularities that could arise from an SST pattern that decreases in width as the poles are approached, as would be the case for a zonally varying SST. Since one of the core principles of RCEMIP is to be able to compare limited-area CRMs and GCMs, we chose this configuration in which the warm SSTs are confined in only one direction in both model types.
3.3.3 GCRMs
For GCRMs on a sphere with a reduced Earth radius of , where RE is the observed radius of Earth,
where 〈SST〉 is the mean SST, ΔSST is the difference between the maximum SST and the minimum SST, and ϕ is latitude in degrees. For km, which yields a radius of , λ=180° corresponds to distance of 6000 km to match the CRM configuration. If a smaller Earth radius of is used, as was used by some GCRMs in RCEMIP-I, λ=180° corresponds to a distance of approximately 5000 km. Smaller Earth radii than this are not recommended.
3.4 Radiative processes
The shortwave and longwave radiative heating rates are to be calculated interactively from the modeled state using a radiative transfer model. Trace gases are to be fixed and spatially uniform, according to Table 2 in Wing et al. (2018). The ozone profile is an analytic approximation of the horizontally uniform equatorial profile derived from the Aqua-Planet Experiment ozone climatology and given by Eq. (1) in Wing et al. (2018). Aerosol effects are to be ignored. The incoming solar radiation is to be spatially uniform and constant in time; there is to be no diurnal or seasonal cycle, and every model grid point should receive the same incident radiation. Following Wing et al. (2018), a reduced solar constant of 551.58 W m−2, a fixed zenith angle of 42.04°, and a fixed surface albedo of 0.07 should be used. The values of the zenith angle and surface albedo are equal to the Equator to 20° and global average insolation-weighted values, respectively (Cronin, 2014).
3.5 Required simulations
The five required RCEMIP-II experiments are listed in Table 1, as are several optional experiments. Figure 1 shows the SST pattern in each experiment.
In selecting the suite of ΔSST values, we took inspiration from, but did not attempt to exactly reproduce, observed SST gradients in the equatorial Pacific and Atlantic. Based on the 1950–2022 HadISST (Rayner et al., 2023) climatological SSTs, averaged between 5.5° N and 5.5° S, there is a difference between the warmest and coldest SSTs of ∼4 K over a distance of 14 344 km in the Pacific (between longitudes of 140.5 and −90.5°) and ∼2 K over a distance of 6783 km in the Atlantic (between longitudes of −50.5 and 10.5°). In both basins, the climatological mean SST is ∼300 K. We desire to use an identical domain configuration as in RCEMIP-I so that the simulations can be compared. This sets λ=6000 km to accommodate a full wavelength with the warmest SSTs at the center of the domain, doubly periodic lateral boundary conditions, and the RCE_large
CRM domain length of Lx∼6000 km. Thus, a difference between the warmest and coldest SSTs of ΔSST=0.83–0.89 K over 3000 km (half a wavelength) would be comparable to the observed SST gradients.
We considered three different values of 〈SST〉 (295, 300, and 305 K) and nine different values of ΔSST (0.625, 0.75, 1.0, 1.25, 1.5, 2.0, 2.5, 3.0, and 5.0 K). Our central value of 〈SST〉=300 K is consistent with current observed mean SSTs in the equatorial Pacific and Atlantic. For the three values of ΔSST at 〈SST〉=300 K (Table 1; Fig. 1), we selected values that resulted in distinctly different weak, moderate, and strong SST gradients, balancing the similarity to observations with choices that have distinct spatial structures of convection and circulation. Our control experiment (MW_300dT1p25
) thus has an SST gradient that is ∼45 % stronger than observed. MW_300dT0p65
and MW_300dT2p5
consider SST gradients that are half and twice as strong as the control; MW_300dT0p65
has an SST gradient that is ∼37 % weaker than observed, while MW_300dT2p5
has an SST gradient that is ∼190 % stronger than observed.
The MW_295dT1p25
and MW_305dT1p25
simulations, along with MW_300dT1p25
, will reveal the effect of uniform warming on the large-scale circulation, convective aggregation, and other features of the simulated climate. This set of simulations can be compared to the corresponding RCE_large
simulations with the same mean SSTs from RCEMIP-I to evaluate the impact of the SST gradient and forced circulation on both the model mean state and response to warming. One slight caveat is that, assuming the weak temperature gradient condition, the domain mean free-tropospheric temperature in the RCEMIP-II simulations will be set not by the mean SST but by the warmest SST (where deep convection is expected to occur). In the RCEMIP-I simulations, the SST is uniform. Thus, there will be a slight difference between the expected mean free-tropospheric temperatures in the RCEMIP-II simulations and the RCEMIP-I simulations at the same mean SST.
![https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f01](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f01-thumb.png)
Figure 1SST pattern for the five required RCEMIP-II experiments for (a) a CRM with domain length km and (b) a GCM with a real Earth radius and λ=54°, following Eqs. (1)–(2). Panel (c) depicts the doubly periodic CRM domain as rotated and tiled 3×24 times to conceptually match the GCM domain in panel (d), where shading in yellow indicates warm SSTs.
The set of required RCEMIP-II simulations will also be used to study the response to warming with the SST gradient held constant, decreased, and increased, respectively. Changing ΔSST under mean warming to amplify or dampen the SST pattern will change the strength of the large-scale circulation and the forced component of aggregation, whereas when ΔSST is kept fixed then only the changes in unforced aggregation will occur. This will facilitate an improved understanding of the modulation of climate and hydrological sensitivity by convective aggregation by having ΔSST as an external parameter that controls the forced component of aggregation. Varying ΔSST under mean warming will also permit us to assess how different models represent the pattern effect on climate feedbacks, which has not previously been examined in a model ensemble that includes models with explicit convection. RCEMIP-II will facilitate investigating the response of tropospheric stability and clouds to different patterns of warming.
We performed test simulations using the System for Atmospheric Modeling (SAM; a CRM), version 6.11.2 (Khairoutdinov and Emanuel, 2013), and the Community Atmosphere Model (CAM; a GCM), version 6 (Danabasoglu et al., 2020). Similar to the paper describing the RCEMIP-I protocol (Wing et al., 2018), we show here sample results from some of these test simulations with SAM and CAM to motivate our choice of required simulations for the RCEMIP-II protocol and as an example of what those simulations might look like; this is not intended as a comprehensive comparison.
The SAM simulations shown in Figs. 2, 7, and 10 use Eq. (1) to set the SST pattern as required for RCEMIP-II with km (Lx=6000 km is not possible due to numerical and code optimization limitations). The SAM simulations used to test the value of ΔSST in Figs. 4 and 6 and to test the domain size in Figs. 11, 12, and 13 were performed before Eq. (1) was set and thus instead used with λ=6000 km and Lx=6144 km. While the details of an individual simulation with λ=6000 km versus λ=6144 km vary, the conclusions regarding the dependence on ΔSST and domain size do not qualitatively depend on this difference. Other than the SST boundary condition, the SAM simulations are identical to their configuration for RCEMIP-I.
The CAM simulations shown all use Eq. (2) with λ=54°, as required by RCEMIP-II. The CAM simulations use the version 6 physics package as in RCEMIP-I (Wing et al., 2020a; Reed et al., 2021), but the model code has been updated to improve some flaws in the cloud microphysics and ice nucleation (Zhu et al., 2022). Comprehensive climate simulations with the updated version of the Community Earth System Model, version 2 (CESM2), result in a lower equilibrium climate sensitivity, a more realistic Last Glacial Maximum, and a weaker shortwave cloud feedback than the version of CESM2 submitted to CMIP6 (Zhu et al., 2022) from which the version of CAM6 used in RCEMIP-I was drawn (Reed et al., 2021). Therefore, we rerun the RCEMIP-I simulations with this updated version of CAM6 as well for comparison.
![https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f02](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f02-thumb.png)
Figure 2Daily mean SST (K), precipitation (mm d−1), vertical motion at 500 hPa (m s−1), outgoing longwave radiation (OLR; W m−2) at day 150, and the moist static energy (MSE; J kg−1) anomaly from the spatial mean and streamfunction averaged over days 100–200 in the MW_300dT1p25
simulation with SAM. Streamfunction is calculated by integrating y-averaged ρw across the x axis of the simulation with contours every 500 kg m−1 s−1.
![https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f03](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f03-thumb.png)
Figure 3SST (K), precipitation (mm d−1), vertical pressure velocity on the 500 hPa pressure level (hPa d−1), and the outgoing longwave radiation (W m−2) averaged over 1 d in the third year of the MW_300dT1p25
simulation with CAM. Note that the color bars are saturated for visualization purposes; the maximum and minimum values are 93.94 and mm d−1 for precipitation, −1263.09 and 118.18 hPa d−1 for vertical velocity, and 102.57 and 315.85 W m−2 for OLR, respectively.
4.1 Convective structure and evolution
Figures 2–3 show the spatial structure in the MW_300dT1p25
simulation. Deep convection and precipitation are generally absent from the regions of coldest SSTs. In SAM, where there is one peak in SST at the center of the domain, the precipitation and low values of outgoing longwave radiation (OLR) associated with cold cloud tops are located in the region of the warmest SSTs and their periphery. There are numerous convective systems within the envelope of warm SST, but there are also convective systems on the flanks of the warm SSTs. If one averages over a longer timescale, however, the strongest mean rising motion occurs over the warmest SSTs where the moist static energy (MSE) is largest, which is in association with an overturning circulation that spans the length of the domain, as shown in the bottom row in Fig. 2. The overturning circulation has both shallow and deep components. This is consistent with prior work, which has shown that mock-Walker simulations may develop stacked overturning circulations (Grabowski et al., 2000; Yano et al., 2002; Larson and Hartmann, 2003; Liu and Moncrieff, 2008; Silvers and Robinson, 2021) due to interactions between radiation and detrained condensate and water vapor (Nuijens and Emanuel, 2018; Sokol and Hartmann, 2022; Lutsko and Cronin, 2023).
In CAM, an SST pattern with three peaks in latitude is used to match the spatial scale of the SST pattern used by SAM. Accordingly, the CAM results show three latitudinal bands of precipitation. However within these latitudinal bands, there are certain longitudes with more precipitation than others. This varies with time (only a single time is shown in the figure) and is similar to the zonal contraction of convection in the simulations of Müller and Hohenegger (2020). While the long-term average circulation in CAM exhibits the ascent over the peaks of SST, there is substantial spatial and short-term variability among the CAM simulations. The longitudinal variations in circulation, including both the ascent and descent within regions of peak SST, are apparent in the precipitation, pressure velocity, and OLR (Fig. 3).
![https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f04](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f04-thumb.png)
Figure 4Hovmöller diagrams of y-averaged column relative humidity in the simulations with 〈SST〉=300 and ΔSST=0.625, 0.75, 1, 1.25, 1.5, 2, 2.5, 3, and 5 K, with SAM. The values of ΔSST that are required for RCEMIP-II are highlighted in red. Note that these SAM simulations, used to test the value of ΔSST, were performed with , with λ=6000 km and Lx=6144 km rather than Eq. (1) and λ=Lx, as required for RCEMIP-II.
![https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f05](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f05-thumb.png)
Figure 5Hovmöller diagrams of column relative humidity averaged over 4° of longitude over the first 200 d of the simulations with 〈SST〉=300 and ΔSST=0.625, 0.75, 1, 1.25, 1.5, 2, 2.5, 3, and 5 K, with CAM. The values of ΔSST that are required for RCEMIP-II are highlighted in red.
Figures 4–5 show the temporal evolution of the column relative humidity in a set of simulations with 〈SST〉=300 K and different values of ΔSST. Column relative humidity is defined as the ratio of the water vapor path to the saturated water vapor path. Localization of convection associated with the development of dry and moist patches emerges within the first 20 d of the simulations, though it takes more than 50 d for the simulations to reach a statistical equilibrium. The localized dry patches are more spatiotemporally intermittent in CAM (Fig. 5) than SAM (Fig. 4), especially at the lower values of 〈SST〉. However, note that Fig. 5 depicts the average across a random 4° longitude in order to average over a region roughly consistent with the SAM domain. Other moist and dry patches are found at other longitudes (and latitudes) in the CAM simulation, though Hovmöller diagrams of the zonal mean column relative humidity are qualitatively similar to Fig. 5, albeit smoother. While only the first 200 d are plotted, the CAM simulation is run for over 1000 d. The structure of the column relative humidity is remarkably consistent with what is shown in Fig. 5 over the full simulation. The only exception is that in the MW_300dT5
, the moist band near the Equator breaks down around day 350 before re-emerging by day 450 (not shown). This could indicate the moist band moving to a different longitude, or it could suggest a temporary reorganization of atmospheric moisture content due to unpredictable convective interactions, as seen in similar simulations by Silvers and Robinson (2021), and is speculated to be a result of parameterized convection.
![https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f06](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f06-thumb.png)
Figure 6Domain mean net shortwave flux at the top of the atmosphere (W m−2) in the simulations with (top row; panels a–c) SAM and (bottom row; panels d–f) CAM at 〈SST〉=300 K. (a, d) Weak (ΔSST=0.625, 0.75, 1 K), (b, e) moderate (ΔSST=1.25, 1.5, 2 K), and (c, f) strong (ΔSST=2.5, 3, 5 K) SST gradients. In each panel, a 5 d running mean is shown with the line thickness increasing with the value of ΔSST. Note that these SAM simulations, used to test the value of ΔSST, were performed with , with λ=6000 km and Lx=6144 km rather than Eq. (1) and λ=Lx, as required for RCEMIP-II.
The differences across different values of ΔSST are discussed below in Sect. 4.2, but here we note that all simulations display rich variability in the structure of the dry and moist regions (and convection). The SAM simulations in particular (Fig. 4) depict disturbances that propagate along the long axis of the domain in both directions, with propagation speeds consistent with convectively coupled gravity waves. There are instances in which a dry band emerges from within a broader moist band, splitting the latter into two moist bands, as well as instances in which two moist bands merge together. There is also lower-frequency variability; for example, the location of the moist and dry regions in the SAM simulations with ΔSST=0.625 K and ΔSST=0.75 K oscillates back and forth with time. There are also fluctuations in the size of the moist, convective region which contracts and expands along the long axis of the domain with a period of ∼30–40 d. This is most prominent in the SAM simulations with ΔSST=2.5 K and ΔSST=3 K but is also seen to a lesser extent in the CAM simulation at those values of ΔSST, as well as both the SAM and CAM simulations at other values of ΔSST.
These low-frequency oscillations can have a substantial influence on the domain mean quantities such as the top-of-atmosphere radiative fluxes (Figs. 6–7). In the SAM simulations with strong SST gradients (Fig. 6c), the OLR and net shortwave flux at the top of the atmosphere vary on ∼30 d timescales with amplitudes of O(10) W m−2. This behavior has also been seen in similar mock-Walker simulations with other CRMs in which the amplitude of the fluctuations can be larger (tens of W m−2) depending on the model, the microphysics settings, and other parameters (Guy Dagan, Andrew Williams, Peter Hill, and Nick Lutsko, personal communication, 2023). Amongst the test simulations examined here with SAM and CAM, it seems that the oscillations are less extreme and less regular in CAM than they are in SAM (Fig. 6). It is not clear whether this is due to the coarser grid spacing and parameterization of convection in CAM or its larger domain and spherical geometry.
![https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f07](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f07-thumb.png)
Figure 7Domain mean (left column; panels a and c) outgoing longwave radiation (OLR; W m−2) and (right column; panels b and d) net shortwave flux at the top of the atmosphere (W m−2) in the simulations required for RCEMIP-II (MW_295dT1p25
, MW_300dT0p625
, MW_300dT1p25
, MW_300dT2p5
, and MW_305dT1p25
) with (top row; panels a and b) SAM and (bottom row; panels c and d) CAM over the first 200 d of the simulation. A 5 d running mean has been applied to all data; note the different axes in the left and right column.
We do not attempt to explain the origins of these low-frequency oscillations here or their dependence on 〈SST〉 and ΔSST. Their physical mechanisms and model dependence are a likely target of analysis across the full suite of RCEMIP-II simulations, once complete, and are thus beyond the scope of this protocol paper. In terms of defining the protocol, however, the presence of these low-frequency oscillations suggests that longer simulations than in RCEMIP-I (which consisted of 100 d simulations for the CRMs) will be needed to achieve robust equilibrium statistics in the RCEMIP-II mock-Walker configuration. Our test simulations indicate that simulations of at least 200 d in length for the CRMs, in which the first 75 d are excluded from the equilibrium averaging period, should be sufficient to average over multiple cycles of the ∼30 d oscillations and achieve robust statistics while still being computationally feasible. If an individual model develops oscillations of an even lower frequency, then longer simulations with that model may be needed.
4.2 Sensitivity to ΔSST
Figures 4–5 indicate that the moist, convective region narrows and that the range of column relative humidity (CRH) values in the domain increases as the SST gradient increases. SAM systematically exhibits a larger range of CRH values than CAM. The simulations tend to divide into three groups, with weak (ΔSST=0.625, 0.75, and 1 K), moderate (ΔSST=1.25, 1.5, and 2 K), and strong (ΔSST=2.5, 3, and 5 K) SST gradients, though these groups are more distinct for SAM than CAM. In the group with weak SST gradients (ΔSST=0.625, 0.75, and 1 K), there are alternating moist and dry regions across the entire domain in the SAM simulation (Fig. 4) reminiscent of the self-aggregation of convection seen in RCE simulations with the same geometry but uniform SST (Wing et al., 2020a). This suggests that ΔSST≤1 K provides such a weak SST gradient that the spatial structure of convection is not influenced by the SST pattern. The circulations set up by the intrinsic self-aggregation are thus stronger than those forced by the SST gradient. This is also apparent in Fig. 8 in which the domain-scale overturning circulation is weak in the simulations with ΔSST≤1 K. Compared to the RCE_large
simulations with SAM from RCEMIP-I in which ΔSST=0 K and self-aggregation generates multiple circulation cells across the domain, the streamfunction in the SAM simulations with ΔSST≤1 K is weaker because the multiple aggregated regions have more variability in their spatial location, and thus their circulations are damped in the time average.
In the group with moderate SST gradients (ΔSST=1.25, 1.5, and 2 K), there is typically one moist, convecting region centered over the warmest SSTs with some variability in its spatial extent (Figs. 4–5). In these simulations with SAM, the SST gradient has forced a domain-spanning overturning circulation (left column in Fig. 8). In the group with strong SST gradients (ΔSST=2.5, 3, and 5 K), stronger low-frequency oscillations emerge (particularly in SAM), and the convective region is increasingly narrow. In SAM, the overturning circulation and MSE anomalies are strongest in these simulations (left column in Fig. 8). A tight connection between rising branches of the overturning circulations and positive MSE anomalies is apparent across all SAM simulations. In simulations with strong and moderate SST gradients, the positive MSE anomalies are generally co-located with the warmest SSTs at the center of the domain, though off-centered positive anomalies are found in a few simulations. The different circulation structures across different values of ΔSST collapse if one instead considers a moisture-sorted circulation (Bretherton et al., 2005; Muller and Held, 2012). To compute the moisture-sorted circulation, CRH is smoothed to a grid scale of 60 km, w and MSE are averaged over each percentile of CRH (in 1 % bins), and ρw is integrated from the 0th to 100th percentile of CRH. All simulations exhibit a moisture-sorted circulation from dry to moist regions with both shallow and deep components (right column in Fig. 8). The shallow overturning circulation has been shown to contribute to convective self-aggregation (Muller and Held, 2012; Muller and Bony, 2015; Coppin and Bony, 2015). Though the differences are muted, the strength of the moisture-sorted circulation tends to be slightly weaker with larger ΔSST or larger 〈SST〉 (not shown), which is consistent with Silvers et al. (2023). This is in contrast with the increase in the strength of the physical circulation in SAM with ΔSST in the left column in Fig. 8. The time mean physical circulation represents the circulation forced by the prescribed SST gradient, while the moisture-sorted circulation emphasizes transport between dry and moist regions regardless of how they are organized in space.
![https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f08](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f08-thumb.png)
Figure 8Streamfunction and moist static energy (MSE; shading; J kg−1) in the SAM simulations. MSE is plotted as anomalies from the spatial mean. The left column shows the streamfunction in physical space, as in Fig. 2, with contours every 500 kg m−1 s−1. The right column shows the streamfunction in moisture space with contours every 0.05 kg % m−2 s−1. The top row shows the RCE_large300
simulations from RCEMIP-I with uniform SST (ΔSST=0 K), averaged over days 75–100, while the other rows show mock-Walker simulations with various ΔSST, averaged over days 100–200. Note that these SAM simulations, used to test the value of ΔSST, were performed with , with λ=6000 km and Lx=6144 km rather than Eq. (1) and λ=Lx, as required for RCEMIP-II.
![https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f09](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f09-thumb.png)
Figure 9Streamfunction and moist static energy (MSE; shading; J kg−1) in the CAM simulations. MSE is plotted as anomalies from the spatial mean. The streamfunction is plotted with contours every 4×1010 kg s−1, following the typical convention with global models. Positive and negative streamfunction values are represented by solid and dashed lines with the zero contour omitted. The top row shows the RCE_large300
simulations from RCEMIP-I with uniform SST (ΔSST=0 K). The remaining rows show mock-Walker simulations with various ΔSST. All panels show data averaged over 25 d.
As in the SAM simulations, the MSE anomalies grow stronger with increasing ΔSST, and the regions of vertical ascent are centered over the peaks of SST for all of the CAM simulations (Fig. 9). All of the mock-Walker simulations have stronger overturning circulations than the RCE_large
simulation from RCEMIP-I (Fig. 9; top panel). Pairs of opposite-signed streamfunction centered over the SST maxima are most well-defined in the simulations with the strongest SST gradients, but systematic variations in the streamfunction as ΔSST changes are less clear in CAM relative to SAM. We think this is in part due to the large regions of descending motion that can form within a region of peak SST (as in Fig. 3). It is also possible that the overturning circulations from each of the three peak SST regions interact with each other to either strengthen or weaken the mean circulation. It should be noted that the stacked circulations seen in SAM and previous literature with CRMs do not appear in the CAM circulations.
The low-frequency oscillations discussed above are quite regular in the SAM simulations with ΔSST=2.5 and 3 K but are more extreme and more irregular with ΔSST=5 K. The MW_300dT5
SAM simulation exhibits oscillations in the domain mean net shortwave flux at the top of the atmosphere of up to ∼20 W m−2, but the magnitude varies (Fig. 6c). This in part motivates us to avoid the use of such strong SST gradients in our selection of the required RCEMIP-II simulations. In the CAM simulations, the low-frequency oscillations are less prevalent and have less dependence on ΔSST (Fig. 6). Another reason to exclude simulations with ΔSST=5 K is that they, along with other simulations that have very warm SSTs in at least part of the domain (such as those with 〈SST〉=305 K), take longer to reach equilibrium – at least in SAM (Figs. 6–7). MW_305dT1p25
is a required simulation, but at least in this particular model, simulations longer than 200 d may be required for it to reach equilibrium.
Since the response of clouds to warming is one of the themes of RCEMIP, we also examine the sensitivity of the cloud amount to the prescribed ΔSST. Figure 10 shows cloud fraction profiles for ΔSST=0.625, 1.25, and 2.5 K, as well as the corresponding RCE_large
simulations from RCEMIP-I, with uniform SST (ΔSST=0 K) for reference. The dependencies of the cloud fraction on ΔSST generally hold across the entire set of ΔSST values examined and not just those shown in Fig. 10. In SAM, increasing ΔSST increases the low-cloud fraction due to the presence of colder SSTs and stronger subsidence in the cold region. A more nuanced examination of clouds using the International Satellite Cloud Climatology Project (ISCCP) histograms demonstrated that the SAM mock-Walker simulations have more numerous, thicker, and deeper low clouds (perhaps indicative of stratocumulus) than their RCEMIP-I counterparts with uniform SST (Stauffer, 2023). Increasing ΔSST also tends to decrease the high-cloud fraction. In CAM, ΔSST has little effect on low clouds, except perhaps in the simulations at 〈SST〉=295 K (Fig. 10d). As in SAM, increasing ΔSST decreases the high-cloud fraction in CAM (Fig. 10), except in those simulations with 〈SST〉=305 K (Fig. 10f), which have much larger high-cloud fractions (∼0.8) that tend to increase with increasing ΔSST (not shown). CAM6 is known to generate a lot of thin, high clouds over warm SSTs in RCEMIP-I (Reed et al., 2021), which is also seen in RCEMIP-II.
![https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f10](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f10-thumb.png)
Figure 10Domain and time mean profiles of cloud fraction, excluding the first 100 d of simulation, for simulations with 〈SST〉 of (left) 295, (middle) 300, and (right) 305 K. The corresponding RCE_large
simulations from RCEMIP-I with uniform SST (ΔSST=0 K) are plotted as the solid black line. The line style indicates the value of ΔSST (see the legend).
The response of clouds to warming is similar across all values of ΔSST (Fig. 10). The upward shift and decrease in high-cloud fraction are robust and consistent with the response in the RCEMIP-I simulations with uniform SST (Stauffer and Wing, 2022). There is some suggestion of an increase in the low-cloud fraction with warming, particularly in SAM. A detailed calculation of the cloud feedback using cloud-radiative kernels and its decomposition into contributions from changes in cloud amount, altitude, and optical depth indicated that the SAM mock-Walker simulations have cloud feedbacks of the same sign but larger magnitudes than their RCEMIP-I counterparts (Stauffer, 2023).
4.3 Sensitivity to domain size
The test simulations with SAM described above follow the domain configuration for CRMs specified in Sect. 3.2.1. The domain is the same as that used for the SAM RCE_large
simulations in RCEMIP-I; there are 2048×128 grid points in the horizontal with a grid spacing of 3 km, resulting in a horizontal domain that is 6144 km × 384 km or an aspect ratio of 16:1. There are 74 vertical levels, following Table 3 of Wing et al. (2018). The experimental design of RCEMIP-II calls for the domain configuration to be the same as in RCEMIP-I to allow for a clean comparison, as well as keep the simulations computationally inexpensive. However, the results could be sensitive to both the long and short dimensions of the domain. For example, prior mock-Walker studies have found that the structure of precipitation and vertical profiles of cloudiness are sensitive to the long dimension of the domain (Silvers and Robinson, 2021; Bretherton et al., 2006). In order to provide some context for the results in our chosen domain configuration, we performed a few test simulations with SAM with different domain sizes: (1) a wide
domain, which is twice as big in the short dimension and has an aspect ratio of 8:1 (2048×256 grid points; 6144 km × 768 km); (2) a long
domain, which is twice as big in the long dimension and has an aspect ratio of 32:1 (4096×128 grid points; 12 288 km × 384 km); and (3) a long and wide (longwide
) domain, which is twice as big in both dimensions and has the same 16:1 aspect ratio as the control simulation (4096×256 grid points; 12 288 km × 768 km). All simulations used to test the domain size utilize , with λ=6000 km and Lx=6144 km rather than Eq. (1) and λ=Lx, as required for RCEMIP-II. The long
and longwide
simulations include two wavelengths of the SST pattern but maintain the same SST gradient as the control domain. We conduct these domain size sensitivity tests for the MW_300dT1p25
and MW_300dT2p5
simulations (longwide
is only performed for the MW_300dT1p25
simulation).
![https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f11](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f11-thumb.png)
Figure 11Hovmöller diagrams of y-averaged column relative humidity in the MW_300dT1p25
simulation with SAM for (a) the standard domain size, (b) a domain that is twice as wide in the y dimension, (c) a domain that is twice as long in the x dimension, and (d) a domain that is twice as wide and twice as long. Note that these SAM simulations, used to test the domain size, were performed with , with λ=6000 km and Lx=6144 km rather than Eq. (1) and λ=Lx, as required for RCEMIP-II. The top panel shows the SST pattern. Note that the output for the first 5 simulation days in panels (c) and (d) was corrupted and thus is not plotted.
![https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f12](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f12-thumb.png)
Figure 12Hovmöller diagrams of y-averaged column relative humidity in the MW_300dT2p5
simulation with SAM for (a) the standard domain size, (b) a domain that is twice as wide in the y dimension, and (c) a domain that is twice as long in the x dimension. Note that these SAM simulations, used to test the domain size, were performed with , with λ=6000 km and Lx=6144 km rather than Eq. (1) and λ=Lx, as required for RCEMIP-II. The top panel shows the SST pattern.
The MW_300dT1p25wide
and MW_300dT2p5wide
simulations exhibit similar spatial structures of convection to that in their narrower counterparts (MW_300dT1p25
and MW_300dT2p5
), as seen in the CRH field in Figs. 11–12, as well as the precipitation, OLR, and vertical velocity fields (not shown). The wide simulations are able to fit more convective cells across their short dimension, but otherwise the convection is similarly confined to the warmest SSTs and exhibits a similar temporal evolution (Figs. 11a–b and 12a–b). This indicates that, at least in this situation with doubly periodic boundary conditions, mock-Walker simulations with different domain widths but the same domain length exhibit qualitatively similar behaviors. When comparing the MW_300dT1p25longwide
simulation to the MW_300dT1p25long
simulation, the spatial structures are also generally similar, though it takes longer for MW_300dT1p25longwide
to evolve to match the MW_300dT1p25long
simulation.
The MW_300dT2p5long
simulation exhibits similar spatial structures and temporal evolution of convection to that in its shorter counterpart MW_300dT2p5
, except that the pattern is repeated twice (Fig. 12a, c). That is, the moist convecting regions are confined to the regions with warmest SSTs at the middle and left and right edges of the domain. There are a few differences; the central warm region is not quite as uniformly moist in the long
simulation, and the low-frequency oscillations in the size of the moist region are less notable. However, in the first-order simulation, the MW_300dT2p5long
simulation appears to be an extension of the MW_300dT2p5
simulation for another wavelength.
The MW_300dT1p25long
simulation, however, has quite different behavior compared to its shorter counterpart (MW_300dT1p25
). Initially, the moist convecting regions are confined to the regions with the warmest SSTs, but unlike the control MW_300dT1p25
simulation, this is not maintained in the long
simulation (Fig. 11a, c). Instead, convection and precipitation begin to also develop on the periphery of the warm SSTs, and eventually, narrower bands emerge with no preference for occurring over the regions of the domain where the SST is maximized. In this regard, the MW_300dT1p25long
simulation is similar to the simulations with weaker SST gradients (i.e., MW_300dT0p625
) and reminiscent of self-aggregation of the convection with uniform SST. This domain length dependence is similar to that found in the mock-Walker simulations of Silvers and Robinson (2021), though in their longer simulations they utilized the same ΔSST between the center of the domain and the edges, resulting in a weaker SST gradient. Thus, it is somewhat surprising that the breakdown of convection into narrow bands spanning the whole domain occurs in our long simulation in which the SST gradient is unchanged. It suggests that ΔSST=1.25 K is near the threshold for when the circulation forced by the SST gradient is strong enough to overcome that of intrinsic self-aggregation, and when provided with a longer domain with multiple wavelengths of SST, intrinsic self-aggregation outweighs the forced circulation. In the MW_300dT1p25
simulation with the control domain size, there is also some suggestion that convection is beginning to break down into narrower bands right at the end of the simulation (Fig. 11a), so, if given enough time, simulation may also evolve to be dominated by intrinsic self-aggregation. In the MW_300dT1p25longwide
simulation, the moist convecting regions are confined to the regions with warmest SSTs for the first ∼100 d of the simulation, but they eventually break down into narrow bands spanning the whole domain as in MW_300dT1p25long
(Fig. 11c, d).
![https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f13](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-f13-thumb.png)
Figure 13Domain mean (left column; panels a and d) outgoing longwave radiation (OLR; W m−2), (b, e) net shortwave flux at the top of the atmosphere (W m−2), and (right column; panels c and f) precipitable water (mm) in simulations with SAM with 〈SST〉=300 K and (top row; panels a–c) ΔSST=1.25 K and (bottom row; panels d–f) ΔSST=2.5 at different domain sizes. The domain sizes include the standard domain size (thin solid teal line), a domain that is twice as wide (medium solid blue line), a domain that is twice as long (thick solid purple line), and a domain that is both twice as wide and twice as long (thick pink dashed–dotted line). A 5 d running mean has been applied to all data. Note that these SAM simulations, used to test the domain size, were performed with , with λ=6000 km and Lx=6144 km rather than Eq. (1) and λ=Lx, as required for RCEMIP-II.
Despite the differences in convective structures discussed above, the simulations with different domain widths generally have similar mean statistics (Fig. 13). The wide
simulations at both ΔSST=1.25 K and ΔSST=2.5 K have similar mean values and temporal variability to their narrower counterparts. The simulations with different domain length exhibit more differences. The MW_300dT2p5long
simulation has a lower amplitude and somewhat less regular lower-amplitude and low-frequency oscillations than its shorter counterpart though roughly similar mean values. The precipitable water in the MW_300dT1p25long
and MW_300dT1p25longwide
simulations diverges to much lower values than their shorter counterparts, as higher fractions of the domain are occupied by dry areas. Consistent with the lower precipitable water, the OLR also tends to be higher in these simulations.
In summary, a mock-Walker configuration with a specified sinusoidal SST boundary condition is proposed for RCEMIP-II. This is intended to provide a partial constraint on the structure of convection while still allowing for rich interactions between convection, clouds, and circulations. RCEMIP-II will build on the success of RCEMIP-I to facilitate deeper understanding of cloud–circulation coupling and convective aggregation and its role in climate.
After testing several equations for the SST boundary condition and nine different values of ΔSST with both SAM, a CRM, and CAM, a GCM, we selected five experiments to be required for RCEMIP-II (Table 1; Fig. 1), with the experimental design as described in Sect. 3. The protocol is designed to allow for comparison between CRMs on a limited-area Cartesian domain and GCMs on the global sphere. The selection of the required ΔSST values balances the similarity to the observed SST gradients in the equatorial Pacific and Atlantic with choices that, in our test simulations, exhibit different spatial structures of convection and circulation. As was the philosophy for RCEMIP-I, the set of five required simulations for RCEMIP-II was carefully selected to facilitate addressing the scientific objectives while minimizing the computational expense (and thus maximizing participation). The SST boundary condition is the only difference in the set-up between RCEMIP-II and RCEMIP-I. While this enables comparison between RCEMIP-I and RCEMIP-II, we welcome participation in RCEMIP-II from models that did not participate in RCEMIP-I.
In addition to providing example results for the five required RCEMIP-II simulations, we also test the sensitivity to ΔSST and the domain geometry. As ΔSST increases, the moist, convective region narrows and becomes increasingly confined to the warmest SSTs. With small values of ΔSST, there are alternating moist and dry regions across the entire domain reminiscent of the self-aggregation of convection seen in RCE simulations with uniform SST. At larger values of ΔSST, low-frequency oscillations in the extent of the moist, convective region and associated variability in domain mean quantities emerge. While the low-cloud fraction tends to increase and the high-cloud fraction tends to decrease with increasing ΔSST, the response of clouds to warming is similar across all values of ΔSST. Simulations with different domain sizes generally have similar mean statistics. The convective structures are relatively insensitive to the width (short horizontal dimension) of the CRM domain but can be sensitive to the length (long horizontal dimension), depending on the value of ΔSST.
The breakdown to self-aggregated convection that spans the entire domain with longer time or larger length simulations (Sect. 4.3) is a potential concern for the choice of ΔSST=1.25 K as representing the moderate SST gradient regime in the set of required RCEMIP-II simulations. It is possible that in some models, the behavior with ΔSST=1.25 K could look similar to that in ΔSST=0.625 K, as the transition between regimes is likely model-dependent. While we considered instead selecting ΔSST=0.75, 1.5, and 3 K as the set of required simulations, since ΔSST=1.5 K might be more securely in the moderate SST gradient range than ΔSST=1.25 K, differences across models are likely larger than those due to a 0.25 K ΔSST difference. Our chosen ΔSST values of 0.625, 1.25, and 2.5 K result in SST gradients that are more symmetrically weaker and stronger than observed than the alternatives. Furthermore, in RCEMIP-I, SAM tended to be one of the most strongly self-aggregated models (Wing et al., 2020a), which could indicate that the dominance of intrinsic self-aggregation over the SST-forced circulation could hold to higher ΔSST values in SAM than in other models. Indeed, there is no guarantee that ΔSST=0.625 will result in a self-aggregation-like regime in models that have weaker self-aggregation tendencies. While intermodel differences in which regime each ΔSST value belongs to may complicate analysis, identifying where these regime transitions occur across the spectrum of models is a goal of RCEMIP-II. To aid in this, while only the top five simulations listed in Table 1 are required, we encourage participants to conduct simulations at additional ΔSST values near regime transitions if possible, as well as with different 〈SST〉 values, as listed as optional simulations at the bottom of Table 1.
(Heus et al., 2010)(Romps, 2008)(Zhou et al., 2019)(Hohenegger et al., 2023)(Lac et al., 2018)(Stratton et al., 2018)(Cotton et al., 2003)(Khairoutdinov and Randall, 2003)(Morrison et al., 2005)(Morrison and Milbrandt, 2015; Gasparini et al., 2022)(Nishizawa et al., 2015; Sato et al., 2015)(Caldwell et al., 2021; Bogenschutz et al., 2023)(Li and Chen, 2019)(Wu et al., 2019)(Hohenegger et al., 2023)(Satoh et al., 2014)(Khairoutdinov et al., 2022)(Neale et al., 2012)(Danabasoglu et al., 2020)(Roehrig et al., 2020; Voldoire et al., 2019)(Chen and Li, 2024)(Hannah et al., 2020)(E3SM Project, 2024)(Zhao et al., 2018a, b)(Hourdin et al., 2020)(Tatebe et al., 2019)(Randall et al., 2016)(Walters et al., 2019)In addition to the objectives related to the simulated mean state, the response of clouds to warming, the role of convective aggregation in climate described in Sect. 2, and the analysis plans described in Sect. 3.5, there are numerous other avenues of investigation that could be explored in the full RCEMIP-II ensemble. By including both models with explicit convection and those with parameterized convection, RCEMIP-II maintains the ability to determine how behavior depends on the representation of convection. Other possible lines of inquiry include determining the physical mechanisms leading to low-frequency oscillations and their dependence on ΔSST, investigating the development of stacked overturning circulations, and investigating what controls the transition between weak, moderate, and strong SST gradient regimes.
While we selected mock-Walker simulations for RCEMIP-II, this does not preclude other types of experiments from being performed by individuals or small groups of models or being led as full intercomparisons as an offshoot of RCEMIP. Indeed, other community efforts utilizing the RCEMIP set-up to investigate other questions have been proposed. One effort involves repeating the RCEMIP-I simulations with aerosol–cloud interactions (Lorian and Dagan, 2023) to compare the response of clouds to aerosol perturbations at equilibrium under a wide range of SST values (RCEMIP-ACI). Such simulations could also be performed using the RCEMIP-II protocol to investigate aerosol–cloud interactions in the presence of a forced circulation. There has also been some interest in performing RCE (uniform thermal forcing) simulations with global models with rotation (e.g., Silvers et al., 2024), which would generate convectively coupled equatorial waves and tropical cyclones. In this way, the RCEMIP-I and RCEMIP-II protocols are a foundation upon which auxiliary investigations with modified experimental designs can be built. A potential future phase III of RCEMIP could focus on adding a slab mixed layer and interactive SSTs to the RCEMIP-I and RCEMIP-II configurations. While adding substantial complexity, this would allow the convective aggregation and cloud–circulation coupling to be influenced by ocean–atmosphere interactions, which is one of the primary physical processes currently missing from the RCEMIP set-up.
One potential limitation of the mock-Walker RCEMIP-II simulations is that the test simulations, as well as prior work, demonstrate the potential for complex behaviors, for example, the emergence of intrinsic unforced self-aggregation that can overwhelm the forced circulation, low-frequency variability, and stacked overturning circulations. While these are all targets of investigation for RCEMIP-II, it is likely that these behaviors will differ across models and thus may complicate interpretation. For this reason, we considered pairing the required mock-Walker simulations with fully interactive radiation with additional simulations with fixed and horizontally and vertically uniform radiative cooling. Prior work has suggested that this can suppress stacked overturning cells (Grabowski et al., 2000; Wofsy and Kuang, 2012; Sokol and Hartmann, 2022; Lutsko and Cronin, 2023). Horizontally uniform radiative cooling has also been shown to suppress self-aggregation (Wing et al., 2017). However, there are numerous nuances regarding how to best implement this, and the design of such simulations requires testing that is beyond the scope of this current paper. Performing both sets of simulations may also be prohibitively difficult for some models, and we desire to maintain the ability for many groups to participate. Therefore, we do not request or provide a protocol for fixed radiation simulations at this time but will continue to investigate their feasibility and optimal design for inclusion in future work in which such additional simulations would be a valuable complement to the required RCEMIP-II simulations described here.
Table 2 shows the list of models that intend to contribute to RCEMIP-II. It includes 18 models that participated in RCEMIP-I and 9 models that did not. This list may grow with participation from additional modeling groups and scientists across the world. We also welcome multiple configurations of a given model (i.e., the same model with various microphysics schemes). Appendix A details the output specification for RCEMIP-II, which closely follows that of RCEMIP-I with a few changes and additions. The additions are requested to facilitate analysis that was not possible with the RCEMIP-I data and are divided into required and optional output requests.
The RCEMIP-II output request closely follows that of RCEMIP-I, as described by Wing et al. (2018) and its corresponding corrigendum. We highlight here a few variables to pay careful attention to and a few changes and additions.
A1 Output variables
In all tables, the italicized variables are non-standard outputs; all others are standard CMIP6 output. Bolded variables are new compared to RCEMIP-I. The variables with a (–)! symbol are outputs for GCMs only. The output should be “CMOR-ized” such that the output variable names and units are the same as in CMIP6, as listed in the tables below.
In addition to the listed variables, the horizontal coordinates, vertical coordinate, and time coordinate should also be output. The time coordinate should be in units of days since the beginning of the simulation. Output should be submitted on an x–y (for CRMs) or latitude–longitude (for global models) grid. In the vertical, the variables should be on model levels, and the necessary information to compute pressure on model levels should be provided. If your model does not employ pressure levels (i.e., it uses height levels or a type of hybrid level), please also output the domain and time mean values of pressure on your model levels for approximate plotting purposes. Ideally this would be included in all files, but it is especially useful in the 1D files.
CRMs should output all variables, including 3D variables, over the full simulation. GCMs should output 0D, 1D, and 2D variables over the full simulation and 3D variables over last 200 d of simulation. GCRMs should output all variables over the full simulation but only upload 0D, 1D, and 2D variables to the RCEMIP data repository (3D variables should be archived locally). Note that the 3D output request is different than RCEMIP-I.
Table A10D hourly averaged variables (t). Italicized variables are not standard CMIP output. TOA stands for top of atmosphere.
![](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-t03-thumb.png)
Table A32D hourly averaged variables (x, y, and t). Italicized variables are not standard CMIP output. Bolded variables are a new request compared to RCEMIP-I. Variables with a “!” symbol are required only for models with parameterized convection.
![](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-t05-thumb.png)
Table A43D instantaneous 6 h variables (x, y, z, and t). Italicized variables are not standard CMIP output. Bolded variables are a new request compared to RCEMIP-I. Variables with a “!” symbol are required only for models with parameterized convection.
![](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-t06-thumb.png)
Table A5Optional output variables. Italicized variables are not standard CMIP output. Bolded variables are a new request compared to RCEMIP-I. FMSE stands for frozen moist static energy.
![](https://gmd.copernicus.org/articles/17/6195/2024/gmd-17-6195-2024-t07-thumb.png)
Tables A1–A4 list the required output. Table A1 indicates the list of zero-dimensional domain-averaged variables (functions of t only) that are to be computed and output as hourly averages. Table A2 indicates the list of one-dimensional domain-averaged profiles (functions of z and t) that are to be computed and output as hourly averages. Table A3 indicates the list of the two-dimensional variables (functions of x, y, and t) to output as hourly averages. All models should output tabot, uabot, and vabot (the air temperature, eastward wind, and northward wind, respectively, at the lowest model level). Those models which routinely estimate and output 2 m air temperature, 10 m eastward wind, and 10 m northward wind (tas, uas, and vas, respectively) should also output those variables. Models that use height coordinates should output vertical velocity, wa500, whereas models that use pressure-based coordinates should output omega, wap500. Table A4 indicates the list of three-dimensional variables to output as instantaneous 6 h snapshots. Models that use height coordinates should output vertical velocity, wa, and pressure, pa, whereas models that use pressure-based coordinates should output omega, wap, and geopotential height zg.
Table A5 lists optional, but recommended, output which consists of additional 3D variables to be output as 6 h instantaneous snapshots, two 2D variables to be output as 30 min instantaneous snapshots (in contrast to the hourly averages requested in Table A3), and additional 2D variables, including those associated with the frozen moist static energy (FMSE) budget. The 30 min instantaneous output will facilitate the tracking of mesoscale convective systems, and the FMSE budget output will facilitate diagnosing physical mechanisms of convective aggregation. We also add as optional output several microphysical variables which would facilitate the analysis of microphysically driven changes in deep convection, clouds, and climate. Since different microphysics schemes might predict a different set of variables, models should only output those quantities that are applicable to them. We also request tntm, the tendency of air temperature due to microphysical latent heating, as an optional 3D output to facilitate heat budget analysis. This is the term L(c−e), where c is the condensation rate, and e is the evaporation rate. The crvi variable is the 2D mass-weighted vertical integral of the gross condensation rate which should be computed as the mass-weighted vertical integral of the negative part of the microphysical tendency of water vapor. We note that if water vapor (or equivalent) is not a prognostic variable, this might not be feasible. The optional requests for the microphysical variables and the 30 min instantaneous 2D variables are new compared to RCEMIP-I.
A2 Cloud fraction
We request the diagnosis of a global cloud fraction profile (cldfrac_avg) that includes all clouds and is the fraction of the entire domain covered by cloud at a given height (Table A2). This 1D variable should be a function of the vertical level and time. Following Stauffer and Wing (2022), the presence of a cloud in CRMs (or models without a cloud scheme) should be defined, where the mixing ratio of the total cloud condensate (cloud liquid water + cloud ice) is greater than g g−1. Note that this is different from the original RCEMIP-I definition in Wing et al. (2018). GCMs or other models with cloud schemes should continue to provide the cloud fraction as output from the cloud scheme. For GCMs, we also request the output of a total cloud fraction for each grid column as the 2D variable cl, which is a function of x, y, and t (Table A3), as well as the output of cloud fraction for each grid column as a function of height as the 3D variable cldfrac, which is a function of x, y, t, and the vertical level (Table A4). Cloud fraction should be output as a fraction (between 0 and 1) and not a percentage out of 100. The request of the 3D variable cldfrac for GCMs is new compared to RCEMIP-I.
A3 Cloud water variables
Following the conventions of CMIP6 (see http://clipc-services.ceda.ac.uk/dreq/index.html (last access: 1 December 2023) for variable descriptions), several of the cloud water variables have confusing variable names. The 0D variable clwvi_avg (Table A1) and the 2D variable clwvi (Table A3) represent the condensed water path which includes both cloud ice and cloud liquid water. The 1D variable clw_avg (Table A2) and the 3D variable clw (Table A4) represent only the cloud liquid water.
A4 Relative humidity
When computing relative humidity (the 1D variable hur_avg; Table A2), care should be taken to compute the saturation with respect to liquid for temperatures above freezing and with respect to ice for temperatures below freezing. The formulas for saturation vapor pressure should follow those used in a given model's thermodynamics. If the model's thermodynamics interpolates between saturation over liquid and saturation over ice for temperatures near freezing, then that should also be followed.
A5 Radiative heating rates and fluxes
The tendency of air temperature due to shortwave and longwave radiative heating, assuming that clear-sky conditions should be output as 3D variables (tntrscs and tntrlcs), in addition to the total radiative heating tendencies (tntrs and tntrl); see Table A4. The request of the 3D variables tntrscs and tntrlcs is new compared to RCEMIP-I. The domain average shortwave and longwave radiative heating rate profiles should also be output as 1D variables for all-sky and clear-sky conditions (tntrs_avg, tntrl_avg, tntrscs_avg, and tntrlcs_avg; see Table A2). The downwelling and upwelling shortwave and longwave radiative fluxes at the surface and top of atmosphere, and for all-sky and clear-sky conditions, are requested as 2D variables (Table A3).
A6 CFMIP Observational Simulator Package (COSP)
COSP simulator outputs are requested for the ISCCP simulator for GCMs that have COSP available as a diagnostic package. This will facilitate comparison with observations, as well as the calculation and decomposition of cloud feedbacks (Zelinka et al., 2012a, b, 2013). More information about COSP can be found on the COSP website (https://cfmip.github.io, last access: 1 December 2023). The ISCCP simulator provides pseudo-retrievals of cloud-top pressure (CTP) and cloud optical thickness (tau) (Klein and Jakob, 1999; Webb et al., 2001). The following ISCCP simulator outputs are requested: hourly averages of ISCCP 2D diagnostics (cltisccp is for cloud fraction, albisccp is for cloud albedo, and pctisccp is for cloud-top pressure; see Table A3) and 6 h diagnostics of the ISCCP CTP–tau histograms (clisccp is for cloud area percentage in the seven pressure and seven optical depth bins). The clisccp variabe is listed in Table A4 because it is requested every 6 h, but to clarify, we expect one histogram for each horizontal grid point. CRMs are not required to provide ISCCP simulator output, but if they have the capability to run COSP, they may provide the domain mean ISCCP histogram as a comparison to planned offline calculations either as a function of time or averaged over the last 100 d. The request for ISCCP simulator output for GCMs is new compared to RCEMIP-I.
A7 Frozen moist static energy budget
Four optional 2D variables associated with the frozen moist static energy (FMSE) budget are requested in Table A5. If these variables (functions of x, y, and t) are diagnosed online in the model, their values may be output as hourly averages. If they are diagnosed offline from the instantaneous 3D output, they may be provided as instantaneous 6 h snapshots. For example, if tnfmse, the tendency of FMSE, is diagnosed offline, then it should be diagnosed from instantaneous fmse output. If it is diagnosed online from instantaneous variables, then its hourly average can then be output.
FMSE is defined as . The values of , and Lf used by the model formulation should be used to compute h. qice is the mass fraction of all ice-phase condensates (cloud ice, snow, etc.). The mass-weighted vertical integral of FMSE is given by
or, in pressure coordinates,
Care should be taken to make sure the same limits of integration are used at all times/locations. The mass-weighted vertical integral of the advective tendency of FMSE (advfmse) is given by
Ideally, FMSE would be diagnosed online, and each model's advection scheme would be used to advect it. For instance, one could diagnose FMSE from the prognostic variables just before and just after they are advected, and then the difference could be taken as a measure of the FMSE advective tendency. If this is not possible, we ask that groups make their best effort to estimate these terms. The spatial variance in the mass-weighted vertical integral of frozen moist static energy is computed using the squared anomalies from the horizontal mean of the mass-weighted vertical integral of moist static energy (). Its tendency (tnfmsevar) is given by
where ′ indicates an anomaly from the horizontal mean.
Analysis scripts used to generate the figures are available in a Zenodo archive at https://doi.org/10.5281/zenodo.11043720 (Wing, 2024a). The version of the System for Atmospheric Modeling (Khairoutdinov and Randall, 2003) used here is publicly available in a Zenodo archive at https://doi.org/10.5281/zenodo.12626604 (Wing, 2024b). CAM6 is publicly available as part of CESM2 at https://github.com/ESCOMP/CESM/releases/tag/release-cesm2.1.3 (last access: 1 December 2023), and the version used here is publicly available in a Zenodo archive at https://doi.org/10.5281/zenodo.12668076 (Wing et al., 2024). Model configuration files and a subset of the model data needed to reproduce the figures are available in a Zenodo archive at https://doi.org/10.5281/zenodo.10137266 (Wing and Silvers, 2023). We thank the German Climate Computing Center (DKRZ) for hosting the standardized RCEMIP-I data which are publicly available online at http://hdl.handle.net/21.14101/d4beee8e-6996-453e-bbd1-ff53b6874c0e (Wing et al., 2020b). The standardized RCEMIP-II data will also be hosted and made publicly available in the DKRZ archive once all the RCEMIP-II simulations have been completed and submitted.
AAW led the writing of the text. AAW and LGS performed and analyzed simulations. AAW, LGS, and KAR contributed to discussing the goals and specifications of RCEMIP-II and editing the text.
The contact author has declared that none of the authors has any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
The authors thank Tim Cronin, Guy Dagan, Peter Hill, Nick Lutsko, Andrew Williams, and the participants of the RCEMIP breakout session at the 2023 CFMIP–GASS meeting in Paris for helpful discussions. The authors are grateful to Catherine Stauffer for assistance with initial analysis of the test simulations and to Graham O'Donnell for calculating streamfunction in the SAM simulations and generating Figs. 2 and 8. We also thank Peter Bogenschutz, Jean-Pierre Chaboureau, Guy Dagan, Stephan de Roode, Blaž Gasparini, Walter Hannah, Keiichi Hashimoto, Peter Hill, Cathy Hohenegger, Chris Holloway, Fredrik Jansson, Todd Jones, Marat Khairoutdinov, Gabrielle Leung, Cheng Li, Shuhei Matsugishi, Hiroaki Miura, David Randall, Camille Risi, Romain Roehrig, David Romps, Masaki Satoh, Lorenzo Tomassini, Sue van den Heever, Chien-Ming Wu, Tomoro Yanase, and Ming Zhao for agreeing to perform RCEMIP-II simulations with the models listed in Table 2. We also thank Marat Khairoutdinov for maintaining and providing the SAM cloud-resolving model which was used to perform simulations in this paper and two anonymous reviewers for in-depth reviews that led to improvements in the article. We would like to acknowledge high-performance computing support from Cheyenne (https://doi.org/10.5065/D6RX99HX) provided by the National Center for Atmospheric Research's (NCAR) Computational and Information Systems Laboratory and sponsored by the National Science Foundation (CISL, 2019).
Allison A. Wing has been supported by the National Science Foundation (NSF; grant nos. AGS-2140419 and AGS-1830724). Levi G. Silvers and Kevin A. Reed have been supported by the NSF (grant no. AGS-2327958). Kevin A. Reed has also been partially supported by the Department of Energy, Office of Science (grant no. DE-SC0016605; “A Framework for Improving Analysis and Modeling of Earth System and Intersectoral Dynamics at Regional Scales”).
This paper was edited by Travis O'Brien and reviewed by two anonymous referees.
Andrews, T. and Webb, M. J.: The Dependence of Global Cloud and Lapse Rate Feedbacks on the Spatial Structure of Tropical Pacific Warming, J. Climate, 31, 641–654, https://doi.org/10.1175/JCLI-D-17-0087.1, 2018. a
Andrews, T., Bodas-Salcedo, A., Gregory, J. M., Dong, Y., Armour, K. C., Paynter, D., Lin, P., Modak, A., Mauritsen, T., Cole, J. N. S., Medeiros, B., Benedict, J. J., Douville, H., Roehrig, R., Koshiro, T., Kawai, H., Ogura, T., Dufresne, J.-L., Allan, R. P., and Liu, C.: On the Effect of Historical SST Patterns on Radiative Feedback, J. Geophys. Res., 127, e2022JD036675, https://doi.org/10.1029/2022JD036675, 2022. a
Back, L. E. and Bretherton, C. S.: On the relationship between SST gradients, boundary layer winds, and convergence over the tropical oceans, J. Climate, 22, 4182–4196, https://doi.org/10.1175/2009JCLI2392.1, 2009. a, b
Becker, T. and Wing, A. A.: Understanding the extreme spread in climate sensitivity within the Radiative-Convective Equilibrium Model Intercomparison Project, J. Adv. Model. Earth Sy., 12, e2020MS002165, https://doi.org/10.1029/2020MS002165, 2020. a, b
Bogenschutz, P. A., Eldred, C., and Caldwell, P. M.: Horizontal Resolution Sensitivity of the Simple Convection-Permitting E3SM Atmosphere Model in a Doubly-Periodic Configuration, J. Adv. Model. Earth Sy., 15, e2022MS003466, https://doi.org/10.1029/2022MS003466, 2023. a
Bony, S., Stevens, B., Frierson, D. M. W., Jakob, C., Kageyam, M., Pincus, R., Shepherd, T. G., Sherwood, S. C., Siebesma, A. P., Sobel, A. H., Watanabe, M., and Webb, M. J.: Clouds, circulation and climate sensitivity, Nat. Geosci., 8, 261–268, https://doi.org/10.1038/ngeo2398, 2015. a
Bretherton, C. S. and Sobel, A. H.: A simple model of a convectively-coupled Walker circulation using the weak temperature gradient approximation, J. Climate, 15, 2907–2920, 2002. a, b, c, d, e
Bretherton, C. S., Blossey, P. N., and Khairoutdinov, M.: An energy-balance analysis of deep convective self-aggregation above uniform SST, J. Atmos. Sci., 62, 4237–4292, https://doi.org/10.1175/JAS3614.1, 2005. a
Bretherton, C. S., Blossey, P. N., and Peters, M. E.: Interpretation of simple and cloud-resolving simulations of moist-convection-radiation interaction with a mock-Walker circulation, Theor. Comput. Fluid Dyn., 20, 421–442, https://doi.org/10.1007/s00162-006-0029-7, 2006. a, b, c
Caldwell, P. M., Terai, C. R., Hillman, B., Keen, N. D., Bogenschutz, P., Lin, W., Beydoun, H., Taylor, M., Bertagna, L., Bradley, A. M., Clevenger, T. C., Donahue, A. S., Eldred, C., Foucar, J., Golaz, J.-C., Guba, O., Jacob, R., Johnson, J., Krishna, J., Liu, W., Pressel, K., Salinger, A. G., Singh, B., Steyer, A., Ullrich, P., Wu, D., Yuan, X., Shpund, J., Ma, H.-Y., and Zender, C. S.: Convection-Permitting Simulations With the E3SM Global Atmosphere Model, J. Adv. Model. Earth Sy., 13, e2021MS002544, https://doi.org/10.1029/2021MS002544, 2021. a
Ceppi, P. and Gregory, J. M.: Relationship of tropospheric stability to climate sensitivity and Earth's observed radiation budget, P. Natl. Acad. Sci. USA, 114, 13126–13131, https://doi.org/10.1073/pnas.1714308114, 2017. a
Chen, S. and Li, C.: ExoCubed: A Riemann-Solver based Cubed-Sphere Dynamic Core for Planetary Atmospheres, arXiv [preprint], arXiv:2403.06844, 2024. a
CISL: Cheyenne: HPE/SGI ICE XA System (University Community Computing), Boulder, CO, National Center for Atmospheric Research, https://doi.org/10.5065/D6RX99HX, 2019. a
Coppin, D. and Bony, S.: Physical mechanisms controlling the initiation of convective self-aggregation in a General Circulation Model, J. Adv. Model. Earth Sy., 7, 2060–2078, https://doi.org/10.1002/2015MS000571, 2015. a
Coppin, D. and Bony, S.: Internal variability in a coupled general circulation model in radiative-convective equilibrium, Geophys. Res. Lett., 44, 5142–5149, https://doi.org/10.1002/2017gl073658, 2017. a
Coppin, D. and Bony, S.: On the Interplay Between Convective Aggregation, Surface Temperature Gradients, and Climate Sensitivity, J. Adv. Model. Earth Sy., 10, 3123–3138, https://doi.org/10.1029/2018MS001406, 2018. a
Cotton, W. R., Pielke Sr., R. A., Walko, R. L., Liston, G. E., Tremback, C. J., Jiang, H., McAnelly, R. L., Harrington, J. Y., Nicholls, M. E., Carrio, G. G., and McFadden, J. P.: RAMS 2001: Current status and future directions, Meteorol. Atmos. Phys., 82, 5–29, https://doi.org/10.1007/s00703-001-0584-9, 2003. a
Cronin, T. W.: On the choice of average solar zenith angle, J. Atmos. Sci., 71, 2994–3003, https://doi.org/10.1175/JAS-D-13-0392.1, 2014. a
Cronin, T. W. and Emanuel, K. A.: The climate time scale in the approach to radiative-convective equilibrium, J. Adv. Model. Earth Sy., 5, 843–849, https://doi.org/10.1002/jame.20049, 2013. a
Danabasoglu, G., Lamarque, J.-F., Bacmeister, J., Bailey, D. A., DuVivier, A. K., Edwards, J., Emmons, L. K., Fasullo, J., Garcia, R., Gettelman, A., Hannay, C., Holland, M. M., Large, W. G., Lauritzen, P. H., Lawrence, D. M., Lenaerts, J. T. M., Lindsay, K., Lipscomb, W. H., Mills, M. J., Neale, R., Oleson, K. W., Otto-Bliesner, B., Phillips, A. S., Sacks, W., Tilmes, S., van Kampenhout, L., Vertenstein, M., Bertini, A., Dennis, J., Deser, C., Fischer, C., Fox-Kemper, B., Kay, J. E., Kinnison, D., Kushner, P. J., Larson, V. E., Long, M. C., Mickelson, S., Moore, J. K., Nienhouse, E., Polvani, L., Rasch, P. J., and Strand, W. G.: The Community Earth Ssytem Model Version 2 (CESM2), J. Adv. Model. Earth Sy., 12, e2019MS001916, https://doi.org/10.1029/2019MS001916, 2020. a, b
Dong, Y., Proistosescu, C., Armour, K. C., and Battisti, D. S.: Attributing Historical and Future Evolution of Radiative Feedbacks to Regional Warming Patterns using a Green's Function Approach: The Preeminence of the Western Pacific, J. Climate, 32, 5471–5491, https://doi.org/10.1175/JCLI-D-18-0843.1, 2019. a
Drotos, G., Becker, T., Mauritsen, T., and Stevens, B.: Global variability in radiative-convective equilibrium with a slab ocean under a wide range of CO2 concentrations, Tellus A, 72, 1–19, https://doi.org/10.1080/16000870.2019.1699387, 2020. a
Dygert, B. D. and Hartmann, D. L.: The Cycle of Large-Scale Aggregation in Tropical Radiative–Convective Equilibrium, J. Geophys. Res., 128, e2022JD037302, https://doi.org/10.1029/2022JD037302, 2023. a
E3SM Project: Energy Exascale Earth System Model (E3SM), E3SM [computer software], https://doi.org/10.11578/E3SM/dc.20240301.3, 2024. a
Fueglistaler, S. and Silvers, L.: The Peculiar Trajectory of Global Warming, J. Geophys. Res., 126, e2020JD033629, https://doi.org/10.1029/2020JD033629, 2021. a
Gasparini, B., Sokol, A. B., Wall, C. J., Hartmann, D. L., and Blossey, P. N.: Diurnal Differences in Tropical Maritime Anvil Cloud Evolution, J. Climate, 35, 1655–1677, 2022. a
Grabowski, W. W., Yano, J.-I., and Moncrieff, M. W.: Cloud resolving modeling of tropical circulations driven by large-scale SST gradients, J. Atmos. Sci., 57, 2022–2040, https://doi.org/10.1175/1520-0469(2000)057<2022:CRMOTC>2.0.CO;2, 2000. a, b, c, d, e, f
Hannah, W. M., Jones, C. R., Hillman, B. R., Norman, M. R., Bader, D. C., Taylor, M. A., Leung, L. R., Pritchard, M. S., Branson, M. D., Lin, G., Pressel, K. G., and Lee, J. M.: Initial Results From the Super-Parameterized E3SM, J. Adv. Model. Earth Sy., 12, e2019MS001863, https://doi.org/10.1029/2019MS001863, 2020. a
Hartmann, D. L. and Dygert, B. D.: Global Radiative Convective Equilibrium With a Slab Ocean: SST Contrast, Sensitivity and Circulation, J. Geophys. Res., 127, e2021JD036400, https://doi.org/10.1029/2021JD036400, 2022. a
Held, I. M., Hemler, R. S., and Ramaswamy, V.: Radiative-convective equilibrium with explicit two-dimensional moist convection, J. Atmos. Sci., 50, 3909–3927, 1993. a
Held, I. M., Zhao, M., and Wyman, B.: Dynamic radiative-convective equilibria using GCM column physics, J. Atmos. Sci., 64, 228–238, https://doi.org/10.1175/JAS3825.11, 2007. a
Heus, T., van Heerwaarden, C. C., Jonker, H. J. J., Pier Siebesma, A., Axelsen, S., van den Dries, K., Geoffroy, O., Moene, A. F., Pino, D., de Roode, S. R., and Vilà-Guerau de Arellano, J.: Formulation of the Dutch Atmospheric Large-Eddy Simulation (DALES) and overview of its applications, Geosci. Model Dev., 3, 415–444, https://doi.org/10.5194/gmd-3-415-2010, 2010. a
Hohenegger, C., Korn, P., Linardakis, L., Redler, R., Schnur, R., Adamidis, P., Bao, J., Bastin, S., Behravesh, M., Bergemann, M., Biercamp, J., Bockelmann, H., Brokopf, R., Brüggemann, N., Casaroli, L., Chegini, F., Datseris, G., Esch, M., George, G., Giorgetta, M., Gutjahr, O., Haak, H., Hanke, M., Ilyina, T., Jahns, T., Jungclaus, J., Kern, M., Klocke, D., Kluft, L., Kölling, T., Kornblueh, L., Kosukhin, S., Kroll, C., Lee, J., Mauritsen, T., Mehlmann, C., Mieslinger, T., Naumann, A. K., Paccini, L., Peinado, A., Praturi, D. S., Putrasahan, D., Rast, S., Riddick, T., Roeber, N., Schmidt, H., Schulzweida, U., Schütte, F., Segura, H., Shevchenko, R., Singh, V., Specht, M., Stephan, C. C., von Storch, J.-S., Vogel, R., Wengel, C., Winkler, M., Ziemen, F., Marotzke, J., and Stevens, B.: ICON-Sapphire: simulating the components of the Earth system and their interactions at kilometer and subkilometer scales, Geosci. Model Dev., 16, 779–811, https://doi.org/10.5194/gmd-16-779-2023, 2023. a, b
Hourdin, F., Rio, C., Grandpeix, J.-Y., Madeleine, J.-B., Cheruy, F., Rochetin, N., Jam, A., Musat, I., Idelkadi, A., Fairhead, L., Foujols, M.-A., Mellul, L., Traore, A.-K., Dufresne, J.-L., Boucher, O., Lefebvre, M.-P., Millour, E., Vignon, E., Jouhaud, J., Bint Diallo, F., Lott, F., Gastineau, G., Caubel, A., Meurdesoif, Y., and Ghattas, J.: LMDZ6A: The atmospheric component of the ISPL climate model with improved and better tuned physics, J. Adv. Model. Earth Sy., 12, e2019MS001892, https://doi.org/10.1029/2019MS001892, 2020. a
Jeevanjee, N., Hassanzadeh, P., Hill, S. A., and Sheshadri, A.: A perspective on climate model hierarchies, J. Adv. Model. Earth Sy., https://doi.org/10.1002/2017MS001038, 2017. a
Khairoutdinov, M. F. and Emanuel, K.: Rotating radiative-convective equilibrium simulated by a cloud-resolving model, J. Adv. Model. Earth Sy., 5, 816–825, https://doi.org/10.1002/2013MS000253, 2013. a
Khairoutdinov, M. F. and Randall, D.: Cloud resolving modeling of the ARM Summer 1997 IOP: Model formulation, results, uncertainties, and sensitivities, J. Atmos. Sci., 60, 607–625, 2003. a, b
Khairoutdinov, M. F., Blossey, P. N., and Bretherton, C. S.: Global System for Atmospheric Modeling: Model Description and Preliminary Results, J. Adv. Model. Earth Sy., 14, e2021MS002968, https://doi.org/10.1029/2021MS002968, 2022. a
Klein, S. A. and Jakob, C.: Validation and senstiivities of frontal clouds simulated by the ECMWF model, Mon. Weather Rev., 127, 2514–2531, 1999. a
Lac, C., Chaboureau, J.-P., Masson, V., Pinty, J.-P., Tulet, P., Escobar, J., Leriche, M., Barthe, C., Aouizerats, B., Augros, C., Aumond, P., Auguste, F., Bechtold, P., Berthet, S., Bielli, S., Bosseur, F., Caumont, O., Cohard, J.-M., Colin, J., Couvreux, F., Cuxart, J., Delautier, G., Dauhut, T., Ducrocq, V., Filippi, J.-B., Gazen, D., Geoffroy, O., Gheusi, F., Honnert, R., Lafore, J.-P., Lebeaupin Brossier, C., Libois, Q., Lunet, T., Mari, C., Maric, T., Mascart, P., Mogé, M., Molinié, G., Nuissier, O., Pantillon, F., Peyrillé, P., Pergaud, J., Perraud, E., Pianezze, J., Redelsperger, J.-L., Ricard, D., Richard, E., Riette, S., Rodier, Q., Schoetter, R., Seyfried, L., Stein, J., Suhre, K., Taufour, M., Thouron, O., Turner, S., Verrelle, A., Vié, B., Visentin, F., Vionnet, V., and Wautelet, P.: Overview of the Meso-NH model version 5.4 and its applications, Geosci. Model Dev., 11, 1929–1969, https://doi.org/10.5194/gmd-11-1929-2018, 2018. a
Larson, K. and Hartmann, D. L.: Interactions among cloud, water vapor, radiation, and large-scale circulation in the tropical climate. part ii: Sensitivity to spatial gradients of sea surface temperature, J. Climate, 16, 1441–1455, 2003. a, b
Li, C. and Chen, X.: Simulating nonhydrostatic atmospheres on planets (SNAP): Formulation, validation, and application to the Jovian atmosphere, The Astrophysical Journal Supplement Series, 240, 37, https://doi.org/10.3847/1538-4365/aafdaa, 2019. a
Lindzen, R. S. and Nigam, S.: On the role of sea surface temperature gradients in forcing low-level winds and convergence in the tropics, J. Atmos. Sci., 44, 2418–2436, 1987. a, b, c
Liu, C. and Moncrieff, M. W.: Explicitly simulated tropical convection over idealized warm pools, J. Geophys. Res., 113, D21121, https://doi.org/10.1029/2008JD010206, 2008. a, b, c, d
Lorian, S. and Dagan, G.: On the sensitivity of aerosol-cloud interactions to changes in sea surface temperature in radiative-convective equilibrium, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2023-2096, 2023. a
Lutsko, N. J. and Cronin, T.: Mock-Walker Simulations: Mean Climates, Responses to Warming and Transition to Double-Cell Circulations, ESS Open Archive, https://doi.org/10.22541/essoar.167591087.76248356/v1, 2023. a, b, c
Lutsko, N. J. and Cronin, T. W.: Increase in precipitation efficiency with surface warming in radiative-convective equilibrium, J. Adv. Model. Earth Sy., 10, 2992–3010, https://doi.org/10.1029/2018MS001482, 2018. a, b
Manabe, S. and Strickler, R. F.: Thermal equilibriation of the atmosphere with a convective adjustment, J. Atmos. Sci., 21, 361–385, 1964. a
Merlis, T. M. and Held, I. M.: Aquaplanet simulations of tropical cyclones, Curr. Clim. Change Rep., 5, 185–195, https://doi.org/10.1007/s40641-019-00133-y, 2019. a
Morrison, H. and Milbrandt, J. A.: Parameterization of cloud microphysics based on the prediction of bulk ice particle properties. Part I: Scheme description and idealized tests, J. Atmos. Sci., 72, 287–311, 2015. a
Morrison, H., Curry, J., and Khvorostyanov, V.: A new double-moment microphysics parameterization for application in cloud and climate models. Part I: Description, J. Atmos. Sci., 62, 1665–1677, 2005. a
Muller, C. J. and Bony, S.: What favors convective aggregation and why?, Geophys. Res. Lett., 42, 5626–5643, https://doi.org/10.1002/2015GL064260, 2015. a
Muller, C. J. and Held, I. M.: Detailed investigation of the self-aggregation of convection in cloud resolving simulations, J. Atmos. Sci., 69, 2551–2565, https://doi.org/10.1175/JAS-D-11-0257.1, 2012. a, b
Müller, S. K. M. and Hohenegger, C.: Self-aggregation of convection in spatially-varying sea surface temperatures, J. Adv. Model. Earth Sy., 12, e2019MS001698, https://doi.org/10.1029/2019MS001698, 2020. a, b, c, d, e
Neale, R., Chen, C.-C., Gettelman, A., Lauritzen, P., Park, S., WIlliamson, D., Conley, A., Garcia, R., Kinnison, D., Lamarque, J.-F., Marsh, D., Mills, M., Smith, A., Tilmes, S., Vitt, F., Morrison, H., Cameron-Smith, P., Collins, W., Iacono, M., Easter, R., Ghan, S., Liu, X., Rasch, P., and Taylor, M.: Description of the NCAR Community Atmosphere Model (CAM 5.0), NCAR Tech. Note TN-486, http://www.cesm.ucar.edu/models/cesm1.0/cam/docs/description/cam5_desc.pdf (last access: 1 December 2023), 2012. a
Nishizawa, S., Yashiro, H., Sato, Y., Miyamoto, Y., and Tomita, H.: Influence of grid aspect ratio on planetary boundary layer turbulence in large-eddy simulations, Geosci. Model Dev., 8, 3393–3419, https://doi.org/10.5194/gmd-8-3393-2015, 2015. a
Nuijens, L. and Emanuel, K.: Congestus modes in circulating equilibria of the tropical atmosphere in a two-column model, Q. J. Roy. Meteor. Soc., 144, 2676–2692, https://doi.org/10.1002/qj.3385, 2018. a
Pope, K. N., Holloway, C. E., Jones, T. R., and Stein, T. H. M.: Radiation, Clouds, and Self-Aggregation in RCEMIP Simulations, J. Adv. Model. Earth Sy., 15, e2022MS003317, https://doi.org/10.1029/2022MS003317, 2023. a
Popke, D., Stevens, B., and Voigt, A.: Climate and climate change in a radiative-convective equilibrium version of ECHAM6, J. Adv. Model. Earth Sy., 5, 1–14, https://doi.org/10.1029/2012MS000191, 2013. a
Randall, D. A., DeMott, C., Stan, C., Khairoutdinov, M., Benedict, J., McCrary, R., Thayer-Calder, K., and Branson, M.: Simulations of the Tropical General Circulation with a Multiscale Global Model, Meteorological Monographs, 56, 15.1–15.15, https://doi.org/10.1175/amsmonographs-d-15-0016.1, 2016. a
Raymond, D. J.: Convective processes and tropical atmospheric circulations, Q. J. Roy. Meteor. Soc., 120, 1431–1455, https://doi.org/10.1002/qj.49712052002, 1994. a
Rayner, N. A., Parker, D. E., Horton, E. B., Folland, C. K., Alexander, L. V., Powell, D. P., Kent, E. C., and Kaplan, A.: Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century, J. Geophys. Res., 108, 4407, https://doi.org/10.1029/2002JD002670, 2023. a
Reed, K. A., Silvers, L. G., Wing, A. A., Hu, I.-K., and Medeiros, B.: Using Radiative Convective Equilibrium to Explore Clouds and Climate in the Community Atmosphere Model, J. Adv. Model. Earth Sys., 13, e2021MS002539, https://doi.org/10.1029/2021MS002539, 2021. a, b, c
Roehrig, R., Beau, I., Saint-Martin, D., Alias, A., Decharme, B., Guérémy, J.-F., Voldoire, A., Abdel-Lathif, A. Y., Bazile, E., Belamari, S., Blein, S., Bouniol, D., Bouteloup, Y., Cattiaux, J., Chauvin, F., Chevallier, M., Colin, J., Douville, H., Marquet, P., Michou, M., Nabat, P., Oudar, T., Peyrillé, P., Piriou, J.-M., Salas y Mélia, D., Séférian, R., and Sénési, S.: The CNRM global atmosphere model ARPEGE-Climate 6.3: description and evaluation, J. Adv. Model. Earth Sy., 12, e2020MS002075, https://doi.org/10.1029/2020MS002075, 2020. a
Romps, D. M.: The dry-entropy budget of a moist atmosphere, J. Atmos. Sci., 65, 3779–3799, 2008. a
Sato, Y., Nishizawa, S., Yashiro, H., Miyamoto, Y., Kajikawa, Y., and Tomita, H.: Impacts of cloud microphysics on trade wind cumulus: which cloud microphysics processes contribute to the diversity in a large eddy simulation?, Prog. Earth Planet. Sci., 2, 23, https://doi.org/10.1186/s40645-015-0053-6, 2015. a
Satoh, M., Tomita, H., Yashiro, H., Miura, H., Kodama, C., Seiki, T., Noda, A. T., Yamada, Y., Goto, D., Sawada, M., Miyoshi, T., Niwa, Y., Hara, M., Ohno, T., Iga, S.-I., Inoue, T. A. T., and Kubokawa, H.: The Non-hydrostatic ICosahedral Atmospheric Model: Description and development, Prog. Earth Planet. Sci., 1, 18, https://doi.org/10.1186/s40645-014-0018-1, 2014. a
Satoh, M., Arakami, K., and Sawada, M.: Structure of tropical convective systems in aqua-planet experiments: Radiative-convective equilibrium versus the Earth-like experiments, SOLA, 12, 220–224, 2016. a
Silvers, L., Reed, K., and Wing, A.: The Response of the Large-Scale Tropical Circulation to Warming, J. Adv. Model. Earth Sy., 15, e2021MS002966, https://doi.org/10.1029/2021MS002966, 2023. a, b, c, d
Silvers, L. G. and Robinson, T.: Clouds and radiation in a Mock-Walker circulation, J. Adv. Model. Earth Sy., 13, e2020MS002196, https://doi.org/10.1029/2020MS002196, 2021. a, b, c, d, e, f, g
Silvers, L. G., Stansfield, A. M., and Reed, K. A.: The Impact of Rotation on Tropical Climate, the Hydrologic Cycle, and Climate Sensitivty, Geophys. Res. Lett., 51, e2023GL105850, https://doi.org/10.1029/2023GL105850, 2024. a, b
Sobel, A. H. and Bretherton, C. S.: Modeling Tropical Precipitation in a Single Column, J. Climate, 13, 4378–4392, https://doi.org/10.1175/1520-0442(2000)013<4378:MTPIAS>2.0.CO;2, 2000. a
Sobel, A. H., Wing, A. A., Camargo, S. J., Patricola, C. M., Vecchi, G. A., Lee, C.-Y., and Tippett, M. K.: Tropical cyclone frequency, Earth's Future, 9, e2021EF002275, https://doi.org/10.1029/2021EF002275, 2021. a
Sokol, A. B. and Hartmann, D. L.: Congestus Mode Invigoration by Convective Aggregation in Simulations of Radiative-Convective Equilibrium, J. Adv. Model. Earth Sy., 14, e2022MS003045, https://doi.org/10.1029/2022MS003045, 2022. a, b
Stauffer, C. and Wing, A.: Properties, changes, and controls of deep-convecting clouds in Radiative-Convective Equilibrium, J. Adv. Model. Earth Sy., 14, e2021MS002917, https://doi.org/10.1029/2021MS002917, 2022. a, b, c, d
Stauffer, C. L.: Cloud Feedbacks and Convective Self-Aggregation in the Radiative-Convective Equilibrium Model Intercomparison Project, Ph. D. thesis, Florida State University, Tallahassee, FL, https://fsu-flvc.primo.exlibrisgroup.com/permalink/01FALSC_FSU/1pc67ru/cdi_proquest_journals_2868559250 (last access: 1 December 2023), 219 pp., 2023. a, b, c
Stauffer, C. L. and Wing, A. A.: Explicitly Resolved Cloud Feedbacks in the Radiative-Convective Equilibrium Model Intercomparison Project, J. Adv. Model. Earth Sy., 15, e2023MS003738, https://doi.org/10.1029/2023MS003738, 2023. a
Stratton, R., Senior, C., Vosper, S., Folwell, S., Boutle, I., Earnshaw, P., Kendon, E., Lock, A., Malcolm, A., Manners, J., Morcrette, C., Short, C., Stirling, A. J., Taylor, C. M., Tucker, S., Webster, S., and Wilkinson, J. M.: A Pan-African convection-permitting regional climate simulation with the Met Office Unified Model: CP4-Africa, J. Climate, 31, 3485–3508, 2018. a
Tatebe, H., Ogura, T., Nitta, T., Komuro, Y., Ogochi, K., Takemura, T., Sudo, K., Sekiguchi, M., Abe, M., Saito, F., Chikira, M., Watanabe, S., Mori, M., Hirota, N., Kawatani, Y., Mochizuki, T., Yoshimura, K., Takata, K., O'ishi, R., Yamazaki, D., Suzuki, T., Kurogi, M., Kataoka, T., Watanabe, M., and Kimoto, M.: Description and basic evaluation of simulated mean state, internal variability, and climate sensitivity in MIROC6, Geosci. Model Dev., 12, 2727–2765, https://doi.org/10.5194/gmd-12-2727-2019, 2019. a
Tompkins, A. M.: On the relationship between tropical convection and sea surface temperature, J. Climate, 14, 633–637, https://doi.org/10.1175/1520-0442(2001)014<0633:OTRBTC>2.0.CO;2, 2001. a, b
Tompkins, A. M. and Craig, G. C.: Radiative-convective equilibrium in a three-dimensional cloud-ensemble model, Q. J. Roy. Meteor. Soc., 124, 2073–2097, 1998. a
Voldoire, A., Saint-Martin, D., Sénési., S., Decharme, B., Alias, A., Chevallier, M., Colin, J., Guérémy, J.-F., Michou, M., Moine, M.-P., Nabat, P., Roehrig, R., y Mélia, D. S., Séférian, R., Valcke, S., Beau, I., Belamari, S., Berthet, S., Cassou, C., Cattiaux, J., Deshayes, J., Douville, H., Ethé, C., Franchistéguy, L., Geoffroy, O., Lévy, C., Madec, G., Meurdesoif, Y., Msadek, R., Ribes, A., Sanchez-Gomez, E., Terray, L., and Waldman, R.: Evaluation of CMIP6 DECK expeirments with CNRM-CM6.1, J. Adv. Model. Earth Sy., 11, 2177–2213, https://doi.org/10.1029/2019MS001683, 2019. a
Walters, D., Baran, A. J., Boutle, I., Brooks, M., Earnshaw, P., Edwards, J., Furtado, K., Hill, P., Lock, A., Manners, J., Morcrette, C., Mulcahy, J., Sanchez, C., Smith, C., Stratton, R., Tennant, W., Tomassini, L., Van Weverberg, K., Vosper, S., Willett, M., Browse, J., Bushell, A., Carslaw, K., Dalvi, M., Essery, R., Gedney, N., Hardiman, S., Johnson, B., Johnson, C., Jones, A., Jones, C., Mann, G., Milton, S., Rumbold, H., Sellar, A., Ujiie, M., Whitall, M., Williams, K., and Zerroukat, M.: The Met Office Unified Model Global Atmosphere 7.0/7.1 and JULES Global Land 7.0 configurations, Geosci. Model Dev., 12, 1909–1963, https://doi.org/10.5194/gmd-12-1909-2019, 2019. a
Webb, M., Senior, C., Bony, S., and Morcrette, J. J.: Combining ERBE and ISCCP data to assess clouds in the Hadley Centre, ECWMF and LMD atmospheric climate models, Clim. Dynam., 17, 905–922, 2001. a
Wing, A. A.: allison-wing/RCEMIP-II: v1.2 (v1.2), Zenodo [code], https://doi.org/10.5281/zenodo.11043720, 2024a. a
Wing, A. A.: SAM6.11.2-AWing, Zenodo [code], https://doi.org/10.5281/zenodo.12626604, 2024b. a
Wing, A. A. and Silvers, L.: Code and data for RCEMIP-II: Mock-Walker Simulations as Phase II of the Radiative-Convective Equilibrium Model Intercomparison Project (v1.0), Zenodo [data set], https://doi.org/10.5281/zenodo.10137266, 2023. a
Wing, A. A. and Singh, M. S.: Control of Stability and Relative Humidity in the Radiative-Convective Equilibrium Model Intercomparison Project, J. Adv. Model. Earth Sy., 16, e2023MS003914, https://doi.org/10.1029/2023MS003914, 2024. a, b
Wing, A. A., Emanuel, K., Holloway, C. E., and Muller, C. J.: Convective self-aggregation in numerical simulations: A review, Surv. Geophys., 38, 1173–1197, https://doi.org/10.1007/s10712-017-9408-4, 2017. a
Wing, A. A., Reed, K. A., Satoh, M., Stevens, B., Bony, S., and Ohno, T.: Radiative–convective equilibrium model intercomparison project, Geosci. Model Dev., 11, 793–813, https://doi.org/10.5194/gmd-11-793-2018, 2018. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p
Wing, A. A., Stauffer, C., Becker, T., Reed, K., Ahn, M.-S., Arnold, N., Bony, S., Branson, M., Bryan, G., Chaboureau, J.-P., de Roode, S., Gayatri, K., Hohenegger, C., Hu, I.-K., Jansson, F., Jones, T., Khairoutdionv, M., Kim, D., Martin, Z., Matsugishi, S., Medeiros, B., Miura, H., Moon, Y., Möller, S., Ohno, T., Popp, M., Prabhakaran, T., Randall, D., Rios-Berrios, R., Rochetin, N., Roehrig, R., Romps Jr., D., J. R., Satoh, M., Silvers, L., Singh, M., Stevens, B., Tomassini, L., van Heerwaarden, C., Wang, S., and Zhao, M.: Clouds and convective self-aggregation in a multi-model ensemble of radiative-convective equilibrium simulations, J. Adv. Model. Earth Sy., 12, e2020MS002138, https://doi.org/10.1029/2020MS002138, 2020a. a, b, c, d, e, f, g
Wing, A. A., Stauffer, C., Becker, T., Reed, K., Ahn, M.-S., Arnold, N., Bony, S., Branson, M., Bryan, G., Chaboureau, J.-P., de Roode, S., Gayatri, K., Hohenegger, C., Hu, I.-K., Jansson, F., Jones, T., Khairoutdionv, M., Kim, D., Martin, Z., Matsugishi, S., Medeiros, B., Miura, H., Moon, Y., Möller, S., Ohno, T., Popp, M., Prabhakaran, T., Randall, D., Rios-Berrios, R., Rochetin, N., Roehrig, R., Romps Jr., D., J. R., Satoh, M., Silvers, L., Singh, M., Stevens, B., Tomassini, L., van Heerwaarden, C., Wang, S., and Zhao, M.: Radiative-Convective Equilibrium Model Intercomparison Project (RCEMIP) Simulation Dataset, WDC Climate [data set], http://hdl.handle.net/21.14101/d4beee8e-6996-453e-bbd1-ff53b6874c0e (last access: 15 August 2024), 2020b. a
Wing, A. A., Silvers, L., and Reed, K.: CESM2p1p3rcemipII, Zenodo [code], https://doi.org/10.5281/zenodo.12668076, 2024. a
Wofsy, J. and Kuang, Z.: Cloud-resolving model simulations and a simple model of an idealized Walker cell, J. Climate, 25, 8090–8107, https://doi.org/10.1175/JCLI-D-11-00692.1, 2012. a, b
Wu, C.-M., Lin, H.-C., Cheng, F.-Y., and Chien, M.-H.: Implementation of the Land Surface Processes into a Vector Vorticity Equation Model (VVM) to Study its Impact on Afternoon Thunderstorms over Complex Topography in Taiwan, Asia-Pacific Journal of Atmospheric Sciences, 55, 701–717, https://doi.org/10.1007/s13143-019-00116-x, 2019. a
Yano, J.-I., Moncrieff, M. W., and Grabowski, W. W.: Walker-type mean circulations and convectively coupled tropical waves as an interacting System, J. Atmos. Sci., 59, 1566–1577, 2002. a
Zelinka, M. D., Klein, S. A., and Hartmann, D. L.: Computing and partitioning cloud feedbacks using cloud property histograms. Part I: Cloud radiative kernels, J. Climate, 25, 3715–3735, https://doi.org/10.1175/JCLI-D-11-00248.1, 2012a. a
Zelinka, M. D., Klein, S. A., and Hartmann, D. L.: Computing and partitioning cloud feedbacks using cloud property hitsograms. Part II: Attribution to changes in cloud amount, altitude, and optical depth, J. Climate, 25, 3736–3754, https://doi.org/10.1175/JCLI-D-11-00249.1, 2012b. a
Zelinka, M. D., Klein, S. A., Taylor, K. E., Andrews, T., Webb, M. J., Gregory, J. M., and Forster, P. M.: Contributions of different cloud types to feedbacks and rapid adjustments in CMIP5, J. Climate, 26, 5007–5027, https://doi.org/10.1175/JCLI-D-12-00555.1, 2013. a
Zhao, M., Golaz, J.-C., Held, I. M., Guo, H., Balaji, V., Benson, R., Chen, J.-H., Chen, X., Donner, L. J., Dunne, J. P., Dunne, K., Durachta, J., Fan, S.-M., Freidenreich, S. M., Garner, S. T., Ginoux, P., Harris, L. M., Horowitz, L. W., Krasting, J. P., Langenhorst, A. R., Liang, Z., Lin, P., Lin, S.-J., Malyshev, S. L., Mason, E., Milly, P. C. D., Ming, Y., Naik, V., Paulot, F., Paynter, D., Philipps, P., Radhakrishnan, A., Ramaswamy, V., Robinson, T., Schwarzkopf, D., Seman, C. J., Shevliakova, E., Shen, Z., Shin, H., Silvers, L. G., Wilson, J. R., Winton, M., Wittenberg, A. T., Wyman, B., and Xiang, B.: The GFDL global atmosphere and land model AM4.0/LM4.0: 1. Simulation characteristics with prescribed SSTs, J. Adv. Model. Earth Sy., 10, 691–734, https://doi.org/10.1002/2017MS001208, 2018a. a
Zhao, M., Golaz, J.-C., Held, I. M., Guo, H., Balaji, V., Benson, R., Chen, J.-H., Chen, X., Donner, L. J., Dunne, J. P., Dunne, K., Durachta, J., Fan, S.-M., Freidenreich, S. M., Garner, S. T., Ginoux, P., Harris, L. M., Horowitz, L. W., Krasting, J. P., Langenhorst, A. R., Liang, Z., Lin, P., Lin, S.-J., Malyshev, S. L., Mason, E., Milly, P. C. D., Ming, Y., Naik, V., Paulot, F., Paynter, D., Philipps, P., Radhakrishnan, A., Ramaswamy, V., Robinson, T., Schwarzkopf, D., Seman, C. J., Shevliakova, E., Shen, Z., Shin, H., Silvers, L. G., Wilson, J. R., Winton, M., Wittenberg, A. T., Wyman, B., and Xiang, B.: The GFDL global atmosphere and land model AM4.0/LM4.0: 2. Model description, sensitivity studies, and tuning strategies, J. Adv. Model. Earth Sy., 10, 735–769, https://doi.org/10.1002/2017MS001209, 2018b. a
Zhou, C., Zelinka, M. D., and Klein, S. A.: Impact of decadal cloud variations on the Earth's energy budget, Nat. Geosci, 9, 871–874, https://doi.org/10.1038/ngeo2828, 2016. a
Zhou, L., Lin, S.-J., Chen, J.-H., Harris, L. M., Chen, X., and Rees, S. L.: Toward Convective-Scale Prediction within the Next Generation Global Prediction System, B. Am. Meteorol. Soc., 100, 1225–1243, https://doi.org/10.1175/bams-d-17-0246.1, 2019. a
Zhu, J., Otto-Bliesner, B. L., Brady, E. C., Gettelman, A., Bacmeister, J. T., Neale, R. B., Poulsen, C. J., Shaw, J. K., McGraw, Z. S., and Kay, J. E.: LGM Paleoclimate Constraints Inform Cloud Parameterizations and Equilibrium Climate Sensitivity in CESM2, J. Adv. Model. Earth Sy., 14, e2021MS002776, https://doi.org/10.1029/2021MS002776, 2022. a, b
The requested paper has a corresponding corrigendum published. Please read the corrigendum first before downloading the article.
- Article
(15006 KB) - Full-text XML