- Open Access
Effusive crises at Piton de la Fournaise 2014–2015: a review of a multi-national response model
© The Author(s). 2017
- Received: 27 January 2017
- Accepted: 29 May 2017
- Published: 12 June 2017
Many active European volcanoes and volcano observatories are island-based and located far from their administrative “mainland”. Consequently, Governments have developed multisite approaches, in which monitoring is performed by a network of individuals distributed across several national research centers. At a transnational level, multinational networks are also progressively emerging. Piton de la Fournaise (La Réunion Island, France) is one such example. Piton de la Fournaise is one of the most active volcanoes of the World, and is located at the greatest distance from its “mainland” than any other vulnerable “overseas” site, the observatory being 9365 km from its governing body in Paris. Effusive risk is high, so that a well-coordinated and rapid response involving near-real time delivery of trusted, validated and operational product for hazard assessment is critical. Here we review how near-real time assessments of lava flow propagation were developed using rapid provision, and update, of key source terms through a dynamic and open integration of near-real time remote sensing, modeling and measurement capabilities on both the national and international level. The multi-national system evolved during the five effusive crises of 2014–2015, and is now mature for Piton de la Fournaise. This review allows us to identify strong and weak points in an extended observatory system, and demonstrates that enhanced multi-national integration can have fundamental implications in scientific hazard assessment and response during an on-going effusive crisis.
- Effusive crisis
- Volcano observatory
- Piton de la Fournaise
- Time averaged discharge rates
- Lava flow model
- Inundation forecasts
- Hazard response
Sub-areal island volcanoes under European governance with historic activity, populations >100 within 30 km of the active center, and/or eruptions since 1900, as listed by the Smithsonian Institution Global Volcanism Program data base (http://volcano.si.edu/search_volcano.cfm). Distance from each state capital (Paris, Athens, Rome, The Hague, Oslo, Lisbon and Madrid) was obtained using the world distance calculator of GlobalFeed.com. Last eruption is as of 31/12/2016
Piton de la Fournaise
Ile de La Réunion
Saint Eustatius Island
La Palma Island
Agua de Pau
Picos volcanic system
Fossa di Vulcano
However, the same active volcanic islands tend to be something of an exotic notion to the mainland population, often being a popular sun-and-beach, tropical vegetation-and-temperature or exotic food-and-rum (or wine) holiday option, with the island being a “sun-mass tourist” destination or popular cruise-ship stop over (e.g., Bardolet and Sheldon 2008; Etcheverria 2014; Garau-Vadell et al. 2014; Silvestre et al. 2008). In some cases, recognizing the island as potentially active may even be deemed unwanted due to potential damage to the same tourism (Dominey-Howes and Minos-Minopoulos 2004). Montserrat is a well-known recent example. During the 1980s Montserrat was an exotic island in the Caribbean, not well-known for its volcanic activity. However, beginning in 1995 renewed activity covered a large part of the southern half of the island, including the principle town (Plymouth), in pyroclastic deposits; necessitating evacuation (Brown 2010). Vulcano (Italy) may be argued to have suffered a similar fate during the 1888–90 eruption. Having purchased and developed the northern part of the island for sulfur-and-alum mining, as well as grape cultivation, in 1870 for £ 8000, the Tyneside-based (UK) entrepreneur, James Stevenson, sold-up and left having seen his beautiful villa and prosperous enterprise destroyed by air fall and ballistics, Vulcano having become ‘an awful place’ (Stevenson 2009). Thus, for the local populations of such “exotic” locations, with well-established and tight-knit local communities, the hazard, risk, impact and loss, both tangible and intangible, due to a volcanic event is very real (e.g., Payet 2007), as was witnessed during the loss of Kalapana on Kilauea (Hawaii) to lava flow inundation during the 1990s (Weisel and Stapleton 1992).
Worse, on-site observatories on active European volcanic islands are either: (i) non-existent, (ii) offshore and/or (iii) lacking in numbers and resource. Thus, during a crisis, the local staff (if there are any) may become spread so thin in collecting and interpreting data, while maintaining equipment, reporting, forecasting and situation-advisory duties, that there is no time to put out calls for help to expand the monitoring and response network. Such a beleaguered staff will need all the help that they can get; but, help needs to involve the implementation of tested and trusted techniques that provide product that is useful and that can be merged seamlessly into their response, forecasting and reporting duties (RED SEED Working Group 2016). That is, in the terminology of the remote sensing community, tools need to have been validated against reference data or ‘ground-truth’ (Lillesand and Kiefer 1987), so that they have gone from experimental to operational (Rudd 1974), thereby being known, trusted, valid and useable. Such help needs to be tested before, not during, a crisis (RED SEED Working Group 2016). In this case, the best way-forward is an ensemble approach whereby external (to the observatory) partners with pertinent expertise are invited to contribute to the response effort. In such a case, all partners need to integrate fully, and openly, with the group so that each partner adapts their strengths, weaknesses and roles as situations and data dictate, with all partners being open to communication and data sharing across the entire group. We here review just such a response model by focusing on a multinational and multidisciplinary group active during the five recent eruptive crises of Piton de la Fournaise (La Réunion Island, France) that occurred between 2014 and 2015; Piton de la Fournaise being the furthest active European island volcano from its mainland administration center, Paris (France).
French Polynesia (Polynesie Française) is a “Collectivité d’Outre Mer” (COM) or an “Overseas Collective”. Here the seismic network is run by CEA (“Commissariat à l’énergie atomique et aux énergies Alternatives”) and real time data are locally transmitted to the “Laboratoire de Géophysique in Papeete” on Tahiti.
The “Terres Australes et Antarctiques Françaises” (TAAF) or the “French Southern and Antarctic Territories” includes those French overseas islands in the Indian and Antarctic Oceans, apart from La Réunion and Mayotte. The “Institut Polaire Française Paul-Émile Victor” in collaboration with the “Institut de Physique du Globe de Strasbourg” (and the Geoscope Observatory) are in charge of monitoring the Antarctic part of TAAF, which are labelled TOM (Territoires d’outre mer). The Austral part of TAAF is not permanently monitored.
Guadeloupe, Martinique, La Réunion and Mayotte are “Departements d’Outre Mer” (DOM). Active DOM volcanoes are monitored by “Institut de Physique du Globe de Paris” (IPGP) via a network of local volcano observatories. These are, respectively, Observatoire Volcanologique et Simologique de la Guadeloupe (OVSG), Observatoire Volcanologique et Sismologique de la Martinique (OVSM) and Observatoire Volcanologique Piton de la Fournaise (OVPF). The National Observation Service for Volcanology (SNOV) operated by the “Institut National des Sciences de l’Univers” (INSU) of the “Centre National de la Recherche Scientifique” (CNRS) is in charge of scientific duties, as well as collection and distribution of geological and geophysical data. Although Mayotte was built by volcanic activity (as was all of the Comoros Archipelago), it does not have a permanent seismic network. The most recent volcanic activity on Mayotte was 6.5 kyr BP (Zinke et al. 2001).
As part of this monitoring system, OVPF was built in 1980 in La Plaine des Cafres (15 km away from Piton de la Fournaise) to monitor volcanic activity on Piton de la Fournaise and Piton des Neiges (Fig. 1b), as well as to track seismic activity on and around La Réunion island. OVPF was set-up in the aftermath of the eccentric 1977 eruption whose lavas inundated the village of Piton Sainte Rose (Fig. 1b). Led by IPGP, whose headquarters are in Paris, OVPF has (as of December 2016) just 12 permanent staff who are based in La Plaine des Cafres (Fig. 1b). Five of these staff are scientists charged with data monitoring (checking data, derived parameters, trends, etc.), five others are engineers charged with instrument and network maintenance and monitoring. All ten have reporting duties both during eruptive and non-eruptive periods (checking activity and data availability, situation updates, and reports to the head of the observatory who prepare official bulletins, etc.). Within the framework of SNOV, OVPF staff collaborate closely with IPGP staff in Paris and other French groups, or National Partners (NP), mainly at La Réunion University, at the Observatoire de Physique du Globe in Clermont-Ferrand (OPGC), and OPGC’s academic companion department, Laboratoire Magmas et Volcans (LMV). The time difference between the OVPF and French sites is 2 or 3 h depending on the season, where Paris and Clermont Ferrand is UTC + 1 and La Plaine des Cafres is UTC + 4. Specific agreements between OVPF and other international agencies, or International partners (INP), such as the Hawaiian Volcano Observatory (USGS-HVO, Hawaii, USA), Istituto Nazionale di Geofisica e Vulcanologia (INGV – Pisa, Palermo and Catania, Italy) and Chalmers University of Technology (Gothenberg, Sweden), as well as informal arrangements with INPs such as Università di Torino (Turin, Italy), have also permitted data sharing, technology upgrade and knowledge transfer beyond the national framework. Here we note that, informal agreements based exclusively on mutual and collaborative efforts, as carried out here with Università di Torino, are developed during an eruptive crisis. However, as discussed later, continuity of service provision, data validation efforts, transparency and efficiency then benefit from development of more robust, formal agreements developed during non-eruptive periods.
During and between crises the role of OVPF is to communicate with the national, regional and local responding agencies. Following the French national plan for response protocols during a crisis at Piton de la Fournaise, the call-down procedure is laid out in “Organisation de la Réponse de SEcurité Civile” (ORSEC), i.e., the “ORSEC-Piton de la Fournaise” plan. Within this plan OVPF communicates only with the prefecture via the “Etat Major de Zone et de Protection Civile de Ocean Indien” (EMZPCOI). The prefecture (decentralized administrative service of the French state) then communicates with other actors in the response chain. For each eruption, and for any change in volcanic activity, a Volcano Observatory Notice for Aviation (VONA) is also sent to the Volcanic Ash Advisory Center (VAAC) in Toulouse (France). OVPF reporting duties also include communication with the air quality office in Saint Denis (La Réunion). The period between 2014 and 2015 was particularly challenging because five eruptions occurred in short succession, with one during June 2014, and then four during 2015 (in February, May, July–August and August–October), the last of which ranked as the fifth largest eruption since records began in 1700 (Michon et al. 2015; Peltier et al. 2016). During these events round the clock service was maintained by OVPF. During non-eruptive periods, when the observatory is not staffed out of working hours, full 24 h service (during all days of the year) is maintained by instrumental monitoring and an alarm, which triggers in the case of a change in seismic activity, where the alert is sent to a scientist-on-duty. Duties at the observatory between eruptions thus continue on a 24/7 basis, and include daily checking of activity and proper functioning of all networks, plus creation of daily reports.
OVPF monitoring network
On site, OVPF-IPGP maintains four types of ground-based real-time monitoring networks, namely (in order of network size): (i) seismic, (ii) geodetic (deformation), (iii) geochemical (gases), and (iv) imagery, including permanently installed visible and infrared cameras (Fig. 1c). The OVPF team is also charged with performing detailed syn-eruptive sampling and mapping of eruptive products (Additional file 1). At the national level, OPGC and LMV are charged with satellite remote sensing (OI2: “Observatoire InSAR de l’Océan Indien”, https://wwwobs.univ-bpclermont.fr/SO/televolc/volinsar/; HotVolc: HotVolc Observing System, https://wwwobs.univ-bpclermont.fr/SO/televolc/hotvolc/), plus petrochemical and volcanological analysis of the eruptive products (Dynvolc: Dynamics of Volcanoes, http://wwwobs.univ-bpclermont.fr/SO/televolc/dynvolc/; GazVolc: Observation des gaz volcaniques, http://wwwobs.univ-bpclermont.fr/SO/televolc/gazvolc/).
La Réunion: Volcanic hazard, risk and perception
The relatively late creation of a continuous monitoring system, where the observatory was established in 1979, together with the recent age of permanent human settlement on the island, where the first people arrived from Brittany (France) in the seventeenth century (Vaxelaire 2012a), have meant that the observatory and the population have had to deal with a growing awareness of the variability, in time and space, of the eruptive behavior of Piton de la Fournaise (see Morandi et al. 2016 for review). Since the creation of the observatory, 66 eruptions (1981–2016) have occurred with durations of between 0.3 and 196 days emitting, on average, a bulk volume of 9 × 106 m3 (Peltier et al. 2009; Roult et al. 2012). With the exception of the 1986 and 1998 eruptions, all lava emissions have been confined inside the uninhabited ‘Enclos Fouqué’ caldera (Fig. 1b), and have involved vents opening inside the summit craters or on the flanks of the central cone, or further away to the east on the floor of the Enclos Fouqué caldera. Eight of these eruptions (March 1986, March 1998, June 2001, January 2002, November 2002, August–October 2004, February 2005, and April 2007) have cut the island belt road – the only link between the southern and northern part of the island on the eastern flank (RN2, Fig. 1b). Large-volume eruptions have been documented, with that of April 2007 (Staudacher et al. 2009), being of quite short-duration (less than 1 month). Effusive events can also be long-lasting (several years), where decade-long phases of continuous activity, punctuated by short-lived explosive events, have been observed in the geological record (Peltier et al. 2012; Michon et al. 2013). In 2007, the most voluminous eruption of historical times occurred at an altitude of 570 m above sea level, 400 m north of the southern wall of the Enclos Fouqué caldera. The eruption buried 1.5 km of the belt road under 60 m of lava, caused gas exposure problems in towns on the east and west coasts, and prompted the evacuation of the nearby village of Le Tremblet. Likewise the eruptions of 1986 and 2002 also required evacuation of villages north or south of the Enclos Fouqué caldera.
There are around 840,000 habitants on La Réunion island, and over 245,000 permanent residents on the volcano flanks. Tourism is a major industry, accounting (in 2014) for 7.8% of the Gross Domestic Product, 8500 jobs (3.2% of total employment), 40.7% of total exports and 3.1% of total investment (WTTC 2015), where the warm waters, exotic marine life, white-sand beaches, surfing and tropical coastline are big draws (e.g., The Lonely Planet 2015; Michelin Green Guide 2015). However, tourism has been affected by two water-related hazards, the “shark crisis”, where there have been 42 attacks since 1990, and Chikungunya epidemics, both of which have been widely disseminated by global media (e.g., Santora 2015, Stewart 2015, Surfer Today 2016). Chikungunya is a mosquito-borne virus characterized by arthralgia or arthritis, where an outbreak between March 2005 and April 2006 resulted in 255,000 cases (i.e., it affected 30% of the population), was responsible for 87% of the deaths on the island during the same period (Josseran et al. 2006), and was followed by a sharp down turn in visitor numbers (INSEE 2016). As a result, the tourism strategy has turned to the attractiveness of the terrestrial environment. This includes the draw of an active volcano (Gaudru 2010), which is classed as a World Heritage Site by UNESCO. This policy re-orientation has allowed La Réunion to maintain visitor numbers at a level of around 420,000 per year since 2007, this being the same as the pre-2006 level; where numbers dipped to 300,000 in 2006 following the Chikungunya epidemic (INSEE 2016). Although the tourist office and local media promotes Piton de la Fournaise and its activity as a tourist attraction, portraying the volcano as dangerous or hazardous would have a negative effect on tourism. In terms of the resident population, a recent survey of Réunion’s resident population on volcanic risk perception thus revealed a relatively poor knowledge of the volcano and its activity; although, the same people had a high level of trust in scientists to provide accurate and reliable information (Nave et al. 2016).
The tourist draw and hazard to tourists
An on-going eruption is a positive draw to the Parc national de La Réunion, which covers 105,000 ha (or 42% of La Réunion island), and the access town (for Piton de la Fournaise) of La Plaine des Cafres and the visitors center (La Cité du Volcan) dedicated to the volcano and its activity. The first tourist activity at Piton de la Fournaise was in the form of scientific expeditions during the nineteenth century. These expeditions were usually composed of a foreign naturalist on a visit to the island, some local volcano experts, guides, porters, and governors or other senior administrative officials. At that time, a guide was quite rare and difficult to find on the island and the porters, initially slaves, were often frightened at the idea of going to the “territory of the devil”. In 1863, for example, Baron Carl Claus Von Der Decken (Kersten 1871; Kersten et al. 2016) spent 6 days exploring the volcano, departing from Saint Denis and pausing to find a guide at La Plaine des Palmistes. Such expeditions were increasingly facilitated by the opening of a railway, which was constructed along the coast from Saint Denis to Saint Pierre and Saint Benoît, in 1882 (Vaxelaire 2012b). In 1925 a first Gîte was built on the site of the present one near Pas de Bellecombe (Fig. 1b), and in 1933 construction of Madame Brunel’s Hotel at La Plaine des Cafres (near the current OVPF buildings – Fig. 1b) gave explorers a base camp. In 1957, the Office National des Forêts (National Forestry Office) initiated the construction of the “Route du Volcan”, which was completed as far as Pas de Bellecombe in 1968.
On a normal day 2000–3000 people use the main visitor access point for the Dolomieu crater (Bello 2010), this being Pas de Bellecombe (Fig. 1b). This load increases enormously during activity (Fig. 3b, c). Given a traffic density of 270 cars per kilometer at peak flow (estimated from Fig. 3b on the basis of nose-to-tail traffic in both directions at peak flow), this gives 1160 cars parked on the final 4 km of the road to Pas de Bellecombe. If we add the capacity of the Pas de Bellecombe car park (1200 cars), then this amounts to around 2400 cars. Given an average number of four people per car, this is a visitor load on Pas de Bellecombe during the night of 31 July - 1 August 2015 of almost 10,000. We round up because this does not take into account emergency parking opened at La Plaine des Sables, as well as the car park at, and 1.3 km road, to the Gite du Volcan (Fig. 1b), a restaurant/lodge 900 m north of Pas de Bellecombe into which around 100 stranded tourists broke-in, so as to shelter for the same (cold) night. However, this increase in traffic overwhelms parking and access facilities, causing severe congestion in towns lying along access roads to the Enclos Fouqué caldera and Dolomieu crater, including the towns of Saint Philippe and Sainte Rose on the RN2 (Bello 2010). A high tourist load can also be damaging to the flora of the park itself (Bello 2010), and will (i) be associated with an increased number of accidents and illness among park visitors (Heggie and Heggie 2004; Heggie 2005), (ii) result in increased need for emergency search and rescue operations (Heggie 2008; Heggie and Heggie 2008; 2009), and (iii) cause fatalities if tourists stray into dangerous environments (e.g. Heggie 2009) or unstable areas subject to collapse (e.g. Perkins 2006). Following the 31 July - 2 August 2015 eruption, the local newspaper reported 41 hikers evacuated from closed zones, 14 victims of minor injuries and illnesses, and 11 cases of hypothermia at Pas de Bellecombe, plus four cases of illness, some minor injuries and cases of fatigue on the road between La Plaine des Cafres and Pas de Bellecombe (L.R. 2015). The eruption was just 50 h in duration.
Our focus here is on running and delivering, in as timely-a-fashion as possible, lava flow simulations during an effusive crisis to assess and update likely inundation zones. We here use the merged output from two lava flow models: FLOWGO and DOWNFLOW, the FLOWGO model already having been initialized and validated for channel-fed flow at Piton de la Fournaise (Harris et al. 2015). Merged application of FLOWGO-DOWNFLOW has also been tested for near-real time lava flow simulation using a feed of lava time-averaged discharge rate (TADR) derived from 1 km spatial resolution thermal data from satellite-based sensors such as MODIS (Wright et al. 2008).
Any lava flow simulation model first requires initialization with vent location, which is provided by OVPF as part of their monitoring and response procedures. The key source term then becomes TADR, which can potentially be obtained in near-real time, and multiple times per day, from satellite-based sensors imaging in the thermal infrared. The Spinning Enhanced Visible and Infrared Imager (SEVIRI) sensor on the MeteoSsat Second Generation (MSG) satellite series provides thermal infrared data for Piton de la Fournaise at a spatial resolution of 3 km and nominal temporal resolution of 15 min (Gouhier et al. 2016). These data are potentially capable of providing TADR in a timely fashion. However, tests on GOES data have shown that, although variations in volcanic radiance can be trusted to provide arrival time of magma at the surface with a precision of 7 ± 7.5 min (Harris et al. 1997a) and to track variations in effusive activity at 15 min time steps (e.g., Mouginis-Mark et al. 2000), TADR derived from 4-km pixel data were not reliable (Dawn Pirie, Hawaii Institute of Geophysics and Planetology, unpublished data, 1999). TADR were not, thus, included as part of “hot spot” products delivered to recipient observatories by the Hawaii Institute of Geophysics and Planetology (HIGP, University of Hawaii, Honolulu, USA) because they were not deemed valid or trustworthy (Harris et al. 2001). This is a result of the large pixel size and mixed pixel problems, as well as pixel deformation, where pixels become increasingly large, ovoid, overlapping and rotated with scan angle (Harris 2013). Worse at high scan angles, or with extreme Earth curvature (which is the case for SEVIRI observations of Piton de la Fournaise), unreliable spectral radiances have been recorded (Holben and Fraser 1984; Singh 1988; Coppola et al. 2010), so that spurious data at high scan angles tend to be filtered out of quantitative analyses (e.g., Tucker et al. 1984; Goward et al. 1991; Harris et al. 1997b). Following Frulla et al. (1995), radiances for volcanic hot spots are deemed unreliable at scan angles of greater than 50°, where SEVIRI views Piton de la Fournaise at an angle of 63.4°. At such high scan angles, while over-estimates of spectral radiance, and hence TADR, will result from smearing of the anomaly due to extreme pixel overlap effects and point-spread-function problems (e.g., Markham 1985; Breaker 1990; Schowengerdt 2007), underestimates will result from atmospheric effects (Coppola et al. 2010) and topographic shadowing of all or part of the thermal anomaly (Dehn et al. 2002). Even local topographic features, such as cones, levees and skylights have been shown to play a role in shadowing the anomaly at high scan angles (Mouginis-Mark et al. 1994), causing detection problems even at quite low scan angles for active lava surrounded by topographic highs (e.g., Wooster et al. 1998; Harris et al. 1999; Calder et al. 2004).
Other source terms involve chemical and petrological data, to set and check rheological models used, as well as physical volcanological measurements (at-vent temperature, crystallinity and vesicularity). These data can, in turn, also be used to track the effusive event (Coppola et al. 2017). Finally, data processed post-event can then be used for de-briefing purposes and supplementary validation checks. While we use lava unit area and length derived from observations, photogrammetric surveys and InSAR data (Bato et al. 2016) for checking model-derived lava flow run-outs, SO2 and thermal camera data are used to check the validity of TADR obtained from MODIS data. We here review the methodologies (as well as their application in near-real time, problems and delivery delays) that formed this chain from source term provision, through model execution to output validation during the five effusive eruptions that comprised the 2015–2016 cycle of activity (Peltier et al. 2016; Coppola et al. 2017) at Piton de la Fournaise.
At the onset of an eruption, dyke location plus propagation direction and velocity is tracked on monitors in the OVPF operations room, in real-time, using the live data stream from the permanent seismic and geodetic networks. Upon eruption, vent location, eruptive fissure length and geometry, plus the number and activity of eruptive sources, are assessed using the OVPF permanent camera network, whose images are also streamed live to the operations room, and field reconnaissance. Remote surveillance is followed by in-situ inspection and/or civil protection helicopter-based over-flight including one OVPF agent. In-situ GPS-location of effusive vents, sampling of eruptive products (solids and gases), and infrared and visible image surveys are typically performed during the first few hours of an eruptive event, when weather conditions are favorable. Additional file 1 provides an overview of the frequency of: (i) sampling of solid products (pyroclasts, lavas, sublimates), (ii) in situ analyses of gas composition and fluxes, and (iii) thermal and visible camera surveys during the 2014–2015 events. A representative set of solid samples is then sent during, and/or immediately after, each eruptive event to LMV and IPGP for textural, chemical and petrological analyses.
Lava flow volume estimates are currently based on field mapping, photogrammetry and InSAR analysis (Peltier et al. 2016). Precise volume estimations are generally performed post-event or late during a long-lasting event due to (i) the cost of, and preparation time required for, satellite and/or aerial photography acquisitions, (ii) the need for Ground Control Point (GCP) measurement, and (iii) the complexities of data processing. However, due to new research funding from the Agence National de Recherche (ANR), a new satellite image-purchase program (Kalidéos 2 by CNES) and implementation of crowd sourcing techniques, volume estimation and lava flow mapping are becoming increasingly available during eruptive events and with low degrees of latency. Crowd sourcing involves use of high definition drone-derived images available on YouTube, collaboration with a network of professional photographers and journalists, and networking with drone pilots who all provide OVPF with images for photogrammetry free of charge.
TADR derivation using MIROVA
Satellite-based thermal data have been used operationally to estimate the lava discharge rates during effusive eruptions since first application in 1997 (Harris et al. 1997a). This approach relies on the observed relationship between lava discharge rate, lava flow area and thermal flux (e.g. Pieri and Baloga 1986; Wright et al. 2001; Harris and Baloga 2009, Garel et al. 2012). For any given eruptive condition, this relationship allows VRP to be set as proportional to the time-averaged lava discharge rate (TADR) using the coefficient of proportionality, crad = TADR/VRP (Coppola et al. 2013). Validation by Coppola et al. (2009; 2010; 2013) indicates that the eruptions of Piton de la Fournaise are characterized by a best-fit coefficient of between 1.4 × 108 and 2.9 × 108 J m−3. This range likely reflects variation in eruptive conditions, such as different emplacement styles (i.e., channel- versus tube-fed) or underlying topography (steep versus gentle slopes). Comparison with post-event lava flow volumes indicated that short-lived (<15 days), low-volume (< 3 × 106 m3) eruptions are best described by the upper bound of the coefficient range (Coppola et al. 2017), with the lower bounds being more applicable to eruptions lasting more than 2 weeks and emitting more than 5 × 106 m3 of lava. However, in the absence of syn-event validation, upper, median and lower bounds on MIROVA-derived TADR are given to take into account this uncertainty (Fig. 5). During the opening phases of each eruption (i.e., during the first 48 h), as well as during periods of major changes in output, TADR time-series were updated at-least four times per day and were delivered, via email, to OVPF for ingestion into the on-going hazard assessment and response. They were also sent to INGV-PI and LMV for initialization and updating of model-based lava flow run-out assessments. This service was maintained during the June 2014, May 2015, July–August 2015 and August–October 2015 eruptions.
Initialization and execution of DOWNFLOW and FLOWGO
Lava flow model product and TADR were delivered to OVPF with a delay of up to 24 h. For a lava flow front advancing at a few tens of meters per hour, and several kilometers distance from vulnerable infrastructure, then a 24 h delay may be reasonable. But for faster moving flows, closer to vulnerable sites, this may need to be reduced; and delivery delay can be reduced to an hour or so.
TADR used for initialization of FLOWGO needs to wait for the first cloud-free MODIS overpass. In the best case, this wait-time was just 23 min, but was typically 3-to-4 h.
Due to the lack of GSM (Global System for Mobile Communications) coverage over a large part of the volcano, the ability to communicate between on-site observers and the observatory is extremely limited, and messages need –literally – to be carried back by hand resulting in delivery delays for vent location if observers are out of range.
The management priorities at OVPF, where a very small team needs to deal with all scientific, media and civil defense reporting duties, mean that it may take some time to communicate vent coordinates, especially at the beginning of an eruption. For example, at the beginning of the July 2015 eruption, the director took calls from at-least seven journalists during the first 2 h of activity, while also having to organize field crews and meet civil protection call down duties. Remaining staff were spread thin keeping up with real-time geophysical and field-based surveillance duties and reporting.
However, the resulting delay in product delivery of typically 3–4 h was acceptable for the cases tracked here, where lava flows had their sources high on the volcano flanks, at-least 5 km from vulnerable infrastructure.
Real time lava flux estimation based on SO2 flux measurements
At the eruptive vent, lava effusion rate is proportional to the gas flux, provided that volatiles are dominantly released by melt upon ascent, decompression and degassing, and assuming that external sources (or magmatic mineral phases) do not comprise a significant fraction of the gas emission. At Piton de la Fournaise, SO2 flux has been demonstrated to scale linearly with effusion rate during small-to-intermediate volume eruptions (Hibert et al. 2015; Coppola et al. 2017). However, large volume, intense eruptions potentially degas a large volume of magma with respect to the volume erupted (Gouhier and Coppola 2011; Di Muro et al. 2014).
During the 2014–2015 eruptions, SO2 fluxes were quantified in real-time by OVPF’s permanent DOAS network and through completion of walked traverses along the caldera rim. The three scanning DOAS have been installed as permanent stations, these being Partage – nord, Enclos – west, and Bert – south (Galle et al. 2010). These locations are all close to the rim of the caldera, and perform continuous scanning of the sky above the Dolomieu cone during day light hours. Real-time integration of the plume cross section is performed using a set of standard and constant values for wind speed and plume height. Daily post-processing allows the spectral analysis to be refined by using the actual plume height and direction through triangulation of simultaneous scans from the three stations, and by taking into account wind speed data. During the 2014–2015 eruptions, wind speed data were provided by Meteo-France and were recorded at a station located at Bellecombe, i.e., between the DOAS stations at Partage and Enclos (see Fig. 1 for locations). Wind speed data are acquired hourly by an anemometer installed on a mast 10 m above the ground. Post-processing is carried out in collaboration with Chalmers University and allows SO2 flux to be correlated with daily rain fall data, acquired by the OVPF stations, to farther constrain the environmental effects on the gas flux estimates.
While the short-lived June 2014 and July–August 2015 eruptions occurred during good weather conditions, periods of rain interrupted the longer eruptions of May 2015 and August–October 2015. Bad weather conditions were dominant during the February 2015 eruption, making any real time assessment of gas emissions unreliable during most of the event. Precise assessment of gas fluxes is also challenging during very short-lived eruptions (lasting just a few hours), especially if a significant part of the eruption occurs at night when the UV-reliant DOAS acquisition cannot be performed. DOAS sessions are acquired at a high sampling rate (one complete sky-scan every 13 min), but still only cover one third of a day, being limited to daylight hours (which is <8 h during the winter).
Texture, geochemistry and petrology
Since 1981, lava samples have been collected during or shortly after all eruptions at Piton de la Fournaise. Efforts are made to collect early-erupted lavas and tephra, and to quench molten lava in water. For geochemical measurements, samples need to be as un-contaminated as possible. During the events considered here, 110 samples were collected during or shortly after each eruption by the OVPF team. As listed in Additional file 1, they included lavas, scoria and lapilli which were both water and air quenched.
After preliminary inspection and characterization at OVPF, a representative subset of samples are mailed to LMV for textural analysis, and to LMV and IPGP for geochemical and petrological analysis. As soon as the samples arrive they are macroscopically and microscopically described, and then each sample is divided according to the needs of textural, petrological, and geochemical measurements. Time for delivery of analyses can be up to 20 days. This is mostly due to the delay imposed by sample shipping to the mainland, where sample preparation and lab-time can be prioritised during an eruptive crisis reducing delays due to booking of the preparation and measurement facilities by other projects.
Samples for textural analysis are dried in an oven for 24 h and then used for grain size, componentry, connectivity, density, porosity and permeability measurements. Vesicle and crystal contents, as well as their size distributions, are also derived from all pyroclasts and lava samples (see Gurioli et al. (2015), Latutrie et al. (2017) and Colombier et al. (2017) for the detail regarding standard procedures, plus the meaning and application of the measurements). These measurements are performed to check variation in space (down a fissure or vent system) and in time, both within single eruptions and between different eruptions. Results are also used to check, and update if necessary, the validity of FLOWGO source terms, such as the chemistry-based rheological model, and vesicularity, density and crystal content values used by FLOWGO (Harris et al. 2016), as well as to allow SO2 emission conversions. All textural measurements are performed at the LMV textural laboratory as part of the DYNVOLC “service d’observation” (SO) or observation service (wwwobs.univ-bpclermont.fr/SO/televolc/dynvolc/).
For textural purposes, the first objective of sampling in active lava is to try to quench the sample so as to preserve the texture of the active flow and its chemistry (e.g., Cashman et al. 1999; Robert et al. 2014). The second objective is to be as representative of the flow source conditions, which for modelling means sampling as close to the vent as possible to allow source term validation, and then – if there is the luxury – to sample farther down channel so-as to provide ground-truth for the model in terms of cooling and crystallization rates (Harris et al. 2016). The third requirement is to be representative of the eruption itself. That is, to be sure that the at-vent thermal, chemical and textural conditions are not changing. To perform systematic observations, it is best to sample always at the flow front at the very beginning, during, and at the end of the eruption. In this way we always sample the same population. However, flow fronts are not always accessible, and to sample an active channel is not trivial; so that the reality of the situation is that we have to collect those samples that we can given difficult and challenging situations (e.g., www.youtube.com/embed/iwwV4hGVEcQ).
Whole rock major and trace element concentrations are analysed by ICP-OES and ICP-MS, respectively. Major element composition of minerals and glass are analysed by EPMA on a subset of quenched samples (naturally quenched in air or water quenched). For geochemical and petrological analysis, samples are cut into centimetre-sized chips, before being crushed into millimetre-sized chips using a set of thermally hardened steel jaws (which were not chemically doped). Finally the sample is powdered in a motorised agate mortar. To reduce cross-contamination, the first powder fraction is discarded, and the second and third powders are used for major and trace element/isotope analysis, respectively. Major elements are analysed by Inductively Coupled Plasma Optical Emission Spectrometry (ICP-OES, HORIBA Jobin Yvon Ultima C) following a Lithium metaborate (LiBO2) fusion method, and trace elements are analysed using a Quadrupole Mass Spectrometer (ICP-MS, Agilent 7500) following acid dissolution (HF-HNO3) of the sample in teflon vials. This method allows routine analysis of 47 elements (Li, Be, Sc, Ti, V, Cr, Co, Ni, Cu, Ga, Ge, As, Rb, Sr, Rb, Sr, Y, Zr, Nb, Cd, In, Sn, Sb, Cs, Ba, Rare Earth Elements, Hf, Ta, W, Tl, Pb, Bi, Th and U), but does not dissolve resistant minerals such as olivine-hosted Cr spinel. High-temperature (220 °C) dissolution of samples with ammonium bifluoride (NH4HF2) (Zhang et al. 2012) is currently being tested to overcome this issue. Routinely measured magma source tracers include Sr, Nd and Pb long-lived radiogenic isotopes. Strontium and Nd are purified using Eichrom specific resins (Sr.Spec and Tru.Spec) and their isotopic compositions are measured by thermal ionisation mass spectrometry (TIMS Thermo Triton). Lead is separated using Biorad AG-1X8 anionic resin and its isotopic composition is measured using a Neptune plus Multi-collector ICP-MS. A detailed description of trace and isotope analytical methods is given in Vlastelic et al. (2009).
Trace element and Pb isotopic composition of the tools used to collect and quench molten lava since 1998
In total 127 TADR values were derived and delivered to OVPF by the MIROVA system during the 2014–2015 crises, cloud cover meaning that reporting varied between zero and four TADRs per day; with an average of two per day during the August–October 2015 eruption. During the five 2014–2015 eruptions, the results of 31 sample analyses, plus all DOWNFLOW and FLOWGO model runs were delivered directly to OVPF. All measurements and products were cross-checked, and merged with data fed back into the communication loop by the monitoring system operated by OVPF. This allowed source term update, uncertainty control and validation, as well as fully constrained event tracking.
MIROVA-derived TADR and cumulative volume
Statistics for the MIROVA-derived eruption parameters for, and total SO2 mass emitted during, each eruption
No. cloud-free MODIS images
(× 106 m3)
Mean Output Rate (m3/s)
Viewing of the flow field on 23–24 October was complicated by a wild fire that was burning on the caldera wall and in the same pixel (s) as the cooling flow field. The fire was ignited by lava contact at the base of the caldera wall around 13:00 on 23 October and quickly spread up the wall, requiring rapid evacuation of footpaths and viewing points (Le Journal de L’Ile de la Réunion, 24 October 2015, p. 4). Landsat ETM+ images acquired after the event revealed a 0.2 km2 (1000 × 200 m) burn scar. However, given an estimate of 100–300 MW for the fire radiative power, the contribution to the total radiative power (2000–6000 MW) was not significant, and of the order of the uncertainty due to the contribution from cooling lava flows emplaced during the main phase of activity (400–500 MW).
Textural characterization and geochemical evolution
July 2015 lava samples were aphyric basalts that mainly contained clinopyroxene and plagioclase microphenocrysts (< 500 μm) with rare olivine microphenocrysts set in a glassy or fine-grained matrix. Lavas emitted at the beginning of the August–October 2015 eruption had the same modal composition as those of July 2015, with microphenocrysts mainly of clinopyroxene and plagioclase set in a glassy-to-microlitic groundmass. After 15 September (the date on which the porosity trend turned around from its minimum value, Fig. 13a), a change in magma composition became evident as olivine mesocrysts became more frequent, and plagioclase microphenocrysts disappeared between 15 and 27 September (an exception being the lavas of 9 October). Thus, from the end of September until the beginning of October, lavas were aphyric basalts with clinopyroxene and olivine microphenocrysts and mesocrysts. From 9 October onwards, clinopyroxene was no longer observed as microphenocrysts, and only olivine microphenocrysts and mesocrysts (with Cr-spinels in inclusion) were observed. From mid-October to the end of the eruption, lavas were olivine basalts that contained 5–10% of olivine crystals (>500 μm, up to 6 mm in size) set in a matrix containing microlites of clinopyroxene, olivine and plagioclase (+ glass).
These changes in chemistry had little effect on the modelled lava flow run-outs, with TADR being the main determinate on controlling FLOWGO-derived run-out. Likewise, textural changes were within the error of source parameters set using the 2010 eruption conditions (Harris et al. 2015). However, following the time evolution of bulk magma composition was critical for the interpretation of time evolution of TADR, which increased during the second half of August eruption, concomitantly with the emission of more magnesian basalts. What the evolution of the geochemical and textural parameters did show was that the system was evolving towards an unloading scenario to result in a terminating effusive “paroxysm” which, in hind sight, signalled the end of this particular cycle (Coppola et al. 2017).
MIROVA-derived time-averaged discharge rate (TADR) for dates on which InSAR data are available with FLOWGO lava flow run-outs that each TADR gives and the InSAR-derived lava flow length for the same day. Δ gives the difference between the InSAR-derived flow length and FLOWGO run-out
InSAR-derived flow length (km)
29 August 2015
8.1 ± 2.8
3.1 ± 0.4
3.9 ± 0.2
0.8 ± 0.6
14 September 2015
6.8 ± 2.3
3.0 ± 0.1
3.7 ± 0.1
0.7 ± 0.2
30 September 2015
7.1 ± 2.5
3.0 ± 0.2
3.7 ± 0.1
0.7 ± 0.3
04 October 2015
10.5 ± 3.6
3.5 ± 0.5
3.6 ± 0.1
0.1 ± 0.6
08 October 2015
9.1 ± 3.1
3.1 ± 0.6
3.7 ± 0.1
0.6 ± 0.7
13 October 2015
10.5 ± 3.6
3.5 ± 0.5
3.6 ± 0.1
0.1 ± 0.6
17 October 2015
24.8 ± 8.5
9.7 ± 1.6
5.6 ± 0.1
4.1 ± 1.7
24 October 2015
22.2 ± 3.8
8.7 ± 1.2
5.7 ± 0.5
3.0 ± 1.7
01 November 2015
0.6 ± 0.2
0.4 ± 0.1
5.7 ± 0.5
5.3 ± 0.6
Measurements to support validation of MIROVA-derived TADR were also made during the July–August 2015 eruption. Observations were made of the main active vent and its outlet channel between 10:00 and 15:00 on 1 August 2015. In addition, thermal video was taken at the head of the master channel where it exited the eruptive fissure for 5 min at 11:40 on 1 August, and a water-quenched sample was collected from the same location at 13:15. The channel was 2 m wide and contained a 2 m deep flow with its surface 1 m below the levee rim. Velocities obtained from the thermal video were 0.05–0.1 m/s, for an effusion rate of 0.2–0.4 m3/s. This agrees with MIROVA-derived TADRs of 0.8 ± 0.4 m3/s obtained from the two evening MODIS overpasses at 19:25 and 20:55 on 1 August, and with the value of 0.18 ± 0.6 m3/s obtained on 2 August at 06:20, 5 h before the eruption ended. Because our observations were made during a period of waning activity, where flow levels and velocities in the master channel underwent a noticeable decline after 11:30 on 1 August, we used Jeffreys (1925) to estimate flow velocities under peak-flow conditions. Viscosity was calculated on the basis of thermal-camera-derived flow temperature (1114 °C), plus sample crystal (19–20%) and/or vesicle (50–58%) content, using Villeneuve et al. (2008), Roscoe (1952), Manga and Loewenberg (2001), Pal (2003) and Llewellin and Manga (2005). Results were in the range 370–700 Pa s. Using this viscosity range in Jeffreys (1925), with the sample density of 1510 kg/m3 and underlying slope of 5°, yields a peak-flow velocity of 0.9–1.7 m/s, for an effusion rate of 5.35 ± 1.65 m3/s. This matches the MIROVA-derived TADR of 6.25 ± 2.25 m3/s for 10:00 on 1 August. Aerial photographs taken by journalists during the opening day of the eruption (H. Douris, Le Journal de L’Ile de la Réunion, 02/08/15, p. 7) reveals the master channel to have been brim full (i.e., flow depth = 3 m) at that point. This higher flow level yields velocities of 2.1–3.9 m/s, which convert to effusion rates of 12 ± 3.7 m3/s. The same photographs indicated that two channel systems of similar dimensions were active during 31 July, so the total effusion rate could have been as high as 24 ± 8 m3/s during the opening hours of the eruption. These values again compare well with those derived from MIROVA which gave 22 ± 8 m3/s from MODIS images acquired at 18:40 and 21:50 on 31 July.
The response: An example from the may 2015 eruption
The May 2015 eruption began at 13:45 (local time) on 17 May at three en-echelon fissures. The shortest (30 m long) and most western fissure was located at 2285 m asl. The second fissure was located 200 m farther east, between 2250 and 2100 m asl, and was 500 m long. The third, and most eastern, fissure was located a farther 1100 m downslope. Located between 2060 and 1980 m asl it was 360 m long. It is the lava flow that spread from this third fissure that is the subject of this case-study.
On 18 May around 08:00, the flows entered a zone of vegetation that extended between 1450 and 1150 m asl covering an area of 42,500 m2 at a distance between 5 and 4.3 km west of the road. At this location, the flow front appeared to slow somewhat (despite the mean slope increasing to 40%). This was probably due to the effect of the vegetation. In addition, average SO2 fluxes quickly decreased from 2700 tons/day on 18 May to 390 tons/day on 19 May, decreasing more progressively thereafter until 23 May. After 23 May SO2 fluxes were <80 tons/day until the end of the eruption. By 21 May, TADRs had also dropped to 2.5 m3/s (Fig. 5) so that predicted run-outs reduced to 2 km, 6 km above the road. The eruption ended at 20:50 on 31 May, with the longest-reached flows being 4.05 km long, having extended to within 3.9 km of the road.
We have reviewed a response model for a crisis at an island volcano that is distant from its administrative center. The model involves synergy between multiple distant nodes so as to create an extended, virtual observatory. At the onset of each eruption the system was triggered by an email from the OVPF director, declaring the date and time of eruption. This was distributed to the email distribution list linking the five institutes involved in this exercise (IPGP, OPGC, LMV, Università di Torino, INGV, Chalmers). Shortly thereafter, the coordinates of the vent location were emailed to all partners and DOWNFLOW launched at INGV-Pisa. The output line of steepest descent was then handed on to LMV for initialization of FLOWGO. FLOWGO was then run, and output delivered to OVPF, as soon as the first MIROVA-derived TADR arrived (again via email from the Università di Torino to the email distribution list). Initially FLOWGO was run using the chemical model validated for the 2010 eruption, and the magma composition, crystal and vesicle properties of the same eruption. These source terms were then updated upon receipt of first chemical and textural analyses.
Source terms, ground truth and products, with source, frequency of provision and delivery delay, developed as part the effusive crisis response network of Fig. 13
Frequency of provision
In-situ sampling by OVPF staff
(delivery time to mainland)
(flow length, area)
Field measurements by OVPF staff
Scanning DOAS system
Once every 30 min
Vent location & DEM
GPS & drone-based photogrammetry (OVPF)
Vent location = hours
DEM = days
Chemistry, crystallinity & vesicularity
Source Term & Product
Geochemistry and textural laboratory (OPGC/LMV)
Same as sampling
Source Term & Product
(Università di Torino)
Four times per day
Source Term & Product
On-demand as vent locations & DEM changes
Four times per day (as TADR-arrives)
On the basis of working together during these five effusive crises we can identify five crucial components to ensure smooth information flow: (1) a need for TADR validation; (2) exclusion of fires from satellite signals to isolate the volcanic component; (3) an understanding of lava – vegetation interactions; (4) timely provision and validation of model source terms, and (5) a clear statement of, and if possible a reduction in, uncertainty (including cleaning of data sets of unreliable or untrustworthy data points). One unexpected uncertainty which we encountered was due to sample contamination by the sampling device. Chemical and textural samples critical for source term checking need to be uncontaminated, so we need to find a way to sample without (i) introducing chemical artifacts due to the sampling device, (ii) changing the vesicle structure due to shearing on withdrawal of the sampling device from the fluid, or (iii) from vesiculation during quenching and boiling in water. A realistic sampling protocol also needs to be defined, with no redundancy. That is, if a small, over-worked field- crew charged with monitoring an eruption several kilometers from the nearest vehicle access is to be efficient, then the minimum number of most useful measurements need to be made in a limited amount of time at the most viable (accessible and safe) sites. That is: a realistic and viable sampling protocol needs to be put in place before the crisis.
The best approach to TADR validation is an ensemble-based one that looks for consistency across approaches. In our case the test was:
satellite-based and ground SO2- and/or IR-derived TADR are in agreement,
when fed into an appropriately initialized and tested flow model the simulation runs out to field-observed flow lengths
the TADR can be trusted.
This is an important need because TADR is intimately related to flow length and area (e.g., Walker 1973; Malin 1980; Pinkerton and Wilson 1994; Calvari and Pinkerton 1998; Murray and Stevens 2000; Harris and Rowland 2009), and hence defines the potential for a lava flow to enter vulnerable areas to cause damage. Collection of field based TADRs are viable, and can be obtained on the basis of thermal emission or flow dynamics from IR cameras. The thermal approach requires a flow-wide thermal mosaic, which requires a platform, such as ultra-light aircraft, helicopter, or drone which may be beyond capacity in terms of funds, weather or technology. The ground-based approach requires set-up of a camera within sight of the master channel to acquire flow dimensions and velocity. Although removing the requirement of a platform, it still means that viewing conditions and geometries have to be opportune, and that a field crew can reach the observation site. At this point, satellite-based methods (validated by occasional field spot checks) appear best suited for rapid response, although the satellite-based measurement has to consistently cross-check with ground truth if it is to be trusted. Inter-validation between satellite and SO2-flux derived TADR is valid only when the pre-eruptive sulfur content is constant. If SO2-derived TADR does not agree with the satellite-based values, but when the latter values provide model-based flow lengths consistent with observed lengths, then we may be able to infer the sulfur output of the erupted lavas.
Problems with fires
If, as is likely on a heavily vegetated, tropical island volcano, lava ingress into vegetation ignites a fire, then the fire has to be removed from the heat budget if the TADR is not to be over-estimated using the thermal approach. This is an old problem, but solutions are possible using a combination of spectral and field observations, as was the case here. Indeed, separating the intensity, location and spreading direction of a fire from that of the lava is an extremely useful exercise as the fire poses a hazard in its own right. A-priori knowledge of the pre-burn fuel load (kg m−2) will help to estimate the “fire radiative power” (e.g., Van Wagner 1967; Viskanta and Mengüç 1987; Mell et al. 2009), once the fire-affected area is estimated. A map of the potential biomass energy released per unit mass of fully burnt, dry fuel, may thus be useful to correct the radiant power measured by satellite if lava ignites a fire.
Trees may, or may not, affect lava spreading in terms of both a thermal and mechanical effect. Thermally, a heavily vegetated zone may cause flows to excessively cool and crystalize due to the need to dry and then ignite trees (Van Wagner 1967). Open tree molds may then serve as skylights that allow heat to radiatively escape from the flow interior as is the case for lava tubes (Witter and Harris 2007). The solidification of lava around trees, as well as the trees themselves, then cause mechanical obstructions. Both effects may serve to slow flow advance, but the problem is totally unconstrained.
Provision and validation of model source terms
Any TADR-derivation model or lava flow simulation model requires input of source terms and validation of output. Crucially this requires data for: chemistry (for the flow viscosity model); eruption temperature (for the flow and TADR model); SO2 gas flux (for TADR validation); flow crystallinity and vesicularity (for the flow and TADR model); plus vent location and up-to-date DEM, with a horizontal and vertical resolution of less than one meter, for flow direction runs. We find use of InSAR data to be extremely promising in this regard.
The response model described here integrates external partners who need to enter the communication network seamlessly, and provide product that is: (i) trusted and validated, (ii) in a format that is immediately useable, (iii) useful for monitoring and execution of assessment and reporting duties, and (iv) trusted. Within this network, the operation of any methodology, and sources of uncertainty, need to be well-known and spurious (or un-necessary) information removed. Transparency, efficiency and full documentation is thus key. Raw data will not be used, neither will product which is difficult to interpret or whose source is unknown. What are needed are answers, where the observatory will have defined the questions; and it will be up to the partners to iterate on their answers until product is seamlessly integrated into the workflow.
A good state-of-the-art example of a similar, but internally-developed, model is that of the HOTSAT system at INGV’s Osservatorio Etnea (Ganci et al. 2016). HOTSAT is designed to allow near-real time assessment and simulation of effusive crises at Mount Etna (Sicily, Italy). The system uses satellite (SEVIRI and MODIS) data which are fed into an in-house physics-based lava flow propagation model (MAGFLOW), and is based on detection algorithms (Ganci et al. 2012a, 2012b), TADR conversions (Ganci et al. 2012a, 2012b), and lava flow simulations (Del Negro et al. 2008; Hérault et al. 2009) that have been developed, tested and validated in-house (Vicari et al. 2009), before being launched operationally to allow improved crisis response (Vicari et al. 2011a). The workflow is almost identical to that applied here, where (i) hot pixels are located using satellite data, (ii) the flow model is initialized with a DEM and appropriate chemical and physical volcanological parameters, (iii) the model is executed with vent coordinates and satellite-derived TADR, allowing (iv) flow coverage assessments to be delivered for observatory-based hazard assessment purposes (Vicari et al. 2011b). In the case of HOTSAT, the workflow is completed by an observatory-based remote sensing and modeling group which has five members: the same number as the total staff at OVPF charged with monitoring. In our case, the group currently working on the same task (TADR derivation and model execution) does actually number four, but that grouping comes from four institutions in two countries.
Periods of major unrest or high magnitude (explosive) eruptive events have, in the past, prompted assembly of multidisciplinary response teams to support small groups charged with tracking the event. A recent example is the interaction of the UK’s Meteorlogical Office and British Geological Survey with the Icelandic Meteorological Office and the Institute of Earth Sciences at the University of Iceland during the 2010 eruption of Eyjafjällajökull (Donovan and Oppenheimer 2012). Such groupings have also been developed, for example, during the eruption of Nevado del Ruiz in 1985 (Hall 1990), at Nyiragongo during and following the 2002 eruption (Allard et al. 2002; Tedesco et al. 2007), and for Merapi in 2010 (Jousset et al. 2012). Another example is the response to eruptive events provided by the USGS Volcano Disaster Assistance Program, which has collaborated with volcano observatories in 12 countries in connection with at least 30 eruptive crises since 1985 (https://volcanoes.usgs.gov/vdap/). In these cases, large international collaborations were configured to assess volcanic activity and their impacts during and following a single event. The case presented in our study differs from these examples in that (i) the collaboration was coordinated from the local volcano observatory, instead of by an external, international (e.g., USGS) or transnational (e.g., UN, VAAC) organization, (ii) our response model is permanent (that is, we activate the model each time there is an eruption), and (iii) the response targets frequent, but low-magnitude eruptive events.
The cases tested here have, instead, been typical, low magnitude effusive events at Piton de la Fournaise with limited (in a geographical sense) impact, with lava flow activity being confined to the Enclos Fouqué caldera. In this case the main causes for concern are (i) burial of the island belt road (RN2), (ii) destruction of observatory instrumentation, (iii) evacuation of the Enclos Fouqué caldera, (iv) security of the footpaths, (v) injury to (and evacuation of) tourists, (vi) fatalities among tourists entering the closed zone, and (vii) forest fires ignited by the active lava. For an eruption in inhabited areas outside of the Enclos Fouqué caldera, including densely populated areas such as around Le Tampon (Fig. 1a), the impact and response model would have to be scaled up, other external participants called, and the component models adjusted to suit, and tested and validated on, the new case. We are currently preparing for such an eventuality through an initiative entitled ANR-LAVA (Lava Advance in Vulnerable Areas). This initiative, funded by the French ANR (Agence National de Recherche), supports the group to develop and test the response model and its component parts, including the simulation model which will be based on that of Bernabeu et al. (2016), for effusive events that enter heavily vegetated and/or populated areas. In such a sensitive case there is even less room for error or mis-communication of uncertainty.
In France, volcano observatories are dedicated to observations and measurements, plus recording, archiving, communicating and distributing data. During an effusive crisis the observatory director needs to provide local civil protection, and therefore the local municipality, with factual elements that are often quantitative – but always based on trusted measurements. For example, the director needs to support statements such as: “as of 09:45 this-morning discharge rate and flow length was increasing, and the flow front was 5 km above the road”. To answer legitimate questions regarding risk, the director will give responses only based on statistics that have been validated and published with appropriate error bars along with providing possible scenarios based on knowledge of past activity at the volcano. A small staff thus needs all the measurement and model based support that it can obtain, as well as base-line data, in order to support such communications. That support needs to be trusted, timely and, above all, validated.
We thank two anonymous reviews for their input, and the patience of the Applied Volcanology Editors Office for allowing us two extensions as we were distracted by two eruptive crisis during the submission (31 January – 5 February 2017) and correction (18-19 May 2017) of this manuscript, respectively. We thank the STRAP project funded by the Agence Nationale de la Recherche (ANR-14-CE03-0004-04) and the OMNCG OSU program from la Reunion University for supporting DOAS data acquisition and processing. This work was funded by the Agence National de la Recherche (ANR: www.agence-nationale-recherche.fr) through project ANR-LAVA (ANR Program: DS0902 2016; Project: ANR-16 CE39-0009, PI: A. Harris, Link: www.agence-nationale-recherche.fr/?Project=ANR-16-CE39-0009). This is ANR-LAVA contribution no. 1.
All authors read and approved the final manuscript.
The authors declare that they have no competing interests
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
- Allard P, Baxter P, Halbwachs M, Komorowski J-C. Final report of the French-British scientific team: submitted to the Ministry for Foreign Affairs, Paris, France, foreign office, London, United Kingdom and respective embassies in Democratic Republic of Congo and Republic of Rwanda; 2002. p. 24.Google Scholar
- Bardolet E, Sheldon PJ. Tourism in archipelagos: Hawai'i and the Balearics. Ann Tour Res. 2008;35(4):900–23.View ArticleGoogle Scholar
- Bato M, Froger J, Harris A, Villeneuve N. Monitoring an effusive eruption at Piton de la Fournaise using radar and thermal infrared remote sensing data: insights into the October 2010 eruption and its lava flows. In: Harris AJL, De Groeve T, Garel F, Carn SA, editors. Detecting, Modelling and responding to effusive eruptions, vol. 426. London: Geological Society, London, Special Publications; 2016. p. 533–52.Google Scholar
- Bello A. Valorisation des eruptions du Piton de la Fournaise, Ile de la Réunion. Report to the Parc National de La Réunion (Secteur Est); 2010. p. 25.Google Scholar
- Bernabeu N, Saramito C, Smutek C. Modelling lava flow advance using a shallow-depth approximation for three-dimensional cooling of viscoplastic flows. In: Harris AJL, De Groeve T, Garel F, Carn SA, editors. Detecting, Modelling and responding to effusive eruptions, vol. 426. London: Geological Society, London, Special Publications; 2016. p. 409–23.Google Scholar
- Boivin P, Bachèlery P. Petrology of 1977 to 1998 eruptions of Piton de la Fournaise, la Réunion Island. J Volcanol Geotherm Res. 2009;184:109–5.View ArticleGoogle Scholar
- Breaker LC. Estimating and removing sensor-induced correlation from advanced very high resolution radiometer satellite data. J Geophys Res. 1990;95(C6):9601–711.View ArticleGoogle Scholar
- Brown L. Birth of a mountain: Montserrat’s volcano – an eyewitness account. Blackwater: Sargeant Press; 2010. p. 278.Google Scholar
- Calder E, Harris A, Peña P, Pilger E, Flynn L, Fuentealba G, et al. Combined thermal and seismic analysis of the Villarrica volcano lava lake, Chile. Rev Geol Chile. 2004;31(2):259–72.Google Scholar
- Calvari S, Pinkerton H. Formation of lava tubes and extensive flow field during the 1991–1993 eruption of Mount Etna. J Geophys Res. 1998;B103:27291–301.View ArticleGoogle Scholar
- Cashman KV, Thornber C, Kauahikaua. Cooling and crystallization of lava in open channels, and the transition of Pāhoehoe lavat o ‘A’ā. Bull Volcanol. 1999;61:306–23.View ArticleGoogle Scholar
- Chester DK, Duncan AM, Dibben CJL. The importance of religion in shaping volcanic risk perception in Italy, with special reference to Vesuvius and Etna. J Volcanol Geotherm Res. 2008;172:216–28.View ArticleGoogle Scholar
- Colombier M, Gurioli L, Druitt TH, Shea T, Boivin P, Miallier D, et al. Textural evolution ofmagma during the 9.4-ka trachytic explosive eruption at Kilian volcano, Chaîne des Puys, France. Bull Volcanol. 2017;79:17. doi:10.1007/s00445-017-1099-7.
- Coppola D, Piscopo D, Staudacher T, Cigolini C. Lava discharge rate and effusive pattern at Piton de la Fournaise from MODIS data. J Volcanol Geotherm Res. 2009;184(1–2):174–92. doi:10.1016/j.jvolgeores.2008.11.031.View ArticleGoogle Scholar
- Coppola D, James MR, Staudacher T, Cigolini C. A comparison of field- and satellite-derived thermal flux at Piton de la Fournaise: implications for the calculation of lava discharge rate. Bull Volcanol. 2010;72(3):341–56. doi:10.1007/s00445-009-0320-8.View ArticleGoogle Scholar
- Coppola D, Laiolo M, Piscopo D, Cigolini C. Rheological control on the radiant density of active lava flows and domes. J Volcanol Geotherm Res. 2013;249:39–48. doi:10.1016/j.jvolgeores.2012.09.005.View ArticleGoogle Scholar
- Coppola D, Laiolo M, Cigolini C, Delle Donne D, Ripepe M. Enhanced volcanic hot-spot detection using MODIS IR data: results from the MIROVA system. In: Harris AJL, De Groeve T, Garel F, Carn SA, editors. Detecting, Modelling and responding to effusive eruptions, vol. 426. London: Geological Society, London, Special Publications; 2016. p. 181–205. First published online 14 May 2015. http://doi.org/10.1144/SP426.5.Google Scholar
- Coppola D, Di Muro A, Peltier A, Villeneuve N, Ferrazzini V, Favalli M, et al. Shallow system rejuvenation and magma discharge trends at Piton de la Fournaise volcano (la Réunion Island). Earth Planet Sci Lett. 2017;463:13–24.Google Scholar
- Dehn J, Dean KG, Engle K, Izbekov P. Thermal precursors in satellite images of the 1999 eruption of Shishaldin volcano. Bull Volcanol. 2002;64:525–45.View ArticleGoogle Scholar
- Del Negro C, Fortuna L, Herault A, Vicari A. Simulations of the 2004 lava flow at Etna volcano by the MAGFLOWcellular automata model. Bull Volcanol. 2008;70:805–12. http://doi.org/10.1007/s00445-007-0168-8 View ArticleGoogle Scholar
- Di Muro A, Métrich N, Vergani D, Rosi M, Armienti P, Fougeroux T, et al. The shallow plumbing system of Piton de la Fournaise volcano (la Réunion Island, Indian Ocean) revealed by the major 2007 caldera forming eruption. J Petrol. 2014;55:1287–315.Google Scholar
- Di Muro A, Métrich N, Allard P, Aiuppa A, Burton M, Galle B, et al. Magma degassing at Piton de la Fournaise volcano. In: Bachelery P, Lenat J-F, Di Muro A, Michon L, editors. Active volcanoes of the Southwest Indian Ocean. Berlin: Springer; 2016. p. 203–22.Google Scholar
- Dietterich HR, Poland MP, Schmidt DA, Cashman KV, Sherrod DR, Espinosa AT. Tracking lava flow emplacement on the east rift zone of Kīlauea, Hawai‘i, with synthetic aperture radar coherence, Geochem. Geophys Geosyst. 2012;13:Q05001. doi:10.1029/2011GC004016.
- Dominey-Howes D, Minos-Minopoulos D. Perceptions of hazard and risk on Santorini. J Volcanol Geotherm Res. 2004;137:285–310.View ArticleGoogle Scholar
- Donovan A, Oppenheimer C. Governing the lithosphere: Insights from Eyjafjallajökull concerning the role of scientists in supporting decision-making on active volcanoes. J Geophys Res. 2012;117(B03214). doi:10.1029/2011JB009080.
- Etcheverria O. Du vignoble à la destination oenotouristique. L’exemple de l’Ile de Santorin. Cultur – Revista de Cultura e Turismo. 2014;8(3):188–210.Google Scholar
- Favalli M, Pareschi MT, Neri A, Isola I. Forecasting lava flow paths by a stochastic approach. Geophys Res Lett. 2005;32(L03305). doi: 10.1029/2004GL021718.
- Frulla LA, Milovich JA, Gagliardini DA. Illumination and observation geometry for NOAA-AVHRR images. Int J Remote Sens. 1995;16(12):2233–53.View ArticleGoogle Scholar
- Galle B, Johansson M, Rivera C, Zhang Y, Kihlman M, Kern C, Lehmann T, Platt U, Arellano S, Hidalgo S. Network for Observation of Volcanic and Atmospheric Change (NOVAC)—A global network for volcanic gas monitoring: Network layout and instrument description. J Geophys Res. 2010;115(D05304). doi:10.1029/2009JD011823.
- Ganci G, Vicari A, Cappello A, Del Negro C. An emergent strategy for volcano hazard assessment: from thermal satellite monitoring to lava flow modeling, remote Sens. Environment. 2012a;119:197–207. doi:10.1016/j.rse.2011.12.021.Google Scholar
- Ganci G, Harris AJL, Del Negro C, Guehenneux Y, Cappello A, Labazuy P, et al. A year of lava fountaining at Etna: volumes from SEVIRI. Geophys Res Lett. 2012b;39:L06305. doi:10.1029/2012GL051026.
- Ganci G, Bilotta G, Capello A, Herault A, Del Negro C. HOTSAT: a multiplatform system for the thermal monitoring of volcanic activity using satellite data. In: Harris AJL, De Groeve T, Garel F, Carn SA, editors. Detecting, Modelling and responding to effusive eruptions, vol. 426. London: Geological Society, London, Special Publications; 2016. p. 207–21.Google Scholar
- Garau-Vadell JB, Díaz-Armas R, Gutierrez-Taño D. Residents’ perceptions of tourism impacts in island destinations: a comparative analysis. Int J Tour Res. 2014;16:578–85.View ArticleGoogle Scholar
- Garel F, Kaminski E, Tait S, Limare A. An experimental study of the surface thermal signature of hot subaerial isoviscous gravity currents: implications for thermal monitoring of lava flows and domes. J Geophys Res. 2012;117(B02205). http://dx.doi.org/10.1029/2011JB008698.
- Gaudru H. Case study 1: Reunion Island, France – Piton de la Fournaise volcano. In: Erfurt-Cooper P, Cooper M, editors. Volcano and geothermal tourism. London: earthscan; 2010. p. 54–5.Google Scholar
- Germanaz C. Du pont des navires au bord des cratères: regards croisés sur le Piton de la Fournaise (1653–1964). Itinéraires iconographiques et essai d’iconologie du volcan actif de La Réunion. Paris: Université Paris-Sorbonne, thèse de doctorat; 2005.Google Scholar
- Germanaz C. Le haut lieu touristique comme objet spatial linéaire : le somin Volcan (île de Réunion) Fabrication, banalisation et patrimonialisation. Cahiers de géographie du Québec. 2013;V57(N162):379–405. doi:10.7202/1026525ar.View ArticleGoogle Scholar
- Gouhier M, Coppola D. Satellite-based evidence for a large hydrothermal system at Piton de la Fournaise volcano. Geophys Res Lett. 2011;38(2):L02302. doi:10.1029/2010GL046183.View ArticleGoogle Scholar
- Gouhier M, Guéhenneux Y, Labazuy P, Cacault P, Decriem J, Rivet S. HOTVOLC: a web-based monitoring system for volcanic hot spots. In: Harris AJL, De Groeve T, Garel F, Carn SA, editors. Detecting, modelling and responding to effusive eruptions, vol. 426, The Geological Society of London. London: Geological Society, London, Special Publications; 2016. p. 223–42. doi:10.1144/SP426.31.Google Scholar
- Goward SN, Markham B, Dye DG, Dulaney W, Yang J. Normalized difference vegetation index measurements from the advanced very high resolution radiometer. Remote Sens Environ. 1991;35:257–77.View ArticleGoogle Scholar
- Gurioli L, Andronico D, Bachèlery P, Balcone-Boissard H, Battaglia J, Boudon G, et al. MeMoVolc consensual document: a review of cross-disciplinary approaches to characterizing small explosive magmatic eruptions. Bull Volcanol. 2015;77:49. doi:10.1007/s00445-015-0935-x.
- Hall M. Chronology of the principal scientific and governmental actions leading up to the November 13, 1985 eruption of Nevado del Ruiz, Colombia. J Volcanol Geotherm Res. 1990;42(1):101–15. doi:10.1016/0377-0273(90)90072-N.View ArticleGoogle Scholar
- Harris AJL, Baloga S. Lava discharge rates from satellite-measured heat flux. Geophys Res Lett. 2009;36(L19302). doi:10.1029/2009GL039717.
- Harris AJL, Rowland SK. FLOWGO: a kinematic thermo-rheological model for lava flowing in a channel. Bull Volcanol. 2001;63:20–44. doi:10.1007/s004450000120.View ArticleGoogle Scholar
- Harris AJL, Rowland SK. Effusion rate controls on lava flow length and the role of heat loss: a review. In: Thordarson T, Self S, Larsen G, Rowland S K & Hoskuldsson A, editors. Studies in Volcanology: The Legacy of George Walker. Special Publications of IAVCEI 2; 2009. p. 33–51.Google Scholar
- Harris AJL, Keszthelyi L, Flynn LP, Mouginis-Mark PJ, Thornber C, Kauahikaua J, et al. Chronology of the episode 54 eruption at Kilauea volcano, Hawaii, from GOES-9 satellite data. Geophys Res Lett. 1997a;24(24):3281–4.Google Scholar
- Harris AJL, Blake S, Rothery DA, Stevens NF. A chronology of the 1991 to 1993 Etna eruption using AVHRR data: implications for real time thermal volcano monitoring. J Geophys Res. 1997b;102(B4):7985–8003.View ArticleGoogle Scholar
- Harris AJL, Wright R, Flynn LP. Remote monitoring of Mount Erebus Volcano, Antarctica, using polar orbiters: progress and prospects. Int J Remote Sens. 1999;20(15&16):3051–71.View ArticleGoogle Scholar
- Harris AJL, Murray JB, Aries SE, Davies MA, Flynn LP, Wooster MJ, et al. Effusion rate trends at Etna and Krafla and their implications for eruptive mechanisms. J Volcanol Geotherm Res. 2000;102(3–4):237–69.Google Scholar
- Harris AJL, Pilger E, Flynn LP, Garbeil H, Mouginis-Mark PJ, Kauahikaua J, et al. Automated, high temporal resolution, thermal analysis of Kilauea volcano, Hawaii, using GOES-9 satellite data. Int J Remote Sens. 2001;22(6):945–67.Google Scholar
- Harris A, Rhéty M, Gurioli L, Villeneuve N, Paris R. Simulating the thermorheological evolution of channel-contained lava: FLOWGO and its implementation in EXCEL. In: Harris AJL, De Groeve T, Garel F, Carn SA, editors. Detecting, modelling and responding to effusive eruptions, vol. 426. London: Geological Society, London Special Publication; 2016. p. 313–36. doi:10.1144/SP426.9.Google Scholar
- Heggie TW. Reported fatal and non-fatal incidents involving tourists in Hawaii volcanoes National Park, 1992–2002. Travel med Infect Dis. 2005;3:123–31.View ArticleGoogle Scholar
- Heggie TW. Search and Rescue in Alaska’s National Parks. Travel med Infect Dis. 2008;6:355–61.View ArticleGoogle Scholar
- Heggie TW. Death by volcanic laze. Wilderness Environ Med. 2009;20(1):101–3.View ArticleGoogle Scholar
- Heggie TW, Heggie TM. Viewing lava safely: an epidemiology of hiker injury and illness in Hawaii volcanoes national park. Wilderness Environ Med. 2004;15:77–81.View ArticleGoogle Scholar
- Heggie TW, Heggie TM. Search and rescue trends and the emergency medical service workload in Utah’s national parks. Wilderness Environ Med. 2008;19:164–71.View ArticleGoogle Scholar
- Heggie TW, Heggie TM. Search and rescue trends associated with recreational travel in US national parks. J Transl Med. 2009;16(1):23–7.Google Scholar
- Hérault A, Vicari A, Ciraudo A, Del Negro C. Forecasting lava flow hazards during the 2006 Etna eruption: using the MAGFLOW cellular automata model. Comput Geosci. 2009;35:1050–60. http://doi.org/10.1016/j.cageo.2007.10.008 View ArticleGoogle Scholar
- Hibert C, Mangeney A, Polacci M, Di Muro A, Vergniolle S, Ferrazzini V, Taisne B, Burton M, Dewez T, Grandjean G, Dupont A, Staudacher T, Brenguier F, Shapiro N.M, Kowalski P, Boissier P, Catherine P, Lauret F (2015) Multidisciplinary monitoring of the January 2010 eruption of Piton de la Fournaise volcano, la Réunion island: J Geophys Res 120: 3026. doi:10.1002/2014JB011769.
- Holben B, Fraser RS. Red and near-infrared sensor response to off-nadir viewing. Int J Remote Sens. 1984;5(1):145–60.View ArticleGoogle Scholar
- INSEE. Réunion: Fréquentation touristique 2015 – La frequentation touristique repart à la hausse. La Réunion, St. Denis), No. 16 (Mai 2016): INSEE Analyses; 2016. p. 4.Google Scholar
- Jeffreys H. The flow of water in an inclined channel of rectangular section. Philos Mag. 1925;49:793–807.View ArticleGoogle Scholar
- Josseran L, Paquet C, Zehgnoun A, Caillere N, Le Tertre A, Solet J-L, et al. Chikungunya disease outrbreak, Reunion Island. Emerg Infect Dis. 2006;12(12):1994–5. doi:10.3201/eid1212.060710.
- Jousset S, Pallister J, Boichu M, Buongiorno M, Budisantoso A, Costa F, et al. The 2010 explosive eruption of Java's Merapi volcano—a ‘100-year’ event. J Volcanol Geotherm Res. 2012;241(242):121–35. doi:10.1016/j.jvolgeores.2012.06.018.
- Kersten O. Baron Carl Claus von der Decken's Reisen in Ost-Afrika in den jahren 1859 bis 1865 (Vol. 2). Ed Leipzig/Heidelberg: Winterliche Verlagshandlung; 1871.Google Scholar
- Kersten O, Tolède M, Fois-Kaschel G (2016). Les Voyages en Afrique orientale du baron Carl Claus von der Decken. : La Réunion (28 mai - 7 août 1863). Les éditions de Villèle; Cercle des Muséophiles de Villèle, 108 p. DOI: 978-2-905861-28-3. <hal-01367548>.Google Scholar
- Kilburn CRJ, Lopes RMC. The growth of aa lava fields on Mount Etna, Sicily. J Geophys Res. 1988;93:14759–72.View ArticleGoogle Scholar
- L. R. Une trentaine d'interventions pour le SDIS et le PGHM. Le Journal de la Ile de La Réunion: Monday 3 August 2015. 2015;21 364:7.Google Scholar
- Latutrie B, Harris A, Médard E, Gurioli L. Eruption and emplacement dynamics of a thick trachytic lava flow of the Sancy volcano (France). Bull Volcanol. 2017;79:4. doi:10.1007/s00445-016-1084-6.View ArticleGoogle Scholar
- Lillesand TM, Kiefer RW. Remote sensing and image interpretation. New York: Wiley; 1987. p. 24–7.Google Scholar
- Llewellin EW, Manga M. Bubble suspension rheology and implications for conduit flow. J Volcanol Geotherm Res. 2005;143:205–17. doi:10.1016/j.jvolgeores.2004.09.018.View ArticleGoogle Scholar
- Malin MC. Lengths of Hawaiian lava flows. Geology. 1980;8:306–8.View ArticleGoogle Scholar
- Manga M, Loewenberg M. Viscosity of magmas containing highly deformable bubbles. J Volcanol Geotherm Res. 2001;105:19–24. doi:10.1016/S0377-0273(00)00239-0.View ArticleGoogle Scholar
- Markham BL. The Landsat sensors’ spatial response. IEEE Trans Geosci Remote Sens. 1985;GE-23:864–75.View ArticleGoogle Scholar
- Mell W, Maranghides A, McDermott R, Manzello SL. Numerical simulation and experiments of burning douglas fir trees. Combustion and Flame. 2009;156:2023–41.View ArticleGoogle Scholar
- Michelin Green Guide. La Réunion. Boulonge Billancourt: Le Guide Verte, Michelin Propriétaires-éditeurs; 2015. p. 258.Google Scholar
- Michon L, Di Muro A, Villeneuve N, Saint-Marc C, Fadda P, Manta F. Explosive activity of the summit cone of Piton de la Fournaise volcano (la Réunion island): a historical and geological review. J Volcanol Geotherm res. 2013;263:117–33.View ArticleGoogle Scholar
- Mouginis-Mark PJ, Garbeil H, Flament P. Effects of viewing geometry on AVHRR observations of volcanic thermal anomalies. Remote Sens Environ. 1994;48:51–60.View ArticleGoogle Scholar
- Mouginis-Mark PJ, Snell H, Ellisor R. GOES satellite and field observations of the 1998 eruption of Volcan Cerro Azul, Galápagos. Bull Volcanol. 2000;62:188–98.View ArticleGoogle Scholar
- Murray JB, Stevens N. New formulae for estimating lava flow volumes at Mt. Etna volcano, Sicily. Bull Volcanol. 2000;61:515–26.View ArticleGoogle Scholar
- Nave R, Ricci T, Pacilli MG. Perception of risk for volcanic hazard in Indian Ocean: la Réunion Island case study. In: Bachelery P, Lenat J-F, Di Muro A, Michon L, editors. Active volcanoes of the Southwest Indian Ocean. Berlin: Springer; 2016. p. 315–26.View ArticleGoogle Scholar
- Pal R. Rheological behavior of bubble-bearing magmas, Earth Planet. Sci Lett. 2003;207:165–79. doi:10.1016/S0012-821X(02)01104-4.Google Scholar
- Payet G. Les Réunionnais et leur Volcan. Antenne Reunionnaise de l’Institute de Victimologie (St Denis, Réunion); 2007. p. 147.Google Scholar
- Peltier A, Bachèlery P, Staudacher T. Magma transport and storage at Piton de la Fournaise (la Réunion) between 1972 and 2007: a review of geophysical and geochemical data. J Volcanol Geotherm Res. 2009;184:93–108.View ArticleGoogle Scholar
- Peltier A, Massin F, Bachèlery P, Finizola A. Internal structures and building of basaltic shield volcanoes : the example of the Piton de la Fournaise terminal cone (la Réunion). Bull Volcanol. 2012;74:1881–97.View ArticleGoogle Scholar
- Peltier A, Beauducel F, Villeneuve N, Ferrazzini V, Di Muro A, Aiuppa A, et al. Deep fluid transfer evidenced by surface deformation during the 2014-2015 unrest at Piton de la Fournaise volcano. J Volcanol Geotherm Res. 2016;321:140–8. doi:10.1016/j.jvolgeores.2016.04.031.
- Perkins MC. Surviving Paradise. USA: Quidnunc Press; 2006. p. 338.Google Scholar
- Pieri DC, Baloga SM. Eruption rate, area, and length relationships for some Hawaiian lava flows. J Volcanol Geotherm Res. 1986;30:29–45.View ArticleGoogle Scholar
- Pinkerton H, Wilson L. Factors effecting the lengths of channel-fed lava flows. Bull Volcanol. 1994;56:108–20.View ArticleGoogle Scholar
- RED SEED Working Group. Conclusion: recommendations and findings of the RED SEED working group. In: AJL H, De Groeve T, Garel F, Carn SA, editors. Detecting Modelling and responding the effusive eruptions, vol. 426. London: Geological Society London Special Publications; 2016. p. 567–648.Google Scholar
- Robert B, Harris A, Gurioli L, Médard E, Sehlke A, Whittington A. Textural and rheological evolution of basalt flowing down a lava channel. Bull Volcanol. 2014;76:824. doi:10.1007/s00445-014-0824-8.View ArticleGoogle Scholar
- Roscoe R. The viscosity of suspensions of rigid spheres. Br J Appl Phys. 1952;3:267–9. doi:10.1088/0508-3443/3/8/306.View ArticleGoogle Scholar
- Rudd DR. Remote Sensing: A better view. Belmont: Duxbury Press; 1974. p. 12–3.Google Scholar
- Santora M. Réunion, once a surfer’s paradise, finds only sharks in its waters. NY: The New York Times (New York Edition), 12 August 2015; 2015. p. A5.Google Scholar
- Schmidt A. Volcanic gas and aerosol hazards from a future Laki-type eruption in Iceland. In: Papale P, editor. Volcanic hazards, risks and disasters, hazards and disasters series. Amsterdam: Elsevier; 2015. p. 377–97.View ArticleGoogle Scholar
- Schowengerdt RA. Remote sensing: models and methods for image processing. Burlington: Academic Press; 2007. p. 515.Google Scholar
- Silvestre AL, Santos CM, Ramalho C. Satisfaction and behavioural intentions of cruise passengers visiting the Azores. Tour Econ. 2008;14(1):169–84.View ArticleGoogle Scholar
- Singh SM. Simulation of solar zenith angle effect on global vegetation index (GVI) data. Int J Remote Sens. 1988;9(2):237–48.View ArticleGoogle Scholar
- Staudacher T, Ferrazzini V, Peltier A, Kowalski P, Boissier P, Catherine P, et al. The April 2007 eruption and the Dolomieu crater collapse, two major events at Piton de la Fournaise (la Réunion Island, Indian Ocean). J Volcanol Geotherm Res. 2009;184:126–37. doi:10.1016/j.jvolgeores.2008.11.005.
- Stevenson H. Jobs for the boys: the story of a family in Britain’s imperial heyday. Ipswich: Dove Books; 2009. p. 377.Google Scholar
- Stewart R. Reunion shark attacks scare surfers and tourists from beaches. The Daily Telegraph; 2015. www.telegraph.co.uk/travel/news/. Downloaded 20/01/17.
- Surfer Today. Is there a solution for the shark attack drama in Reunion Island? SurferToday; 2016. www.surfertoday.com/environment/12880. Downloaded 20/01/17.
- Tedesco D, Badiali L, Boschi E, Papale P, Tassi F, Vaselli O, et al. Cooperation on Congo volcanic and environmental risks. Eos, Transactions American Geophysical Union. 2007;88(16):2324–9250. http://dx.doi.org/10.1029/2007EO160001
- Tucker CJ, Gatlin JA, Schneider SR. Monitoring vegetation in the Nile delta with NOAA-6 and NOAA-7 AVHRR imagery. Photogramm eng Remote Sens. 1984;50(1):53–61.Google Scholar
- Van Wagner CE. Calculations on forest fire spread by flame radiation. Ottowa: Forestry Branch Departmental Publication 1185, Queen’s Printer and Controller of Stationary; 1967. p. 14.Google Scholar
- Vaxelaire D. L’histoire de La Réunion 1. Des origines à 1848. Saint-Denis: Édtions Orphie; 2012. p. 350.Google Scholar
- Vaxelaire D. L’histoire de La Réunion 2. De 1848 à 2012. Saint-Denis: Édtions Orphie; 2012. p. 703.Google Scholar
- Vicari A, Ciraudo A, Del Negro C, Herault A, Fortuna L. Lava flow simulations using discharge rates from thermal infrared satellite imagery during the 2006 Etna eruption. Nat Hazards. 2009;50:539–50. http://doi.org/10.1007/s11069-008-9306-7 View ArticleGoogle Scholar
- Vicari A, Ganci G, Behncke B, Cappello A, Neri M, Del Negro C. Near-real-time forecasting of lava flow hazards during the 12–13 January 2011 Etna eruption. Geophys Res Lett. 2011;38:L13317. http://doi.org/10.1029/2011GL047545 View ArticleGoogle Scholar
- Vicari A, Bilotta G, et al. LAV@HAZARD: web-Gis interface for volcanic hazard assessment. Ann Geophys. 2011;54:662–70. http://doi.org/10.4401/ag-5347 Google Scholar
- Villeneuve N, Neuville DR, Boivin P, Bachèlery P, Richet P. Magma crystallization and viscosity: a study of molten basalts from the Piton de la Fournaise volcano (la Réunion island). Chem Geol. 2008;256:242–51.View ArticleGoogle Scholar
- Viskanta R, Mengüç MP. Radiation heat transfer in combustion systems. Prog Energy Combust Sci. 1987;13:97–160.View ArticleGoogle Scholar
- Vlastélic I, Staudacher T, Semet M. Rapid change of lava composition from 1998 to 2002 at Piton de la Fournaise (Réunion) inferred from Pb isotopes and trace elements: evidence for variable crustal contamination. J Petrol. 2005;46:79–107.View ArticleGoogle Scholar
- Vlastélic I, Deniel C, Bosq C, Télouk P, Boivin P, Bachèlery P, et al. Pb isotope geochemistry of Piton de la Fournaise historical lavas. J Volcanol Geotherm Res. 2009;184:63–78.Google Scholar
- Vlastélic I, Gannoun A, Di Muro A, Gurioli L, Bachèlery P, Henot JM. Origin and fate of sulfide liquids in hotspot volcanism (la Réunion): Pb isotope constraints from residual Fe-cu oxides. Geochim Cosmochim Acta. 2016;194:179–92.View ArticleGoogle Scholar
- Wadge G. The variation of magma discharge during basaltic eruptions. J Volcanol Geotherm Res. 1981;11:139–68.View ArticleGoogle Scholar
- Walker GPL. Lengths of lava flows. Philos Trans R Soc Lond. 1973;274:107–18.View ArticleGoogle Scholar
- Weisel D, Stapleton F. Aloha O Kalapana. Honolulu: Bishop Museum Press; 1992. p. 153.Google Scholar
- Witter JB, Harris AJL. Field measurements of heat loss from skylights and lava tube systems. J Geophys Res. 2007;112 (B01203). doi:10.1029/2005JB003800.
- Wooster MJ, Rothery DA. Time-series analysis of effusive volcanic activity using the ERS along track scanning radiometer: the 1995 eruption of Fernandina volcano, Galápagos Islands. Remote Sens Environ. 1997;62:109–17.View ArticleGoogle Scholar
- Wooster MJ, Rothery DA, Kaneko T. Geometric considerations for the remote monitoring of volcanoes: studies of lava domes using ATSR and the implications for MODIS. Int J Remote Sens. 1998;19(13):2585–91.View ArticleGoogle Scholar
- Wright R. MODVOLC: 14 years of autonomous observations of effusive volcanism from space. In: Harris AJL, De Groeve T, Garel F, Carn SA, editors. Detecting, Modelling and responding to effusive eruptions, vol. 426. London: Geological Society, London, Special Publications; 2016. p. 23–53.Google Scholar
- Wright R, Blake S, Harris AJL, Rothery DA. A simple explanation for the space-based calculation of lava eruption rates. Earth Planet Sci Lett. 2001;192:223–33. doi:10.1016/S0012-821X(01)00443-5.View ArticleGoogle Scholar
- Wright R, Garbeil H, Harris AJL. Using infrared satellite data to drive a thermo-rheological/stochastic lava flow emplacement model: A method for near-real-time volcanic hazard assessment. Geophys Res Lett. 2008;35(L19307): doi:10.1029/2008GL035228.
- WTTC. Travel & Tourism Economic Impact 2015 Reunion. London: Annual Report of the World Travel and Tourism Council; 2015. p. 24.Google Scholar
- Zebker HA, Rosen P, Hensley S, Mouginis-Mark PJ. Analysis of active lava flows on Kilauea volcano, Hawaii, using SIR-C radar correlation measurements. Geology. 1996;24:495.View ArticleGoogle Scholar
- Zhang W, Hu Z, Liu Y, Chen H, Gao S, Gasching RM. Total rock dissolution using ammonium bifluoride (NH4HF2) in screw-top Teflon vials: a new development in open-vessel digestion. Anal Chem. 2012;84:10686–93.View ArticleGoogle Scholar
- Zinke J, Reijmer JJG, Thomassin BA. Seismic architecture and sediment distribution within the Holocene barrier reef-lagoon complex of Mayotte (Comoro archipelago, SW Indian Ocean). Palaeogeogr Palaeoclimatol Palaeoecol. 2001;175:343–68.View ArticleGoogle Scholar
- Zuccaro G, Cacace F, Spence RJS, Baxter PJ. Impact of explosive eruption scenarios at Vesuvius. J Volcanol Geotherm Res. 2008;178:416–53.View ArticleGoogle Scholar