Towards Global Volcano Monitoring Using Multisensor Sentinel Missions and Artiﬁcial Intelligence: The MOUNTS Monitoring System

: Most of the world’s 1500 active volcanoes are not instrumentally monitored, resulting in deadly eruptions which can occur without observation of precursory activity. The new Sentinel missions are now providing freely available imagery with unprecedented spatial and temporal resolutions, with payloads allowing for a comprehensive monitoring of volcanic hazards. We here present the volcano monitoring platform MOUNTS (Monitoring Unrest from Space), which aims for global monitoring, using multisensor satellite-based imagery (Sentinel-1 Synthetic Aperture Radar SAR, Sentinel-2 Short-Wave InfraRed SWIR, Sentinel-5P TROPOMI), ground-based seismic data (GEOFON and USGS global earthquake catalogues), and artiﬁcial intelligence (AI) to assist monitoring tasks. It provides near-real-time access to surface deformation, heat anomalies, SO 2 gas emissions, and local seismicity at a number of volcanoes around the globe, providing support to both scientiﬁc and operational communities for volcanic risk assessment. Results are visualized on an open-access website where both geocoded images and time series of relevant parameters are provided, allowing for a comprehensive understanding of the temporal evolution of volcanic activity and eruptive products. We further demonstrate that AI can play a key role in such monitoring frameworks. Here we design and train a Convolutional Neural Network (CNN) on synthetically generated interferograms, to operationally detect strong deformation (e.g., related to dyke intrusions), in the real interferograms produced by MOUNTS. The utility of this interdisciplinary approach is illustrated through a number of recent eruptions (Erta Ale 2017, Fuego 2018, Kilauea 2018, Anak Krakatau 2018, Ambrym 2018, and Piton de la Fournaise 2018–2019). We show how exploiting multiple sensors allows for assessment of a variety of volcanic processes in various climatic settings, ranging from subsurface magma intrusion, to surface eruptive deposit emplacement, pre / syn-eruptive morphological changes, and gas propagation into the atmosphere. The data processed by MOUNTS is providing insights into eruptive precursors and eruptive dynamics of these volcanoes, and is sharpening our understanding of how the integration of multiparametric datasets can help better monitor volcanic hazards.


Introduction
About 1500 volcanoes are considered active worldwide [1], with about 50-85 erupting volcanoes each year [2]. Unfortunately, due to the cost and difficulty to maintain instrumentation in volcanic designed to easily incorporate new targets to monitor, thanks in particular to the global-coverage and free-access of the data. We use freely-available processing toolboxes to perform routine data processing, and allow the user to investigate the multi-parametric results related to a given volcano through an open-access website (www.mounts-project.com). We demonstrate the utility of such an interdisciplinary approach through a number of recent eruptions, describing a range of volcanic processes which can be tracked, from subsurface magma migration, to surface eruptive deposit emplacement, pre/syn-eruptive morphological changes, and aerosol propagation into the atmosphere. In addition, we show how artificial intelligence can play a key role in monitoring tasks. We show how a pre-trained Convolutional Neural Network (CNN) can be incorporated into the processing pipeline to detect large deformation in InSAR interferograms (e.g., related to dyke intrusions), in an automated, timely, and robust fashion.

Background: Existing Space-Based Monitoring Systems
The majority of existing volcano monitoring systems are set up and operated by national monitoring agencies, and consequently focus on volcanoes located within their own political boundaries or limited regions of interest. The U.S. Geological Survey (USGS) for example, provides earthquake information on a global scale, but their Volcano Hazards Program monitors volcanoes on U.S. territory only (https://volcanoes.usgs.gov). Each monitored volcano has a dedicated webpage, where a standardized red-yellow-green Volcano Alert Level System (VALS) informs on the state of volcanic activity [39]. Similar systems are operated by several other countries in charge of monitoring active volcanoes on their territory. However, many volcanoes are located on political boundaries between neighboring countries and states, and therefore leave uncertainties concerning the responsibilities during eruptive crisis. In this regard, only a few satellite-based platforms have the scope to monitor volcanoes on a global scale. Such platforms are often operated by scientists at research institutions and universities, and have proven to be particularly important for volcano observatories which do not have the resources to operate a monitoring network. These however, usually focus on the analysis of one type of data: spectroradiometric imagery (IR and UV), or SAR. We hereafter briefly describe a selection of relevant operating satellite-based platforms which provide open data access on a global scale (Table 1), and expose afterwards the rationale behind MOUNTS.

Spectroradiometry-Based Systems (IR, UV)
The development of thermal remote sensing volcanic activity has seen a long evolution since the 1960s, with a significant expansion and development of autonomous monitoring systems after the beginning of the new millennium [14].
MODVOLC, developed in 2002 by the Hawai'i Institute of Geophysics and Planetology [40][41][42], is a widely used volcano thermal monitoring system which can be considered as the archetype of open-access satellite monitoring services. It is based on the autonomous analysis of MODIS (Moderate Resolution Imaging Spectroradiometer) IR data, which has 1 km spatial resolution. The system provides the location of high-temperature thermal anomalies, spectral radiance data, and heat flux estimates in NRT over all volcanoes around the globe, through an open interactive Web-GIS interface (http://modis.higp.hawaii.edu/). MIROVA (Middle InfraRed Observation of Volcanic Activity, www.mirovaweb.it) is an enhanced NRT volcanic hot-spot detection system, based on the analysis of the MODIS Middle InfraRed (MIR) bands, developed by the University of Torino and the University of Firenze [43]. It currently monitors 215 volcanoes, and displays time series of various measurements, including the Volcanic Radiative Power (VRP) which measures the heat radiated by the volcanic activity, and higher level post-processed products such as lava effusion rates [44,45]. MOUNTS is strongly inspired by MIROVA, whose website template was used.
HOTVOLC is a monitoring system designed to achieve NRT monitoring of volcanic activity using geostationary satellite data (e.g., MSG-0 SEVIRI). It uses both IR and UV spectroradiometry analysis to monitor volcanic thermal activity, lava effusion rate [46], SO 2 gas plume [47], and volcanic ash plume [48]. It is developed at the Observatoire de Physique du Globe de Clermont-Ferrand (France), and monitors~50 volcanoes worldwide. HOTVOLC provides an interactive Web-GIS interface on which geocoded raster images and time series of the above mentioned parameters are visualized (http://hotvolc.opgc.fr).
The Global Sulfur Dioxide Monitoring Home Page, developed and maintained by NASA's Atmospheric Chemistry and Dynamics Laboratory (www. https://so2.gsfc.nasa.gov/), is a reference portal that brings together volcanic and anthropogenic SO 2 emission data from various monitoring systems. The results are based on the analysis of IR and UV data from different sensors that allow the development of the first inventory of volcanic SO 2 emission at global scale [49]. The website provides an archive of daily SO 2 emissions detected over different regions of volcanic interest that can be freely downloaded. NRT maps are also available through the various systems linked to the Global Sulfur Dioxide Monitoring Home Page.

SAR-Based Systems
SARVIEWS, developed by the University of Alaska Fairbanks (http://sarviews-hazards.alaska.edu), is a monitoring service based on the analysis of SAR imagery. The processing flow is triggered by an external alert system informing on an ongoing hazardous event, including volcanic eruptions, earthquakes, floods, or fires. By default it provides wrapped interferograms (6-day and 12-day) and Radiometric Terrain Correction (RTC) images computed over large spatial scales (~250 km wide). A Hybrid Pluggable Processing Pipeline (HyP3) allows on-demand higher level analysis, such as change detection algorithms. The derived products, although not available in open-access, have proven to be useful for operational volcano monitoring [50,51].
LiCS (Looking inside the Continents from Space) is developed by COMET (Centre for Observation and Modelling of Earthquakes, Volcanoes and Tectonics). The InSAR portal (https://comet.nerc.ac.uk/ COMET-LiCS-portal/) is an interactive map which provides links to download interferograms and coherence maps computed from Sentinel-1 images on large spatial scales. Although the initial focus was the Alpine-Himalayan tectonic belt, the system has recently expanded to process volcanic areas globally. A volcano deformation database provides access to interferograms and coherence maps at several hundreds of volcanoes (https://comet.nerc.ac.uk/volcanoes/), with the intention to be used for monitoring purposes in the future [26].
The University of Miami Geodesy Laboratory developed a system which provides InSAR displacement time-series at a smaller number of selected volcano (https://insarmaps.miami.edu/). The system relies on automatic InSAR time series analysis of Sentinel-1 products, which allows visualization of the temporal evolution of deformation on a large number of points across the volcano. Such approaches can be useful to detect pre-eruptive inflation, which is often observed in closed volcanic systems [52]. Table 1. Selected operational space-based volcano monitoring systems having a global scope. Although both MOUNTS (Monitoring Unrest from Space) and the Insar Viewer (Univ. Miami) incorporate ground-based open-access data in their platform, we here only describe satellite-based data. We refer the reader to the websites and publications for more in-depth detail on the various systems. Revisit and repeat frequencies are indicative, as these vary depending on the latitude and acquisition plan of the space agencies. The revisit frequency indicates the rate at which a volcano is imaged without considering repetitiveness of relative orbits, whereas the repeat frequency considers coverage ensured from the same repetitive relative orbit (i.e., necessary for InSAR analysis in particular). Acronyms used in the

Commercial Platforms
In recent years, commercial platforms have emerged, offering free visualization of open-access satellite imagery. Even though these are not specifically dedicated to volcano monitoring applications, they can be useful to volcano observatories which do not have the capability to download and process large volumes of data.
EO-Browser and Sentinel-hub Playground in particular (https://apps.sentinel-hub.com/sentinelplayground/ and https://apps.sentinel-hub.com/eo-browser/, respectively), make it possible to browse within a Web-GIS the complete archive of ESA Sentinel missions (Sentinel-1, Sentinel-2, Sentinel-3, Sentinel-5P), NASA Landsat missions (Landsat-5, 7 and 8), NASA MODIS imagery, and more. It gives the possibility to export image products in various formats, create time-lapse animation clips, and provides simple statistical analysis tools. Additional features are offered in commercial packages.
Google Earth Engine (GEE) gives access to several catalogs of open-access satellite imagery and geospatial datasets, and provides a platform to develop and run algorithms on Google's cloud infrastructure at no cost. While the GEE Explorer provides a simple web interface to visualize satellite imagery, the GEE Code Editor (https://code.earthengine.google.com/) provides a web-based IDE (Integrated Development Environment) for the GEE JavaScript API (Application User Interface). This allows users to develop algorithms and run them on the cloud without having to download the data. Although useful for processing large amounts of data, such service is not at the moment an immediate resource for volcano monitoring.

Rationale of the MOUNTS System
Global volcano monitoring has the aim to assemble and provide timely information of physical and geochemical parameters at a large number of (or all) volcanoes worldwide. The information can in turn be useful to both the scientific and operational communities, to understand volcanic processes and mitigate volcanic hazards. The rationale behind MOUNTS is based on the following three requisites: (1) To analyze the multi-sensor dataset offered by the Sentinel missions to provide the most complete and comprehensive understanding of volcanic processes ( Figure 1). Specifically, the datasets considered are: (i) SAR imagery from Sentinel-1 (S1) to monitor surface deformation, morphological changes, and reflectivity changes, (ii) SWIR imagery from Sentinel-2 (S2) to monitor surface thermal anomalies, and (iii) TROPOMI SO 2 data from Sentinel-5P (S5P) to monitor volcanic SO 2 gas plumes. Ground-based data is also incorporated, by querying global earthquake catalogues (GEOFON and USGS) for earthquakes recorded in the vicinity of the volcano. The volcanic processes which can be apprehended from the integration of these datasets is summarized in Figure 1, and those flagged with an asterisk are illustrated in the paper.
(2) To disseminate the results of the analysis as: (i) geocoded images, and (ii) time series of key parameters extracted from these, so as to have a comprehensive understanding of their evolution through time. The results are published on an open-access website (www.mounts-project.com), with a design similar to MIROVA, whereby each monitored volcano is assigned a webpage following a standard template. Dissemination of these results is achieved in NRT, typically within <1-6 h following the availability of the source product on ESA's Sentinel Data Hub. The availability of such products after sensing varies depending on the data type and processing level: typically <24 h from sensing for Sentinel-1 SLC (Single Look Complex) products, between 2 and 12 h from sensing for Sentinel-2 L1C (Level-1 C) products, and <3 h from sensing for Sentinel-5P NRTI (Near Real Time) SO 2 products.
(3) To implement a modular architecture, allowing easy incorporation of (i) new data types to analyze, (ii) new targets to monitor (i.e., new volcanoes, or new regions of interests for a given volcano depending on the activity, e.g., Kilauea 10 × 10 km extent during regular volcanic activity, and/or Kilauea East Rift Zone 40 × 60 km extent during flank eruption), and (iii) new algorithms to process the data (i.e., new processing chains, or trained neural networks to solve specific tasks).
MOUNTS is the first system to integrate all these components on a unique platform, with open-access, and global coverage capability. Figure 1. Utility of different sensor types, both spaceborne and ground-based, used by MOUNTS to monitor various volcanic processes. The occurrence of these may vary strongly from one volcano to another [11,26], and therefore the timings reported are given in an indicative manner. Processes flagged with an asterisk are those illustrated in this paper. The approximate revisit frequency of each satellite product used is indicated with a colored-marker (see Table 1 for details). Apart from Sentinel-3, all data are available on MOUNTS website. USGS: U.S. Geological Survey; SAR: Synthetic Aperture Radar; SWIR: short-wave infrared (SWIR).

Workflow of the MOUNTS System
MOUNTS is based on freely available remote sensing data, open source toolboxes, and a modular architecture. The workflow is outlined in Figure 2, and the basic processing information is described hereafter. We refer the reader to the Supplementary Material S1 for more in-depth technical details on various aspects of the processing chain. The entire framework is managed with Python, and is composed of the following elements: (1) Data query: the ESA Sentinel Data Hub is queried for new products at regular time intervals (typically every hour) using the Open Search API. Query options include the product type (i.e., platform name, file type), the product sensing time, and the footprint of the area of interest (AOI). Products types used by MOUNTS are the following: Sentinel-1 Level 1 (L1) Single Look Complex (SLC) Interferometric Wide (IW) products, Sentinel-2 Level 1C (L1C) products, and Sentinel-5P Level 2 (L2) Near Real Time (NRTI) SO 2 products. The standard AOI for a specific volcano is defined as a 10 × 10 km mask centered around the volcano summit, but can be adapted based on specific requirements (e.g., Kilauea Rift Zone 40 × 60 km to monitor flank eruptions). Volcano name, coordinates, and identification number are taken from the Global Volcanism Program (GVP) database [1], maintained by the Smithsonian Institution. The volcanoes monitored by MOUNTS and the associated relevant information for data query and processing are stored in a local database.
(2) Data download: when a new product is found on the Sentinel data hub, the data is automatically downloaded, and its metadata is stored in a local database.
(3) Data processing: as soon as a product is downloaded it goes through the dedicated processing chain, which intends to generate both image files (geocoded) and extract key parameters which help visualize volcanic activity. Products are opened and processed using ESA's open-source SNAP software (Sentinel Application Platform) through its Python API. The processing techniques applied to Sentinel-1, Sentinel-2, and Sentinel-5P products are described in Section 3.3.
(4) Data dissemination: the outputs of the processing chain are displayed on an open-access website (www.mounts-project.com). The web-design is based upon that of MIROVA, whereby all monitored volcanoes have a dedicated webpage displaying the same set of information. A range of tools have been incorporated to conveniently visualize the multi-parametric data, and thereby evaluate the volcanic activity ( Figure 3).

Sentinel-1
The Sentinel-1 mission comprises two polar-orbiting satellites (Sentinel-1A and Sentinel-1B, launched in April 2014 and April 2016, respectively), each carrying a SAR operating at C-band [57]. As it is an active remote sensing instrument sending its own electromagnetic radiation, and as radar wavelengths are mostly unaffected by weather clouds, it provides data in all-weather, day and night conditions, making it particularly advantageous for volcano monitoring. The two-satellite constellation offers a 6 day exact repeat cycle at the equator. However, since the orbit track spacing varies with latitude, the revisit rate is significantly greater at higher latitudes than at the equator. MOUNTS uses Interferometric Wide (IW), Single Look Complex (SLC) products, where each pixel is encoded as a complex number containing information on both the intensity and the phase of the signal backscattered to the radar sensor. Both information can be exploited to retrieve parameters relevant to volcano monitoring [12,26]. In turn, MOUNTS operates two processing pipelines: a SAR pipeline to monitor reflectivity changes due to morphological changes (incoherent change detection), and a DInSAR (Differential SAR Interferometry) pipeline to monitor surface deformation and emplacement of eruptive deposits (coherent change detection).

SAR Processing
The amplitude component of a SAR image measures the energy backscattered to the radar antenna after the electromagnetic wave hits the surface. It depends on the surface slope (i.e., incidence angle of the wave front onto the surface), the surface roughness (i.e., specular or diffuse reflector), and the surface material (i.e., the dielectric constant, particularly affected by moisture content). For instance, a terrain with a slope facing towards the radar sensor will backscatterer a lot of energy to the radar (high amplitude signal), which will appear as bright pixels on the SAR image. Conversely, a slope facing away from the radar will appear as dark pixels (low amplitude). Volcanic ash will tend to smooth the terrain on which it is deposited, thereby acting like a specular reflector, which will tend to backscatter more energy if the radar viewing angle is orthogonal to the slope, but less energy if looking at an oblique angle [23].
As soon as an S1 product is downloaded, it goes through the SAR processing chain which creates a geocoded SAR intensity image (see Supplementary Material S1 for details). Visual inspection of this image can be useful for volcano monitoring even for non-experts, especially when optical images are not available because of cloud coverage. In particular, large surface topography changes and eruptive deposits are usually easily identifiable (see Section 4.3).
Retrieving parameters which can help characterize volcanic activity over time in a robust and automated fashion is however a difficult task. Various approaches have been used to detect and visualize changes in SAR imagery. "Low-cost" approaches usually involve "incoherent change detection", which identifies surface intensity changes in successive SAR images. These changes can be visualized in a number of ways: (i) as difference maps, whereby the arithmetic difference, logarithmic ratio, or normalized change index (NCI) between two images is plotted [58][59][60], or (ii) as composite RGB image, where two bands are assigned to the first and second SAR image, and the third band to one of the difference map described above (e.g., [23,24]). More advanced techniques track features in the radar image prior to geocoding (i.e., radar coordinates in range and azimuth) to quantify morphological changes. Tracking radar shadows, for example, has proven useful to detect crater deepening [23], infilling of valleys by pyroclastic flows [24], and lava lake level variations [61]. Tracking the position of specific reflectors on the other hand, has proven capable of quantifying volcanic dome growths [62]. Such techniques unfortunately require a lot of manual tuning to adapt to each volcano, and are consequently a difficult approach for global monitoring. In this paper we instead favor the logarithmic ratio approach, in which we color code areas where the intensity has increased and decreased to red and blue respectively (see Section 4.2.2). DInSAR Processing SAR interferometry exploits the phase difference between two SAR images having similar viewing geometries, in order to measure both surface topography and displacements. More specifically, InSAR processing is commonly used to generate digital elevation maps, whereas DInSAR processing (Differential InSAR) successively subtracts the terrain elevation to retrieve surface deformation between the two acquisitions [63][64][65][66]. While these techniques typically exploit only two images, more advanced techniques known as "time series processing" exploit multiple images to retrieve more accurate deformation time series (i.e., Permanent Scatters PS and Small Baseline SBAS, see [67] for a review). The accuracy difference between DInSAR and time series approaches is about one order of magnitude, i.e., from centimeter to millimeter accuracy [68]. However, time series techniques require a lot of computing power and are therefore not performed by MOUNTS, which instead focuses on the two-image DInSAR approach.
Each time a new SAR image (S1 IW SLC product) is downloaded by MOUNTS, it goes through a DInSAR processing chain with the preceding acquisition (shortest temporal baseline, typically 6-12 days) from the same relative orbit. Processing of the image pair is performed with SNAP, freely distributed by ESA, and managed through a Python API. Back-geocoding and terrain correction are by default achieved with the SRTM-3 DEM. No multilooking is performed, so that the ground pixel resolution of the output products is 14 × 14 m. We refer the reader to the Supplementary Material S1 for more details on the processing chain. The outputs are typically an interferogram and a coherence map. The interferogram depicts the deformation of the ground along the sensor's line-of-sight (LOS) direction, and is displayed with values ranging between 0 and 2π which appear as a series of fringes (i.e., "wrapped" interferogram). The coherence is a measure of the interferometric quality, as it is related to the level of noise affecting the phase difference between the two SAR images. It ranges between 0 and 1, and reflects the degree of surface change, where values towards 0 indicate a loss of InSAR coherence (i.e., decorrelation). The main contributions to decorrelation are phase noise due to (i) the temporal change of the scatterers (e.g., vegetation, water, sand-covered areas will appear highly incoherent as the scatterers change continuously), (ii) geometric decorrelation (i.e., images have slightly different look-angles), and (iii) variation in atmospheric water vapor content. Although decorrelation usually renders interferograms useless for measuring ground deformation, it can be used to detect changes in the ground properties (e.g., emplacement of lava flows and eruptive deposits, delineation of eruptive fissures, etc.).
Both the interferogram and coherence image are successively analyzed to derive quantitative parameters, which are meant to be plotted as time series and consequently inform on the evolution of relevant volcanic processes through time. The interferogram analysis is intended to detect strong ground deformation occurring in short time intervals (i.e., between two successive S1 acquisitions, usually 6 to 12 days apart), which typically imprint as multiple narrow fringes (one fringes in Sentinel-1 interferograms represents 2.8 cm displacement in the radar line-of-sight). Processes such as magma intrusion (dykes), caldera collapse, or sudden flank movements are typically detected. The deformation detection is performed using artificial intelligence, which recovers filtered phase gradients in wrapped geocoded interferograms. Unwrapping of the phase gradients is successively performed to recover a deformation map, from which a deformation score (DEF) informing on the amount of ground displacement is computed. The approach is described in detail in Section 3.4 and in Supplementary Material S2. The coherence map, on the other hand, is analyzed to derive a decorrelation score (COH), which counts the number of pixels below a given coherence threshold. Although this parameter can fluctuate due to a number of non-volcanic processes (see phase noise sources described above), it can be useful to track eruptive deposits (e.g., fresh lava flows, pyroclastic deposits, etc.).
Hence, for each new SAR acquisition, an interferogram and coherence map are generated (shortest temporal baseline possible, typically 6 or 12 days), and disseminated on the website. From these, a deformation and decorrelation score are respectively computed, which are plotted as time series on the website.

Sentinel-2
The Sentinel-2 mission comprises two polar-orbiting satellites, Sentinel-2A and Sentinel-2B (launched in June 2015 and March 2017, respectively), placed at 180 • from each other in the same sun-synchronous orbit. They carry a Multispectral Instrument (MSI) providing multispectral data in 13 bands spanning from visible (VIS) to short-wave infrared (SWIR) wavelengths. The SWIR channels in particular, allow for detection of thermal emission produced by hot bodies. Considering a revisit frequency of five days under cloud-free conditions (reduced to 2-3 days at mid-latitudes), and a spatial resolution of 20 m/pixel in the SWIR bands, Sentinel-2 provide important measurements for volcano thermal monitoring. Specifically, analysis of thermal and morphometric features of heat sources facilitates tracking of a series of volcanic processes, including lava flow advancements, lava lake pulses, extrusion phases of lava domes, fumarolic activity, thermal activity at multiple active craters, and rise of magma column in open-vent volcanic systems (e.g., [69,70]).
MOUNTS downloads and analyzes Sentinel-2 Level 1C products. We use an approach presented in [55], which is an enhancement of the HOTMAP detection algorithm developed by [71]. It is based on fixed ratios between SWIR bands, and on a contextual threshold derived from a statistical distribution of the thermal anomaly clusters. Images are analyzed considering the TOA (Top of the Atmosphere) reflectance of the 12-11-8A bands (R: 2190 nm, 12; G: 1610 nm, 11; B: 865 nm, 8a). The algorithm detects the number of "hot" pixels in each downloaded image, and stores this parameter (SWIR) which can be displayed in time series. This thermal algorithm was already successfully verified on different volcanic cases worldwide, and compared with MIROVA data, showing a strong correlation between the detected number of "hot" pixels in Sentinel-2, and the VRP [55].

Sentinel-5P
SO 2 emissions are recovered from the imaging spectrometer TROPOMI on-board the Sentinel-5P satellite (launched in 2017, data available since December 2018). It offers data with a revisit frequency of 1 day, at exceptionally high spatial resolution (7 × 3.5 km 2 ) and detection limit (respectively 13 and 4 times better than the heritage Ozone Monitoring Instrument (OMI)), thereby opening the possibility to detect volcanic SO 2 with an unprecedented precision [16]. SO 2 slant column densities are retrieved by TROPOMI using differential optical absorption spectroscopy (DOAS) in the UV spectra (312-326 nm). Final gas densities are provided as vertical column densities (VCDs), for 1 km thick boxes at three different altitudes: 0-1 km (planet boundary layer, PBL), 6.5-7.5 km (mid-troposphere), and 14.5-15.5 km (upper troposphere). ESA disseminates the data as near-real time (NRTI) products, available during 1 month after acquisition, after which they are replaced by Offline (OFFL) products. OFFL products are based upon consolidated calibration and auxiliary data and thus provide better data quality data compared to NRTI products.
MOUNTS downloads Level 2 NRTI products, and processes them to retrieve SO 2 mass found in a 500 × 500 km box centered around the volcano. Vertical column SO 2 densities are first converted from mol·m −2 to the more commonly used Dobson Unit (DU), using a multiplication factor of 2241,15. SO 2 mass is then calculated following [72]: where MSO 2 is the mass of SO 2 (in tons), and A i and SO2 i are respectively the area (7 × 3.5 km 2 ) and the density (in DU) at each pixel i of the 0-1 km VCD (PBL). Pixels contaminated with SO 2 are successively isolated by creating a mask where DU > 1, to which is successively applied a morphological filter (i.e., erosion + dilatation operators, structuring element of 5 × 5 pixels), which allows for removal of noisy pixels. The remaining pixel clusters are assumed to be of volcanic origin, and the resulting SO 2 mass is assigned to the volcano for the time series display. The 500 × 500 km box size is chosen arbitrarily to capture large SO 2 plumes. However, several monitored volcanoes may be located within this area, potentially resulting in redundant SO 2 anomalies.
Although visual inspection of the image helps identify the emitting volcano, future developments will focus on using wind models to discriminate the emission source, and therefore recover more accurate SO 2 emissions for each volcano.

Machine Learning in Support of Deformation Detection
Machine learning algorithms can help solve specific tasks without being explicitly programmed, relying instead on patterns learned from a training dataset. We here show how trained Convolutional Neural Networks (CNN) can be plugged into the workflow to automatically detect large ground deformation in interferograms (Figure 4).
Synthetic data was first generated to train the network. The goal was to produce synthetic interferograms, wrapped and orthorectified, depicting realistic deformation fringes and phase decorrelation related to various artifacts (i.e., atmospheric phase delays, geometric distortions due to orthorectification, phase decorrelation due to vegetated areas, etc.). We synthesized the data from a combination of procedural noise and empirical rules, which although not physically motivated, can express a wide range of patterns similar to those found in real data (at least locally). The network was designed as a fully convolutional auto encoder built out of residual learning blocks (inspired by ResNet, [73]). Details on the synthetic data generation and network architecture can be found in Supplementary Material S2, and the trained CNN can be accessed on a Github repository (https://github.com/Andreas-Ley/SAR-InterfPhaseFilter).
The network was then trained to recover from the synthetic interferograms (provided as input), the associated phase gradients and phase decorrelation mask (expected outputs), Figure 4a. Once trained, the network is plugged into the processing pipeline ( Figure 4b): newly generated interferograms are fed through this pre-trained network, which yields both the phase gradients and phase decorrelation (i.e., noise mask). Unwrapping of these phase gradients is successively performed to recover approximate ground displacements (Figure 4c), where robustness against artifacts that might trigger false alarms is preferred over precise phase differences. The unwrapping simply solves for a displacement map whose gradients match the gradients estimated by the CNN. Since the CNN performs the filtering of the gradients and detection of noise, treatment of decorrellated areas in the unwrapping is rather simple. From this displacement map a deformation score is ultimately computed (DEF), as the standard deviation of the displacements expressed in meters in the radar line-of-sight (LOS). This data is then disseminated (Figure 4d) and can be used to plot time series that allow for identification of strong deformation events. In addition, automated email alerts are sent to interested users when a threshold is exceeded. Based on empiric evaluation of all 17 volcanoes currently monitored by MOUNTS, representing a total of >1360 interferograms, the threshold was set to 0.001 m. This value can however be tuned for each volcano, to adapt to specific volcanic activity or monitoring requirements.
A similar approach involving CNN has been presented by [34] to detect large deformation signals in short-duration, wrapped and geocoded interferograms. The main differences to our approach are the following: (i) training is performed on real interferograms, combined with a data augmentation approach to increase the number of interferograms where deformation is recorded; we instead create synthetic training data, allowing the generation of an unlimited number of interferograms, and avoiding the time-consuming task of labeling interferograms where deformation is identified; (ii) edge detection is applied to the interferograms (during both training and prediction processes), to identify phase jumps between fringes, which are subsequently considered as noise-free regions where deformation can be detected; our network instead is trained to output a decorrelation mask informing on the degree of phase noise, successively used to filter noisy regions in the interferogram where erroneous gradients could be detected, (iii) the network output is a deformation probability; our network outputs clean phase gradients (and decorrelation mask for additional cleaning), which can be successively unwrapped to compute ground deformation maps, from which a physically meaningful score can be derived (i.e., standard deviation expressed in meters), (iv) the authors fine-tuned a popular network (AlexNet, [74]), whereas the architecture of our network was designed from scratch and trained on our synthetic dataset, thereby allowing more flexibility in the desired outputs. Figure 4. Workflow to detect ground deformation using a pre-trained convolutional network: (a) Convolutional Neural Network (CNN) training using synthetic data (see Supplementary Material S2), (b) operational usage of the pre-trained CNN to detect to detect phase gradients and decorrelation mask in real interferograms (Piton de la Fournaise in this example), (c) phase unwrapping to recover ground displacements, and derive a deformation score defined as the standard deviation of the displacement map, (d) data dissemination, as images and time series on a public website (www.mounts-project.com), and email alerts to interested users when a threshold is overcome.
We compare the performances of the two approaches using the same evaluation metrics used in their study, which are derived from the confusion matrix: accuracy (ratio between correctly predicted samples and all testing samples), true positive rate (TPR, ratio between correctly identified positive samples and all positive samples), and true negative rate (TNR, ratio between correctly identified negative samples and all negative samples). Table 2 shows the performances of both networks, tested on the same volcanoes (Erta Ale, Ethiopia, and Etna, Italy), and highlights the good performances of our CNN. Table 2. Evaluation of the performances of the proposed CNN, and comparison with the CNN presented by [34] (AlexNet). The metrics are derived from the confusion matrix, where P is the number of real positives in the dataset (i.e., number of interferograms where deformation is recorded), N is the number of real negatives, TPR is the true positive rate, and TNR is the true negative rate. For each volcano we evaluate on the entire interferogram dataset (since training is performed on an independent synthetic dataset), whereas [34] report the average over multiple random splits into training and testing groups of equal size. Interferograms where deformation is expected were marked manually for ground truth (i.e., real positives P), and during evaluation deformation detection threshold was DEF = 0.001, the same threshold as for the email notifications. Synthetic data can only seldomly cover the full range of variations and all the corner cases that can be observed in the real world, and is thus usually inferior to real data. However, there are also benefits that make it an attractive choice: with synthetic data the ground truth values of what is to be predicted are known and can be used for the training. Also, since arbitrary amounts of data can be generated, overfitting due to insufficient data is not an issue, and the generalization capability of the final model depends primarily on the realism of the synthetic data. While the present work was under review, a study using synthetic interferograms to train neural networks was published [35]. Unlike our study, the authors generate deformation patterns based on analytic models simulating realistic deformation sources in volcanic settings (i.e., Mogi and Okada sources in particular), and model stratified atmospheric effects from weather models (i.e., GACOS). This approach incorporates expert knowledge into the network via the training data, and opens pathways to potentially recover deformation sources properties from the interferogram in an automatic fashion.

Results
We hereafter show how the products derived from S1, S2, and S5P analysis can help monitor various aspects of volcanic activity. The usefulness of each of these parameters varies from one volcano to another, due to both varying volcanic activity (i.e., effusive vs. explosive, open-vent vs. close-vent) and climatic setting (i.e., desert vs. tropical environment). For this reason, we here take various eruptive case examples, where the utility of each of these parameters is best illustrated. Below we successively describe processes related to magma migration towards the surface, effusive and explosive eruptive deposit emplacement, as well as SO 2 gas plume propagation in the atmosphere.

Detection of Surface Deformation (DInSAR, AI)
A number of processes can cause ground deformation at volcanoes: dyke intrusions, reservoir pressurization, caldera subsidence, cooling of eruptive deposits (lava/pyroclastic flows), landslides, etc. These can to some extent be distinguished based on distinctive patterns of ground displacement imaged in the interferograms [75]. As described previously, MOUNTS focuses on the automatic detection of processes generating large ground deformations in short time intervals (i.e., between two successive S1 acquisitions, usually 6 to 12 days apart), which typically imprint on interferograms as many colored fringes. In particular, dyke intrusions (whereby magma intrudes into a fissure and pushes the surrounding rock aside), are commonly characterized by two lobes with opposite displacements directions, that imprint as a "butterfly" shape. If the dyke reaches the surface, it will result in eruptive fissures/eruptive vents from which lava flows are emplaced. Figure 5 illustrates the efficiency of this automated deformation detection system, taking the case of Piton de la Fournaise (Réunion Island) as an example. The volcano experienced 5 intrusive episodes in the past year (April 2018-April 2019), all of which reached the surface and resulted in lava flows that lasted between 0.7-47 days (episode 1 and 4 respectively). The system successfully detected the interferograms where ground deformation is recorded (Figure 5a-c), illustrating its robustness against various artifacts in the geocoded interferogram (related to strong topography gradients, strong atmospheric phase noise, and decorrelation in vegetated areas). Furthermore, thermal anomalies detected by S2 SWIR analysis (i.e., number of hot pixels) are compared with the Volcanic Radiative Power (VRP) provided by the MIROVA system (Figure 5d) to testify the good correspondence between the active flow area and the heat radiated by the flow surface.
The systematic analysis of both ground deformation (amplitude, pattern, orientation) related to the intrusion of magma in the shallow portions of the edifice, and lava effusion rate once the magma breaches the surface, should help better understand the mechanisms controlling effusive eruptions in closed-vent volcano systems such as Piton de la Fournaise. . The eruptive episodes are numbered from 1 to 5, and highlighted by gray areas based on the eruption timing provided by the OVPF (Observatoire Volcanologique du Piton de la Fournaise): 0.7, 34.6, 0.8, 47, and 19.9 days respectively, publically available at http://www.ipgp.fr/fr/ovpf/activite-recente-piton-de-fournaise. Eruptive episode 5 is shown in detail in Figure 9, in which SO 2 data is also displayed.

Detection of Eruptive Deposits (SWIR, DInSAR, SAR)
Volcanic eruptions can generate various types of eruptive deposits, including lava flows, pyroclastic flows, mud flows, pyroclastic deposits, ash fall, etc. The best-suited parameter to detect these from space will depend on both the volcanic product type and the surface on which it is emplaced (i.e., arid surface, vegetated surface). We here describe two extreme eruptive case scenarios: the Erta Ale eruption (Ethiopia, ongoing since January 2017), characterized by an effusive activity where lava flows are emplaced in a desert environment, and the Fuego eruption (Guatemala, June 2018), characterized by a violent explosive activity which generated pyroclastic flows (i.e., avalanche of hot blocks and ash which flowed down the surrounding valleys in a densely vegetated environment).

Using DInSAR Coherence and SWIR
Detecting lava flows in non-vegetated regions is particularly efficient using the interferometric coherence [76]. Active lava flows will appear as highly incoherent areas (coherence values close to zero), contrasting with the coherent surface on which they are emplaced (coherence values typically >0.5). Moreover, because the summits of active volcanoes are usually vegetation-free, the coherence can also be utilized to monitor summital volcanic activity, as eruptive deposits will appear incoherent.
Both applications are illustrated in Figure 6, which depicts the ongoing Erta Ale eruption and its precursory activity. Indeed, in the months preceding the eruption onset, intense summit activity is indicated from the coherence map analyzed in a 2 × 2 km box around the summit (Figure 6a black curve). In particular, from June 2016 onwards the northwestern lake becomes active (in addition to the permanently active southeastern lake), and in the following months the activity at both lakes progressively intensifies, leading to multiple lava overflows (e.g., small overflow on 20 August 2016 visible on coherence image in Figure 6(i), and large overflow on 19 January 2017 visible on Sentinel-2 image). On 28 January 2017, the S1 interferogram shows a strong deformation signal with a distinctive pattern characteristic of a dyke intrusion (Figures 6c and 6(ii)), which marks the onset of a >2-year eruption which is still ongoing today [77]. The eruptive fissure opened multiple new vents located SE of the summit lava lakes (clearly identifiable from S1 coherence and amplitude images), which progressively focused on a single new eruptive vent from which lava flowed during the following months. During this time, the activity at the summit lava lakes significantly decreases (i.e., coherence analysis on 2 × 2 km extent, Figure 6a black curve), suggesting a possible drainage of the summit lava lakes. Both the coherence and the SWIR analyses computed on a large >20 km spatial scale (Figure 6a blue curve and Figure 6b respectively), depict the lava flow emplacement. The good agreement between the two highlights the fact that coherence can be used to track the active flow front, even when S2 images may not be usable due to cloud coverage. From this basic analysis, more elaborate parameters can be extracted offline, such as the lava flow front position through time, or the lava effusion rate, both key for hazard mitigation issues.
The coherence threshold can be adapted to each volcano, in order to account for specific environments (i.e., how incoherent is the surrounding land surface). More elaborate strategies can also be implemented to enhance the decorrelation sensibility to volcano-related processes, such as using NDVI masks (computed from S2 images) to exclude vegetated areas from the analysis. Nevertheless, in certain climatic settings coherence-based detection is simply not the best suited, and intensity-based detection can prove more useful. Figure 6. Example of lava flow emplacement detection at Erta Ale (Ethiopia) using both S1 interferometric coherence and S2 SWIR analysis. (a) Decorrelation area related to volcanic activity, calculated as the number of pixels where coherence <0.5, multiplied by the pixel area (14 × 14 m). Black curve: decorrelation on 2 × 2 km mask centered at the summit to image the lava lake activity. Blue curve: decorrelation on 50 × 50 km mask to image the lava flow emplacement. Because on such a large scale the coherence map also includes sandy areas where coherence is <0.5, these regions where manually excluded for this specific figure so that the curve only depicts the area of the active lava flow. (b) Number of hot pixels (×10 6 ) detected in the S2 SWIR image (50 × 50 km mask). (c) Deformation score DEF computed from the deformation map (unwrapped interferogram). Coherence maps (top image row) and SWIR images (bottom image row) are displayed at selected dates, with spatial scales varying according to the eruptive phase: (i) pre-eruptive phase with intense activity at the summit lava lakes, 2 × 2 km mask; (ii) eruptive onset, manifested by surface deformation and aperture of a new eruptive vent,~12 × 12 km mask; (iii) early stages of the effusive activity (highlighted by a gray box), showing the emplacement of lava flows on both NE and SW flanks of the volcano,~12 × 12 km mask; (iv) advanced stages of the effusive activity, when the lava flow front reaches its maximum distance from the vent, 11.5 × 19.5 km mask. The red box displayed in the coherence images refer to the 2 × 2 km mask centered on the active lava lakes.

Using SAR Intensity
The interferometric coherence is hardly exploitable in regions with very dense vegetation, or where the surface is likely to change rapidly due to various environmental factors (e.g., snow or sand). In such context, changes in the intensity of SAR images can help identify eruptive deposits. Figure 7 illustrates this, taking as example the 2018 eruption of Fuego (Guatemala) which killed over 200 persons due to pyroclastic flows [5]. The interferometric coherence image (Figure 7c) computed between SAR image 1 and 2, respectively acquired before and after the eruption, is entirely incoherent and therefore unusable. However, computing the log ratio between the two intensity images (Figure 7a) reveals substantial changes: areas in blue are those where the intensity has decreased, whereas areas in red are those where the intensity has increased. The latter are mainly confined in the valleys, and likely correspond to the rougher block-and-ash deposit [24]. Conversely, areas where intensity decreases are concentrated around the summit vent, which could be associated to the deposition of ash.

Detection of Morphological Changes (SAR)
The intensity of SAR images is strongly dependent on the terrain slope, and is therefore useful to monitor morphological changes affecting the volcano. Because radar wavelengths penetrate through clouds, SAR intensity images provide crucial insights into the volcanic activity when optical imagery is obstructed by atmospheric and/or volcanic gas clouds. We here give two examples taken from very different volcanological settings: the summit crater collapse of Kīlauea (Hawai'i) during the 2018 effusive eruption (Figure 8a), and the Anak Krakatau (Indonesia) island growth and destruction during the 2018 explosive eruption (Figure 8b).  Kilauea is well known for its persistent active lava lake. In 2018, it experienced its largest flank eruption and caldera collapse in the last 200 years [79,80]. During spring 2018, the lava lake activity was high, which was clearly detected as a hotspot in the S2 SWIR images (Figure 8a, 13 April 2018). On 30 April 2018, seismicity indicated the intrusion of a dyke along the East Rift Zone, which generated a~38 km long deformation zone, and multiple eruptive fissures with lava flows rapidly reaching the sea (see Supplementary Material S5 for analysis of S1 and S2 over the entire rift zone). During this time the summit underwent significant changes: lava lake withdrawal (i.e., hotspot disappears in SWIR images, Figure 8a), accompanied by summit subsidence (i.e., deflation detected in interferogram), progressively evolving in a~3 km wide caldera collapse (see LIDAR digital elevation model in [80]), as the shallow magma reservoir was being drained. The progression of the caldera collapse is clearly imaged with SAR intensity images, which reveal the progressive formation of fractures and the profound summit morphological changes accompanying the flank eruption (Figure 8a and video in Supplementary Material S3). Ash deposits following the eruption onset is also captured, identified by a decrease in the SAR backscattered intensity on the SE flank of the volcano, also visible in the SWIR images. This decrease can be explained by the fact that fresh ash is less reflective than bare rock owing to its loose structure and high porosity, and that ash deposits smooth the surface, resulting in a more specular reflector which backscatters less energy towards when the slope is facing away from the sensor [23].
Krakatau is well-known for its volcano-induced tsunamis. Just over 135 years after the famous 1883 event, the volcano triggered on 22 December 2018 another deadly wave. Analysis of the SAR intensity images clearly shows the progressive island growth in the months preceding the tsunami, due to multiple lava flows reaching the sea and extending the island's coast line (Figure 8b

Detection of SO 2 Gas Flux (UV)
On 18 February 2019 a new eruption of the Piton de la Fournaise began. According to OVPF Reports, 14 ± 5 Mm 3 of lava were erupted during the 18 days of activity, fed by several eruptive fissures located on the upper east flank of the volcanic cone. Following a phase of gradual increase in volcanic tremor and the intensification of surface activity, the eruption ended abruptly on March 10 (OVPF Reports). Figure 9 shows the time series of SO 2 mass burden recovered from S5P (Figure 9b), complemented by thermal anomalies recorded from S2 SWIR data (Figure 9c). The SO 2 mass obtained from MOUNTS is compared with the SO 2 mass computed by NASA as a measure of correlation between the two datasets. In addition, thermal anomalies detected by S2 are overlaid with the VRP data provided by MIROVA (Figure 9c) to show correspondence between the active flow area and the heat radiated by the flow surface. Once suitably calibrated, the combination of SO 2 and thermal data provides a synoptic view of the gas and magma fluxes during the course of the February-March 2019 eruption which can be used to track eruptive trends and patterns in real time [82]. Notably for this specific case, the two datasets show consistent trends, indicating a gradual intensification of the effusive and degassing activity during the final phases of the eruption. Thermal anomalies recorded after 10 March and in the absence of gas emission are attributed to the cooling of the lava field. The last SO 2 detection (10 March 2019 09:38 UTC) is attributed to the gas plume, no longer fed from the eruptive vent but still inside the 500 × 500 km AOI as it slowly drifts away from the island.

Combining Ground-Based and Space-Based Sensors
Magma migration within the crust generates stresses, which can result in earthquakes as the surrounding rocks are displaced or fractured. This seismicity, commonly known as volcano-tectonic (VT) seismicity, is often recorded both prior and during volcanic eruptions, within and around the volcanic edifice [83]. High magnitude VT earthquakes can be recorded and located by global seismological networks, even when the nearest seismic stations are installed several hundreds of kilometers away. In turn, their timing, location, magnitude, and sometimes focal mechanism, are stored in open access global earthquake databases, particularly GEOFON and USGS catalogues. MOUNTS facilitates the interrogation of such catalogs, recovering potential earthquakes recorded in a region centered around the monitored volcano. This data can support the analysis of the volcanic phenomena, especially when the volcano is not equipped with ground-based monitoring instrumentation. Figure 10 shows the recent eruption of Ambrym (Vanuatu), and illustrates how combining ground-based and space-based sensors helps understand the eruptive dynamics of this volcano located in a very remote and cloud-prone region. On 15 December 2018 and in the days that followed, a swarm of volcano-tectonic earthquakes were recorded in the vicinity of the volcano, with magnitudes ranging between~4.5 and 5.5 (Figure 10d). The volcano was known until then for its persistent activity characterized by two active volcanic lakes (Figure 10(b.1)), responsible for high heat and gas fluxes (Figure 10b,c respectively), [84,85]. Analysis of the SAR intensity images immediately before and after this swarm reveal profound morphological changes ( Figure 10(d.1,d.2)), in particular the collapse and enlargement of the summit crater. DInSAR analyses indicate very strong ground deformation during this period (Figure 10a red curve, Figure 10(a.1)), related to dyke intrusion and caldera subsidence [86]. Simultaneously, the decorrelation in the coherence map increases (Figure 10a blue curve), due to both the ground deformation and perhaps also pyroclastic deposits. Once stabilized, the coherence map reveals the presence of a new eruptive vent (Figure 10(a.2)), from which lava was most likely emitted, as suggested by the SWIR image acquired on 15 December 2018 (Figure 10(b.2)). Interestingly, following this event the volcano completely changed dynamics: the summit lava lakes were most likely drained, as suggested by the absence of thermal anomalies and the cessation of SO 2 gas emissions.

Discussion
The key to detecting volcano unrest and understanding the underlying mechanisms is to be able to recognize when a volcano is deviating from its background level of activity. Once the eruption starts on the other hand, the key to decipher the eruptive dynamics and to mitigate the related hazards, is to integrate multiparametric dataset streaming from both space-and ground-based sensors, in order to provide the most comprehensive view of the eruptive phenomena. Both require "monitoring", i.e., observing the volcanic activity over long periods of time, during both quiescent and eruptive phases. As such, the aim of monitoring platforms such as MOUNTS is twofold: (1) a scientific one, aiming at deepening our understanding of volcanic processes and patterns at stake at active volcanoes, by processing in a systematic way large amounts of data in an effort to construct global databases, and (2) a societal one, aiming at producing more successful eruption forecasts, and providing additional information to the operational community (e.g., local volcano observatories and civil protection) in order to mitigate the risks related to volcanic hazards.
The results presented in this paper intend to demonstrate how the monitoring platform MOUNTS can contribute to both scientific and operational aims. We here discuss the benefits, limitations, and future developments of the system. ,2) S1 intensity images immediately before and after the eruption onset, revealing profound morphological changes. Spatial extent of images a.1,2, b.1,2 and d.1,2 = 10 × 10 km, extent of images c.1,2 = 500 × 500 km.

Benefits of MOUNTS
The benefits of the developed system are the following: • Automated processing of free multisensor dataset which provide key parameters for volcano monitoring: surface deformation and reflectivity changes (Sentinel-1), heat anomalies (Sentinel-2), SO 2 gas emission (Sentinel-5P), and seismic activity (USGS and GEOFON earthquake catalogues). This interdisciplinary approach allows for assessment of a variety of volcanic phenomena in various volcanological contexts. Moreover, exploiting multiple sensors spanning across the electromagnetic spectrum reduces the dependency to sensing conditions (e.g., night, clouds).

•
Flexible design allowing fast implementation of new targets to monitor, with freedom regarding the size and shape of the region of interest. This allows to rapidly respond to new eruptive crisis, and adapt to the specific scenarios (e.g., monitoring of summit activity on small~1-2 km spatial extent, and/or monitoring of effusive activity on large >50 km extent). The system currently monitors 17 volcanoes in various volcanological and climatic settings across the globe, many of which recently experienced large eruptive crisis.

•
Visualization through an open-access website (www.mounts-project.com) of both geocoded images (i.e., DInSAR interferograms wrapped/unwrapped, DInSAR coherence map, SAR VV intensity image, SWIR B12-B11-B8A image, and SO 2 PBL concentration), and time series of parameters extracted from each image type (i.e., deformation score, number of decorrelated pixels, number of hot pixels, and SO 2 mass, respectively). This allows to apprehend the evolution through time of the volcanic activity and eruptive products. Download of full resolution images and time series graphs (PNG format) is readily possible from the website; more specific data download based on user-defined queries is planned, but for the moment possible on-demand only.

•
Based on the free SNAP toolboxes, providing a unique framework to manipulate data from various satellites, with state-of-the-art processing algorithms (e.g., DInSAR). MOUNTS is open-source, with a Github repository (https://github.com/sebastienValade/mounts) storing both the source code and a changelog informing on all the notable changes made to the system and website. • Modular architecture, allowing the implementation of new processing algorithms to extract relevant volcanological parameters, or solve specific tasks. As a matter of example, a pre-trained CNN was plugged to detect strong deformation in the interferograms generated by the system.

•
Automated email alert messaging to dedicated users when specific thresholds are overcome. Interaction with other monitoring systems such as MIROVA is achieved by facilitating access to volcano-dedicated webpages ( Figure 3). Strengthening the interactivity between the systems is planned, in particular by sharing database access in order to confront datasets more easily.
The operational community such as volcano observatories can use MOUNTS and contribute to its development in a number of ways. The IGEPN (Instituto Geofísico de la Escuela Politécnica Nacional) for example, responsible for volcano monitoring in Ecuador, suggested to add Sangay to the list of monitored volcanoes in order to contribute to the surveillance of this remote edifice. The data available on the platform was used freely, and a collaborative exchange was initiated upon request to provide more specific data processing. The resulting material was further analyzed by IGEPN staff according to their needs, and was used in the activity reports describing the ongoing crisis for public information [87]. (Disclaimers on the data usage and appropriate acknowledgements can be found on the website). Scientific collaborations to investigate specific volcanic processes, or to develop specific methods (based on either the dataset available on the website, or on datasets resulting from more complex analysis) are also welcomed.

Limitations of MOUNTS and Future Developments
The quantity of data available for volcano monitoring is increasing exponentially, but so is the difficulty to transform it into knowledge. Indeed a number of limitations arise, related to the extraction of meaningful parameters (are we looking at the right variables?), resolution issues (is the sampling in time and space accurate enough?), and data handling issues (how do we deal with the growing mass of data?). MOUNTS is at this stage still a proof-of-concept, which has large potential for improvement. We hereafter discuss the main limitations and development directions.

1.
Improve MOUNTS' capability to recover parameters informing on the state of volcanic activity, eruptive precursors in particular. While IR and UV spectroradiometry is able to provide rather straightforward parameters (i.e., heat and gas flux respectively), recovering parameters from SAR in a robust and automated fashion is more challenging. In this paper, we show how trained neural networks can achieve complex tasks in a timely and reliable manner, and can be easily implemented in operational processing chains. In particular, strong deformation typically imprint on interferograms as many colored fringes, which are successfully detected. Further development however is needed to detect slow deformation mechanisms, which do not generate deformation patterns with numerous fringes. Future developments should also focus on designing and training neural networks to recover from SAR data other relevant parameters that can inform on volcanic activity. For example, efficient change detection able to exclude changes non related to volcanic activity (e.g., snow fall, vegetation growths, etc.) would prove extremely useful during both pre-eruptive and syn-eruptive phases.

2.
Incorporate additional data types in the processing chain to provide further insights into the volcanic activity. A priority is to analyze Sentinel-3 (S3) TIR data routinely, crucial to monitor ground thermal anomalies at high sampling rate (sensor characteristics similar to MODIS), but also to detect ash plumes in the atmosphere. Ash detection is commonly achieved using the brightness temperature difference (BTD) procedure [88], applied to two channels centered around 11 and 12 µm. This approach is easily applicable, but also prone to generate false alarms. A number of methods have been developed to overcome this issue [20], including 3-band algorithms [48], the BTD algorithm with water vapor correction (BTD-WVC), the Robust Satellite Technique (RST) specifically configured for volcanic ash, and shallow neural networks [89,90]. Future developments should therefore implement automated S3 processing to monitor volcanic ash propagation in the atmosphere, which poses a major threat for air traffic in particular.

3.
Incorporate modeling tools to predict the propagation of volcanic eruptive products, using the recovered multiparametric data as input source terms for the models. Such strategies are not new, and are now, thanks to increasing computing power and data availability, becoming achievable in NRT to forecast the propagation of lava flows [91] or the dispersion of ash plumes into the atmosphere (e.g., [92,93]), as well as to predict the geometry and depth of magma bodies responsible for volcano deformation (e.g., [94][95][96]). Monitoring platforms such as MOUNTS should not necessarily include such modeling routines, but should at least strive to provide parameters than can be fed to such models.

4.
Migrate processing tasks on cloud platforms where data is archived. The development of MOUNTS was done by automatizing data download and processing, using free and open-source data and software. In doing so, we were able to keep the costs of this proof-of-concept platform very low (i.e., a single desktop computer manages data download, processing, and hosting of the web server). Nevertheless, this architecture limits the ability to process larger amounts of data, which would require the analysis of the entire available Sentinel dataset over hundreds of volcanoes worldwide. To achieve this, cloud computing strategies are preferable, whereby algorithms would run on a platform where the data is hosted, thereby preventing data download, and at the same time offering higher computing power. Commercial platforms offering such services exist: the Copernicus DIAS platforms (Data and Information Access Services, which include Sobloo, Onda, Creodias, Mundi, and Wekeo), the Amazon Web Services (AWS), or the Google Cloud Platform.

5.
Analyze the recovered multiparametric volcanic time-trends. Monitoring efforts such as the one presented here, allow the development of consistent multiparametric databases on a variety of volcanic settings (i.e., various volcano types, tectonic settings and magma compositions). Such databases are crucial to decipher eruptive patterns, and potentially better estimate future activity [25][26][27]29]. As a matter of example, decadal heat and gas emission time-trends help decrypt slow mechanisms of magma/gas accumulation and release at active volcanoes [97]. Clustering of these trends in categorically similar patterns, together with time-series analysis and probabilistic approaches (e.g., [98,99]), should be investigated to help decision making and potentially lead to better eruption forecasting. Moreover, incorporating standardized volcano alert level classifications (whereby color-codes help flag the activity of volcanoes [39,100]), will help better communicate the level of volcanic unrest and eruption likelihood to local populations and governmental authorities [101,102], and potentially lead to better early warning systems.

Conclusions
We present an operational volcano monitoring system, based on the automated download and processing of multisensor satellite-based data (Sentinel-1 SAR, Sentinel-2 SWIR, Sentinel-5P TROPOMI). The recovered data aim at providing key parameters able to inform on the state of volcanic activity, namely: surface deformation, surface reflectivity changes, surface heat anomalies, and SO 2 gas emissions. The results are disseminated in NRT on a public website (www.mounts-project.com), where both geocoded images and multi-parametric time series help understand the activity. Moreover, we demonstrate how artificial intelligence can be used in such monitoring system to solve complex tasks. In particular, we designed and trained a convolutional neural network to detect large deformation signals in wrapped interferograms with no atmospheric corrections. The training was done on synthetically generated interferograms, and evaluated on >1360 real interferograms produced by MOUNTS. Due to the very good performances of the network, it is now incorporated into the operational processing chain, which delivers automatic email alerts to dedicated users when strong deformation is recorded at the monitored volcanoes.
In addition to the set of parameters recovered from spaceborne sensors, we incorporate information available from global earthquake catalogues (GEOFON and USGS) to inform on the seismicity located in the vicinity of the volcano. The utility of integrating both satellite-based parameters (deformation, heat and gas) and ground-based parameters (seismicity) are demonstrated through a number of recent eruptions: Erta Ale 2017, Piton de la Fournaise 2018-2019, Fuego 2018, Kilauea 2018, Anak Krakatau 2018, and Ambrym 2018. We show how this interdisciplinary approach allows for assessment of a variety of volcanic phenomena, ranging from subsurface magma migration, to surface eruptive deposit emplacement, pre/syn-eruptive morphological changes, and SO 2 gas emission into the atmosphere. The data processed by MOUNTS is providing insights into the eruptive dynamics of these volcanoes, and is sharpening our understanding of how the integration of such multiparametric datasets can help better monitor volcanic hazards.