This study presents the version of the LMDZ global atmospheric model used as the atmospheric component of the Institut Pierre Simon Laplace coupled model (IPSL-CM6A-LR) to contribute to the 6th phase of the international Coupled Model Intercomparison Project (CMIP6). This LMDZ6A version includes original convective parameterizations that define the LMDZ "New Physics": a mass flux parameterization of the organized structures of the convective boundary layer, the "thermal plume model," and a parameterization of the cold pools created by reevaporation of convective rainfall. The vertical velocity associated with thermal plumes and gust fronts of cold pools are used to control the triggering and intensity of deep convection. Because of several shortcomings, the early version 5B of this New Physics was worse than the previous "Standard Physics" version 5A regarding several classical climate metrics. To overcome these deficiencies, version 6A includes new developments: a stochastic triggering of deep convection, a modification of the thermal plume model that allows the representation of stratocumulus and cumulus clouds in a unified framework, an improved parameterization of very stable boundary layers, and the modification of the gravity waves scheme targeting the quasi-biennal oscillation in the stratosphere. These improvements to the physical content and a more well-defined tuning strategy led to major improvements in the LMDZ6A version model climatology. Beyond the presentation of this particular model version and documentation of its climatology, the present paper underlines possible methodological pathways toward model improvement that can be shared across modeling groups.Plain Language Summary The improvement of global numerical models is essential for the anticipation of future climate changes. We present significant advances in the physical content of a particular atmospheric model which contributes to the simulations of the Coupled Model Intercomparison project CMIP that feed reports from the IPCC. We document in particular the improvements of the representation through "parameterizations" of convective and cloudy processes. The article emphasizes the importance of strengthening the formalization of the methodology of development and tuning of models, so that new physical ideas can be translated into effective improvement of the climate representation.

Abstract. Water stable isotopes in central Antarctic ice cores are critical to quantify past temperature changes. Accurate temperature reconstructions require one to understand the processes controlling surface snow isotopic composition. Isotopic fractionation processes occurring in the atmosphere and controlling snowfall isotopic composition are well understood theoretically and implemented in atmospheric models. However, post-deposition processes are poorly documented and understood. To quantitatively interpret the isotopic composition of water archived in ice cores, it is thus essential to study the continuum between surface water vapour, precipitation, surface snow and buried snow. Here, we target the isotopic composition of water vapour at Concordia Station, where the oldest EPICA Dome C ice cores have been retrieved. While snowfall and surface snow sampling is routinely performed, accurate measurements of surface water vapour are challenging in such cold and dry conditions. New developments in infrared spectroscopy enable now the measurement of isotopic composition in water vapour traces. Two infrared spectrometers have been deployed at Concordia, allowing continuous, in situ measurements for 1 month in December 2014–January 2015. Comparison of the results from infrared spectroscopy with laboratory measurements of discrete samples trapped using cryogenic sampling validates the relevance of the method to measure isotopic composition in dry conditions. We observe very large diurnal cycles in isotopic composition well correlated with temperature diurnal cycles. Identification of different behaviours of isotopic composition in the water vapour associated with turbulent or stratified regime indicates a strong impact of meteorological processes in local vapour/snow interaction. Even if the vapour isotopic composition seems to be, at least part of the time, at equilibrium with the local snow, the slope of δD against δ18O prevents us from identifying a unique origin leading to this isotopic composition.

Investigation of meteorological measurements along a 45 m tower at Dome C on the high East Antarctic Plateau revealed two distinct stable boundary layer (SBL) regimes at this location. The first regime is characterized by strong winds and continuous turbulence. It results in full vertical coupling of temperature, wind magnitude and wind direction in the SBL. The second regime is characterized by weak winds, associated with weak turbulent activity and very strong temperature inversions reaching up to 25 K in the lowest 10 m. Vertical temperature profiles are generally exponentially shaped (convex) in the first regime and ‘convex–concave–convex’ in the second. The transition between the two regimes is particularly abrupt when looking at the near‐surface temperature inversion and it can be identified by a 10 m wind‐speed threshold. With winds under this threshold, the turbulent heat supply toward the surface becomes significantly lower than the net surface radiative cooling. The threshold value (including its range of uncertainty) appears to agree with recent theoretical predictions from the so‐called ‘minimum wind speed for sustainable turbulence’ (MWST) theory. For the quasi‐steady, clear‐sky winter cases, the relation between the near‐surface inversion amplitude and the wind speed takes a characteristic ‘S’ shape. Closer analysis suggests that this relation corresponds to a ‘critical transition’ between a steady turbulent and a steady ‘radiative’ regime, with a dynamically unstable branch in the transition zone. These fascinating characteristics of the Antarctic boundary layer challenge present and future numerical models to represent this region in a physically correct manner.

Abstract. Regional climate model MAR (Modèle Atmosphérique Régional) was run for the region of Dome C located on the East Antarctic plateau, during Antarctic summer 2011-2012, in order to refine our understanding of meteorological conditions during the OPALE tropospheric chemistry campaign. A very high vertical resolution is set up in the lower troposphere, with a grid spacing of roughly 2 m. Model output is compared with temperatures and winds observed near the surface and from a 45 m high tower as well as sodar and radiation data. MAR is generally in very good agreement with the observations, but sometimes underestimates cloud formation, leading to an underestimation of the simulated downward long-wave radiation. Absorbed short-wave radiation may also be slightly overestimated due to an underestimation of the snow albedo, and this influences the surface energy budget and atmospheric turbulence. Nevertheless, the model provides sufficiently reliable information about surface turbulent fluxes, vertical profiles of vertical diffusion coefficients and boundary layer height when discussing the representativeness of chemical measurements made nearby the ground surface during field campaigns conducted at Concordia station located at Dome C (3233 m above sea level).

No abstract

Observations evidence extremely stable boundary layers (SBL) over the Antarctic Plateau and sharp regime transitions between weakly and very stable conditions. Representing such features is a challenge for climate models. This study assesses the modeling of the dynamics of the boundary layer over the Antarctic Plateau in the LMDZ general circulation model. It uses 1 year simulations with a stretched‐grid over Dome C. The model is nudged with reanalyses outside of the Dome C region such as simulations can be directly compared to in situ observations. We underline the critical role of the downward longwave radiation for modeling the surface temperature. LMDZ reasonably represents the near‐surface seasonal profiles of wind and temperature but strong temperature inversions are degraded by enhanced turbulent mixing formulations. Unlike ERA‐Interim reanalyses, LMDZ reproduces two SBL regimes and the regime transition, with a sudden increase in the near‐surface inversion with decreasing wind speed. The sharpness of the transition depends on the stability function used for calculating the surface drag coefficient. Moreover, using a refined vertical grid leads to a better reversed “S‐shaped” relationship between the inversion and the wind. Sudden warming events associated to synoptic advections of warm and moist air are also well reproduced. Near‐surface supersaturation with respect to ice is not allowed in LMDZ but the impact on the SBL structure is moderate. Finally, climate simulations with the free model show that the recommended configuration leads to stronger inversions and winds over the ice‐sheet. However, the near‐surface wind remains underestimated over the slopes of East‐Antarctica.

A conceptual model is used in combination with observational analysis to understand regime transitions of near-surface temperature inversions at night as well as in Arctic conditions. The model combines a surface energy budget with a bulk parameterization for turbulent heat transport. Energy fluxes or feedbacks due to soil and radiative heat transfer are accounted for by a ''lumped parameter closure,'' which represents the ''coupling strength'' of the system.Observations from Cabauw, Netherlands, and Dome C, Antarctica, are analyzed. As expected, inversions are weak for strong winds, whereas large inversions are found under weak-wind conditions. However, a sharp transition is found between those regimes, as it occurs within a narrow wind range. This results in a typical S-shaped dependency. The conceptual model explains why this characteristic must be a robust feature. Differences between the Cabauw and Dome C cases are explained from differences in coupling strength (being weaker in the Antarctic). For comparison, a realistic column model is run. As findings are similar to the simple model and the observational analysis, it suggests generality of the results.Theoretical analysis reveals that, in the transition zone near the critical wind speed, the response time of the system to perturbations becomes large. As resilience to perturbations becomes weaker, it may explain why, within this wind regime, an increase of scatter is found. Finally, the so-called heat flux duality paradox is analyzed. It is explained why numerical simulations with prescribed surface fluxes show a dynamical response different from more realistic surface-coupled systems.

In this work we study the dynamics of the surface‐based temperature inversion over the Antarctic Plateau during the polar winter. Using 6 years of observations from the French–Italian Antarctic station Concordia at Dome C, we investigate sudden regime transitions in the strength of the near‐surface temperature inversion. Here we define “near‐surface” as being within the domain of the 45‐m measuring tower. In particular, we consider the strongly nonlinear relation between the 10‐m inversion strength ( T 10m – T s ) and the 10‐m wind speed. To this end, all individual events for which the 10‐m inversion strength increases or decreases continuously by more than 15 K in time are considered. Composite time series and vertical profiles of wind and temperature reveal specific characteristics of the transition from weak to very strong inversions and vice versa. In contrast to midlatitudes, the largest variations in temperature are not found at the surface but at a height of 10 m. A similar analysis was performed on results from an atmospheric single‐column model (SCM). Overall, the SCM results reproduce the observed characteristics of the transitions in the near‐surface inversion remarkably well. Using model output, the underlying mechanisms of the regime transitions are identified. The nonlinear relation between inversion strength and wind speed at a given level is explained by variations in the geostrophic wind speed, changes in the depth of the turbulent layer and the vertical divergence of turbulent fluxes. Moreover, the transitions between different boundary layer regimes cannot be explained without considering the contribution of subsidence heating.