The Radiative‐Convective Equilibrium Model Intercomparison Project (RCEMIP) is an intercomparison of multiple types of numerical models configured in radiative‐convective equilibrium (RCE). RCE is an idealization of the tropical atmosphere that has long been used to study basic questions in climate science. Here, we employ RCE to investigate the role that clouds and convective activity play in determining cloud feedbacks, climate sensitivity, the state of convective aggregation, and the equilibrium climate. RCEMIP is unique among intercomparisons in its inclusion of a wide range of model types, including atmospheric general circulation models (GCMs), single column models (SCMs), cloud‐resolving models (CRMs), large eddy simulations (LES), and global cloud‐resolving models (GCRMs). The first results are presented from the RCEMIP ensemble of more than 30 models. While there are large differences across the RCEMIP ensemble in the representation of mean profiles of temperature, humidity, and cloudiness, in a majority of models anvil clouds rise, warm, and decrease in area coverage in response to an increase in sea surface temperature (SST). Nearly all models exhibit self‐aggregation in large domains and agree that self‐aggregation acts to dry and warm the troposphere, reduce high cloudiness, and increase cooling to space. The degree of self‐aggregation exhibits no clear tendency with warming. There is a wide range of climate sensitivities, but models with parameterized convection tend to have lower climate sensitivities than models with explicit convection. In models with parameterized convection, aggregated simulations have lower climate sensitivities than unaggregated simulations.

[1] Subtropical marine low cloud sensitivity to an idealized climate change is compared in six large-eddy simulation (LES) models as part of CGILS. July cloud cover is simulated at three locations over the subtropical northeast Pacific Ocean, which are typified by cold sea surface temperatures (SSTs) under well-mixed stratocumulus, cool SSTs under decoupled stratocumulus, and shallow cumulus clouds overlying warmer SSTs. The idealized climate change includes a uniform 2 K SST increase with corresponding moist-adiabatic warming aloft and subsidence changes, but no change in free-tropospheric relative humidity, surface wind speed, or CO 2 . For each case, realistic advective forcings and boundary conditions are generated for the control and perturbed states which each LES runs for 10 days into a quasi-steady state. For the control climate, the LESs correctly produce the expected cloud type at all three locations. With the perturbed forcings, all models simulate boundary-layer deepening due to reduced subsidence in the warmer climate, with less deepening at the warm-SST location due to regulation by precipitation. The models do not show a consistent response of liquid water path and albedo in the perturbed climate, though the majority predict cloud thickening (negative cloud feedback) at the cold-SST location and slight cloud thinning (positive cloud feedback) at the cool-SST and warm-SST locations. In perturbed climate simulations at the cold-SST location without the subsidence decrease, cloud albedo consistently decreases across the models. Thus, boundary-layer cloud feedback on climate change involves compensating thermodynamic and dynamic effects of warming and may interact with patterns of subsidence change.

Abstract. The current version of the Dutch AtmosphericLarge-Eddy Simulation (DALES) is presented. DALES is a large-eddy simulation code designed for studies of the physics of the atmospheric boundary layer, including convective and stable boundary layers as well as cloudy boundary layers. In addition, DALES can be used for studies of more specific cases, such as flow over sloping or heterogeneous terrain, and dispersion of inert and chemically active species. This paper contains an extensive description of the physical and numerical formulation of the code, and gives an overview of its applications and accomplishments in recent years.

[1] CGILS-the CFMIP-GASS Intercomparison of Large Eddy Models (LESs) and single column models (SCMs)-investigates the mechanisms of cloud feedback in SCMs and LESs under idealized climate change perturbation. This paper describes the CGILS results from 15 SCMs and 8 LES models. Three cloud regimes over the subtropical oceans are studied: shallow cumulus, cumulus under stratocumulus, and wellmixed coastal stratus/stratocumulus. In the stratocumulus and coastal stratus regimes, SCMs without activated shallow convection generally simulated negative cloud feedbacks, while models with active shallow convection generally simulated positive cloud feedbacks. In the shallow cumulus alone regime, this relationship is less clear, likely due to the changes in cloud depth, lateral mixing, and precipitation or a combination of them. The majority of LES models simulated negative cloud feedback in the wellmixed coastal stratus/stratocumulus regime, and positive feedback in the shallow cumulus and stratocumulus regime. A general framework is provided to interpret SCM results: in a warmer climate, the moistening rate of the cloudy layer associated with the surface-based turbulence parameterization is enhanced; together with weaker VOL. 5, 826-842, doi:10.1002/2013MS000246, 2013 large-scale subsidence, it causes negative cloud feedback. In contrast, in the warmer climate, the drying rate associated with the shallow convection scheme is enhanced. This causes positive cloud feedback. These mechanisms are summarized as the ''NESTS'' negative cloud feedback and the ''SCOPE'' positive cloud feedback (Negative feedback from Surface Turbulence under weaker Subsidence-Shallow Convection PositivE feedback) with the net cloud feedback depending on how the two opposing effects counteract each other. The LES results are consistent with these interpretations.Citation: Zhang, M., et al. (2013), CGILS: Results from the first phase of an international project to understand the physical mechanisms of low cloud feedbacks in single column models, J. Adv. Model. Earth Syst., 5, 826-842,

The length scale evolution of various quantities in a clear convective boundary layer (CBL), a stratocumulustopped boundary layer, and three radiatively cooled (''smoke cloud'') convective boundary layers are studied by means of large-eddy simulations on a large horizontal domain (25.6 ϫ 25.6 km 2 ). In the CBL the virtual potential temperature and the vertical velocity fields are dominated by horizontal scales on the order of the boundary layer depth. In contrast, the potential temperature and the specific humidity fields become gradually dominated by mesoscale fluctuations. However, at the mesoscales their effects on the virtual potential temperature fluctuations nearly compensate. It is found that mesoscale fluctuations are negligibly small only for conserved variables that have an entrainment to surface flux ratio close to Ϫ0.25, which is about the flux ratio for the buoyancy. In the CBL the moisture and potential temperature flux ratios can have values that significantly deviate from this number.The geometry of the buoyancy flux was manipulated by cooling the clear convective boundary layer from the top, in addition to a positive buoyancy flux at the surface. For these radiatively cooled cases it is found that both the vertical velocity as well as the virtual potential temperature spectra tend to broaden. The role of the buoyancy flux in their respective prognostic variance equations is discussed. It is argued that in the upper part of the clear CBL, where the mean vertical stratification is stable, vertical velocity variance and virtual potential temperature variance cannot be produced simultaneously. For the stratocumulus case, in which latent heat release effects in the cloud layer play an important role in its dynamics, the field of any quantity, except for the vertical velocity, becomes dominated by mesoscale fluctuations.In general, the location of the spectral peak of any quantity becoming constrained by the domain size should be avoided. The answer to the question of how large the LES horizontal domain size should be in order to include mesoscale fluctuations will, on the one hand, depend on the type of convection to be simulated and the kind of physical question one aims to address, and, on the other hand, the time duration of the simulation. Only if one aims to study the dynamics of a dry CBL that excludes moisture, a rather small domain size suffices. In case one aims to examine either the spatial evolution of the fields of any arbitrary conserved scalar in the CBL, or any quantity in stratocumulus clouds except for the vertical velocity, a larger domain size that allows the development of mesoscale fluctuations will be necessary.

Previously observed twice-Clausius-Clapeyron (2CC) scaling for extreme precipitation at hourly time scales has led to discussions about its origin. The robustness of this scaling is assessed by analyzing a subhourly dataset of 10-min resolution over the Netherlands. The results confirm the validity of the previously found 2CC scaling for extreme convective precipitation.Using a simple entraining plume model, an idealized deep convective environmental temperature profile is perturbed to analyze extreme precipitation scaling from a frequently used relation based on the column condensation rate. The plume model simulates a steady precipitation increase that is greater than ClausiusClapeyron scaling (super-CC scaling). Precipitation intensity increase is shown to be controlled by a flux of moisture through the cloud base and in-cloud lateral moisture convergence. Decomposition of this scaling relation into a dominant thermodynamic and additional dynamic component allows for better understanding of the scaling and demonstrates the importance of vertical velocity in both dynamic and thermodynamic scaling. Furthermore, systematically increasing the environmental stability by adjusting the temperature perturbations from constant to moist adiabatic increase reveals a dependence of the scaling on the change in environmental stability. As the perturbations become increasingly close to moist adiabatic, the scaling found by the entraining plume model decreases to CC scaling. Thus, atmospheric stability changes, which are expected to be dependent on the latitude, may well play a key role in the behavior of precipitation extremes in the future climate.

As part of the European Project on Cloud Systems in Climate Models, the diurnal cycle of stratocumulus has been simulated with Large-Eddy Simulation (LES) models and Single Column Models (SCMs). The models were initialized and compared with observations collected in marine stratocumulus in July 1987 during the First International Satellite Cloud Climatology Project Regional Experiment. The results of the six LES models are found to be in a fair agreement with the observations. They all capture the distinct diurnal variation in the cloud liquid-water path, the turbulence profiles and clearly show a decoupled boundary layer during daytime and a vertically well-mixed boundary layer during the night. Entrainment of relatively dry and warm air from just above the inversion into the boundary layer is the major process modifying the thermodynamic structure of the boundary layer during the night. The differences that arise in the liquid-water path evolution can therefore be attributed mainly to differences in the entrainment rate. The mean entrainment rates computed from the LES model results are 0.58 +/- 0.08 cm s(-1) and 0.36 +/- 0.03 cm s(-1) for the night-time and daytime periods, respectively. If the horizontal domain size in a LES model is enlarged, mesoscale fluctuations develop. This leads to a broader liquid-water path distribution and a reduction of the cloud albedo. To assess the quality of the representation of stratocumulus in general-circulation models, results from ten SCMs are compared with observations and LES results. The SCM latent and sensible heat fluxes at the surface agree fairly well with the LES results. Many of the SCMs predict a liquid-water path which is much too low, a cloud cover smaller than unity, and cloud tops that are lower than the observations and the LES results. This results in a much larger amount of downwelling short-wave radiation absorbed at the sea surface. Improvement of entrainment parametrizations is needed for a better representation of stratocumulus in SCMs. Observations and LES results of entrainment rates for different stratocumulus cases are compared. The observed entrainment rates in Atlantic stratocumulus clouds during the Atlantic Stratocumulus Transition Experiment (ASTEX) are larger than for the ones over the Pacific Ocean off the coast of California. Results from LES models corroborate these findings. The differences in the entrainment rate can likely be attributed to the smaller inversion jumps of the liquid-water potential temperature for the ASTEX stratocumulus case

The application of a steady-state vertical velocity equation for parameterized moist convective updrafts in climate and weather prediction models is currently common practice. This equation usually contains an advection, a buoyancy, and a lateral entrainment term, whereas the effects of pressure gradient and subplume contributions are typically incorporated as proportionality constants a and b for the buoyancy and the entrainment terms, respectively. A summary of proposed values of these proportionality constants a and b in the literature demonstrates that there is a large uncertainty in their most appropriate values. To shed new light on this situation an analysis is presented of the full vertical budget equation for shallow cumulus clouds obtained from large eddy simulations of three different Global Energy and Water Cycle Experiment (GEWEX) Cloud System Study (GCSS) intercomparison cases. It is found that the pressure gradient term is the dominant sink term in the vertical velocity budget, whereas the entrainment term only gives a small contribution. This result is at odds with the parameterized vertical velocity equation in the literature as it employs the entrainment term as the major sink term. As a practical solution the damping effect of the pressure term may be parameterized in terms of the lateral entrainment rates as used for thermodynamic quantities like the total specific humidity. By using a least squares method, case-dependent optimal values are obtained for the proportionality constants a and b, which are linearly related with each other. This relation can be explained from a linear relationship between the lateral entrainment rate and the buoyancy.