ALBERT

Description: As knowledge about the cirrus clouds in the lower stratosphere is limited, reliable long-term measurements are needed to assess their characteristics, radiative impact and important role in upper troposphere and lower stratosphere (UTLS) chemistry. We used 6 years (2006–2012) of Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) measurements to investigate the global and seasonal distribution of stratospheric cirrus clouds and compared the MIPAS results with results derived from the latest version (V4.x) of the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) data. For the identification of stratospheric cirrus clouds, precise information on both the cloud top height (CTH) and the tropopause height is crucial. Here, we used lapse rate tropopause heights estimated from the ERA-Interim global reanalysis. Considering the uncertainties of the tropopause heights and the vertical sampling grid, we define CTHs more than 0.5 km above the tropopause as stratospheric for CALIPSO data. For MIPAS data, we took into account the coarser vertical sampling grid and the broad field of view so that we considered cirrus CTHs detected more than 0.75 km above the tropopause as stratospheric. Further sensitivity tests were conducted to rule out sampling artefacts in MIPAS data. The global distribution of stratospheric cirrus clouds was derived from night-time measurements because of the higher detection sensitivity of CALIPSO. In both data sets, MIPAS and CALIPSO, the stratospheric cirrus cloud occurrence frequencies are significantly higher in the tropics than in the extra-tropics. Tropical hotspots of stratospheric cirrus clouds associated with deep convection are located over equatorial Africa, South and Southeast Asia, the western Pacific, and South America. Stratospheric cirrus clouds were more often detected in December–February (15 %) than June–August (8 %) in the tropics (±20∘). At northern and southern middle latitudes (40–60∘), MIPAS observed about twice as many stratospheric cirrus clouds (occurrence frequencies of 4 %–5 % for MIPAS rather than about 2 % for CALIPSO). We attribute more frequent observations of stratospheric cirrus clouds with MIPAS to the higher detection sensitivity of the instrument to optically thin clouds. In contrast to the difference between daytime and night-time occurrence frequencies of stratospheric cirrus clouds by a factor of about 2 in zonal means in the tropics (4 % and 10 %, respectively) and at middle latitudes for CALIPSO data, there is little diurnal cycle in MIPAS data, in which the difference of occurrence frequencies in the tropics is about 1 percentage point in zonal mean and about 0.5 percentage point at middle latitudes. The difference between CALIPSO day and night measurements can also be attributed to their differences in detection sensitivity. Future work should focus on better understanding the origin of the stratospheric cirrus clouds and their impact on radiative forcing and climate.

Description: Polar stratospheric clouds (PSCs) play a key role in polar ozone depletion in the stratosphere. Improved observations and continuous monitoring of PSCs can help to validate and improve chemistry–climate models that are used to predict the evolution of the polar ozone hole. In this paper, we explore the potential of applying machine learning (ML) methods to classify PSC observations of infrared limb sounders. Two datasets were considered in this study. The first dataset is a collection of infrared spectra captured in Northern Hemisphere winter 2006/2007 and Southern Hemisphere winter 2009 by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) instrument on board the European Space Agency's (ESA) Envisat satellite. The second dataset is the cloud scenario database (CSDB) of simulated MIPAS spectra. We first performed an initial analysis to assess the basic characteristics of the CSDB and to decide which features to extract from it. Here, we focused on an approach using brightness temperature differences (BTDs). From both the measured and the simulated infrared spectra, more than 10 000 BTD features were generated. Next, we assessed the use of ML methods for the reduction of the dimensionality of this large feature space using principal component analysis (PCA) and kernel principal component analysis (KPCA) followed by a classification with the support vector machine (SVM). The random forest (RF) technique, which embeds the feature selection step, has also been used as a classifier. All methods were found to be suitable to retrieve information on the composition of PSCs. Of these, RF seems to be the most promising method, being less prone to overfitting and producing results that agree well with established results based on conventional classification methods.

Description: Tropical volcanic eruptions have been widely studied for their significant contribution to stratospheric aerosol loading and global climate impacts, but the impact of high-latitude volcanic eruptions on the stratospheric aerosol layer is not clear and the pathway of transporting aerosol from high latitudes to the tropical stratosphere is not well understood. In this work, we focus on the high-latitude volcano Sarychev (48.1° N, 153.2° E), which erupted in June 2009, and the influence of the Asian summer monsoon (ASM) on the equatorward dispersion of the volcanic plume. First, the sulfur dioxide (SO2) emission time series and plume height of the Sarychev eruption are estimated with SO2 observations of the Atmospheric Infrared Sounder (AIRS) and a backward trajectory approach using the Lagrangian particle dispersion model Massive–Parallel Trajectory Calculations (MPTRAC). Then, the transport and dispersion of the plume are simulated using the derived SO2 emission time series. The transport simulations are compared with SO2 observations from AIRS and validated with aerosol observations from the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS). The MPTRAC simulations show that about 4 % of the sulfur emissions were transported to the tropical stratosphere within 50 days after the beginning of the eruption, and the plume dispersed towards the tropical tropopause layer (TTL) through isentropic transport above the subtropical jet. The MPTRAC simulations and MIPAS aerosol data both show that between the potential temperature levels of 360 and 400 K, the equatorward transport was primarily driven by anticyclonic Rossby wave breaking enhanced by the ASM in boreal summer. The volcanic plume was entrained along the anticyclone flows and reached the TTL as it was transported southwestwards into the deep tropics downstream of the anticyclone. Further, the ASM anticyclone influenced the pathway of aerosols by isolating an aerosol hole inside of the ASM, which was surrounded by aerosol-rich air outside. This transport barrier was best indicated using the potential vorticity gradient approach. Long-term MIPAS aerosol detections show that after entering the TTL, aerosol from the Sarychev eruption remained in the tropical stratosphere for about 10 months and ascended slowly. The ascent speed agreed well with the ascent speed of the water vapor tape recorder. Furthermore, a hypothetical MPTRAC simulation for a wintertime eruption was carried out. It is shown that under winter atmospheric circulations, the equatorward transport of the plume would be suppressed by the strong subtropical jet and weak wave breaking events. In this hypothetical scenario, a high-latitude volcanic eruption would not be able to contribute to the tropical stratospheric aerosol layer.

Description: Volcanic sulfate aerosol is an important source of sulfur for Antarctica, where other local sources of sulfur are rare. Midlatitude and high-latitude volcanic eruptions can directly influence the aerosol budget of the polar stratosphere. However, tropical eruptions can also enhance polar aerosol load following long-range transport. In the present work, we analyze the volcanic plume of a tropical eruption, Mount Merapi in 2010, and investigate the transport pathway of the volcanic aerosol from the tropical tropopause layer (TTL) to the lower stratosphere over Antarctica. We use the Lagrangian particle dispersion model Massive-Parallel Trajectory Calculations (MPTRAC) and Atmospheric Infrared Sounder (AIRS) SO2 measurements to reconstruct the altitude-resolved SO2 injection time series during the explosive eruption period and simulate the transport of the volcanic plume using the MPTRAC model. AIRS SO2 and aerosol measurements, the aerosol cloud index values provided by Michelson Interferometer for Passive Atmospheric Sounding (MIPAS), are used to verify and complement the simulations. The Lagrangian transport simulation of the volcanic plume is compared with MIPAS aerosol measurements and shows good agreement. Both the simulations and the observations presented in this study suggest that volcanic plumes from the Merapi eruption were transported to the south of 60∘ S 1 month after the eruption and even further to Antarctica in the following months. This relatively fast meridional transport of volcanic aerosol was mainly driven by quasi-horizontal mixing from the TTL to the extratropical lower stratosphere, and most of the quasi-horizontal mixing occurred between the isentropic surfaces of 360 to 430 K. When the plume went to Southern Hemisphere high latitudes, the polar vortex was displaced from the South Pole, so that the volcanic plume was carried to the South Pole without penetrating the polar vortex. Although only 4 % of the sulfur injected by the Merapi eruption was transported into the lower stratosphere south of 60∘ S, the Merapi eruption contributed up to 8800 t of sulfur to the Antarctic lower stratosphere. This indicates that the long-range transport under favorable meteorological conditions enables a moderate tropical volcanic eruption to be an important remote source of sulfur for the Antarctic stratosphere.

Description: The European Centre for Medium-Range Weather Forecasts' (ECMWF's) next-generation reanalysis ERA5 provides many improvements, but it also confronts the community with a “big data” challenge. Data storage requirements for ERA5 increase by a factor of ∼80 compared with the ERA-Interim reanalysis, introduced a decade ago. Considering the significant increase in resources required for working with the new ERA5 data set, it is important to assess its impact on Lagrangian transport simulations. To quantify the differences between transport simulations using ERA5 and ERA-Interim data, we analyzed comprehensive global sets of 10-day forward trajectories for the free troposphere and the stratosphere for the year 2017. The new ERA5 data have a considerable impact on the simulations. Spatial transport deviations between ERA5 and ERA-Interim trajectories are up to an order of magnitude larger than those caused by parameterized diffusion and subgrid-scale wind fluctuations after 1 day and still up to a factor of 2–3 larger after 10 days. Depending on the height range, the spatial differences between the trajectories map into deviations as large as 3 K in temperature, 30 % in specific humidity, 1.8 % in potential temperature, and 50 % in potential vorticity after 1 day. Part of the differences between ERA5 and ERA-Interim is attributed to the better spatial and temporal resolution of the ERA5 reanalysis, which allows for a better representation of convective updrafts, gravity waves, tropical cyclones, and other meso- to synoptic-scale features of the atmosphere. Another important finding is that ERA5 trajectories exhibit significantly improved conservation of potential temperature in the stratosphere, pointing to an improved consistency of ECMWF's forecast model and observations that leads to smaller data assimilation increments. We conducted a number of downsampling experiments with the ERA5 data, in which we reduced the numbers of meteorological time steps, vertical levels, and horizontal grid points. Significant differences remain present in the transport simulations, if we downsample the ERA5 data to a resolution similar to ERA-Interim. This points to substantial changes of the forecast model, observations, and assimilation system of ERA5 in addition to improved resolution. A comparison of two Lagrangian trajectory models allowed us to assess the readiness of the codes and workflows to handle the comprehensive ERA5 data and to demonstrate the consistency of the simulation results. Our results will help to guide future Lagrangian transport studies attempting to navigate the increased computational complexity and leverage the considerable benefits and improvements of ECMWF's new ERA5 data set.

Description: Volcanic eruptions can increase the stratospheric sulfur loading by orders of magnitude above the background level and are the most important source of variability in stratospheric sulfur. We present a set of vertical profiles of sulfate aerosol volume densities and derived liquid-phase H2SO4 (sulfuric acid) mole fractions for 2005–2012, retrieved from infrared limb emission measurements performed with the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) on board of the Environmental Satellite (Envisat). Relative to balloon-borne in situ measurements of aerosol at Laramie, Wyoming, the MIPAS aerosol data have a positive bias that has been corrected, based on the observed differences to the in situ data. We investigate the production of stratospheric sulfate aerosol from volcanically emitted SO2 for two case studies: the eruptions of Kasatochi in 2008 and Sarychev in 2009, which both occurred in the Northern Hemisphere midlatitudes during boreal summer. With the help of chemical transport model (CTM) simulations for the two volcanic eruptions we show that the MIPAS sulfate aerosol and SO2 data are qualitatively and quantitatively consistent with each other. Further, we demonstrate that the lifetime of SO2 is explained well by its oxidation by hydroxyl radicals (OH). While the sedimentation of sulfate aerosol plays a role, we find that the long-term decay of stratospheric sulfur after these volcanic eruptions in midlatitudes is mainly controlled by transport via the Brewer–Dobson circulation. Sulfur emitted by the two midlatitude volcanoes resides mostly north of 30∘ N at altitudes of ∼ 10–16 km, while at higher altitudes (∼ 18–22 km) part of the volcanic sulfur is transported towards the Equator where it is lifted into the stratospheric “overworld” and can further be transported into both hemispheres.

Description: Volcanic sulfate aerosol is an important source of sulfur for Antarctica where other local sources of sulfur are rare. Mid- and high latitude volcanic eruptions can directly influence the aerosol budget of the polar stratosphere. However, tropical eruptions can also enhance polar aerosol load following long-range transport. In the present work, we analyze the volcanic plume of a tropical eruption, Mount Merapi in October 2010, using the Lagrangian particle dispersion model Massive-Parallel Trajectory Calculations (MPTRAC), Atmospheric Infrared Sounder (AIRS) SO2 observations and Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) aerosol observations. We investigate the pathway and transport efficiency of the volcanic aerosol from the tropical tropopause layer (TTL) to the lower stratosphere over Antarctica. We first estimated the time- and height-resolved SO2 injection time series over Mount Merapi during the explosive eruption using the AIRS SO2 observations and a backward trajectory approach. Then the SO2 injections were tracked for up to 6 months using the MPTRAC model. The Lagrangian transport simulation of the volcanic plume was compared to MIPAS aerosol observations and showed good agreement. Both of the simulation and the observations presented in this study suggest that a significant amount of aerosols of the volcanic plume from the Merapi eruption was transported from the tropics to the south of 60°S within one month after the eruption and even further to Antarctica in the following two months. This relatively fast meridional transport of volcanic aerosol was mainly driven by quasi-horizontal mixing from the TTL to the extratropical lower stratosphere, which was facilitated by the weakening of the subtropical jet during the seasonal transition from austral spring to summer and linked to the westerly phase of the quasi-biennial oscillation (QBO). When the plume went to southern high latitudes, the polar vortex was displaced from the south pole, so the volcanic plume was carried to the south pole without penetrating the polar vortex. Based on the model results, the most efficient pathway for the quasi-horizontal mixing was in between the isentropic surfaces of 360 and 430K. Although only 4% of the initial SO2 load was transported into the lower stratosphere south of 60°S, the Merapi eruption contributed about 8800 tons of sulfur to the Antarctic lower stratosphere. This indicates that the long-range transport under favorable meteorological conditions enables tropical volcanic eruptions to be an important remote source of sulfur for the Antarctic stratosphere.

Description: Infrared limb emission instruments have a long history in measuring clouds and aerosol. In particular the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) instrument aboard ESA's Envisat provides 10 years of altitude resolved global measurements. Previous studies found systematic over- and underestimations of cloud top heights for cirrus and polar stratospheric clouds. To assess the cloud top height information and to characterize its uncertainty for the MIPAS instrument we performed simulations for ice clouds, volcanic ash, and sulfate aerosol. From the simulation results we found that in addition to the known effects of the field-of-view that can lead to a cloud top height overestimation, and broken cloud conditions that can lead to underestimation, also the cloud extinction plays an important role. While for optically thick clouds the possible cloud top height overestimation for MIPAS reaches up to 1.6 km due to the field-of-view, for optically thin clouds and aerosol the systematic underestimation reaches 5.1 km. For the detection sensitivity and the degree of underestimation of the MIPAS measurements also the cloud layer thickness plays a role. 1 km thick clouds are detectable down to extinctions of 5 x 10-4 km-1 and 6 km thick clouds are detectable down to extinctions of 1 x 10-4 km-1, where the largest underestimations of the cloud top height occur for the optically thinnest clouds with a vertical extent of 6 km. The relation between extinction coefficient, cloud top height estimate, and layer thickness is confirmed by a comparison of MIPAS cloud top heights of the volcanic sulfate aerosol from the Nabro eruption in 2011 with space- and ground-based lidar measurements and twilight measurements between June 2011 and February 2012. In the fresh to two months old plume, where the extinction was between 1 x 10-4 and 7 x 10-4 km-1 and the layer thickness mostly below 4 km, we found for MIPAS an average underestimation of 1.1 km. In the aged plume with extinctions down to 5 x 10-5 km-1 and layer thicknesses of up to 9.5 km the underestimation was higher reaching 7.2 km. The dependency of the cloud top height over- or underestimation on the extinction coefficient can explain seemingly contradictory results of previous studies. In spite of the relatively large uncertainty range of the cloud top height, the comparison of the detection sensitivity towards sulfate aerosol between MIPAS and a suite of widely used UV/VIS limb and IR nadir satellite aerosol measurements shows that MIPAS provides complementary information in terms of detection sensitivity.

Description: We analyse polar stratospheric cloud (PSC) signatures in airborne MIPAS-STR (Michelson Interferometer for Passive Atmospheric Sounding – STRatospheric aircraft) observations under conditions suitable for the existence of nitric acid trihydrate (NAT) above northern Scandinavia on 11 December 2011. The high resolution infrared limb emission spectra of MIPAS-STR show a characteristic "shoulder-like" signature in the spectral region around 820 cm−1, which is attributed to the ν2 deformation mode of NO3− in β-NAT. Using radiative transfer calculations involving Mie and T-Matrix methods, the spectral signatures of spherical and aspherical particles are simulated. The simulations are constrained using collocated in-situ particle measurements. Simulations assuming highly aspherical spheroids with aspect ratios (AR) of 0.1 or 10.0 and a lognormal particle mode with a mode radius of 4.8 µm reproduce the observed spectra to a high degree. A smaller lognormal mode with a mode radius of 2.0 µm, which was also taken into account, plays only a minor role. Best overall agreement is found for elongated spheroids with AR = 0.1. Simulations of spherical particles and spheroids with AR = 0.5 and 2.0 return results very similar to each other and do not allow to reproduce the signature around 820 cm−1. The observed “shoulder-like” signature is explained by the combination of the emission and scattering characteristics of large highly aspherical β-NAT particles. The size distribution supported by our results corresponds with ~9 ppbv of gas-phase equivalent HNO3 at the flight altitude of ~18.5 km. The results are compared with the size distributions derived from the in-situ observations, a corresponding Chemical Lagrangian Model of the Stratosphere (CLaMS) simulation, and excess gas-phase HNO3 observed in a nitrification layer directly below the observed PSC. The presented results suggest that large highly aspherical β-NAT particles involved in denitrification of the polar stratosphere can be identified by means of passive infrared limb emission measurements.