Satellite and Ground Remote Sensing Techniques to Trace the Hidden Growth of a Lava Flow Field : The 2014-15 Effusive Eruption at Fogo Volcano ( Cape Verde )

1Istituto Nazionale di Geofisica e Vulcanologia, Osservatorio Etneo (INGV-OE), Catania, 95125, Italy, 2University of Cabo Verde (UNICV), CP 279 Praia, Santiago, Cape Verde, 3Instituto Tecnológico de Energias Renovábles de Tenerife (ITER), 38611 Granadilla de Abona, Tenerife, Spain, 4Instituto Volcanológico de Canarias (INVOLCAN), 38400 Puerto de la Cruz, Tenerife, Spain, 5National Protection Service of Cabo Verde, Praia, Santiago, Cape Verde


Introduction
When an effusive eruption starts, the maximum distance that a flow can travel can be easily estimated based on the measured peak effusion rate [1][2][3][4].However, when lava tubes develop within a flow, these prevent lava from cooling, increasing its ability to cover longer distances [3,5].Lava tubes are normally hidden below tens to hundreds of meters of lava and thus it is not easy to detect them [3,6].The increasing use of remote sensing techniques allows a more accurate and faster detection of lava tube formation [7], and is essential for hazard assessment and risk mitigation [4,8].
The general shape of a complex lava flow field is defined by a few arterial lava flows generally displaying aa (or rough) texture, with its outline modified by secondary lava flows normally having a pahoehoe (or smooth) surface [2,9,10].Aa is the type of lava that in solidified form is characterized by a rough, jagged, spiny and generally clinkery surface and is erupted at high flow rates.Conversely, pahoehoe is characterized by a smooth, billowy, or ropy surface and is erupted at low flow rates.When lava tubes develop within complex lava flows, their hidden path is revealed by the distribution of ephemeral vents or breakouts [3,11,12].Three types of ephemeral vents have been distinguished at Mt Etna [3], based on age and position with respect to the arterial flows and the tube system: (1) first-order ephemeral vents open at the flow front and normally feed secondary flows with aa morphology, playing a major role in the formation and propagation of lava tubes because they are continuously supplied by lava; (2) second-order ephemeral vents open on top of lava tubes and act as pressure release valves during periods of tube blockage, increase in the supply rate or inflation, and form aa or pahoehoe flow lobes; and (3) third-order ephemeral vents are located at the flow margins around flows that are no longer fed, draining their interior and producing small pahoehoe or toothpaste flows.Thus, pahoehoe and toothpaste morphologies on Etna are common at the margins of aa lava flows and master tubes or in the proximal portion of aa flows that produced extensive lava tube systems [3,6,13].
The process of endogenous growth or inflation in pahoehoe lava flows was described [14] and quantified [15] a long time ago, as also its role in the formation of lava tubes within sheet flows.Since then, other studies have been devoted to recognizing inflation in active and past lava flows [5,16], describing its features [17][18][19], its importance in the growth of aa lava flow fields [3,13] and in the emplacement of flood basalts [20,21].
Inflation of a lava flow is often related to the formation of lava tubes.Preferred pathways develop in the older portions of the liquid-cored flow as the flow advances, and these pathways can evolve into lava tube systems within a few weeks [15].Inflation combined with flow overlapping causes a slow and mostly undetected expansion of a complex lava flow field, and can result in sudden reactivation of previously stagnating lava flow fronts [3,6,15,[22][23][24].
Quantitative measurements of the amount of magma intruded within a lava flow field and causing inflation are scant and normally carried out only on a portion of a compound lava flow field [21,[23][24][25][26][27][28].On pahoehoe flows from Kilauea, a volumetric growth of ten times the initial flow lobe by inflation has been documented [15].However, to the best of our knowledge, there are no published data offering an estimate of the amount of endogenous growth for a whole lava flow field, although a previous study carried out on Kilauea evaluated the partitioning of lava between surface and tube flows but was limited by lava spilling into the ocean [29].In this paper, we present for the first time an estimation of the total volume endogenously intruded within the lava flow field of the Fogo 2014-2015 eruption.We obtained this estimate by comparing results of time averaged discharge rate (TADR) from satellite measurements [30] and of magma supply rate (SR) from the source as retrieved from daily measurements of SO 2 flux, which is related to the amount of degassing magma intruded within the system [31,32].Given that satellite images only detect surface emplacement of lava, and that SO 2 flux instead is related to the total magma volume stored in the source region [32] and feeding the eruption, the difference between the two should give an estimation of the magma intruded within the lava flow field and responsible for its inflation and supplying lava tubes.We compare these results with satellite and field data on the lava flow field growth as well as published erupted volumes, and explain the fast growth of the lava flow field at its northern end that caused the destruction of the villages of Portela and Bangaeira.Provided that routine SO 2 flux measurements and fast calculation of TADR from satellite imagery can be carried out, these results could be applied during future eruptions to detect lava flow field inflation.This would help predict, hopefully well in advance, the formation of lava tubes that extend lava flow fields at this and other basaltic volcanoes.

Geologic Background
Cape Verde Islands comprise a ~50 Ma old volcanic archipelago located ~700 km off West (W) Africa (Figure 1a).The archipelago lies on the African plate, which drifts at a velocity of ~0.9 cm/year above a hot spot mantle plume [33,34].Fogo Island is located at the SW end of the archipelago (Figure 1a), and has a conical shape, with a diameter of ~25 km and an elevation of ~2829 m (Figure 1b).The summit is truncated at ~1700 m a.s.l. by the Cha Caldera (Figure 1c), a 9 km wide depression open to East (E) by a flank collapse scar [35] formed ~73 ka [36].The upper part of this depression has a flat bottom, confined to the North (N), W and South (S) by the 1000 m high vertical cliff known as Bordeira (Figure 1c).The eastern boundary of the depression is occupied by the ~1100-m-high active volcano named Pico do Fogo (Figure 1c), which gave rise between 1500 and 1750 CE to several eruptions from its summit.After this lapse of time, ~9 eruptive events occurred from fissures at its base, with the last episode in 1995 [34] (Table S1 and Figure 1c).Eruptive activity at Fogo is fed directly from the mantle plume at depths greater than 16 km, with an estimated magma SR of 1.7 × 10 6 m 3 year −1 during the last three centuries [34].No shallow magma storage has been found during the last effusive eruptions, and an estimated 16-24 km depth has been inferred for the magma storage that fed the 1995 eruption [37].Erupted lavas range in composition from basanites, tephrites and nephelinites [38], with the 1995 eruption producing cogenetic basanites and phonotephrites [37].
Remote Sens. 2018, 10, x FOR PEER REVIEW 3 of 20 destruction of the villages of Portela and Bangaeira.Provided that routine SO2 flux measurements and fast calculation of TADR from satellite imagery can be carried out, these results could be applied during future eruptions to detect lava flow field inflation.This would help predict, hopefully well in advance, the formation of lava tubes that extend lava flow fields at this and other basaltic volcanoes.

Geologic Background
Cape Verde Islands comprise a ~50 Ma old volcanic archipelago located ~700 km off West (W) Africa (Figure 1a).The archipelago lies on the African plate, which drifts at a velocity of ~0.9 cm/year above a hot spot mantle plume [33,34].Fogo Island is located at the SW end of the archipelago (Figure 1a), and has a conical shape, with a diameter of ~25 km and an elevation of ~2829 m (Figure 1b).The summit is truncated at ~1700 m a.s.l. by the Cha Caldera (Figure 1c), a 9 km wide depression open to East (E) by a flank collapse scar [35] formed ~73 ka [36].The upper part of this depression has a flat bottom, confined to the North (N), W and South (S) by the 1000 m high vertical cliff known as Bordeira (Figure 1c).The eastern boundary of the depression is occupied by the ~1100-m-high active volcano named Pico do Fogo (Figure 1c), which gave rise between 1500 and 1750 CE to several eruptions from its summit.After this lapse of time, ~9 eruptive events occurred from fissures at its base, with the last episode in 1995 [34] (Table S1 and Figure 1c).Eruptive activity at Fogo is fed directly from the mantle plume at depths greater than 16 km, with an estimated magma SR of 1.7 × 10 6 m 3 year −1 during the last three centuries [34].No shallow magma storage has been found during the last effusive eruptions, and an estimated 16-24 km depth has been inferred for the magma storage that fed the 1995 eruption [37].Erupted lavas range in composition from basanites, tephrites and nephelinites [38], with the 1995 eruption producing cogenetic basanites and phonotephrites [37].

Methods
Fogo volcano and the fumaroles located within the summit crater of Pico do Fogo (Figure 1c) have been monitored since 2008 by the Instituto Volcanológico de Canarias (INVOLCAN, Canary Islands, Spain [40][41][42]).The 2014-2015 eruption was jointly monitored by researchers from the University of Cape Verde (UNICV) and from INVOLCAN directly in field, and remotely by the Istituto Nazionale di Geofisica e Vulcanologia-Osservatorio Etneo (INGV-OE) of Italy through satellite images [30,43].The amount and quality of data collected during the eruption allowed us to reconstruct the expansion of the lava flow field even during the phases of inflation and endogenous growth, leading to the identification of a master lava tube that was eventually responsible for the destruction of the villages of Portela and Bangaeira.
The chronology of the eruption described in this paper has been obtained through almost daily field surveys carried out by several authors of this manuscript using photos and thermal images, as well as through satellite images (i.e., Landsat-8 OLI, EO-1 ALI, Pléiades) and time-lapse frames recorded from a fixed monitoring camera.The fixed camera was installed by INVOLCAN on 27 November 2014 on a 4-m-high pole located at the S entrance of the Cha Caldera valley (Figure 1c), allowing a view from the S and from a distance of ~2 km away from the eruptive fissure and proximal lava flow field.It recorded frames every 4-5 s from 30 November 2014 to 27 January 2015, automatically turning to infrared at dusk.Camera failure caused a lack of image acquisition between 1 and 12 December 2014.The morphology of the lava flow field, with the aim of identifying the hidden path of lava tubes, has been analyzed using the images freely available in Google Earth and acquired on 2 March 2016 (Image c 2016 CNES/Astrium).

Magma Supply Rate
Two to twelve (average 5.5) daily plume SO 2 flux measurements were performed along road traverses using a mini-DOAS.The mini-DOAS instrument is based on an Ocean Optics USB2000 UV spectrometer, which collects the UV radiation via an optical fiber coupled to a vertically pointing telescope [44,45].The instrument position was tracked using a handheld Global Positioning System (GPS) receiver.The mini-DOAS measures the SO 2 column density in parts per million per mass at every measurement point along the plume transect.Integrated path values were obtained by adding the products of the SO 2 column density and perpendicular displacement for each segment along the path.The SO 2 emission rate was obtained by multiplying the integrated path values by the average wind speed obtained using the Weather Research and Forecasting (WRF) model (http://www.wrf-model.org/index.php)[43].
Using the daily SO 2 flux estimations, we calculated magma degassing rates and volumes of magma degassed in the period between 28 November 2014 and 7 February 2015.The total volume of degassed magma (V d ) was calculated using [32]: where V S is the volume of elemental sulfur, [S] is the weight fraction of sulfur degassed per unit of magma, x is the crystal fraction and ρ is the density of magma.We used a value of 10% [46] as mean crystal fraction, 0.3 wt % as mean original sulfur content, and 2600 kg/m 3 as density of magma, assimilating Fogo's magma to Etna's basalt [47].Plume SO 2 measurements were undertaken daily, allowing the degassed magma flux (in m 3 s −1 ) to be calculated, and from the daily degassed magma volume we have obtained a time-averaged (daily) magma supply rate (SR) [48].Given SR, we can integrate it through time to calculate the volume of degassed magma over the measurement period.Errors associated with the SO 2 method have been shown to range between 20% and 30% [49].

TADR Estimation
The thermal activity at Fogo volcano was observed and quantified from multispectral-satellite data via the HOTSAT system [50].This volcano monitoring system processes data acquired by the Moderate Resolution Imaging Spectroradiometer (MODIS, pixel nominal resolution: 1 km) as well as Spinning Enhanced Visible and Infrared Imager (SEVIRI, pixel resolution: 3 km at nadir) sensors to: (i) locate thermal anomalies; (ii) compute the associated radiant heat flux; and (iii) provide an estimation of the time averaged discharge rate (TADR [51]).Clouds in satellite data can alter or mask the thermal activity, thus they are retrieved by using a texton-based cloud detection algorithm [52], and a cloud index is provided to interpret the radiant heat flux signal.The hotspot detection algorithm is based on a contextual approach [53] and the radiant heat flux is computed following Wooster et al. [54].Radiant heat flux curves provided by HOTSAT were compared against ground thermal camera data providing comparable results both for a small lava flow [55] and for the Nyiragongo lava lake [56].TADR is computed from radiant heat flux by using: where Q is the total thermal flux obtained summing up the radiant heat flux computed for each hot spot pixel, ρ is the lava density, C p is the specific heat capacity, ∆T is the eruption temperature minus temperature at which flow stops, C L is the latent heat of crystallization, and ∆φ is the volume percent of crystals that form while cooling through ∆T [57].TADR were calculated by considering, in the variability range of each lava parameter, the largest values [30].Moreover, for the percentage of crystals and temperature variations, we used the recent constrained values [46].HOTSAT uses satellite images as soon as they are available, i.e., every 15 min in the case of SEVIRI full disk, up to 5 min in the case of SEVIRI Rapid Scanning for the northern third of Meteosat disc.

Chronology of the Eruption
The first few days of the eruption were characterized by sudden changes in the eruptive activity and fast growth of the lava flows.The main events are summarized in Table S1, and Figure 2a displays the key phases of growth of the lava flow field.Figure 2b-g shows the active portion of the lava flow field (in red) as observed by satellites, compared to the nearest lava flow field map (yellow outline).
The eruption started on 23 November 2014 at 09:45 LT (local time; all times in this paper are LT and given as hh:mm) within the Cha Caldera depression.It began with the opening of an east-northeast-west-southwest (ENE-WSW), ~700-m-long eruptive fissure at the southwest (SW) base of the Pico do Fogo cone between 2200 m and 1800 m above sea level (a.s.l.), close to the previous 1995 fissure (Figures 1 and 2).Four eruptive vents opened along the fissure, giving rise to powerful Strombolian explosive activity that within a few hours became lava fountains, forming a ~6 km high ash plume (estimated by the Capo Verde Express pilots).At 10:00, lava flows erupted from a vent at ~2100 m a.s.l., spreading at first SW along the fissure, then splitting at its base in two directions: south (S) (Flow 1) and north-northwest (NNW, Flow 2), being divided by the hill of Mount Beco (Figure 2a,b), and then following the paved road that crossed the Cha Caldera from S to north (N) (Figure 1c), which was the only access to the villages of Portella and Bangaeira.Flow 1 travelled ~1.2 km during the first 2 h, and by 12:00 covered the Cha Caldera paved road (Figure 2a).Flow 2 was initially flowing parallel to Flow 1 along the fissure, expanding SW towards Mount Beco, and then it moved NNW, bounded by the high topography of Mount Beco and by the 1995 lava flow field (Figure 2a).
On 24 November, the eruptive fissure continued to open upslope, with up to seven active explosive vents along it (Table S1, Figure 2b).The two channel-fed arterial lava flows (Flows 1 and 2, Figure 2b) had aa surface and reached ~2 km in length.Flow 1 slowed down and became stagnated by dusk, while Flow 2 was spreading NNW.On 25 November morning, Flow 1 was again slowly spreading S, whereas Flow 2 was expanding N and approaching the Portela village.By the afternoon, Flow 2 became thicker and faster, and when on 26 November it encountered the N cliff, the flow front split into two branches directed east (E) to Portela and west (W) (Figure 2a,b).The flow front directed E began covering the building of the Fogo Natural Park, completed just a few months earlier.On 27 November, a new vent (Vent 1, Figures 2a and 3a) opened at the base of the eruptive fissure, feeding a small ~30-40 m wide lava flow and expanding above Flow 2. While powerful Strombolian explosions from the upper part of the fissure were still building up two cinder cones (Figure 3a,b), several overflows occurred from the upper S rim of the fissure, breaching and eroding the still growing cinder cones and forming several lava branches overlapping Flow 1 (Figure 3b).This activity also continued the following days (Figure 3c,d), and lava emerging from Vent 1 gradually built a tumulus (Figure 3d).Vent 1 drained lava from the cinder cones on the upper fissure for three days before a significant decline of the explosive activity was observed at the summit vents, accompanied by a reduction of the surface temperature at the cones along the fissure and in the upper lava flow field (Figure 3d,e).On 30 November, most of the cinder cones along the eruptive fissure were already built up, and the uppermost tumulus and lava channel fed by Vent 1 began sealing over, forming several skylights along its path, and thus revealing the presence of a proximal lava tube (Figure 3e).Explosions and pulsating ash emission occurred along the uppermost fissure, indicating a lower magma level when compared to the previous lava fountain activity [58].On 30 November, the time-lapse images from the fixed monitoring camera installed at the S end of the Cha Caldera became available (Figure 1c), showing a fast and continuous sealing of the upper lava tube and occasional overflows from Vent 1.In general, when explosive activity at the summit vents along the fissure was increasing, we observed a decrease in the explosive and degassing activity at the skylights along the upper tube, and vice versa, indicating that the two systems were connected at a shallow depth.by dusk, while Flow 2 was spreading NNW.On 25 November morning, Flow 1 was again slowly spreading S, whereas Flow 2 was expanding N and approaching the Portela village.By the afternoon, Flow 2 became thicker and faster, and when on 26 November it encountered the N cliff, the flow front split into two branches directed east (E) to Portela and west (W) (Figure 2a,b).The flow front directed E began covering the building of the Fogo Natural Park, completed just a few months earlier.On 27 November, a new vent (Vent 1, Figures 2a and 3a) opened at the base of the eruptive fissure, feeding a small ~30-40 m wide lava flow and expanding above Flow 2. While powerful Strombolian explosions from the upper part of the fissure were still building up two cinder cones (Figure 3a,b), several overflows occurred from the upper S rim of the fissure, breaching and eroding the still growing cinder cones and forming several lava branches overlapping Flow 1 (Figure 3b).This activity also continued the following days (Figure 3c,d), and lava emerging from Vent 1 gradually built a tumulus (Figure 3d).Vent 1 drained lava from the cinder cones on the upper fissure for three days before a significant decline of the explosive activity was observed at the summit vents, accompanied by a reduction of the surface temperature at the cones along the fissure and in the upper lava flow field (Figure 3d,e).On 30 November, most of the cinder cones along the eruptive fissure were already built up, and the uppermost tumulus and lava channel fed by Vent 1 began sealing over, forming several skylights along its path, and thus revealing the presence of a proximal lava tube (Figure 3e).Explosions and pulsating ash emission occurred along the uppermost fissure, indicating a lower magma level when compared to the previous lava fountain activity [58].On 30 November, the timelapse images from the fixed monitoring camera installed at the S end of the Cha Caldera became available (Figure 1c), showing a fast and continuous sealing of the upper lava tube and occasional overflows from Vent 1.In general, when explosive activity at the summit vents along the fissure was increasing, we observed a decrease in the explosive and degassing activity at the skylights along the upper tube, and vice versa, indicating that the two systems were connected at a shallow depth.In early December, the front of Flow 2 reached the N cliff of Bordeira, where it slowed down and inflated (Figure 2a).At that time, most of the master tube within Flow 2 had already formed, and only ephemeral vents along the uppermost lava flow field were feeding small lava lobes on the flow surface (Figure 2c).At this stage, the lava flow field began expanding W through a second-order ephemeral vent, giving rise to initial Flow 3, whereas Flow 1 was already inactive, and most of the surface of Flow 2 between Vent 1 and Portela was crusted over (Figure 2c).On 7 December, two first-order ephemeral vents [3] opened at the E margin of the Flow 2 front and at the exit of the lava tube that was already ~5.7 km long.These vents were discharging lava towards Portela and Bangaeira at a high rate, thus expanding the lava flow field a further ~1 km by 16 December (Figure 2c).From this date to the end of December, the lava flow field widened by westward expansion of Flow 3 (Figure 2d,e).In addition, this flow stopped, widened and inflated as soon as it reached the W cliff, but did not expand further because its supply was eventually cut off.From the end of December (Figure 2f) to the end of the eruption (Figure 2g), only small third-order ephemeral vents [3] were observed along the lava flow field, causing thickening of the proximal lava flow field [59], drainage of the main flows, and slight changes in the outline of the lava flow field.The eruption eventually ended on 8 February 2015, after a gradual decline of the output rate [30].

Field Surveys and Observations of Lava Tubes
A survey along the 2014-2015 inactive lava flow field was carried out on 23 November 2015, one year after the start of the eruption.The first lava flow (Flow 1) emerging from the fissure had aa morphology (observed in the field while spreading), but this flow is not exposed.The topographic surface was then covered by a surge of lava that gave rise to a sheet flow of shining blue pahoehoe with a glassy surface, very similar to those observed at Kilauea [60].This was covered by an aa flow that formed a channel but not a lava tube.The levees were ~1 m high outside the channel and above the pahoehoe sheet flow, and ~2 m high inside the channel.The aa surface inside the channel was 2 m lower than the levees, indicating drainage of the channel and possibly erosion of the base.In the lower portion of the channel, there was a break in slope where the inner coating was detached from the walls and collapsed within the channel.This inner coating was a mere 10-cm-thick layer, suggesting that the channel was not active for long (only one coating) and that lava feeding it was rather viscous (indicated by the thickness of the coating [13]).In fact, Flow 1 was only active for six days and intermittently fed by overflows from the S rim of the fissure (Figure 3b-d).Further down along the channel and a few meters below this site, the inner coating was jammed in the middle of the channel forming a significant obstruction, and the channel disappeared below aa flows coming from overflows both from the channel and from the fissure.The intermittent supply and fast drainage of the channel were apparently responsible for its instability and blockage, deactivation, and for the lack of a lava tube within this flow.As already observed [3,5], a steady supply is the key requirement for lava tube formation.
Further N, the eruptive fissure gave rise to several lava benches along the fissure crossing the cone.These benches were formed by lava stagnating at progressively lower levels while the fissure was propagating along the cone and down slope (Figure 4a).During this propagation, at the base of the cones, it formed 3-4 layers of overlapped bubbly grey lava ~30-50 cm thick suggesting a temporary lava pond, just where Vent 1 was located (Figure 4a).This pond eventually broke when the eruptive fissure propagated down slope, confirmed by a fracture crossing it.The dike emptied from the base of the fracture, giving rise to pahoehoe sheet flows spreading like a fan.We have found a lava pond about 20 m down slope from the base of the fissure that eventually drained completely forming a skylight (Figure 4a,b).At the base of this skylight, there was a lava tube ~2-3 m in diameter and located ~20 m below the topographic surface (Figure 4c).About 50 m down slope from this pit, a larger lava pond, ~50 m wide, drained forming a circular collapsed structure ("shatter ring" [5]).Along the path of the lava tube and further W, a 50 m wide tumulus elongated E-W suggested the prolonged feeding to a lava tube in this direction, possibly supplied by lava during the final stages of the eruption [59].The last pahoehoe flows erupted in Portela were very black, compact and poorly vesiculated.

Lava Flow Field Morphology
With the aim of identifying the hidden path of lava tubes within the lava flow field, we analyzed the lava flow morphology from the images freely available in Google Earth. Figure 5 displays the portion of the lava flow field of mainly aa flow surfaces with the red outline, and with blue color the portions that are mainly pahoehoe flows.Flow 1 is essentially made up of aa flows, while Flow 3 is mostly pahoehoe.Flow 2, on the other hand, is the longest and most complex, characterized by pahoehoe surface at the base of the eruptive fissure and at the boundary of the flow, and by aa surface in the middle of the flow and at the flow front.Flow 1 comprises mainly overlapped sheet flows that were intermittently supplied by lava through overflows from the eruptive fissure (Figure 3b-d).It did not show any significant amount of inflation and its outline is not surrounded by pahoehoe flows.The proximal portion of the lava flow field at the base of the eruptive fissure and around Vent 1 is mainly of pahoehoe flows because this portion emplaced during the final stages of the eruption by drainage of degassed lava at low output rate through the skylight that replaced Vent 1 at the base of the eruptive fissure (Figures 2a and 3e).Most of the surface of Flow 2 is made up of aa clinkers.It was emplaced as an arterial flow surrounded by late third-order ephemeral vents causing the partial drainage of the inflated flow interior and the emplacement of pahoehoe flow lobes all around it and especially at the exit of the two first-order ephemeral vents that covered the village of Portela with lava.Flow 3 is mostly formed by pahoehoe surface.It took place from a second-order ephemeral vent draining the middle portion of the inflated flow field when the front of Flow 2 was stopped and inflating against the N Bordeira wall.

Lava Flow Field Morphology
With the aim of identifying the hidden path of lava tubes within the lava flow field, we analyzed the lava flow morphology from the images freely available in Google Earth. Figure 5 displays the portion of the lava flow field of mainly aa flow surfaces with the red outline, and with blue color the portions that are mainly pahoehoe flows.Flow 1 is essentially made up of aa flows, while Flow 3 is mostly pahoehoe.Flow 2, on the other hand, is the longest and most complex, characterized by pahoehoe surface at the base of the eruptive fissure and at the boundary of the flow, and by aa surface in the middle of the flow and at the flow front.Flow 1 comprises mainly overlapped sheet flows that were intermittently supplied by lava through overflows from the eruptive fissure (Figure 3b-d).It did not show any significant amount of inflation and its outline is not surrounded by pahoehoe flows.The proximal portion of the lava flow field at the base of the eruptive fissure and around Vent 1 is mainly of pahoehoe flows because this portion emplaced during the final stages of the eruption by drainage of degassed lava at low output rate through the skylight that replaced Vent 1 at the base of the eruptive fissure (Figures 2a and 3e).Most of the surface of Flow 2 is made up of aa clinkers.It was emplaced as an arterial flow surrounded by late third-order ephemeral vents causing the partial drainage of the inflated flow interior and the emplacement of pahoehoe flow lobes all around it and especially at the exit of the two first-order ephemeral vents that covered the village of Portela with lava.Flow 3 is mostly formed by pahoehoe surface.It took place from a second-order ephemeral vent draining the middle portion of the inflated flow field when the front of Flow 2 was stopped and inflating against the N Bordeira wall.

Pyroclastic Cone
During the first week of the eruption, the explosive activity initially built up a spatter rampart and then above it two cinder cones elongated NE-SW along the eruptive fissure (Figure 2a).They eventually merged, forming a unique cinder cone that extended between 1950 m and 2200 m a.s.l., was ~125 m high and with a ~500 m wide base, and with a crater width of ~165 m (measured from the map [61]).This results in a vesiculated volume of ~11.8 × 10 6 m 3 .Considering ~20% vesiculation for the deposit and ~30% for the magma typical of basaltic compositions [28,51,62,63], a volume of ~5.9 × 10 6 m 3 dense rock equivalent (DRE) is obtained for the pyroclastic cone built up during the explosive phase of the eruption.

Lava Flow Field Growth
We have extracted the parameters of lava flow field growth, taking the maximum length of the three arterial flows measured along their main axes at different stages of the eruption.Figure 6 shows the expansion of Flow 1, Flow 2 and Flow 3 with time.Flow 1 displayed a significant growth only

Pyroclastic Cone
During the first week of the eruption, the explosive activity initially built up a spatter rampart and then above it two cinder cones elongated NE-SW along the eruptive fissure (Figure 2a).They eventually merged, forming a unique cinder cone that extended between 1950 m and 2200 m a.s.l., was ~125 m high and with a ~500 m wide base, and with a crater width of ~165 m (measured from the map [61]).This results in a vesiculated volume of ~11.8 × 10 6 m 3 .Considering ~20% vesiculation for the deposit and ~30% for the magma typical of basaltic compositions [28,51,62,63], a volume of ~5.9 × 10 6 m 3 dense rock equivalent (DRE) is obtained for the pyroclastic cone built up during the explosive phase of the eruption.

Lava Flow Field Growth
We have extracted the parameters of lava flow field growth, taking the maximum length of the three arterial flows measured along their main axes at different stages of the eruption.Figure 6 shows the expansion of Flow 1, Flow 2 and Flow 3 with time.Flow 1 displayed a significant growth only during the first two days of the eruption, expanding more than 2.0 km on the first day, but its length remained stable at 2.9 km after 26 November, when the lava was directed NNW along Flow 2. Flow 2 was the longest, reaching its maximum length of ~7.9 km on 8 December.Flow 3 formed after 26 November by a second-order ephemeral vent opened at the W margin of a partially roofed over lava channel, but its main expansion occurred mostly after 16 December when Flow 2 stopped, and caused a westward widening of the lava flow field.Flow 3 reached its maximum length of 3.4 km on 28 December.From December until the end of the eruption, only small secondary pahoehoe flows changed the outline of the lava flow field, causing a significant thickening and widening of the middle portion, but no more increase in length.November by a second-order ephemeral vent opened at the W margin of a partially roofed over lava channel, but its main expansion occurred mostly after 16 December when Flow 2 stopped, and caused a westward widening of the lava flow field.Flow 3 reached its maximum length of 3.4 km on 28 December.From December until the end of the eruption, only small secondary pahoehoe flows changed the outline of the lava flow field, causing a significant thickening and widening of the middle portion, but no more increase in length.

Lava Flow Field Volume
Figures 6 and 7 show the difference in volume between the TADR measured by satellite (data from [30]) and the SR calculated here from SO2 flux measurements [43].It is worth noting that SO2 flux measurements started on 27 November, thus the SR estimates are lacking for the first three days of eruption.The difference between intruded and erupted magma volume becomes greater than 2 × 10 6 m 3 in the second half of December, when the lava tube within Flow 2 was mostly sealed and started draining, thus feeding Flow 3 with mostly pahoehoe lava.The tube efficiency caused the lack

Lava Flow Field Volume
Figures 6 and 7 show the difference in volume between the TADR measured by satellite (data from [30]) and the SR calculated here from SO 2 flux measurements [43].It is worth noting that SO 2 flux measurements started on 27 November, thus the SR estimates are lacking for the first three days of eruption.The difference between intruded and erupted magma volume becomes greater than 2 × 10 6 m 3 in the second half of December, when the lava tube within Flow 2 was mostly sealed and started draining, thus feeding Flow 3 with mostly pahoehoe lava.The tube efficiency caused the lack of large surface aa flows (Figure 2c-g), thus the lava flow field was dotted with only small pahoehoe lobes.Figure 6 displays the difference between supplied (SR) and erupted (TADR) magma volumes during the whole eruption.The difference between SR and TADR increases from mid-December 2014 until the end of the eruption.The final difference between supplied and erupted magma is 14.7 ± 8.8 × 10 6 m 3 (Figure 7), and we suggest that this volume represents the magma causing inflation and endogenous growth of the lava flow field, undetected by satellite since it did not expand on the surface.

Discussion
Comparing this episode with lava fountains observed at Mt Etna, we assume that the cinder cone was mainly formed by proximal ballistic fallout during the lava fountain phase of the eruption [63][64][65].Given that the cone built up in seven days (between 23 and 30 November 2014), the DRE TADR [49] for its growth is ~9.7 m 3 s −1 .This is much higher than the DRE 0.12 m 3 s −1 measured for the growth of the cinder cone within Etna's Bocca Nuova summit crater in 2012 [28] and than the DRE 1.5 m 3 s −1 of the 2002 Laghetto cone on Etna's S flank [63].It is also lower than the 50 m 3 s −1 that gave rise to the New SE Crater during powerful explosive episodes in 2011-2013 [64].
The direct observation of the eruption compared with thermal images collected from both the ground and satellite, as well as the morphology map distinguishing pahoehoe from aa surfaces, has allowed us to recognize the stages of formation and growth of the lava tubes within the lava flow field, their position and extension (Figure 8).The first sector of lava tube to become sealed was the uppermost portion within the eruptive fissure and up to Vent 1.This was ~1.1 km long and sealed in ~7 days.In fact, as shown in Figure 3e, on 30 November the eruptive fissure was mostly roofed over, there was a tumulus at its base above Vent 1, and several skylights along the upper tube.This caused a significant cooling of the flow surface, as shown by the comparison between Figure 3b-d

Discussion
Comparing this episode with lava fountains observed at Mt Etna, we assume that the cinder cone was mainly formed by proximal ballistic fallout during the lava fountain phase of the eruption [63][64][65].Given that the cone built up in seven days (between 23 and 30 November 2014), the DRE TADR [49] for its growth is ~9.7 m 3 s −1 .This is much higher than the DRE 0.12 m 3 s −1 measured for the growth of the cinder cone within Etna's Bocca Nuova summit crater in 2012 [28] and than the DRE 1.5 m 3 s −1 of the 2002 Laghetto cone on Etna's S flank [63].It is also lower than the 50 m 3 s −1 that gave rise to the New SE Crater during powerful explosive episodes in 2011-2013 [64].
The direct observation of the eruption compared with thermal images collected from both the ground and satellite, as well as the morphology map distinguishing pahoehoe from aa surfaces, has allowed us to recognize the stages of formation and growth of the lava tubes within the lava flow field, their position and extension (Figure 8).The first sector of lava tube to become sealed was the uppermost portion within the eruptive fissure and up to Vent 1.This was ~1.1 km long and sealed in ~7 days.In fact, as shown in Figure 3e, on 30 November the eruptive fissure was mostly roofed over, there was a tumulus at its base above Vent 1, and several skylights along the upper tube.This caused a significant cooling of the flow surface, as shown by the comparison between Figures 3b-d and 3e.Tubes on steep slopes and along eruptive fissures form more rapidly, because local changes in discharge rate can easily trigger inward accretion of levees and roofing over of lava channels [3,5].The second sector of lava tube to form was the longest, extending on the low slope between Vent 1 and the N cliff of Bordeira, in the middle portion of Flow 2. This tube was the down slope continuation of the proximal portion; it was about ~4.8 km long and formed in about 10 days (Figure 2c).It led to the destruction of the two villages of Portela and Bangaeira by fast-spreading lava flows erupted from first-order ephemeral vents at the E margin of Flow 2 front.Indeed, on 16 December (Figure 2c), only small pahoehoe lobes were spreading along the upper flow field, most of the lava being carried within the ~5.9 km master tube to the N end of the lava flow field.
The third sector of lava tube to form was the ~1.1 km branch that fed Flow 3. It formed much more slowly than the previous, because it started on 26 November but then stopped between 29 November and 7 December.It was only after the halting of the Flow 2 front that the middle portion of the lava field inflated, and a second-order ephemeral vent opened on the W margin of the lava flow field, acting as a pressure-release valve [3].It drained the flow field interior feeding Flow 3 with mostly pahoehoe lava.The supply to Flow 3 ended by 23 December, and no further expansion of Flow 3 was observed after this date.No evidence of lava tubes has been found within Flow 1.
The distribution of pahoehoe and aa flow surfaces has important implications for recognizing the hidden path of lava tubes.In fact, Flow 1, which is mainly made up of aa surfaces, did not show any evidence of lava tubes in the field.It was characterized by a drained channel obstructed by its collapsed inner coating, indicating complete drainage of the channel and successive lack of supply.Flow 1 was intermittently fed by overflows from the fissure, a condition that does not enable lava tube development [3,5].Flow 2, which mainly comprises aa surfaces and pahoehoe at the margins of The sector of lava tube to form was the longest, extending on the low slope between Vent 1 and the N cliff of Bordeira, in the middle portion of Flow 2. This tube was the down slope continuation of the proximal portion; it was about ~4.8 km long and formed in about 10 days (Figure 2c).It led to the destruction of the two villages of Portela and Bangaeira by fast-spreading lava flows erupted from first-order ephemeral vents at the E margin of Flow 2 front.Indeed, on 16 December (Figure 2c), only small pahoehoe lobes were spreading along the upper flow field, most of the lava being carried within the ~5.9 km master tube to the N end of the lava flow field.
The third sector of lava tube to form was the ~1.1 km branch that fed Flow 3. It formed much more slowly than the previous, because it started on 26 November but then stopped between 29 November and 7 December.It was only after the halting of the Flow 2 front that the middle portion of the lava field inflated, and a second-order ephemeral vent opened on the W margin of the lava flow field, acting as a pressure-release valve [3].It drained the flow field interior feeding Flow 3 with mostly pahoehoe lava.The supply to Flow 3 ended by 23 December, and no further expansion of Flow 3 was observed after this date.No evidence of lava tubes has been found within Flow 1.
The distribution of pahoehoe and aa flow surfaces has important implications for recognizing the hidden path of lava tubes.In fact, Flow 1, which is mainly made up of aa surfaces, did not show any evidence of lava tubes in the field.It was characterized by a drained channel obstructed by its collapsed inner coating, indicating complete drainage of the channel and successive lack of supply.Flow 1 was intermittently fed by overflows from the fissure, a condition that does not enable lava tube development [3,5].Flow 2, which mainly comprises aa surfaces and pahoehoe at the margins of the flow, has been characterized by a long lava tube extending from the eruptive fissure down to V1 (~1.1 km), and from here all along the flat surface of Cha Caldera up to the N boundary of Bordeira (~4.8 km), for a total of ~5.9 km length.Another small sector of lava tube, ~1.1 km long, formed during the second half of December as a lateral branch of the master tube that fed Flow 2. This fed Flow 3, made up almost entirely by pahoehoe surface.The presence of the tube system was confirmed in the proximal lava flow field by the skylight observed during the field survey (Figure 4c), displaying a ~2-3 m wide lava tube located ~20 m below the flow surface.The skylight inner wall was coated by pahoehoe lava congealed while dribbling inwards and showing final drainage of fluid lava within the upper tube (Figure 4a,b).
Thus, based on the data here described, a master tube developed within the eruptive fissure and along Flow 2 after just one week from the start of the eruption, and its growth completed within about 10 days from the start of the eruption.This tube was responsible for the fast lengthening of Flow 2, and formed a short branch to the W as soon as Flow 2 slowed down when encountering the N cliff of Bordeira.The magma inside the tube then accumulated upslope inflating the lava field and eventually feeding Flow 3. The path of Flow 2 was dictated by the line of maximum steepness [57], but its expansion was probably favored also by the smooth surface of the paved road crossing the Cha Caldera from S to N, given that most of the road had been covered by lava, as already happened also during the 1995 eruption [39].
To estimate the quantity of magma intruding the lava flow field and causing its inflation, we have compared the measurements of SR derived from SO 2 flux measurements [43], and the TADR derived from satellite [30].Although the measurements of SO 2 flux, and then of SR, are lacking for the first three days of the eruption, we can assume that at that time there was no difference between supplied and erupted magma, given that the lava flows were spreading with no significant surface cooling (Figure 2b), and that there was no inflation or endogenous growth of the lava flow field.The difference between SR and TADR increased significantly after 9 December, as soon as Flow 2 stopped advancing (Figures 6 and 7), and this is when Flow 3 increased its length.Once Flow 2 reached its maximum length allowed by surface crust cooling [9], the magma flowing within the master tube accumulated to the rear and upslope, inflating the middle portion of the flow field.This triggered an expansion of the tube towards the W for ~1.1 km, feeding Flow 3. Flow 3 was mostly pahoehoe, and fed by part of the degassed interior causing the inflation of the lava flow field.
Another significant increase in the difference between SR and TADR occurred after 21 December when also Flow 3 stopped advancing (Figures 6 and 7), and only small pahoehoe flow lobes were observed on the proximal lava flow field (Figure 2e).Flow 3 stopped after coming up against the W Bordeira wall and following a significant decline in the output rate [30].At this stage, the supply of lava was not sufficient to keep the ~5.9 km long tube active.Lava thus poured out from the base of the eruptive fissure [59], thickening the proximal lava flow field with a pahoehoe fan of lava from the skylight that eventually drained at the end of the effusive episode (Figure 8).This is the second major inflation stage of the lava flow field, not followed by major lava flow expansion.This is probably due to the confining effect of the lava crustal growth [9].By the end of the eruption, the erupted volume estimated by satellite measurements is 26.5 ± 3.2 × 10 6 m 3 DRE (Figure 7) corresponding only to surface lava, while mass volume obtained from SO 2 flux measurements is 46.9 ± 5.74 × 10 6 m 3 .Considering the respective uncertainty associated to each estimation, a total difference of 14.7 ± 8.8 × 10 6 m 3 DRE was obtained between SR and TADR.We argue that this differential volume, which is about 31% of the total erupted volume, caused inflation of the lava flow field and endogenous growth.
Differences between published bulk and DRE volumes can be also ascribed to voids from the tubes, drained flows and empty spaces between crust and flow.The amount of voids within the flow field is testified by the vertical displacement by cooling and contraction measured using spaceborne synthetic aperture radar interferometry ~9 months after the end of the eruption [64], and observed also at other volcanoes [67].It is noteworthy that the greatest contraction affected the thickest zones of the flow field, these being the eruptive fissure and the middle line of Flow 2, following the path of the master tube that fed the lava flow field (Figure 8).
The comparison between SR and TADR highlighted a growing difference between the two values from early December 2014 onwards.This is consistent with surface cooling and inflation of the lava field, and with the endogenous accumulation of lava.It is worth noting that the main destruction of Portela and Bangaeira occurred at this very stage of the eruption, when there was enough magma available within the inflated flow field to feed fast-spreading secondary flows from the exit of the master tube.
Codes such as MAGFLOW [68,69] and DOWNFLOW [70,71] have been developed to provide near-real-time prediction of lava flow paths.These and other models have generally been shown to provide a good match with the true path of supply-limited lava flows governed by underlying topography.However, the evolution of more complex cooling-limited flow fields may be much harder to predict, especially where far-field breakouts are developing [12,72,73].In these latter cases, inflation of the flow front has been a key precursor to breakout formation and attendant hazards.Thus, an observed mismatch between SR and TADR measurements on an actively emplacing lava flow could be used to inform predictive modeling and hazard identification.If there is no discrepancy between the SR and TADR, then well-established codes could effectively predict the path of the flow.However, if a discrepancy emerges, then inflation is occurring and hazards may be concentrated at stagnant flow fronts.
The development of lava tubes within the Fogo lava flow field was favored by a stable supply of lava along Flow 2, and by the flat topography on which lava was spreading.We suggest that any time lava flow spreads on a flat surface, the early comparison between TADR estimated by satellite and SR obtained from SO 2 flux measurements might allow recognizing if endogenous growth is occurring, leading to lava tube expansion and to potentially much greater length and destructive power of the lava flow field.
The 2014-2015 effusive eruption at Fogo was very similar to the previous 1995 eruption.The effusive vents of both eruptions were located in the same zone, i.e., at the W base of Pico do Fogo cone (Figure 1), sharing the same NE-SW orientation.The duration and erupted volume of both eruptions is similar, with 54 days of eruption and 46 × 10 6 m 3 lava erupted in 1995 [33] and 77 days and ~46 × 10 6 m 3 erupted magma in 2014-2015.Given that the last two lava flow fields affected the same flat ground of Cha Caldera, and that this area is the most prone to lava invasion [59], should a new eruption occur in the same zone, we suggest measuring and comparing the calculated SR from SO 2 flux measurements and the TADR obtained by satellite as soon as possible.This would enable the detection of major phases of inflation of the lava flow field in good time to prevent damage.

Conclusive Remarks
The 2014-2015 eruption at Fogo volcano developed a ~5.9 km long lava tube system that started from the eruptive fissure that became sealed in just six days after the start of the eruption.The master tube soon propagated within the arterial Flow 2, following the line of maximum slope and probably also the smooth paved road of access to the Cha Caldera.When Flow 2 stopped against the N Bordeira cliff, the lava flow field inflated and eventually led to vent opening at the E margin of its front with lava flows that destroyed the villages of Portela and Bangaeira.At this stage, further lengthening of the lava flow field was probably prevented by surface crust cooling.Thus, the lava feeding the eruption accumulated upslope, inflating the proximal flow field and causing a fast spreading of Flow 3 during the second half of December.This is also when the difference between SR and TADR increased significantly (Figures 5 and 6), causing inflation and endogenous growth of the lava flow field.When Flow 3 also halted, we recorded a growing difference between SR and TADR and a significant inflation of the lava flow field.However, the low SR of this stage and the containment offered by the cooling crust did not allow further expansion of the lava flow field, and only proximal minor pahoehoe lobes increased the thickness of the lava flow field.We have estimated that, by the end of the eruption, ~31% of the lava emplaced endogenously.We suggest that early detection of flow field inflation might be invaluable for recognizing the development of potentially devastating flows in time to mitigate their destructive effects.

Figure 1 .
Figure 1.Google Maps with location of: (a) Fogo Island (yellow square/rectangle) at Cape Verde Archipelago, western Africa; (b) Cha Caldera and Pico do Fogo volcano (inside the red rectangle); and (c) Cha Caldera depression and the paved road (yellow line) that crosses it from SE to NW, Pico do Fogo volcano, and the Portela and Bangaeira villages.The blue dot shows the fixed camera location.The 1995 lava flow field is from [39].

Figure 1 .
Figure 1.Google Maps with location of: (a) Fogo Island (yellow square/rectangle) at Cape Verde Archipelago, western Africa; (b) Cha Caldera and Pico do Fogo volcano (inside the red rectangle); and (c) Cha Caldera depression and the paved road (yellow line) that crosses it from SE to NW, Pico do Fogo volcano, and the Portela and Bangaeira villages.The blue dot shows the fixed camera location.The 1995 lava flow field is from [39].

Figure 2 .
Figure 2. (a) Map displaying the most important stages of growth of the lava flow field.The legend shows in different colors the stages of growth of the lava flow field, as well as the eruptive fissure (red dotted line), Mount Beco, Portela, Bangaeira, Vent 1, and the names of the three arterial flows: Flow 1, Flow 2, and Flow 3. (b-g) Satellite images showing in red the active portion of the lava flow field, with the outline of the lava flow field in yellow.The destroyed (dotted black line) and residual portion (black line) of the access road to the Cha Caldera is also shown.

Figure 2 .
Figure 2. (a) Map displaying the most important stages of growth of the lava flow field.The legend shows in different colors the stages of growth of the lava flow field, as well as the eruptive fissure (red dotted line), Mount Beco, Portela, Bangaeira, Vent 1, and the names of the three arterial flows: Flow 1, Flow 2, and Flow 3. (b-g) Satellite images showing in red the active portion of the lava flow field, with the outline of the lava flow field in yellow.The destroyed (dotted black line) and residual portion (black line) of the access road to the Cha Caldera is also shown.

Figure 3 .
Figure 3. Lava tube formation and growth during the early stages of the eruption shown by photo and thermal images.(a) Photo taken from SW (Mt.Beco, see Figure 1) on 27 November 2014 showing the opening of Vent 1 at the base of the eruptive fissure; and (b) thermal image collected on the same day, showing lava overflows from the S crater rim feeding several branches overlapping Flow 1. Powerful Strombolian activity occurred from seven vents along the fissure.(c) Thermal image displaying a decreased explosive activity from two active vents within the cinder cone, several overflows from the S crater rim overlapping Flow 1, and Vent 1 at the base of the eruptive fissure.(d) Composition of two overlapped thermal images, displaying a significant erosion of the S crater rim by lava overflows (right of the image), and a small tumulus formed by lava accumulation around Vent 1 (left of the image).(e) Thermal image collected on 30 November 2014, displaying the inactive overflows from the crater rim, and Vent 1 at the base of the fissure forming a much higher tumulus than the previous day, with several skylights (white dots) displaying the path of the proximal lava tube.Note the much smaller range of the thermal scale (100-400 °C) to the right of this image when compared to the scales of previous thermal images (100-700 °C), suggesting a much cooler upper lava flow field.

Figure 3 .
Figure 3. Lava tube formation and growth during the early stages of the eruption shown by photo and thermal images.(a) Photo taken from SW (Mt.Beco, see Figure 1) on 27 November 2014 showing the opening of Vent 1 at the base of the eruptive fissure; and (b) thermal image collected on the same day, showing lava overflows from the S crater rim feeding several branches overlapping Flow 1. Powerful Strombolian activity occurred from seven vents along the fissure.(c) Thermal image displaying a decreased explosive activity from two active vents within the cinder cone, several overflows from the S crater rim overlapping Flow 1, and Vent 1 at the base of the eruptive fissure.(d) Composition of two overlapped thermal images, displaying a significant erosion of the S crater rim by lava overflows (right of the image), and a small tumulus formed by lava accumulation around Vent 1 (left of the image).(e) Thermal image collected on 30 November 2014, displaying the inactive overflows from the crater rim, and Vent 1 at the base of the fissure forming a much higher tumulus than the previous day, with several skylights (white dots) displaying the path of the proximal lava tube.Note the much smaller range of the thermal scale (100-400 • C) to the right of this image when compared to the scales of previous thermal images (100-700 • C), suggesting a much cooler upper lava flow field.

20 Figure 4 .
Figure 4. Photos of Fogo: (a) The 2014 eruptive fissure and cinder cone, taken from West on 23 November 2015, with at its base a frozen lava pond where Vent 1 was located.In the foreground is a skylight with a tube (visible in c) at its base, about 20 m below the surface.Pico do Fogo is in the background.(b) View of the skylight as taken from east and looking down flow, with the pahoehoe surface coating the inner walls of the skylight, displaying final drainage of lava within the depression.(c) Lava tube, about 2-3 m wide, at the bottom of the skylight shown in a and b.

Figure 4 .
Figure 4. Photos of Fogo: (a) The 2014 eruptive fissure and cinder cone, taken from West on 23 November 2015, with at its base a frozen lava pond where Vent 1 was located.In the foreground is a skylight with a tube (visible in (c)) at its base, about 20 m below the surface.Pico do Fogo is in the background.(b) View of the skylight as taken from east and looking down flow, with the pahoehoe surface coating the inner walls of the skylight, displaying final drainage of lava within the depression.(c) Lava tube, about 2-3 m wide, at the bottom of the skylight shown in a and b.

20 Figure 5 .
Figure 5. Map with outline of the whole lava flow field (red line) mainly comprising aa lava flow surfaces, with the blue area indicating the distribution of pahoehoe lava flows.See text for further explanation.Flow 1 made of aa (arterial flow volume-controlled); Flow 3 made of pahoehoe (secondary flow fed by drainage of Flow 2); and Flow 2 complex flow made of aa (initial arterial flow) + pahoehoe (at the flow margins, mainly by emptying of the flow interior) + proximal pahoehoe by late stage (degassed) flow emplacement.The yellow dotted line indicates the eruptive fissure, the black line the access road to Cha Caldera (dotted where covered by the lava flows).

Figure 5 .
Figure 5. Map with outline of the whole lava flow field (red line) mainly comprising aa lava flow surfaces, with the blue area indicating the distribution of pahoehoe lava flows.See text for further explanation.Flow 1 made of aa (arterial flow volume-controlled); Flow 3 made of pahoehoe (secondary flow fed by drainage of Flow 2); and Flow 2 complex flow made of aa (initial arterial flow) + pahoehoe (at the flow margins, mainly by emptying of the flow interior) + proximal pahoehoe by late stage (degassed) flow emplacement.The yellow dotted line indicates the eruptive fissure, the black line the access road to Cha Caldera (dotted where covered by the lava flows).

Figure 6 .
Figure 6.Graph showing maximum lava flow lengths (in meters) vs. time (dd/mm/yyyy) for the three arterial lava flows that developed during the Fogo 2014-2015 effusive eruption, Flow 1, Flow 2, and Flow 3. Note that Flow 1 increased its length only until the end of November; Flow 2 had two major growing phases, at the end of November and in early December; Flow 3 was mainly active during the second half of December, when the other two flows had largely halted.The difference in volume between the TADR measured by satellite and the Supply Rate calculated based on the SO2 flux measurements (blue bar), with uncertainty bars (gray and pale blue bars), is also shown, with values increasing especially from early December onwards.

Figure 6 .
Figure 6.Graph showing maximum lava flow lengths (in meters) vs. time (dd/mm/yyyy) for the three arterial lava flows that developed during the Fogo 2014-2015 effusive eruption, Flow 1, Flow 2, and Flow 3. Note that Flow 1 increased its length only until the end of November; Flow 2 had two major growing phases, at the end of November and in early December; Flow 3 was mainly active during the second half of December, when the other two flows had largely halted.The difference in volume between the TADR measured by satellite and the Supply Rate calculated based on the SO 2 flux measurements (blue bar), with uncertainty bars (gray and pale blue bars), is also shown, with values increasing especially from early December onwards.

20 Figure 7 .
Figure 7. Graph showing Time Averaged Discharge Rate (TADR, purple diamonds) retrieved by HOTSAT and Supply Rate derived from SO2 data (green diamonds) vs. time (dd/mm/yyyy) during the Fogo 2014-2015 effusive eruption.Red and orange solid lines and relative errors (dashed lines) show volumes obtained by integrating TADR and Supply Rate curves, respectively.
and Figure 3e.Tubes on steep slopes and along eruptive fissures form more rapidly, because local changes in discharge rate can easily trigger inward accretion of levees and roofing over of lava channels [3,5].

Figure 7 .
Figure 7. Graph showing Time Averaged Discharge Rate (TADR, purple diamonds) retrieved by HOTSAT and Supply Rate derived from SO 2 data (green diamonds) vs. time (dd/mm/yyyy) during the Fogo 2014-2015 effusive eruption.Red and orange solid lines and relative errors (dashed lines) show volumes obtained by integrating TADR and Supply Rate curves, respectively.

20 Figure 8 .
Figure 8. Interpretative model.Outline of the final lava flow field (red line) with path of lava tubes (yellow dotted line).The position of V1 at the base of the eruptive fissure is also shown.

Figure 8 .
Figure 8. Interpretative model.Outline of the final lava flow field (red line) with path of lava tubes (yellow dotted line).The position of V1 at the base of the eruptive fissure is also shown.