REtrieval Method for optical and physical Aerosol Properties in the stratosphere (REMAPv1)

Abstract. Stratospheric aerosol is an important climate forcing agent as it scatters some of the incoming solar radiation back to space, thus cooling the Earth's surface and the troposphere. At the same time, it absorbs some of the upwelling terrestrial radiation that heats the stratosphere. It also plays an important role in stratospheric ozone chemistry by hosting heterogeneous reactions. Major volcanic eruptions can cause strong perturbations of stratospheric aerosol, changing its radiative and chemical effects by more than an order of magnitude. Many global climate models require prescribed stratospheric aerosol as input to properly simulate both climate effects in the presence and absence of volcanic eruptions. This paper describes REMAP, a retrieval method and code for aerosol properties that has been used in several model intercomparison projects (under the name Stratospheric Aerosol and Gas Experiment-3λ, SAGE-3λ). The code fits a single-mode lognormal size distribution for a pure aqueous sulfuric acid aerosol to aerosol extinction coefficients from observational or model data sets. From the retrieved size distribution parameters, the code calculates the effective radius; surface area density; and extinction coefficients, single-scattering albedos, and asymmetry factors of the aerosol within the wavelength bands specified for each individual climate model. We validate REMAP using balloon-borne observations after the Mount Pinatubo and Hunga Tonga–Hunga Ha`apai (HTHH) volcanic eruptions, as well as 4 decades of lidar measurements. Within the constraints of a single-mode lognormal distribution, REMAP generates realistic effective radii and surface area densities after volcanic eruptions and generally matches the lidar backscatter time series within measurement uncertainty. Deviations in aerosol backscatter by up to a factor of 2 arise when (non-volcanic) tropospheric intrusions (e.g., from wildfires) are present and the size distribution deviates significantly from the single-mode lognormal type. We describe the products that have been used in CCMI (Chemistry–Climate Model Initiative), CMIP6 (Coupled Model Intercomparison Project Phase 6 ), and other model intercomparison projects and provide practical instructions for use of the code in future applications.