Stefano Corradini

Mt. Etna tropospheric ash retrieval and sensitivity analysis using moderate resolution imaging spectroradiometer measurements

A retrieval of tropospheric volcanic ash from Mt Etna has been carried out, using measurements from the Moderate Resolution Imaging Spectroradiometer (MODIS). The NASA-MODIS satellite instrument acquires images in the 0.4 to 14 μm spectral range with a spatial resolution of 1 km at nadir. The eruption which occurred on 24 November 2006 is considered as a test case in this work. In order to derive the ash plume optical thickness, the particle effective radius and the total mass, the Brightness Temperature Difference procedure has been applied to MODIS channels 31 (centered at 11 μm) and 32 (centered at 12 μm). Channel 5 (centered at 1.24 μm) has been used to refine the cloud discrimination, exploiting the distinct reflectivity of meteorological and volcanic clouds in the near infrared spectral range. The detection of volcanic ash pixels has been significantly improved by applying an atmospheric water vapor correction to MODIS data. This procedure doubles the number of pixels identified as containing volcanic ash compared to the original method. The retrieved mean ash optical thickness at 0.55 μm, mean particle effective radius and the total ash mass in the plume are 0.4, 3.5 μm and 3620 tons, respectively. A detailed sensitivity analysis has been carried out to investigate errors in the retrieval caused by the uncertainty in various parameters: surface temperature and emissivity, plume geometry (altitude and thickness), ash type and atmospheric water vapor. Results show that the largest contributions to retrieval errors are from uncertainty in surface parameters, aerosol type and atmospheric water vapor. The total tropospheric volcanic ash retrieval errors are estimated to be 30%, 30% and 40% for mean AOT, mean effective radius and total mass retrieval, respectively.

A Multi-Sensor Approach for Volcanic Ash Cloud Retrieval and Eruption Characterization: The 23 November 2013 Etna Lava Fountain

Volcanic activity is observed worldwide with a variety of ground and space-based remote sensing instruments, each with advantages and drawbacks. No single system can give a comprehensive description of eruptive activity, and so, a multi-sensor approach is required. This work integrates infrared and microwave volcanic ash retrievals obtained from the geostationary Meteosat Second Generation (MSG)-Spinning Enhanced Visible and Infrared Imager (SEVIRI), the polar-orbiting Aqua-MODIS and ground-based weather radar. The expected outcomes are improvements in satellite volcanic ash cloud retrieval (altitude, mass, aerosol optical depth and effective radius), the generation of new satellite products (ash concentration and particle number density in the thermal infrared) and better characterization of volcanic eruptions (plume altitude, total ash mass erupted and particle number density from thermal infrared to microwave). This approach is the core of the multi-platform volcanic ash cloud estimation procedure being developed within the European FP7-APhoRISM project. The Mt. Etna (Sicily, Italy) volcano lava fountaining event of 23 November 2013 was considered as a test case. The results of the integration show the presence of two volcanic cloud layers at different altitudes. The improvement of the volcanic ash cloud altitude leads to a mean difference between the SEVIRI ash mass estimations, before and after the integration, of about the 30%. Moreover, the percentage of the airborne “fine” ash retrieved from the satellite is estimated to be about 1%–2% of the total ash emitted during the eruption. Finally, all of the estimated parameters (volcanic ash cloud altitude, thickness and total mass) were also validated with ground-based visible camera measurements, HYSPLIT forward trajectories, Infrared Atmospheric Sounding Interferometer (IASI) satellite data and tephra deposits.

Volcanic Ash Cloud Properties: Comparison Between MODIS Satellite Retrievals and FALL3D Transport Model

The moderate Resolution Imaging Spectroradiometer (MODIS) is a multispectral satellite instrument operating from the visible to thermal infrared spectral range. FALL3D is a 3-D time-dependent Eulerian model for the transport and deposition of volcanic particles. In this letter, quantitative comparison between the volcanic cloud ash mass and optical depth retrieved by MODIS and modeled by FALL3D has been performed. Three MODIS images collected on October 28, 29, and 30 on Mt. Etna volcano during the 2002 eruption have been considered as test cases. The results show a general good agreement between the retrieved and the modeled volcanic clouds in the first 300 km from the vents. Even if the modeled volcanic cloud area is systematically wider than the retrieved area, the ash total mass is comparable and varies between 35 and 60 kt and between 20 and 42 kt for FALL3D and MODIS, respectively. The mean aerosol optical depth (AOD) values are in good agreement and approximately equal to 0.8. When the whole volcanic clouds are considered the ash areas, then the total ash masses, computed by FALL3D model, are significantly greater than the same parameters retrieved from the MODIS data, while the mean AOD values remain in very good agreement and equal to about 0.6. The volcanic cloud direction in its distal part is not coincident for the October 29 and 30, 2002 images due to the difference between the real and the modeled local wind fields. Finally, the MODIS maps show regions of high mass and AOD due to volcanic puffs not modeled by FALL3D.

Multisensor Satellite Monitoring of the 2011 Puyehue-Cordon Caulle Eruption

This paper shows the main outcomes of the Puyehue volcano (Chile) eruption monitoring by means of multisensor remote sensing instruments working from thermal infrared (TIR) to microwave (MW) spectral range. Thanks to the use of Synthetic Aperture Radar (SAR) and the Moderate Resolution Imaging Spectroradiometer (MODIS), the eruption evolution was observed, capturing the deformations of volcano edifice, the lava extension, as well as the information on ash and gas emitted. On the one hand, SAR Interferometry applied to ENVISAT-ASAR data allowed the estimation of the deformation occurred just before the beginning of the eruption and the subsequent deflation, with monthly sampling. On the other hand, with the combined use of the very high-resolution (VHR) images taken by COSMO-SkyMed X-band SAR, and ENVISAT-ASAR ones, we were able to follow the lava deposition during the most intense phase of the eruption. Additionally, the joined exploitation of SAR and optical MODIS images allowed ash detection, also in cloudy sky conditions. Finally, the information gathered by both types of sensors allowed to highlight some volcanological features of the eruption and the relationship between surface deformation and the amount of ash and gases emitted by the volcano.

Empirical correction of multifilter rotating shadowband radiometer (MFRSR) aerosol optical depths for the aerosol forward scattering and development of a long-term integrated MFRSR-Cimel dataset at Lampedusa.

Aerosol optical properties have been measured on the island of Lampedusa (35.5°N, 12.6°E) with seven-band multifilter rotating shadowband radiometers (MFRSRs) and a CE 318 Cimel sunphotometer (part of the AERONET network) since 1999. Four different MFRSRs have operated since 1999. The Cimel sunphotometer has been operational for a short period in 2000 and in 2003-2006 and 2010-present. Simultaneous determinations of the aerosol optical depth (AOD) from the two instruments were compared over a period of almost 4 years at several wavelengths between 415 and 870 nm. This is the first long-term comparison at a site strongly influenced by desert dust and marine aerosols and characterized by frequent cases of elevated AOD. The datasets show a good agreement, with MFRSR underestimating the Cimel AOD in cases with low Ångström exponent; the underestimate decreases for increasing wavelength and increases with AOD. This underestimate is attributed to the effect of aerosol forward scattering on the relatively wide field of view of the MFRSR. An empirical correction of the MFRSR data was implemented. After correction, the mean bias (MB) between MFRSR and Cimel simultaneous AOD determinations is always smaller than 0.004, and the root mean square difference is ≤0.031 at all wavelengths. The MB between MFRSR and Cimel monthly averages (for months with at least 20 days with AOD determinations) is 0.0052. Thus, by combining the MFRSR and Cimel observations, an integrated long-term series is obtained, covering the period 1999-present, with almost continuous measurements since early 2002. The long-term data show a small (nonstatistically significant) decreasing trend over the period 2002-2013, in agreement with independent observations in the Mediterranean. The integrated Lampedusa dataset will be used for aerosol climatological studies and for verification of satellite observations and model analyses.

Inspecting MIVIS capability to retrieve chemical-mineralogical information: evaluation and analysis of VNIR-SWIR data acquired on a volcanic area

This work is a contribution to the assessment of MIVIS (Multi-spectral Infrared and Visible Imaging Spectrometer) airborne imaging spectrometer capability in applications of surface characterization. The focus is on the visible and near-infrared–short wave infrared (VNIR–SWIR) spectral region, using a dataset acquired in 1994 on Vulcano Island (Italy), to retrieve chemical–mineralogical information on the altered deposits related to volcanic activity. The main processing steps include data quality evaluation in terms of signal-to-noise ratio, atmospheric and topographic corrections and spectral interpretation of the image. Estimation of surface reflectance is based on atmospheric modelling by MODTRAN3.5 and 6S radiative transfer codes. Representative MIVIS reflectance spectra of the main surface units are compared with spectra measured in the laboratory on field samples, and interpreted to characterize the mineralogy on the basis of their spectral features. A thematic map of the main alteration units is then produced by applying spectral mapping techniques to the surface reflectance image, using a set of channels selected on the basis of their data quality and image-derived end-member spectra.

Stereoscopic Estimation of Volcanic Ash Cloud-Top Height from Two Geostationary Satellites

The characterization of volcanic ash clouds released into the atmosphere during explosive eruptions includes cloud height as a fundamental physical parameter. A novel application is proposed of a method based on parallax data acquired from two geostationary instruments for estimating ash cloud-top height (ACTH). An improved version of the method with a detailed discussion of height retrieval accuracy was applied to estimate ACTH from two datasets acquired by two satellites in favorable positions to fully exploit the parallax effect. A combination of MSG SEVIRI (HRV band; 1000 m nadir spatial resolution, 5 min temporal resolution) and Meteosat-7 MVIRI (VIS band, 2500 m nadir spatial resolution, 30 min temporal resolution) was implemented. Since MVIRI does not acquire data at exactly the same time as SEVIRI, a correction procedure enables compensation for wind advection in the atmosphere. The method was applied to the Mt. Etna, Sicily, Italy, eruption of 23 November 2013. The height of the volcanic cloud was tracked with a top height of ~8.5 km. The ash cloud estimate was applied to the visible channels to show the potential accuracy that will soon be achievable also in the infrared range using the next generation of multispectral imagers. The new constellation of geostationary meteorological satellites will enable full exploitation of this technique for continuous global ACTH monitoring.

Volcanic ash cloud detection from space: a comparison between the RSTASH technique and the water vapour corrected BTD procedure

Volcanic eruptions can inject large amounts (Tg) of gas and particles into the troposphere and, sometimes, into the stratosphere. Besides the main gases (H2O, CO2, SO2 and HCl), volcanic clouds contain a mix of silicate ash particles in the size range from 0.1 μm to 1 mm or larger. The interest in volcanic ash detection is high, particularly because it represents a serious hazard for air traffic. Particles with dimensions of several millimetres can damage the aircraft structure (windows, wings, ailerons), while particles less than 10 μm may be extremely dangerous for the jet engines and are undetectable by the pilots during night or in low visibility conditions. Furthermore, ash detection represents a critical step towards quantitative retrievals of plume parameters. In this paper two different satellite techniques for volcanic cloud detection and tracking are compared, namely a water vapour corrected version of the brightness temperature difference (BTD-WVC) procedure and an implementation of the robust satellite technique, specifically configured for volcanic ash (RSTASH). The BTD method identifies volcanic ash clouds on the basis of the brightness temperature difference measured in two infrared spectral bands at around 11 and 12 μm. To account for the atmospheric water vapour differential absorption in the 11–12 μm spectral range, which tends to reduce (and in some cases completely mask) the BTD signal, a water vapour correction procedure has been developed (BTD-WVC), based on measured or synthetic atmospheric profiles. RSTASH instead, is based on the analysis of a time series of satellite records, aimed at identifying signal anomalies through an automatic unsupervised change detection step. To assess the performance of the BTD-WVC and RSTASH methods in detecting volcanic ash clouds, some eruptive events of Mt Etna, observed by the Advanced Very High Resolution Radiometer (AVHRR) sensor, have been analysed. The obtained results show a good agreement between the BTD-WVC and RSTASH techniques for all the considered images, in terms of pixels detected as ‘ash affected’ (i.e. the ash cloud area). In particular, compared to the traditional BTD procedure, the BTD-WVC and RSTASH techniques significantly improve volcanic ash cloud detection, both in daytime and night-time data, especially in the case of low ash loading.

Theoretical study on SO2 and ash volcanic plume retrievals using ground TIR camera. Sensitivity analysis and retrieval procedure developments

In this paper, a sensitivity analysis and procedure development for volcanic-plume sulfur dioxide and ash retrievals using ground thermal infrared camera have been carried out. The semiconductor device camera, considered as a reference, has a spectral range of 8-14 ¿m with noise equivalent temperature difference that is better than 100 mK at 300 K. The camera will be used to monitor and assess the hazards of Mt. Etna volcano to mitigate the risk and impact of volcanic eruptions on the civil society and transports. A minimum number of filters have been selected for sulfur dioxide (SO2) and volcanic ash retrievals. The sensitivity study has been carried out to determine the SO2 and volcanic ash minimum concentration detectable by the system varying the camera geometry and the atmospheric profiles. Results show a meaningful sensitivity increase considering high instrument altitudes and low camera-elevation angles. For all geometry configurations and monthly profiles, the sensitivity limit varies between 0.5 and 2 g ·m-2 for SO2 columnar abundance and between 0.02 and 1 for ash optical depth. Two procedures to detect SO2 and ash, based on the least square fit method and on the brightness temperature difference (BTD) algorithm, respectively, have also been proposed. Results show that high concentration of atmospheric water vapor columnar content significantly reduces the ash-plume effect on the BTD. A water vapor-correction procedure introduced improves the ash retrievals and the cloud discrimination in every season, considering all the camera geometries.

Theoretical Study on Volcanic Plume

In this paper, a sensitivity analysis and procedure development for volcanic-plume sulfur dioxide and ash retrievals using ground thermal infrared camera have been carried out. The semiconductor device camera, considered as a reference, has a spectral range of 8-14 ¿m with noise equivalent temperature difference that is better than 100 mK at 300 K. The camera will be used to monitor and assess the hazards of Mt. Etna volcano to mitigate the risk and impact of volcanic eruptions on the civil society and transports. A minimum number of filters have been selected for sulfur dioxide (SO2) and volcanic ash retrievals. The sensitivity study has been carried out to determine the SO2 and volcanic ash minimum concentration detectable by the system varying the camera geometry and the atmospheric profiles. Results show a meaningful sensitivity increase considering high instrument altitudes and low camera-elevation angles. For all geometry configurations and monthly profiles, the sensitivity limit varies between 0.5 and 2 g ·m-2 for SO2 columnar abundance and between 0.02 and 1 for ash optical depth. Two procedures to detect SO2 and ash, based on the least square fit method and on the brightness temperature difference (BTD) algorithm, respectively, have also been proposed. Results show that high concentration of atmospheric water vapor columnar content significantly reduces the ash-plume effect on the BTD. A water vapor-correction procedure introduced improves the ash retrievals and the cloud discrimination in every season, considering all the camera geometries.