Isabelle Taylor

Abstract:

Volcanic ash can interact with the earth system on many temporal and spatial scales and is a significant hazard to aircraft. In the event of a volcanic eruption, fast and robust decisions need to be made by aviation authorities about which routes are safe to operate. Such decisions take into account forecasts of ash location issued by Volcanic Ash Advisory Centers (VAACs) which are informed by simulations from Volcanic Ash Transport and Dispersion (VATD) models. The estimation of the time-evolving vertical distribution of ash emissions for use in VATD simulations in real time is difficult which can lead to large uncertainty in these forecasts. This study presents a method for constraining the ash emission estimates by combining an inversion modeling technique with an ensemble of meteorological forecasts, resulting in an ensemble of ash emission estimates. These estimates of ash emissions can be used to produce a robust ash forecast consistent with observations. This new ensemble approach is applied to the 2011 eruption of the Icelandic volcano Grímsvötn. The resulting emission profiles each have a similar temporal evolution but there are differences in the magnitude of ash emitted at different heights. For this eruption, the impact of precipitation uncertainty (and the associated wet deposition of ash) on the estimate of the total amount of ash emitted is larger than the impact of the uncertainty in the wind fields. Despite the differences that are dominated by wet deposition uncertainty, the ensemble inversion provides confidence that the reduction of the unconstrained emissions (a priori), particularly above 4 km, is robust across all members. In this case, the use of posterior emission profiles greatly reduces the magnitude and extent of the forecast ash cloud. The ensemble of posterior emission profiles gives a range of ash column loadings much closer in agreement with a set of independent satellite retrievals in comparison to the a priori emissions. Furthermore, airspace containing volcanic ash concentrations deemed to be associated with the highest risk (likelihood of exceeding a high concentration threshold) to aviation are reduced by over 85%. Such improvements could have large implications in emergency response situations. Future research will focus on quantifying the impact of uncertainty in precipitation forecasts on wet deposition in other eruptions and developing an inversion system that makes use of the state-of-the-art meteorological ensembles which has the potential to be used in an operational setting.

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.