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.

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 δ¹⁸O prevents us from identifying a unique origin leading to this isotopic composition.

This work is motivated by the identification of the land‐atmosphere interactions as one of the key sources of uncertainty in climate change simulations. It documents new developments in related processes, namely, boundary layer/convection/clouds parameterizations and land surface parameterization in the Earth System Model of the Institut Pierre Simon Laplace (IPSL). Simulations forced by prescribed oceanic conditions are produced with different combinations of atmospheric and land surface parameterizations. They are used to explore the sensitivity to the atmospheric physics and/or soil physics of major biases in the near surface variables over continents, the energy and moisture coupling established at the soil/atmosphere interface in not too wet (energy limited) and not too dry (moisture limited) soil moisture regions also known as transition or “hot‐spot” regions, the river runoff at the outlet of major rivers. The package implemented in the IPSL‐Climate Model for the Phase 6 of the Coupled Models Intercomparison Project (CMIP6) allows us to reduce several biases in the surface albedo, the snow cover, and the continental surface air temperature in summer as well as in the temperature profile in the surface layer of the polar regions. The interactions between soil moisture and atmosphere in hotspot regions are in better agreement with the observations. Rainfall is also significantly improved in volume and seasonality in several major river basins leading to an overall improvement in river discharge. However, the lack of consideration of floodplains and human influences in the model, for example, dams and irrigation, impacts the realism of simulated discharge.

While most precipitation in Antarctica falls as snow, little is known about liquid precipitation, although it can have ecological and climatic impacts. This study combines meteorological reports at 10 stations with the ERA5 reanalysis to provide a climatological characterization of rainfall occurrence over Antarctica. Along the East Antarctic coast, liquid precipitation occurs 22 days per year at most and coincides with maritime intrusions and blocking anticyclones. Over the north‐western Antarctic Peninsula, rainfall occurs more than 50 days per year on average and the recent summer cooling was accompanied by a decrease of −35 annual rainy days per decade between 1998 and 2015 at Faraday‐Vernadsky. Projections from seven latest‐generation climate models reveal that Antarctic coasts will experience a warming and more frequent and intense rainfall by the end of the century. Rainfall is expected to impact new regions of the continent, increasing their vulnerability to melting by the preconditioning of surface snow.

Climate models exhibit major radiative biases over the Southern Ocean owing to a poor representation of mixed‐phase clouds. This study uses the remote‐sensing dataset from the Measurements of Aerosols, Radiation and Clouds over the Southern Ocean (MARCUS) campaign to assess the ability of the Weather Research and Forecasting (WRF) model to reproduce frontal clouds off Antarctica. It focuses on the modeling of thin mid‐level supercooled liquid water layers which precipitate ice. The standard version of WRF produces almost fully glaciated clouds and cannot reproduce cloud top turbulence. Our work demonstrates the importance of adapting the ice nucleation parameterization to the pristine austral atmosphere to reproduce the supercooled liquid layers. Once simulated, droplets significantly impact the cloud radiative effect by increasing downwelling longwave fluxes and decreasing downwelling shortwave fluxes at the surface. The net radiative effect is a warming of snow and ice covered surfaces and a cooling of the ocean. Despite improvements in our simulations, the local turbulent circulation related to cloud‐top radiative cooling is not properly reproduced, advocating for the need to develop a parameterization for top‐down convection to capture the turbulence‐microphysics interplay at cloud top.

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.