This paper describes the main characteristics of CNRM‐CM6‐1, the fully coupled atmosphere‐ocean general circulation model of sixth generation jointly developed by Centre National de Recherches Météorologiques (CNRM) and Cerfacs for the sixth phase of the Coupled Model Intercomparison Project 6 (CMIP6). The paper provides a description of each component of CNRM‐CM6‐1, including the coupling method and the new online output software. We emphasize where model's components have been updated with respect to the former model version, CNRM‐CM5.1. In particular, we highlight major improvements in the representation of atmospheric and land processes. A particular attention has also been devoted to mass and energy conservation in the simulated climate system to limit long‐term drifts. The climate simulated by CNRM‐CM6‐1 is then evaluated using CMIP6 historical and Diagnostic, Evaluation and Characterization of Klima (DECK) experiments in comparison with CMIP5 CNRM‐CM5.1 equivalent experiments. Overall, the mean surface biases are of similar magnitude but with different spatial patterns. Deep ocean biases are generally reduced, whereas sea ice is too thin in the Arctic. Although the simulated climate variability remains roughly consistent with CNRM‐CM5.1, its sensitivity to rising CO₂ has increased: the equilibrium climate sensitivity is 4.9 K, which is now close to the upper bound of the range estimated from CMIP5 models.

This study introduces CNRM‐ESM2‐1, the Earth system (ES) model of second generation developed by CNRM‐CERFACS for the sixth phase of the Coupled Model Intercomparison Project (CMIP6). CNRM‐ESM2‐1 offers a higher model complexity than the Atmosphere‐Ocean General Circulation Model CNRM‐CM6‐1 by adding interactive ES components such as carbon cycle, aerosols, and atmospheric chemistry. As both models share the same code, physical parameterizations, and grid resolution, they offer a fully traceable framework to investigate how far the represented ES processes impact the model performance over present‐day, response to external forcing and future climate projections. Using a large variety of CMIP6 experiments, we show that represented ES processes impact more prominently the model response to external forcing than the model performance over present‐day. Both models display comparable performance at replicating modern observations although the mean climate of CNRM‐ESM2‐1 is slightly warmer than that of CNRM‐CM6‐1. This difference arises from land cover‐aerosol interactions where the use of different soil vegetation distributions between both models impacts the rate of dust emissions. This interaction results in a smaller aerosol burden in CNRM‐ESM2‐1 than in CNRM‐CM6‐1, leading to a different surface radiative budget and climate. Greater differences are found when comparing the model response to external forcing and future climate projections. Represented ES processes damp future warming by up to 10% in CNRM‐ESM2‐1 with respect to CNRM‐CM6‐1. The representation of land vegetation and the CO₂‐water‐stomatal feedback between both models explain about 60% of this difference. The remainder is driven by other ES feedbacks such as the natural aerosol feedback.

This paper presents the Meso-NH model version 5.4. Meso-NH is an atmospheric non hydrostatic research model that is applied to a broad range of resolutions, from synoptic to turbulent scales, and is designed for studies of physics and chemistry. It is a limited-area model employing advanced numerical techniques, including monotonic advection schemes for scalar transport and fourth-order centered or odd-order WENO advection schemes for momentum. The model includes state-of-the-art physics parameterization schemes that are important to represent convective-scale phenomena and turbulent eddies, as well as flows at larger scales. In addition, Meso-NH has been expanded to provide capabilities for a range of Earth system prediction applications such as chemistry and aerosols, electricity and lightning, hydrology, wildland fires, volcanic eruptions, and cyclones with ocean coupling. Here, we present the main innovations to the dynamics and physics of the code since the pioneer paper of Lafore et al. (1998) and provide an overview of recent applications and couplings.

This is the first part of a series of two articles analyzing the global thermal properties of atmosphere–ocean coupled general circulation models (AOGCMs) within the framework of a two-layer energy-balance model (EBM). In this part, the general analytical solution of the system is given and two idealized climate change scenarios, one with a step forcing and one with a linear forcing, are discussed. These solutions give a didactic description of the contributions from the equilibrium response and of the fast and slow transient responses during a climate transition. Based on these analytical solutions, a simple and physically based procedure to calibrate the two-layer model parameters using an AOGCM step-forcing experiment is introduced. Using this procedure, the global thermal properties of 16 AOGCMs participating in phase 5 of the Coupled Model Intercomparison Project (CMIP5) are determined. It is shown that, for a given AOGCM, the EBM tuned with only the abrupt 4×CO₂ experiment is able to reproduce with a very good accuracy the temperature evolution in both a step-forcing and a linear-forcing experiment. The role of the upper-ocean and deep-ocean heat uptakes in the fast and slow responses is also discussed. One of the main weaknesses of the simple EBM discussed in this part is its ability to represent the evolution of the top-of-the-atmosphere radiative imbalance in the transient regime. This issue is addressed in Part II by taking into account the efficacy factor of deep-ocean heat uptake.

The current version of the Dutch Atmospheric Large-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.

In this second part of a series of two articles analyzing the global thermal properties of atmosphere–ocean coupled general circulation models (AOGCMs) within the framework of a two-layer energy-balance model (EBM), the role of the efficacy of deep-ocean heat uptake is investigated. Taking into account such an efficacy factor is shown to amount to representing the effect of deep-ocean heat uptake on the local strength of the radiative feedback in the transient regime. It involves an additional term in the formulation of the radiative imbalance at the top of the atmosphere (TOA), which explains the nonlinearity between radiative imbalance and the mean surface temperature observed in some AOGCMs. An analytical solution of this system is given and this simple linear EBM is calibrated for the set of 16 AOGCMs of phase 5 of the Coupled Model Intercomparison Project (CMIP5) studied in Part I. It is shown that both the net radiative fluxes at TOA and the global surface temperature transient response are well represented by the simple EBM over the available period of simulations. Differences between this two-layer EBM and the previous version without an efficacy factor are analyzed and relationships between parameters are discussed. The simple model calibration applied to AOGCMs constitutes a new method for estimating their respective equilibrium climate sensitivity and adjusted radiative forcing amplitude from short-term step-forcing simulations and more generally a method to compute their global thermal properties.