Automated Alert System for Volcanic Unrest

N. Anantrasirichai

Vi-Lab, University of Bristol

Juliet Biggs

School of Earth Sciences

University of Bristol

Fabien Albino

School of Earth Sciences

University of Bristol

Aim

To develop an effective and efficient automated global volcano alert system using Interferometric Synthetic Aperture Radar (InSAR) data, which can detect surface deformation, having a strong statistical link to eruption.

Funder
EPSRC Global Challenges Sponsorship, NERC Impact Acceleration Awards, NERC Innovation Award (NE/S013970/1)

Automated processing and CNN detection of deformation from 592,224 wrapped interferograms of 1,084 volcanoes [May 2021]. (left) Spatial distribution of detections (P>0.5). Size represents the number of detections and colour the percentage of detections. (right) Number of detections (P>0.5) per volcano (log-scale), coloured by maximum probability (Pmax).

Application of Machine Learning to Classification of Volcanic Deformation in Routinely Generated InSAR Data

Recent improvements in the frequency, type, and availability of satellite images mean it is now feasible to routinely study volcanoes in remote and inaccessible regions, including those with no ground-based monitoring. In particular, Interferometric Synthetic Aperture Radar data can detect surface deformation, which has a strong statistical link to eruption. However, the data set produced by the recently launched Sentinel-1 satellite is too large to be manually analyzed on a global basis.

In this study, we systematically process >30,000 short-term interferograms at over 900 volcanoes and apply machine learning algorithms to automatically detect volcanic ground deformation. We use a convolutional neutral network to classify interferometric fringes in wrapped interferograms with no atmospheric corrections. We employ a transfer learning strategy and test a range of pretrained networks, finding that AlexNet is best suited to this task. The positive results are checked by an expert and fed back for model updating. Following training with a combination of both positive and negative examples, this method reduced the number of interferograms to ∼100 which required further inspection, of which at least 39 are considered true positives. We demonstrate that machine learning can efficiently detect large, rapid deformation signals in wrapped interferograms, but further development is required to detect slow or small deformation patterns which do not generate multiple fringes in short duration interferograms. This study is the first to use machine learning approaches for detecting volcanic deformation in large data sets and demonstrates the potential of such techniques for developing alert systems based on satellite imagery.

A Deep Learning Approach to Detecting Volcano Deformation from Satellite Imagery using Synthetic Datasets

Satellites enable widespread, regional or global surveillance of volcanoes and can provide the first indication of volcanic unrest or eruption. Here we consider Interferometric Synthetic Aperture Radar (InSAR), which can be employed to detect surface deformation with a strong statistical link to eruption. Recent developments in technology as well as improved computational power have resulted in unprecedented quantities of monitoring data, which can no longer be inspected manually. The ability of machine learning to automatically identify signals of interest in these large InSAR datasets has already been demonstrated, but data-driven techniques, such as convolutional neutral networks (CNN) require balanced training datasets of positive and negative signals to effectively differentiate between real deformation and noise. As only a small proportion of volcanoes are deforming and atmospheric noise is ubiquitous, the use of machine learning for detecting volcanic unrest is more challenging than many other applications.

In this paper, we address this problem using synthetic interferograms to train the AlexNet CNN. The synthetic interferograms are composed of 3 parts: 1) deformation patterns based on a Monte Carlo selection of parameters for analytic forward models, 2) stratified atmospheric effects derived from weather models and 3) turbulent atmospheric effects based on statistical simulations of correlated noise. The AlexNet architecture trained with synthetic data outperforms that trained using real interferograms alone, based on classification accuracy and positive predictive value (PPV). However, the models used to generate the synthetic signals are a simplification of the natural processes, so we retrain the CNN with a combined dataset consisting of synthetic models and selected real examples, achieving a final PPV of 82%. Although applying atmospheric corrections to the entire dataset is computationally expensive, it is relatively simple to apply them to the small subset of positive results. This further improves the detection performance without a significant increase in computational burden (PPV of 100%). Thus, we demonstrate that training with synthetic examples can improve the ability of CNNs to detect volcano deformation in satellite images, and propose an efficient workflow for the development of automated systems.


The Application of Convolutional Neural Networks to Detect Slow, Sustained Deformation in InSAR Time Series

Automated systems for detecting deformation in satellite interferometric synthetic aperture radar (InSAR) imagery could be used to develop a global monitoring system for volcanic and urban environments. Here, we explore the limits of a convolutional neural networks for detecting slow, sustained deformations in wrapped interferograms. Using synthetic data, we estimate a detection threshold of 3.9 cm for deformation signals alone and 6.3 cm when atmospheric artifacts are considered. Overwrapping reduces this to 1.8 and 5.2 cm, respectively, as more fringes are generated without altering signal to noise ratio. We test the approach on time series of cumulative deformation from Campi Flegrei and Dallol, where overwrapping improves classification performance by up to 15%. We propose a mean-filtering method for combining results of different wrap parameters to flag deformation. At Campi Flegrei, deformation of 8.5 cm/year was detected after 60 days and at Dallol, deformation of 3.5 cm/year was detected after 310 days. This corresponds to cumulative displacements of 3 and 4 cm consistent with estimates based on synthetic data.