Changes in stratospheric aerosol extinction coefficient after the 2018 Ambae eruption as seen by OMPS-LP and MAECHAM5-HAM

Stratospheric aerosols are an important component of the climate system. They not only change the radiative budget of the Earth but also play an essential role in ozone depletion. These impacts are particularly noticeable after volcanic eruptions when SO2 injected with the eruption reaches the stratosphere, oxidizes, and forms stratospheric aerosol. There have been several studies in which a volcanic eruption plume and the associated radiative forcing were analyzed using climate models and/or data from satellite measurements. However, few have compared vertically and temporally resolved volcanic plumes using both measured and modeled data. In this paper, we compared changes in the stratospheric aerosol loading after the 2018 Ambae eruption observed by satellite remote sensing measurements and simulated by a global aerosol model. We use vertical profiles of the aerosol extinction coefficient at 869 nm retrieved at the Institute of Environmental Physics (IUP) in Bremen from OMPS-LP (Ozone Mapping and Profiling Suite – Limb Profiler) observations. Here, we present the retrieval algorithm and a comparison of the obtained profiles with those from SAGE III/ISS (Stratospheric Aerosol and Gas Experiment III on board the International Space Station). The observed differences are within 25 % for most latitude bins, which indicates a reasonable quality of the retrieved limb aerosol extinction product. The volcanic plume evolution is investigated using both monthly mean aerosol extinction coefficients and 10 d averaged data. The measurement results were compared with the model output from MAECHAM5-HAM (ECHAM for short). In order to simulate the eruption accurately, we use SO2 injection estimates from OMPS and OMI (Ozone Monitoring Instrument) for the first phase of eruption and the TROPOspheric Monitoring Instrument (TROPOMI) for the second phase. Generally, the agreement between the vertical and geographical distribution of the aerosol extinction coefficient from OMPS-LP and ECHAM is quite remarkable, in particular, for the second phase. We attribute the good consistency between the model and the measurements to the precise estimation of injected SO2 mass and height, as well as to the nudging to ECMWF ERA5 reanalysis data. Additionally, we compared the radiative forcing (RF) caused by the increase in the aerosol loading in the stratosphere after the eruption. After accounting for the uncertainties from different RF calculation methods, the RFs from ECHAM and OMPS-LP agree quite well. We estimate the tropical (20∘ N to 20∘ S) RF from the second Ambae eruption to be about −0.13 W m−2.