ALBERT

Description: We present version 3 (V3) of the Cloud_cci Along-Track Scanning Radiometer (ATSR) and Advanced ATSR (AATSR) data set. The data set was created for the European Space Agency (ESA) Cloud_cci (Climate Change Initiative) programme. The cloud properties were retrieved from the second ATSR (ATSR-2) on board the second European Remote Sensing Satellite (ERS-2) spanning 1995–2003 and the AATSR on board Envisat, which spanned 2002–2012. The data are comprised of a comprehensive set of cloud properties: cloud top height, temperature, pressure, spectral albedo, cloud effective emissivity, effective radius, and optical thickness, alongside derived liquid and ice water path. Each retrieval is provided with its associated uncertainty. The cloud property retrievals are accompanied by high-resolution top- and bottom-of-atmosphere shortwave and longwave fluxes that have been derived from the retrieved cloud properties using a radiative transfer model. The fluxes were generated for all-sky and clear-sky conditions. V3 differs from the previous version 2 (V2) through development of the retrieval algorithm and attention to the consistency between the ATSR-2 and AATSR instruments. The cloud properties show improved accuracy in validation and better consistency between the two instruments, as demonstrated by a comparison of cloud mask and cloud height with co-located CALIPSO data. The cloud masking has improved significantly, particularly in its ability to detect clear pixels. The Kuiper Skill score has increased from 0.49 to 0.66. The cloud top height accuracy is relatively unchanged. The AATSR liquid water path was compared with the Multisensor Advanced Climatology of Liquid Water Path (MAC-LWP) in regions of stratocumulus cloud and shown to have very good agreement and improved consistency between ATSR-2 and AATSR instruments. The correlation with MAC-LWP increased from 0.4 to over 0.8 for these cloud regions. The flux products are compared with NASA Clouds and the Earth's Radiant Energy System (CERES) data, showing good agreement within the uncertainty. The new data set is well suited to a wide range of climate applications, such as comparison with climate models, investigation of trends in cloud properties, understanding aerosol–cloud interactions, and providing contextual information for co-located ATSR-2/AATSR surface temperature and aerosol products. The following new digital identifier has been issued for the Cloud_cci ATSR-2/AATSRv3 data set: https://doi.org/10.5676/DWD/ESA_Cloud_cci/ATSR2-AATSR/V003 (Poulsen et al., 2019).

Description: We present an optimal-estimation (OE) retrieval scheme for stratospheric sulfur dioxide from the High-Resolution Infrared Radiation Sounder 2 (HIRS/2) instruments on the NOAA and MetOp platforms, an infrared radiometer that has been operational since 1979. This algorithm is an improvement upon a previous method based on channel brightness temperature differences, which demonstrated the potential for monitoring volcanic SO2 using HIRS/2. The Prata method is fast but of limited accuracy. This algorithm uses an optimal-estimation retrieval approach yielding increased accuracy for only moderate computational cost. This is principally achieved by fitting the column water vapour and accounting for its interference in the retrieval of SO2. A cloud and aerosol model is used to evaluate the sensitivity of the scheme to the presence of ash and water/ice cloud. This identifies that cloud or ash above 6 km limits the accuracy of the water vapour fit, increasing the error in the SO2 estimate. Cloud top height is also retrieved. The scheme is applied to a case study event, the 1991 eruption of Cerro Hudson in Chile. The total erupted mass of SO2 is estimated to be 2300 kT ± 600 kT. This confirms it as one of the largest events since the 1991 eruption of Pinatubo, and of comparable scale to the Northern Hemisphere eruption of Kasatochi in 2008. This retrieval method yields a minimum mass per unit area detection limit of 3 DU, which is slightly less than that for the Total Ozone Mapping Spectrometer (TOMS), the only other instrument capable of monitoring SO2 from 1979 to 1996. We show an initial comparison to TOMS for part of this eruption, with broadly consistent results. Operating in the infrared (IR), HIRS has the advantage of being able to measure both during the day and at night, and there have frequently been multiple HIRS instruments operated simultaneously for better than daily sampling. If applied to all data from the series of past and future HIRS instruments, this method presents the opportunity to produce a comprehensive and consistent volcanic SO2 time series spanning over 40 years.

Description: The 6-month-long 2014–2015 Holuhraun eruption was the largest in Iceland for 200 years, emitting huge quantities of sulfur dioxide (SO2) into the troposphere, at times overwhelming European anthropogenic emissions. Weather, terrain and latitude made continuous ground-based or UV satellite sensor measurements challenging. Infrared Atmospheric Sounding Interferometer (IASI) data are used to derive the first time series of daily SO2 mass present in the atmosphere and its vertical distribution over the entire eruption period. A new optimal estimation scheme is used to calculate daily SO2 fluxes and average e-folding time every 12 h. For the 6 months studied, the SO2 flux was observed to be up to 200 kt day−1 and the minimum total SO2 erupted mass was 4.4±0.8 Tg. The average SO2 e-folding time was 2.4±0.6 days. Where comparisons are possible, these results broadly agree with ground-based near-source measurements, independent remote-sensing data and values obtained from model simulations from a previous paper. The results highlight the importance of using high-resolution time series data to accurately estimate volcanic SO2 emissions. The SO2 mass missed due to thermal contrast is estimated to be of the order of 3 % of the total emission when compared to measurements by the Ozone Monitoring Instrument. A statistical correction for cloud based on the AVHRR cloud-CCI data set suggested that the SO2 mass missed due to cloud cover could be significant, up to a factor of 2 for the plume within the first kilometre from the vent. Applying this correction results in a total erupted mass of 6.7±0.4 Tg and little change in average e-folding time. The data set derived can be used for comparisons to other ground- and satellite-based measurements and to petrological estimates of the SO2 flux. It could also be used to initialise climate model simulations, helping to better quantify the environmental and climatic impacts of future Icelandic fissure eruptions and simulations of past large-scale flood lava eruptions.

Description: Increased concentrations of aerosol can enhance the albedo of warm low-level cloud. Accurately quantifying this relationship from space is challenging due in part to contamination of aerosol statistics near clouds. Aerosol retrievals near clouds can be influenced by stray cloud particles in areas assumed to be cloud-free, particle swelling by humidification, shadows and enhanced scattering into the aerosol field from (3-D radiative transfer) clouds. To screen for this contamination we have developed a new cloud–aerosol pairing algorithm (CAPA) to link cloud observations to the nearest aerosol retrieval within the satellite image. The distance between each aerosol retrieval and nearest cloud is also computed in CAPA. Results from two independent satellite imagers, the Advanced Along-Track Scanning Radiometer (AATSR) and Moderate Resolution Imaging Spectroradiometer (MODIS), show a marked reduction in the strength of the intrinsic aerosol indirect radiative forcing when selecting aerosol pairs that are located farther away from the clouds (−0.28±0.26 W m−2) compared to those including pairs that are within 15 km of the nearest cloud (−0.49±0.18 W m−2). The larger aerosol optical depths in closer proximity to cloud artificially enhance the relationship between aerosol-loading, cloud albedo, and cloud fraction. These results suggest that previous satellite-based radiative forcing estimates represented in key climate reports may be exaggerated due to the inclusion of retrieval artefacts in the aerosol located near clouds.

Description: The Community Cloud retrieval for Climate (CC4CL) is a cloud property retrieval system for satellite-based multispectral imagers and is an important component of the Cloud Climate Change Initiative (Cloud_cci) project. In this paper we discuss the optimal estimation retrieval of cloud optical thickness, effective radius and cloud top pressure based on the Optimal Retrieval of Aerosol and Cloud (ORAC) algorithm. Key to this method is the forward model, which includes the clear-sky model, the liquid water and ice cloud models, the surface model including a bidirectional reflectance distribution function (BRDF), and the "fast" radiative transfer solution (which includes a multiple scattering treatment). All of these components and their assumptions and limitations will be discussed in detail. The forward model provides the accuracy appropriate for our retrieval method. The errors are comparable to the instrument noise for cloud optical thicknesses greater than 10. At optical thicknesses less than 10 modeling errors become more significant. The retrieval method is then presented describing optimal estimation in general, the nonlinear inversion method employed, measurement and a priori inputs, the propagation of input uncertainties and the calculation of subsidiary quantities that are derived from the retrieval results. An evaluation of the retrieval was performed using measurements simulated with noise levels appropriate for the MODIS instrument. Results show errors less than 10 % for cloud optical thicknesses greater than 10. Results for clouds of optical thicknesses less than 10 have errors up to 20 %.

Description: This study focuses on two new aspects of inverse modelling of volcanic emissions. First, we derive an observation operator for satellite retrievals of plume height, and second, we solve the inverse problem using an algorithm based on the 4D-Var data assimilation method. The approach is first tested in a twin experiment with simulated observations and further evaluated by assimilating IASI SO2 plume height and total column retrievals in a source term inversion for the 2010 eruption of Eyjafjallajökull. The inversion resulted in temporal and vertical reconstruction of the SO2 emissions during 1–20 May 2010 with formal vertical and temporal resolutions of 500 m and 12 h.The plume height observation operator is based on simultaneous assimilation of the plume height and total column retrievals. The plume height is taken to represent the vertical centre of mass, which is transformed into the first moment of mass (centre of mass times total mass). This makes the observation operator linear and simple to implement. The necessary modifications to the observation error covariance matrix are derived.Regularization by truncated iteration is investigated as a simple and efficient regularization method for the 4D-Var-based inversion. In the twin experiments, the truncated iteration was found to perform similarly to the commonly used Tikhonov regularization, which in turn is equivalent to a Gaussian a priori source term. However, the truncated iteration allows the level of regularization to be determined a posteriori without repeating the inversion.In the twin experiments, assimilating the plume height retrievals resulted in a 5–20 % improvement in root mean squared error but simultaneously introduced a 10–20 % low bias on the total emission depending on assumed emission profile. The results are consistent with those obtained with real data. For Eyjafjallajökull, comparisons with observations showed that assimilating the plume height retrievals reduced the overestimation of injection height during individual periods of 1–3 days, but for most of the simulated 20 days, the injection height was constrained by meteorological conditions, and assimilation of the plume height retrievals had only small impact. The a posteriori source term for Eyjafjallajökull consisted of 0.29 Tg (with total column and plume height retrievals) or 0.33 Tg (with total column retrievals only) erupted SO2 of which 95 % was injected below 11 or 12 km, respectively.

Description: New cloud property datasets based on measurements from the passive imaging satellite sensors AVHRR, MODIS, ATSR2, AATSR and MERIS are presented. Two retrieval systems were developed that include components for cloud detection and cloud typing followed by cloud property retrievals based on the optimal estimation (OE) technique. The OE-based retrievals are applied to simultaneously retrieve cloud-top pressure, cloud particle effective radius and cloud optical thickness using measurements at visible, near-infrared and thermal infrared wavelengths, which ensures spectral consistency. The retrieved cloud properties are further processed to derive cloud-top height, cloud-top temperature, cloud liquid water path, cloud ice water path and spectral cloud albedo. The Cloud_cci products are pixel-based retrievals, daily composites of those on a global equal-angle latitude–longitude grid, and monthly cloud properties such as averages, standard deviations and histograms, also on a global grid. All products include rigorous propagation of the retrieval and sampling uncertainties. Grouping the orbital properties of the sensor families, six datasets have been defined, which are named AVHRR-AM, AVHRR-PM, MODIS-Terra, MODIS-Aqua, ATSR2-AATSR and MERIS+AATSR, each comprising a specific subset of all available sensors. The individual characteristics of the datasets are presented together with a summary of the retrieval systems and measurement records on which the dataset generation were based. Example validation results are given, based on comparisons to well-established reference observations, which demonstrate the good quality of the data. In particular the ensured spectral consistency and the rigorous uncertainty propagation through all processing levels can be considered as new features of the Cloud_cci datasets compared to existing datasets. In addition, the consistency among the individual datasets allows for a potential combination of them as well as facilitates studies on the impact of temporal sampling and spatial resolution on cloud climatologies.For each dataset a digital object identifier has been issued:Cloud_cci AVHRR-AM: https://doi.org/10.5676/DWD/ESA_Cloud_cci/AVHRR-AM/V002Cloud_cci AVHRR-PM: https://doi.org/10.5676/DWD/ESA_Cloud_cci/AVHRR-PM/V002Cloud_cci MODIS-Terra: https://doi.org/10.5676/DWD/ESA_Cloud_cci/MODIS-Terra/V002Cloud_cci MODIS-Aqua: https://doi.org/10.5676/DWD/ESA_Cloud_cci/MODIS-Aqua/V002Cloud_cci ATSR2-AATSR: https://doi.org/10.5676/DWD/ESA_Cloud_cci/ATSR2-AATSR/V002Cloud_cci MERIS+AATSR: https://doi.org/10.5676/DWD/ESA_Cloud_cci/MERIS+AATSR/V002

Description: Increased concentrations of aerosol can enhance the albedo of warm lowlevel cloud. Accurately quantifying this relationship from space is challenging due in part to contamination of aerosol statistics near clouds. Aerosol retrievals near clouds can be influenced by stray cloud particles in areas assumed to be cloud-free, particle swelling by humidification, shadows and enhanced scattering into the aerosol field from (3D radiative transfer) clouds. To screen for this contamination, we have developed a new 5 Cloud-Aerosol Pairing Algorithm (CAPA) to link cloud observations to the nearest aerosol retrieval within the satellite image. The distance between each aerosol retrieval and nearest cloud is also computed in CAPA. Results from two independent satellite imagers, the Advanced Along Track Scanning Radiometer (AATSR) and MODerate Resolution Imaging Spectroradiometer (MODIS) show a marked reduction in the strength of the intrinsic aerosol indirect forcing when selecting aerosol pairs that are located farther away from the clouds (−0.28 ± 0.26 W/m2) compared to those 10 including pairs that are within 15 km of the nearest cloud (−0.49 ± 0.18 W/m2). The larger aerosol optical depths in closer proximity to cloud artificially enhance the relationship between aerosol loading, cloud albedo, and cloud fraction. These results suggest that previous satellite-based radiative forcing estimates represented in key climate reports may be exaggerated due to including retrieval artefacts in the aerosol located near clouds.

Description: We present an optimal estimation (OE) retrieval scheme for stratospheric sulphur dioxide from the High Resolution Infrared Radiation Sounder 2 (HIRS/2) instruments on the NOAA and MetOp platforms, an infrared radiometer that has been operational since 1979. This algorithm is an improvement upon a previous method based on channel brightness temperature differences developed by Prata et al. (2003), which demonstrated the potential for monitoring volcanic SO2 using HIRS/2. The Prata method is fast but of limited accuracy. This algorithm uses an optimal estimation retrieval approach yielding increased accuracy for only moderate computational cost. This is principally achieved by fitting the column water vapour and accounting for its interference in the retrieval of SO2. A cloud and aerosol model is used to evaluate the sensitivity of the scheme to the presence of ash and water/ice cloud. This identifies that cloud or ash above 6 km limits the accuracy of the water vapour fit, increasing the error in the SO2 estimate. Cloud top height is also retrieved. The scheme is applied to a case study event, the 1991 eruption of Cerro Hudson in Chile. The total erupted mass of SO2 is estimated to be 2300 kT ± 600 kT. This confirms it as one of the largest events since the 1991 eruption of Pinatubo, and of comparable scale to the Northern Hemisphere eruption of Kasatochi in 2008. This retrieval method yields a minimum mass per unit area detection limit of 3 DU, which is slightly less than that for the Total Ozone Mapping Spectrometer (TOMS), the only other instrument capable of monitoring SO2 from 1979–1996. We show an initial comparison to TOMS for part of this eruption, with broadly consistent results. Operating in the infrared (IR), HIRS has the advantage of being able to measure both during the day and at night, and there have frequently been multiple HIRS instruments operated simultaneously for better than daily sampling. If applied to all data from the series of past and future HIRS instruments, this method presents the opportunity to produce a comprehensive and consistent volcanic SO2 timeseries spanning over 40 years.

Description: The six-month-long 2014–2015 Holuhraun eruption was the largest in Iceland for 200 years, emitting huge quantities of sulphur dioxide (SO2) into the troposphere, at times overwhelming European anthropogenic emissions. Weather, terrain and latitude, made continuous ground-based or UV satellite sensor measurements challenging. Infrared Atmospheric Sounding Interferometer (IASI) data, is used to derive the first time-series of daily SO2 mass and vertical distribution over the eruption period. A new optimal estimation scheme is used to calculate daily SO2 fluxes and average e-folding time every twelve hours. The algorithm is used to estimate SO2 fluxes of up to 200 kt per day and a minimum total SO2 erupted mass of 4.4±0.8Tg. The average SO2 e-folding time was 2.4±0.6 days. Where comparisons are possible, these results broadly agree with ground-based near-source measurements, independent remote-sensing data and model simulations of the eruption. The results highlight the importance of high-resolution time-series data to accurately estimate volcanic SO2 emissions.