Space-borne Synthetic Aperture Radar (SAR) Interferometry (InSAR) is now a key geophysical tool for surface deformation studies. The European Commission’s Sentinel-1 Constellation began acquiring data systematically in late 2014. The data, which are free and open access, have global coverage at moderate resolution with a 6 or 12-day revisit, enabling researchers to investigate large-scale surface deformation systematically through time. However, full exploitation of the potential of Sentinel-1 requires specific processing approaches as well as the efficient use of modern computing and data storage facilities. Here we present Looking Into Continents from Space with Synthetic Aperture Radar (LiCSAR), an operational system built for large-scale interferometric processing of Sentinel-1 data. LiCSAR is designed to automatically produce geocoded wrapped and unwrapped interferograms and coherence estimates, for large regions, at 0.001° resolution (WGS-84 coordinate system). The products are continuously updated at a frequency depending on prioritised regions (monthly, weekly or live update strategy). The products are open and freely accessible and downloadable through an online portal. We describe the algorithms, processing, and storage solutions implemented in LiCSAR, and show several case studies that use LiCSAR products to measure tectonic and volcanic deformation. We aim to accelerate the uptake of InSAR data by researchers as well as non-expert users by mass producing interferograms and derived products.

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.

There are some 1,500 volcanoes with the potential to erupt, but most are not instrumentally monitored. However, routine acquisition by the Sentinel‐1 satellites now fulfils the requirements needed for interferometric synthetic aperture radar (InSAR) to progress from a retrospective analysis tool to one used for near‐real‐time monitoring globally. However, global monitoring produces vast quantities of data, and consequently, an automatic detection algorithm is therefore required that is able to identify signs of new deformation, or changes in rate, in a time series of interferograms. On the basis that much of the signal contained in a time series of interferograms can be considered as a linear mixture of several latent sources, we explore the use of blind source separation methods to address this issue. We consider principal component analysis and independent component analysis (ICA) which have previously been applied to InSAR data and nonnegative matrix factorization which has not. Our systematic analysis of the three methods shows ICA to be best suited for most applications with InSAR data. However, care must be taken in the dimension reduction step of ICA not to remove important smaller magnitude signals. We apply ICA to the 2015 Wolf Volcano eruption (Galapagos Archipelago, Ecuador) and automatically isolate three signals, which are broadly similar to those manually identified in other studies. Finally, we develop a prototype detection algorithm based on ICA to identify the onset of the eruption.

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.

Pressure influences both magma production and the failure of magma chambers. Changes in pressure interact with the local tectonic settings and can affect magmatic activity. Present-day reduction in ice load on subglacial volcanoes due to global warming is modifying pressure conditions in magmatic systems. The large pulse in volcanic production at the end of the last glaciation in Iceland suggests a link between unloading and volcanism, and models of that process can help to evaluate future scenarios. A viscoelastic model of glacio-isostatic adjustment that considers melt generation demonstrates how surface unloading may lead to a pulse in magmatic activity. Iceland’s ice caps have been thinning since 1890 and glacial rebound at rates exceeding 20 mm yr ⁻¹ is ongoing. Modelling predicts a significant amount of ‘additional’ magma generation under Iceland due to ice retreat. The unloading also influences stress conditions in shallow magma chambers, modifying their failure conditions in a manner that depends critically on ice retreat, the shape and depth of magma chambers as well as the compressibility of the magma. An annual cycle of land elevation in Iceland, due to seasonal variation of ice mass, indicates an annual modulation of failure conditions in subglacial magma chambers.