Earth Observation Data Group

Abstract:

Ash clouds are a geographically far-reaching hazard associated with volcanic eruptions. To minimise the risk that these pose to aircraft and to limit disruption to the aviation industry, it is important to closely monitor the emission and atmospheric dispersion of these plumes. The altitude of the plume is an important consideration and is an essential input into many models of ash cloud propagation. CO2 slicing is an established technique for obtaining the top height of aqueous clouds, and previous studies have demonstrated that there is potential for this method to be used for volcanic ash. In this study, the CO2 slicing technique has been adapted for volcanic ash and applied to spectra obtained from the Infrared Atmospheric Sounding Interferometer (IASI). Simulated ash spectra are first used to select the most appropriate channels and then demonstrate that the technique has merit for determining the altitude of the ash. These results indicate a strong match between the true heights and CO2 slicing output with a root mean square error (RMSE) of less than 800 m. Following this, the technique was applied to spectra obtained with IASI during the Eyjafjallajökull and Grímsvötn eruptions in 2010 and 2011 respectively, both of which emitted ash clouds into the troposphere, and which have been extensively studied with satellite imagery. The CO2 slicing results were compared against those from an optimal estimation scheme, also developed for IASI, and a satellite-borne lidar is used for validation. The CO2 slicing heights returned an RMSE value of 2.2 km when compared against the lidar. This is lower than the RMSE for the optimal estimation scheme (2.8 km). The CO2 slicing technique is a relatively fast tool and the results suggest that this method could be used to get a first approximation of the ash cloud height, potentially for use for hazard mitigation, or as an input for other retrieval techniques or models of ash cloud propagation.

Abstract:

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.

Abstract:

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.

Abstract:

The most widely used satellite products are averages of data onto a regular spatiotemporal grid, known as Level 3 data. Some atmospheric variables can vary rapidly in response to changing conditions. Over the scales of Level 3 averaging, the combination of observations across different conditions may result in data that is not normally distributed, such that a simple mean is not representative. The problem is illustrated by the distribution of aerosol optical depth from different sensors and algorithms. A simple statistical technique is proposed to better convey the diversity of satellite observations to users whereby a multimodal log-normal distribution is fit to the distribution of data observed within each grid cell. Allowing multiple modes within each cell is shown to improve the agreement between satellite products by highlighting regions of significant variability and isolating systematic differences between instruments.

Abstract:

Radiative transfer models used in remote sensing and hazard assessment of volcanic ash require knowledge of ash optical parameters. Here, we characterise the bulk and glass compositions of a representative suite of volcanic ash samples with known complex refractive indices (n + ik: where n is the real and k is the imaginary part). Using a linear regression model, we develop a new parameterization allowing the complex refractive index of volcanic ash to be estimated from ash SiO2 content or ratio of non-bridging oxygens to tetrahedrally-coordinated cations (NBO/T). At visible wavelengths, n correlates better with bulk than glass composition (both SiO2 and NBO/T), and k correlates better with SiO2 content than NBO/T. Over a broader spectral range (0.4–19 μm), bulk correlates better than glass composition, and NBO/T generally correlates better than SiO2 content for both parts of the refractive index. In order to understand the impacts of our new parameterization on satellite retrievals, we compared IASI satellite (wavelengths 3.62–15.5 μm) mass loading retrievals using our new approach with retrievals that assumed a generic (Eyjafjallajökull) ash refractive index. There are significant differences in mass loading using our calculated indices specific to ash type rather than a generic index. Where mass loadings increase, there is often improvement in retrieval quality (corresponding to cost function decrease). This new parameterization of refractive index variation with ash composition will help to improve remote sensing retrievals for the rapid identification of ash and quantitative analysis of mass loadings from satellite data on operational timescales.