The impact of deep convection representation in a global atmospheric model on the warm conveyor belt and jet stream during NAWDEX IOP6

The effect of parameterized deep convection on warm conveyor belt (WCB) activity and the jet stream is investigated by performing simulations of an explosively developing large-scale cyclone that occurred during the North Atlantic Waveguide and Downstream Impact Experiment (NAWDEX) field campaign using the Météo-France global atmospheric model ARPEGE. Three simulations differing only from their deep convection representation are analysed. The first one was performed with the Bougeault (1985) scheme (B85), the second one with the Prognostic Condensates Microphysics and Transport (PCMT) scheme of Piriou et al. (2007), and the third one without any parameterized deep convection. In the latter simulation, the release of convective instability at the resolved scales of the model generates localized cells marked by strong heating with few degrees extent in longitude and latitude along the fronts. In runs with active parameterized deep convection (B85, PCMT), the heating rate is more homogeneously distributed along fronts as the instability release happens at subgrid scales. This difference leads to more rapid and abrupt ascents in the WCB without parameterized deep convection and more moderate but more sustained ascents with parameterized deep convection. While the number of WCB trajectories does not differ much between the three simulations, the averaged heating rates over the WCB trajectories exhibits distinct behaviour. After 1 d of simulations, the upper-level heating rate is on average larger, with the B85 scheme leading to stronger potential vorticity (PV) destruction. The difference comes from the resolved sensible and latent heating and not the parameterized one. A comparison with (re)analyses and a large variety of airborne observations from the NAWDEX field campaign (Doppler radar, Doppler lidar, dropsondes) made during the coordinated flights of two aircraft in the WCB outflow region shows that B85 performs better in the representation of the double jet structure at 1 d lead time than the other two simulations. That can be attributed to the more active WCB at upper levels. However, this effect is too strong and that simulation becomes less realistic than the other ones at forecast ranges beyond 1.5 d. The simulation with the PCMT scheme has an intermediate behaviour between the one with the B85 scheme and without parameterized deep convection, but its impact on the jet stream is closer to the latter one. Finally, additional numerical experiments show that main differences in the impact on the jet between PCMT and B85 largely come from the chosen closure, with the former being based on CAPE and the latter on moisture convergence.

Variations of Tropical Lapse Rates in Climate Models and their Implications for Upper Tropospheric Warming

The vertical temperature structure in the tropics is primarily set by convection and therefore follows a moist adiabat to first order. However, tropical upper tropospheric temperatures differ among climate models and observations, as atmospheric convection remains poorly understood. Here, we quantify the variations in tropical lapse rates in CMIP6 models and explore reasons for these variations. We find that differences in surface temperatures weighted by the regions of strongest convection cannot explain these variations and therefore we hypothesise that the representation of convection itself and associated small scale processes are responsible. We reproduce these variations in perturbed physics experiments with the global atmospheric model ICON-A, in which we vary autoconversion and entrainment parameters. For smaller autoconversion values, additional freezing enthalpy from the cloud water that is not precipitated warms the upper troposphere. Smaller entrainment rates also lead to a warmer upper troposphere, as convection and thus latent heating reaches higher. Furthermore, we show that according to most radiosonde datasets all CMIP6 AMIP simulations overestimate recent upper tropospheric warming. Additionally, all radiosonde datasets agree that climate models on average overestimate the amount of upper tropospheric warming for a given lower tropospheric warming. We demonstrate that increased entrainment rates reduce this overestimation, likely because of the reduction of latent heat release in the upper troposphere. Our results suggest that imperfect convection parameterisations are responsible for a considerable part of the variations in tropical lapse rates and also part of the overestimation of warming compared to the observations.

The impact of spatially varying wetland source signatures on the atmospheric variability of δ D-CH 4

We present the first spatially resolved distribution of the δ D-CH 4 signature of wetland methane emissions and assess its impact on atmospheric δ D-CH 4 . The δ D-CH 4 signature map is derived by relating δ D-H 2 O of precipitation to measured δ D-CH 4 of methane wetland emissions at a variety of wetland types and locations. This results in strong latitudinal variation in the wetland δ D-CH 4 source signature. When δ D-CH 4 is simulated in a global atmospheric model, little difference is found in global mean, inter-hemispheric difference and seasonal cycle if the spatially varying δ D-CH 4 source signature distribution is used instead of a globally uniform value. This is because atmospheric δ D-CH 4 is largely controlled by OH fractionation. However, we show that despite these small differences, using atmospheric records of δ D-CH 4 to infer changes in the wetland emissions distribution requires the use of the more accurate spatially varying δ D-CH 4 source signature. We find that models will only be sensitive to changes in emissions distribution if spatial information can be exploited through the spatially resolved source signatures. In addition, we also find that on a regional scale, at sites measuring excursions of δ D-CH 4 from background levels, substantial differences are simulated in atmospheric δ D-CH 4 if using spatially varying or uniform source signatures. This article is part of a discussion meeting issue ‘Rising methane: is warming feeding warming? (part 1)’.

fv3gfs-wrapper: a Python wrapper of the FV3GFS atmospheric model

Simulation software in geophysics is traditionally written in Fortran or C++ due to the stringent performance requirements these codes have to satisfy. As a result, researchers who use high-productivity languages for exploratory work often find these codes hard to understand, hard to modify, and hard to integrate with their analysis tools. fv3gfs-wrapper is an open-source Python-wrapped version of the NOAA (National Oceanic and Atmospheric Administration) FV3GFS (Finite-Volume Cubed-Sphere Global Forecast System) global atmospheric model, which is coded in Fortran. The wrapper provides simple interfaces to progress the Fortran main loop and get or set variables used by the Fortran model. These interfaces enable a wide range of use cases such as modifying the behavior of the model, introducing online analysis code, or saving model variables and reading forcings directly to and from cloud storage. Model performance is identical to the fully compiled Fortran model, unless routines to copy the state in and out of the model are used. This copy overhead is well within an acceptable range of performance and could be avoided with modifications to the Fortran source code. The wrapping approach is outlined and can be applied similarly in other Fortran models to enable more productive scientific workflows.

Evaluation of coastal Antarctic precipitation in LMDz6 global atmospheric model using ground-based radar observations

In the current context of climate change in the poles, one of the objectives of the APRES3 (Antarctic Precipitation Remote Sensing from Surface and Space) project was to characterize the vertical structure of precipitation in order to better simulate it. Precipitation simulated by models in Antarctica is currently very widespread and it overestimates the data. Sensitivity studies have been conducted using a global climate model and compared to the observations obtained at the Dumont d’Urville coast station, obtained by a Micro Rain Radar (MRR). The LMDz/IPSL general circulation model, with zoomed configuration over Dumont d’Urville, has been considered for this study. A sensitivity study was conducted on the physical and numerical parameters of the LMDz model with the aim of estimating their contribution to the precipitation simulation. Sensitivity experiments revealed that changes in the sedimentation and sublimation parameters do not significantly impact precipitation rate. However, dissipation of the LMDz model, which is a numerical process that dissipates spatially excessive energy and keeps the model stable, impacts precipitation indirectly but very strongly. A suitable adjustment of the dissipation reduces significantly precipitation over Antarctic peripheral area, thus providing a simulated profile in better agreement with the MRR observations.

The impact of deep convection representation in a global atmospheric model on the warm conveyor belt and jet stream during NAWDEX IOP6

The effect of parameterized deep convection on warm conveyor belt (WCB) activity and jet stream is investigated by performing simulations of an explosively-developing large-scale cyclone that occurred during the North Atlantic Waveguide and Downstream Impact Experiment (NAWDEX) field campaign using the Météo-France global atmospheric model ARPEGE. Three simulations differing only from their deep convection representation are analysed. The first one was performed with the Bougeault et al. (1985) scheme (B85), the second one with the Prognostic Condensates Microphysics and Transport (PCMT) scheme of Piriou et al. (2007), and the third one without any parameterized deep convection. In the latter simulation, the release of convective instability at the resolved scales of the model generates localized cells marked by strong heating with few degrees extent in longitude and latitude along the fronts. In runs with active parameterized deep convection (B85, PCMT), the heating rate is more homogeneously distributed along fronts as the instability release happens at sub-grid scales. This difference leads to more rapid and abrupt ascents in the WCB without parameterized deep convection, and more moderate but more sustained ascents with parameterized deep convection. While the number of WCB trajectories does not differ much between the three simulations, the averaged heating rates over the WCB trajectories exhibits distinct behavior. After one day of simulations, the upper-level heating rate is in average larger with B85 scheme leading to stronger potential vorticity (PV) destruction. The difference comes from the large-scale heating and not the parameterized heating.A comparison with (re)analyses and a large variety of airborne observations from the NAWDEX field campaign (Doppler radar, Doppler lidar, dropsondes) made during the coordinated flights of two aircraft in the WCB outflow region shows that B85 performs better in the representation of the double jet structure at 1-day lead time than the other two simulations. That can be attributed to the more active WCB at upper levels. However this effect is too strong and that simulation becomes less realistic at longer forecast range (1.5 to 2 days) than the other ones. The simulation with PCMT scheme has an intermediate behavior between the one with B85 scheme and without parameterized deep convection but its impact on the jet stream is closer to the latter one. Finally, additional numerical experiments show that main differences in the impact on the jet between PCMT and B85 largely come from the chosen closure, the former being based on CAPE and the latter on moisture convergence.

fv3gfs-wrapper: a Python wrapper of the FV3GFS atmospheric model

Simulation software in geophysics is traditionally written in Fortran or C++ due to the stringent performance requirements these codes have to satisfy. As a result, these codes are often hard to understand, hard to modify and hard to interface with high-productivity languages used for exploratory work. fv3gfs-wrapper is an open-source Python-wrapped version of NOAA's FV3GFS global atmospheric model, which is coded in Fortran. The wrapper provides simple interfaces to progress the Fortran main loop and get or set state from the Fortran model. These interfaces enable a wide range of use cases such as modifying the behavior of the model, introducing online analysis code, or saving model state and reading forcings directly to and from cloud storage. Model performance is identical to the fully-compiled Fortran model, unless routines to copy state in and out of the model are used. This copy overhead is well within an acceptable range of performance, and could be avoided with modifications to the Fortran source code. The wrapping approach is outlined and can be applied similarly in other Fortran models to enable more productive scientific workflows.

Topography-correlated atmospheric signal mitigation for InSAR applications in the Tibetan plateau based on global atmospheric models

<p>Atmospheric heterogeneity mainly exposes itself as tropospheric phase delay in satellite interferometric synthetic aperture radar (InSAR) observations, which smears or even overshadows the deformation component of InSAR measurements. In this study, we estimated the performance of four global atmospheric models (GAMs), i.e. ERA5, ERA-Interim (ERA-I), MERRA2 and GACOS, for tropospheric phase delay reduction in InSAR applications in the Tibetan plateau, of which ERA5 is the latest global atmospheric model released by ECMWF. We demonstrated the effectiveness of atmospheric phase screen (APS) correction using the four GAMs for more than 700 Sentinel-1 TOPS interferograms covering two study areas in the southern (R1) and northwest margins (R2) of the Tibetan plateau. Topography-correlated signals have been widely observed in these interferograms, which are most likely due to the APS effects. We calculated the standard deviations (STD) and correlation coefficients between InSAR Line of Sight (LOS) measurements and topography before and after applying APS correction. The results show that the STDs of non-deformation areas from the GAMs decrease to ~4 mm from ~10 mm and ~12 mm originally on average for R1 and R2, respectively, and the correlation coefficients after the APS correction are reduced below 0.4 from ~0.8 for the selected interferometric pairs. In addition, as the newly released GAM, ERA5 has similar performance with GACOS products and outperforms other models generally. This suggests that GAMs, particularly ERA5, have great potentials in the APS correction for InSAR applications in the Tibetan plateau.</p>