Sulphur dioxide (SO₂) fluxes of active degassing volcanoes are routinely measured with ground-based equipment to characterize and monitor volcanic activity. SO₂ of unmonitored volcanoes or from explosive volcanic eruptions, can be measured with satellites. However, remote-sensing methods based on absorption spectroscopy generally provide integrated amounts of already dispersed plumes of SO₂ and satellite derived flux estimates are rarely reported. Here we review a number of different techniques to derive volcanic SO₂ fluxes using satellite measurements of plumes of SO₂ and investigate the temporal evolution of the total emissions of SO₂ for three very different volcanic events in 2011: Puyehue-Cordón Caulle (Chile), Nyamulagira (DR Congo) and Nabro (Eritrea). High spectral resolution satellite instruments operating both in the ultraviolet-visible (OMI/Aura and GOME-2/MetOp-A) and thermal infrared (IASI/MetOp-A) spectral ranges, and multispectral satellite instruments operating in the thermal infrared (MODIS/Terra-Aqua) are used. We show that satellite data can provide fluxes with a sampling of a day or less (few hours in the best case). Generally the flux results from the different methods are consistent, and we discuss the advantages and weaknesses of each technique. Although the primary objective of this study is the calculation of SO₂ fluxes, it also enables us to assess the consistency of the SO₂ products from the different sensors used.

Sensitive and accurate detection of sulfur dioxide (SO₂) from space is important for monitoring and estimating global sulfur emissions. Inspired by detection methods applied in the thermal infrared, we present here a new scheme to retrieve SO₂ columns from satellite observations of ultraviolet back-scattered radiances. The retrieval is based on a measurement error covariance matrix to fully represent the SO₂-free radiance variability, so that the SO₂ slant column density is the only retrieved parameter of the algorithm. We demonstrate this approach, named COBRA, on measurements from the TROPOspheric Monitoring Instrument (TROPOMI) aboard the Sentinel-5 Precursor (S-5P) satellite. We show that the method reduces significantly both the noise and biases present in the current TROPOMI operational DOAS SO₂ retrievals. The performance of this technique is also benchmarked against that of the principal component algorithm (PCA) approach. We find that the quality of the data is similar and even slightly better with the proposed COBRA approach. The ability of the algorithm to retrieve SO₂ accurately is further supported by comparison with ground-based observations. We illustrate the great sensitivity of the method with a high-resolution global SO₂ map, considering 2.5 years of TROPOMI data. In addition to the known sources, we detect many new SO₂ emission hotspots worldwide. For the largest sources, we use the COBRA data to estimate SO₂ emission rates. Results are comparable to other recently published TROPOMI-based SO₂ emissions estimates, but the associated uncertainties are significantly lower than with the operational data. Next, for a limited number of weak sources, we demonstrate the potential of our data for quantifying SO₂ emissions with a detection limit of about 8 kt yr⁻¹, a factor of 4 better than the emissions derived from the Ozone Monitoring Instrument (OMI). We anticipate that the systematic use of our TROPOMI COBRA SO₂ column data set at a global scale will allow missing sources to be identified and quantified and help improve SO₂ emission inventories.

Initial objectives and design of the Benchmark campaign organized within the European COST Action ES1206 (2013–2017) are described in the paper. This campaign has aimed to support the development and validation of advanced Global Navigation Satellite System (GNSS) tropospheric products, in particular high-resolution and ultra-fast zenith total delays (ZTDs) and tropospheric gradients derived from a dense permanent network. A complex data set was collected for the 8-week period when several extreme heavy precipitation episodes occurred in central Europe which caused severe river floods in this area. An initial processing of data sets from GNSS products and numerical weather models (NWMs) provided independently estimated reference parameters – zenith tropospheric delays and tropospheric horizontal gradients. Their provision gave an overview about the product similarities and complementarities, and thus a potential for improvements of a synergy in their optimal exploitations in future. Reference GNSS and NWM results were intercompared and visually analysed using animated maps. ZTDs from two reference GNSS solutions compared to global ERA-Interim reanalysis resulted in accuracy at the 10 mm level in terms of the root mean square (rms) with a negligible overall bias, comparisons to Global Forecast System (GFS) forecasts showed accuracy at the 12 mm level with the overall bias of −5 mm and, finally, comparisons to mesoscale ALADIN-CZ forecast resulted in accuracy at the 8 mm level with a negligible total bias. The comparison of horizontal tropospheric gradients from GNSS and NWM data demonstrated a very good agreement among independent solutions with negligible biases and an accuracy of about 0.5 mm. Visual comparisons of maps of zenith wet delays and tropospheric horizontal gradients showed very promising results for future exploitations of advanced GNSS tropospheric products in meteorological applications, such as severe weather event monitoring and weather nowcasting. The GNSS products revealed a capability of providing more detailed structures in atmosphere than the state-of-the-art numerical weather models are able to capture. In an initial study on the contribution of hydrometeors (e.g. cloud water, ice or snow) to GNSS signal delays during severe weather, the effect reached up to 17 mm, and it was suggested that hydrometeors should be carefully accounted for within the functional model. The reference products will be further exploited in various specific studies using the Benchmark data set. It is thus going to play a key role in these highly interdisciplinary developments towards better mutual benefits from advanced GNSS and meteorological products.

An extensive validation of line-of-sight tropospheric slant total delays (STD) from Global Navigation Satellite Systems (GNSS), ray tracing in numerical weather prediction model (NWM) fields and microwave water vapour radiometer (WVR) is presented. Ten GNSS reference stations, including collocated sites, and almost 2 months of data from 2013, including severe weather events were used for comparison. Seven institutions delivered their STDs based on GNSS observations processed using 5 software programs and 11 strategies enabling to compare rather different solutions and to assess the impact of several aspects of the processing strategy. STDs from NWM ray tracing came from three institutions using three different NWMs and ray-tracing software. Inter-techniques evaluations demonstrated a good mutual agreement of various GNSS STD solutions compared to NWM and WVR STDs. The mean bias among GNSS solutions not considering post-fit residuals in STDs was −0.6 mm for STDs scaled in the zenith direction and the mean standard deviation was 3.7 mm. Standard deviations of comparisons between GNSS and NWM ray-tracing solutions were typically 10 mm ± 2 mm (scaled in the zenith direction), depending on the NWM model and the GNSS station. Comparing GNSS versus WVR STDs reached standard deviations of 12 mm ± 2 mm also scaled in the zenith direction. Impacts of raw GNSS post-fit residuals and cleaned residuals on optimal reconstructing of GNSS STDs were evaluated at inter-technique comparison and for GNSS at collocated sites. The use of raw post-fit residuals is not generally recommended as they might contain strong systematic effects, as demonstrated in the case of station LDB0. Simplified STDs reconstructed only from estimated GNSS tropospheric parameters, i.e. without applying post-fit residuals, performed the best in all the comparisons; however, it obviously missed part of tropospheric signals due to non-linear temporal and spatial variations in the troposphere. Although the post-fit residuals cleaned of visible systematic errors generally showed a slightly worse performance, they contained significant tropospheric signal on top of the simplified model. They are thus recommended for the reconstruction of STDs, particularly during high variability in the troposphere. Cleaned residuals also showed a stable performance during ordinary days while containing promising information about the troposphere at low-elevation angles.

The presence of volcanic clouds in the atmosphere affects air quality, the environment, climate, human health and aviation safety. The importance of the detection and retrieval of volcanic SO₂ lies with risk mitigation as well as with the possibility of providing insights into the mechanisms that cause eruptions. Due to their intrinsic characteristics, satellite measurements have become an essential tool for volcanic monitoring. In recent years, several sensors, with different spectral, spatial and temporal resolutions, have been launched into orbit, significantly increasing the effectiveness of the estimation of the various parameters related to the state of volcanic activity. In this work, the SO₂ total masses and fluxes were obtained from several satellite sounders—the geostationary (GEO) MSG-SEVIRI and the polar (LEO) Aqua/Terra-MODIS, NPP/NOAA20-VIIRS, Sentinel5p-TROPOMI, MetopA/MetopB-IASI and Aqua-AIRS—and compared to one another. As a test case, the Christmas 2018 Etna eruption was considered. The characteristics of the eruption (tropospheric with low ash content), the large amount of (simultaneously) available data and the different instrument types and SO₂ columnar abundance retrieval strategies make this cross-comparison particularly relevant. Results show the higher sensitivity of TROPOMI and IASI and a general good agreement between the SO₂ total masses and fluxes obtained from all the satellite instruments. The differences found are either related to inherent instrumental sensitivity or the assumed and/or calculated SO₂ cloud height considered as input for the satellite retrievals. Results indicate also that, despite their low revisit time, the LEO sensors are able to provide information on SO₂ flux over large time intervals. Finally, a complete error assessment on SO₂ flux retrievals using SEVIRI data was realized by considering uncertainties in wind speed and SO₂ abundance.

We present a multi-sensor archive collecting spatial and temporal information about volcanic SO₂ clouds generated by the 11 largest eruptions of this century. The detection and monitoring of volcanic clouds are an important topic for aviation management, climate issues and weather forecasts. Several studies focusing on single eruptive events exist, but no archive available at the moment combines quantitative data from as many instruments. We archived and collocated the SO₂ vertical column density estimations from three different satellite instruments (AIRS, IASI and GOME-2), atmospheric parameters as vertical profiles from the Global Navigation Satellite Systems (GNSS) Radio Occultations (RO), and the cloud-top height and aerosol type from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP). Additionally, we provide information about the cloud-top height from three different algorithms and the atmospheric anomaly due to the presence of the cloud. The dataset is gathering 206 d of SO₂ data, collocated with 44 180 backscatter profiles and 64 764 radio occultation profiles. The modular structure of the archive allows an easy collocation of the datasets according to the users' needs, and the cross-comparison of the datasets shows different consistency of the parameters estimated with different sensors and algorithms, according to the sensitivity and resolution of the instruments. The data described here are published with a DOI at https://doi.org/10.5880/fidgeo.2020.016 (Tournigand et al., 2020a).