Abstract
The Azores archipelago is composed of nine volcanic islands located at the triple junction between the North American, Eurasian, and Nubian plates. Nowadays the volcanic activity in the archipelago is characterized by the presence of secondary manifestations of volcanism, such as hydrothermal fumaroles, thermal and cold CO2-rich springs as well as soil diffuse degassing areas, and low magnitude seismicity. Soil CO2 degassing (concentration and flux) surveys have been performed at Pico, Faial, and São Jorge islands to identify possible diffuse degassing structures. Since the settlement of the Azores in the fifteenth Century these three islands were affected by seven onshore volcanic eruptions and at least six destructive earthquakes. These islands are crossed by numerous active tectonic structures with dominant WNW-ESE direction, and less abundant conjugate NNW-SSE trending faults. A total of 2,855 soil CO2 concentration measurements have been carried out with values varying from 0 to 20.7 vol.%. Soil CO2 flux measurements, using the accumulation chamber method, have also been performed at Pico and Faial islands in the summer of 2011 and values varied from absence of CO2 to 339 g m−2 d−1. The highest CO2 emissions were recorded at Faial Island and were associated with the Pedro Miguel graben faults, which seem to control the CO2 diffuse degassing and were interpreted as the pathways for the CO2 ascending from deep reservoirs to the surface. At São Jorge Island, four main degassing zones have been identified at the intersection of faults or associated to WNW-ESE tectonic structures. Four diffuse degassing structures were identified at Pico Island essentially where different faults intersect. Pico geomorphology is dominated by a 2,351 m high central volcano that presents several steam emissions at its summit. These emissions are located along a NW-SE fault and the highest measured soil CO2 concentration reached 7.6 vol.% with a maximum temperature of 77°C. The diffuse degassing maps show that anomalous CO2 degassing areas are controlled essentially by the tectonic structures and the lithology of the sites since the youngest volcanic systems are characterized by very low CO2 emissions.
Introduction
Carbon dioxide (CO2) together with water vapor (H2O) and sulfur dioxide (SO2) are usually the most abundant gases released to the atmosphere during volcanic eruptions (Giggenbach, ; Fischer and Chiodini, and references therein). During non-eruptive periods, gas emissions are also frequent in many volcanic systems not only through the presence of permanent fumaroles, but also due to the existence of thermal and cold CO2-rich springs, as well as soil diffuse degassing emissions (Allard et al., ; Chiodini et al., ). In the latter, the most studied gases released through volcanic soils are usually CO2 and the radioactive gas radon (222Rn); in some cases the CO2 diffusely emitted by the soils is of similar magnitude to the CO2 released from fumaroles and crater plumes (e.g., Allard et al., ; Chiodini et al., ; Viveiros et al., 2010; Pedone et al., 2015). Mantle-derived CO2 can also be released away from active volcanoes through deep tectonic structures, as demonstrated by several studies (Chiodini et al., , ; Jung et al., ; Lee et al., 2016). In fact, in the last 30 years, CO2 spatial distribution has been used worldwide to identify hidden tectonic structures, since faults/fractures act as preferential pathways (high permeability zones) for the escape of gases from the deep crust or mantle to the surface (Giammanco et al., , ; Baubron et al., ; Hutchison et al., ; Liuzzo et al., 2015). Based on the release of soil gases in confined areas, Chiodini et al. () named the anomalous CO2 degassing areas, where hydrothermal/volcanic CO2 is released, as diffuse degassing structures (DDS), whose shape depends on morphological, geological, and structural factors, such as the topography, existence of lithological heterogeneities, and presence of faults/fractures (Schöpa et al., 2011; Peltier et al., 2012; Pantaleo and Walter, 2014).
CO2 has been widely used for volcanic monitoring due to its low solubility in silicate melts, being consequently one of the first gases released to the surface in case of replenishment of a magma chamber (e.g., Hernández et al., ; Carapezza et al., ; Giammanco et al., ; Aiuppa et al., ; Liuzzo et al., 2013; de Moor et al., ). Carbon dioxide is also an inert asphyxiant gas if present in high concentrations in the air, and above 10 vol.% can be lethal (Blong, ; Weinstein and Cook, 2005). Several hundreds of deaths have been reported in volcanic and non-volcanic environments due to deep-derived CO2 emissions (Hansell and Oppenheimer, ) and high CO2 concentrations have been reported in buildings located on CO2 anomalous zones (Baxter et al., ; Viveiros et al., 2009, 2016a). This fact highlights the importance of identifying CO2 anomalous areas by producing degassing maps, which should constitute valuable tools for land-use planners (Beaubien et al., ; Viveiros et al., 2009, 2010, 2015, 2016a). Building in high CO2 degassing areas should be avoided, or, depending on the soil gas concentrations measured, buildings should follow few “gas-resistant” construction rules as the ones defined by Viveiros et al. (2016a).
Carbon dioxide is a colorless and odorless gas that is only detected with specific equipment. Several instruments have been developed in the last decades to detect this gas in soils in volcanic/tectonic environments, even if the first studies started to be applied in agricultural fields to measure soil respiration rates (e.g., Kanemasu et al., 1974; Parkinson, 1981). Relatively short-time and low-cost methodologies using infrared CO2 detectors have been implemented to measure soil CO2 emissions since the end of the eighties, which include soil CO2 concentration measurements (Giammanco et al., , ; Evans and Staudacher, ) and soil CO2 fluxes (dynamic concentration and accumulation chamber methods; Gurrieri and Valenza, ; Chiodini et al., ). Studies performed in diffuse degassing environments show that not only CO2 is permanently emitted to the atmosphere from soils, but also significant variations in the amount of gases emitted are observed. Several spike-like and long term variations have been interpreted as resulting from meteorological influences and seasonal effects (Granieri et al., ; Viveiros et al., 2008, 2014; Rinaldi et al., 2012). Gases released from soils in hydrothermal-volcanic environments may also result from organic matter decomposition and fauna respiration, the so-called soil respiration (Luo and Zhou, 2006). Discrimination of the different CO2 sources has been done through the use of statistical methodologies (Chiodini et al., ), as well as the carbon isotopic composition of the CO2 released (e.g., Chiodini et al., ).
Soil CO2 concentration surveys were performed at Faial, Pico, and São Jorge islands in the period between 2001 and 2004, and preliminary results are available in Master theses (Faria, ; Marcos, 2006). In 2011, sporadic soil CO2 flux measurements were also performed in the islands of Faial and Pico, with the main scope of evaluating the stability of the previously recognized anomalous CO2 degassing areas and select an area to install a permanent soil CO2 flux station. This equipment is necessary to integrate with other seismo-volcanic monitoring techniques and to complement the gas geochemistry network already installed in the Azores archipelago (Viveiros et al., 2008). This study reviews the preexisting database on soil CO2 concentration applying more recent statistical tools in order to (1) evaluate the existence of different sources for the CO2 diffusely released from soils, (2) produce CO2 degassing maps, (3) identify the main DDS, and (4) correlate the soil CO2 emitted with the main volcano-tectonic structures identified in the islands.
Geological setting
The Azores archipelago is composed of nine volcanic islands located in the North Atlantic Ocean where the American, Eurasian and the Nubian plates meet at a triple junction (Searle, 1980). The main tectonic features are (1) the Mid-Atlantic Ridge (MAR) that crosses the archipelago between the islands of Flores and Faial, (2) the East Azores Fracture Zone (EAFZ), which extends E-W from the MAR to south of Santa Maria, and the Azores-Gibraltar Fault Zone that includes the E-W trending Gloria Fault, and (3) its western segment, the Terceira Rift (TR), which extends from the MAR to the island of Santa Maria along a general WNW-ESE direction, and corresponds to the present-day EU-NU plate boundary (e.g., Searle, 1980; Madeira and Ribeiro, 1990; Vogt and Jung, 2004; Miranda et al., 2015; Figure 1). The main tectonic structures in the islands are normal dextral faults with a WNW-ESE trend, characteristic of the Terceira rift; the NNW-SSE conjugate fault system exhibits oblique normal left lateral displacement (Madeira and Brum da Silveira, 2003; Madeira et al., 2015). These main fault systems are also present in Faial, Pico, and São Jorge islands, with the dominant structures striking WNW-ESE and dipping 60–90° to the NNE or to the SSW controlling the general shape of some islands (Madeira and Brum da Silveira, 2003).
Figure 1
Most of the identified faults in the Azores archipelago are considered active (Madeira, 1998; Madeira and Brum da Silveira, 2003; Marques et al., 2013; Madeira et al., 2015) and the major faults usually have well-developed scarps; however the volcanic nature of the islands, together with differential erosion, can attenuate or amplify the tectonic slope (Madeira and Brum da Silveira, 2003) and even make the identification of tectonic structures difficult (due to scarp burial by pumice deposits, for instance). In addition, the exuberant vegetation of some of the Azorean islands covers most of the geological structures. Tectonic structures used in the present study were mapped by Madeira and Brum da Silveira (2003) and Madeira et al. (2015) through fieldwork and vertical air-photo interpretation. All the faults in the study islands are considered active by Madeira and Brum da Silveira (2003) as they displace volcanic sequences younger than 100 ka.
The cause for the volcanism in the islands is highly debated in the literature and most studies suggest the presence of a mantle plume (e.g., Schilling, 1975; Cannat et al.,
The complex geodynamic setting of the Azores explains the frequent seismicity and volcanism in the islands. Since the settlement of the islands, in the fifteenth century, at least 28 volcanic eruptions and more than 15 major earthquakes caused a large number of casualties (>6,500 deaths; Gaspar et al.,
Present seismo-volcanic activity in the archipelago is characterized by almost daily record of seismic events (Gaspar et al.,
Faial Island
Tectonic structures and seismicity
The eastern part of Faial Island is dominated by a WNW-ESE trending graben structure (named Pedro Miguel Graben), which is formed by seven normal dextral faults (Madeira and Brum da Silveira, 2003; Figure 2). Pedro Miguel graben formed as a result of the activity of transtensive dextral faults with an extension rate between ~3.4 and ~8.2 mm/y was proposed by Trippanera et al. (2014). Other important WNW-ESE tectonic structures are located south of the summit caldera and in the western part of the island (Figure 2). The conjugate normal left lateral faults system (trending NNW-SSE to NW-SE) shows less developed geomorphic expression (Madeira and Brum da Silveira, 2003).
Figure 2

Faial, Pico, and São Jorge digital elevation models with the identification of the volcanic systems and the main tectonic structures. Legend of the tectonic structures: Faial Island: LG, Lomba Grande; R, Ribeirinha; CC, Chã da Cruz; RR, Ribeira do Rato; RV, Rocha Vermelha; E, Espalamaca; F, Flamengos; Pico Island: LF, Lomba de Fogo-S. João; LC, Lagoa do Capitão; T, Topo; São Jorge Island: RS, Ribeira Seca; USJ, Urze-São João. VS, Volcanic System. The red diamonds represent the approximate locations of the onshore historical volcanic centers (tectonic structures modified from Madeira and Brum da Silveira, 2003).
In the last century Faial Island was affected by four important seismic events in 1926, 1958, 1973, and 1998, with intensities equal or higher than VII (Mercalli Modified Scale). The last destructive earthquake that caused nine casualties occurred offshore in 1998, and had a magnitude of (ML) 5.8 (Dias et al.,
Volcanism
Faial Island is composed by four main volcanic systems: the Ribeirinha (848 to 358 ka) and Caldeira (>440 ka) central volcanoes and two basaltic fissure systems (Horta Platform and Capelo Peninsula fissure systems; Chovelon,
Seismic tomography performed at Faial Island suggests the existence of a low P velocity between 3 and 7 km depth beneath Caldeira Volcano, which was interpreted by Dias et al. (
Two volcanic eruptions affected the island since the settlement and both occurred in the eastern part of the island, in the Capelo Peninsula volcanic system (1672/73 and 1957/58). During this last eruptive episode, a phreatic explosion also occurred inside the caldera of the Caldeira Volcano and the fumarolic emissions associated to this event persisted from May to October 1958 (Machado et al., 1962; Madeira, 1998; Pacheco, 2001). At least three deaths were attributed to the 1672–73 volcanic eruption (Weston, 1964; Pacheco, 2001).
Nowadays, the only visible secondary manifestations of volcanism in Faial correspond to residual steam emissions from the Capelinhos eruption (1957–58) (Ferreira,
Table 1
| Island | Sampling period | Temp. (°C) | Gas composition (molar%) | References | |||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| CO2 | H2S | O2 | Ar | N2 | CH4 | He | H2 | ||||
| Pico | Jun-16 | 77.0 | 40.45 | 0.00 | 12.96 | 0.54 | 46.05 | 0.00 | 0.00 | 0.00 | This study |
| 1994 | 76.1 | 40.50 | 0.00 | 11.84 | 0.55 | 47.11 | 0.00 | n.d. | 0.00 | Nunes, 1999 | |
| Faial | Mar-93 | 91.0 | 1.31 × 10−2 | n.d. | 9.1* | 90.9 | n.d. | n.d. | 6 × 10−5 | Gaspar and Ferreira, | |
Gas composition (molar%) of the fumarolic emissions from Faial and Pico islands.
Represents the sum of the O2 and Ar due to limitations of the analytical procedures.
Pico Island
Tectonic structures and seismicity
Several authors consider that Faial and Pico are the emerged parts of a single main volcanic ridge, the Pico-Faial ridge (Hildenbrand et al.,
The main tectonic structures of Pico Island trend WNW-ESE, similarly to what is observed in Faial Island. The main structures are Lagoa do Capitão and Topo faults that progressively merge to the east and both define a graben structure, smaller than Pedro Miguel graben in Faial Island. Several other WNW-ESE and conjugated NNW-SSE faults are defined by volcanic alignments (Madeira, 1998; Nunes, 1999; Madeira and Brum da Silveira, 2003; Figure 2). The larger cones and Pico Volcano are located in the intersection of WNW-ESE with NNW-SSE conjugate faults (Madeira and Brum da Silveira, 2003).
A magnitude VII earthquake (EMS-98, Silva, 2005) that occurred in 1973 was the strongest earthquake with epicenter in Pico Island. However, and due to the proximity of the studied islands, strong earthquakes that occurred in Faial, São Jorge, and in the surrounding offshore areas also caused damage at Pico. This is the case of the 1757 earthquake that had epicenter close to São Jorge Island and severely damaged Pico causing 11 casualties. Damage on Pico was also reported associated with the seismic events of 1926 and 1998 that mostly affected Faial Island, and the 1980 Terceira earthquake (Silva, 2005; Gaspar et al.,
Volcanism
Pico is the youngest island of the archipelago and is formed by three volcanic systems: the Topo Volcano (186 ± 5 ka), the Planalto da Achada fissure system (270 ± 150 ka to <25 ka), and the Pico central volcano (Féraud et al.,
Three historic volcanic eruptions affected the island (1562–64; 1718 and 1720), and two deaths are reported for the 1718 eruption. Both 1562–64 and 1720 volcanic eruptions were associated with the Pico Island fissural system (Achada Plateau); the 1718 activity occurred at the Pico Volcano. In 1963 a submarine eruption was reported offshore the NW coast of Pico Island (Madeira, 1998; Nunes, 1999; Gaspar et al.,
Fumarolic emissions at Pico Island are located on the summit of the 2,351 m-tall Pico Volcano and a survey performed in 1994 showed that the dry gas is composed of CO2, O2, Ar, and N2 (Table 1; Nunes, 1999). MORB-type values for the 3He/4He ratio (8.5 Ra) were measured by Jean-Baptiste et al. (
São jorge Island
Tectonic structures and seismicity
The main tectonic structures defined for São Jorge are parallel to the WNW-ESE elongation of the island (Figure 2). The younger western half of São Jorge is dominated by the Picos and Pico do Carvão fault zones, which are marked by alignments of cones, craters, scoria ramparts, and short fault scarps (Madeira, 1998; Madeira and Brum da Silveira, 2003; Madeira et al., 2015). In the older eastern part of the island, Madeira and Brum da Silveira (2003) identified a set of faults with WNW-ESE direction. The most important faults in this region are the Urze-S. João Fault that displays a 10 km long continuous scarp (Figure 2). The two sectors of the island are separated by the NNW-SSE Ribeira Seca hidden fault according to Madeira and Brum da Silveira (2003).
The strongest earthquake (M = 7.4) that affected the Azores islands after settlement occurred offshore São Jorge in 1757 and was responsible for the death of ~1,000 persons, corresponding to about 20% of the population of the island at that time (Madeira, 1998; Madeira and Brum da Silveira, 2003; Gaspar et al.,
Volcanism
Three main fissure basaltic volcanic systems have been defined in São Jorge Island: Topo (1.32 ± 0.02 Ma), Rosais (368 ± 6 ka to 270 ka), and Manadas (<6 ka years) volcanic systems (Hildenbrand et al.,
Since settlement two sub-aerial eruptions occurred in the island, one in 1580 and other in 1808; persistent fumarolic activity was reported associated to the main vents (north of Urzelina village) for several years following the 1808 eruption, but nowadays no evidence of these emissions is possible to find out; as previously mentioned a submarine eruption occurred in 1964 associated to intense seismicity (Weston, 1964; Madeira, 1998). Even if the historical volcanic eruptions at São Jorge Island are dominated by Hawaiian/Strombolian eruptive styles, more than 45 deaths were reported as a consequence of the volcanic activity (Gaspar et al.,
No thermal anomalies are reported in the island and two cold CO2-rich springs located in the eastern half of the island are the only secondary manifestations of volcanism presently found at São Jorge (Viveiros, 2003; Ferreira et al.,
Materials and methods
Soil diffuse degassing
Field surveys and instruments
Soil CO2 concentrations were measured at a depth of about 50 cm through the insertion of an iron probe in the soil. The probe is connected to an infrared CO2 detector (Geotechnical Instruments, model Anagas CD95) and the gas is pumped to the instrument until the concentration reaches a stable value. This detector measures the CO2 concentration in the range 0–100 vol.% (resolution of 0.1 and precision of 2%).
The Faial Island survey was carried out on two campaigns: during the summer of 2001 (1712 measurements; Faria,
Figure 3

Sampling location and soil CO2 concentration distribution maps for the studied islands: (A) Faial, (B) Pico, and (C) São Jorge islands.
Soil CO2 fluxes were also measured at Pico and Faial islands during August 2011. These measurements were performed using the accumulation chamber method (Chiodini et al.,
A total of 93 and 107 sites were surveyed, respectively, at Pico and Faial islands (Table 2, Figure 4). Considering the significant impact that environmental conditions can have on the soil gas fluxes (e.g., Granieri et al.,
Table 2
| Variables | Statistics | Faial Island | Pico Island | São Jorge Island |
|---|---|---|---|---|
| Soil CO2 concentration | Minimum (vol.%) | 0.0 | 0.0 | 0.0 |
| Average (vol.%) | 1.3 | 1.5 | 0.9 | |
| Maximum (vol.%) | 20.7 | 13.3 | 13.4 | |
| Standard deviation (vol.%) | 1.7 | 2.4 | 1.4 | |
| Number of points | 2,157 | 316 | 382 | |
| Soil CO2 flux | Minimum (g m−2 d−1) | 0 | 0 | – |
| Average (g m−2 d−1) | 170 | 22 | – | |
| Maximum (g m−2 d−1) | 339 | 40 | – | |
| Standard deviation (g m−2 d−1) | 38 | 9 | – | |
| Number of points | 107 | 93 | – |
Descriptive statistics of the soil CO2 concentration/flux data.
Figure 4

Sampling location and soil CO2 flux distribution maps for Faial (A) and Pico (B) islands.
Field conditions (e.g., geomorphologic features, anthropogenic infrastructures, authorization to access the sites) affected the sampling design resulting in an irregular grid (Figures 3, 4). Nevertheless, sampled sites were distributed throughout the entire islands in order to cover as much as possible the sub-aerial surface of the islands. Locations of the sampling points were determined by a handheld Global Positioning System (GPS, GPSMAP® 76S from Garmin Company) with uncertainty of ±6 m.
Soil temperature measurements have been performed in the Faial (Capelinhos) and Pico (Pico Volcano) fumaroles using a portable thermocouple (Hanna Instruments, model HI93531). Measurements in Faial were carried out in August 2011 and in Pico in July 2003 (Table 3).
Table 3
| Island | UTM M | UTM P | Soil temperature (°C) | Average | Sampling period |
|---|---|---|---|---|---|
| Faial | 340730 | 4274107 | 64.3 | 59.0 | 02/08/2011 |
| 340730 | 4274110 | 66.8 | 02/08/2011 | ||
| 340734 | 4274152 | 79.8 | 02/08/2011 | ||
| 340732 | 4274153 | 77.5 | 02/08/2011 | ||
| 340736 | 4274164 | 48.6 | 02/08/2011 | ||
| 340733 | 4274173 | 50.4 | 02/08/2011 | ||
| 340736 | 4274186 | 49.1 | 02/08/2011 | ||
| 340539 | 4274004 | 35.2 | 02/08/2011 | ||
| Pico | 377894 | 4258704 | 48.6 | 52.5 | 01/07/2003 |
| 377896 | 4258714 | 50.4 | 01/07/2003 | ||
| 377890 | 4258727 | 54.5 | 01/07/2003 | ||
| 377897 | 4258733 | 74.2 | 01/07/2003 | ||
| 377892 | 4258739 | 57.9 | 01/07/2003 | ||
| 377883 | 4258738 | 52.1 | 01/07/2003 | ||
| 377891 | 4258705 | 36.8 | 01/07/2003 | ||
| 377889 | 4258706 | 51.9 | 01/07/2003 | ||
| 377893 | 4258714 | 46.0 | 01/07/2003 |
Soil temperature measured in the fumarolic areas of Pico and Faial islands.
Data processing
CO2 released by soils can be fed by multiple gas sources, such as biogenic and volcanic-hydrothermal origins (Chiodini et al.,
Graphical statistical approach (GSA)
Cumulative probability plots have been widely applied in the literature to distinguish different CO2 populations that can be representative of distinct CO2 sources (biogenic and volcanic-hydrothermal; e.g., Chiodini et al.,
The 95th percentile of the lowest values population, which usually represents the biological CO2, is used as cutt-off for the biological (background) CO2 fluxes. A similar criteria was used in previous studies (e.g., Chiodini et al.,
Interpolated maps
Final CO2 degassing maps have been elaborated based on the sequential Gaussian simulations (sGs) methodology (Deutsch and Journel,
Figure 5

Omnidirectional variograms of the soil CO2 concentration data for Faial (A), Pico (B), and São Jorge (C). Cc: Partial sill; a: range (in m).
Fumarolic emissions
Sampling methodology
Gases released at Pico Volcano fumaroles were collected in June 2016 using the methodology defined by Giggenbach (
Analytical procedures
The analysis of the chemical composition of the gases was carried out at the gas geochemistry laboratories of the University of the Azores, using gas chromatography and titration techniques. Gases in the headspace of the bottle (CH4, N2, O2, Ar, He, and H2) were analyzed with a Perkin Elmer Clarus 580 gas chromatograph. This chromatograph has two channels equipped with two Thermal Conductivity Detectors and both a MS 5A plot column and a MS packed column that allow using He and Ar as carrier gases to separately quantify Ar and O2. Gases dissolved in the alkaline suspension were quantified by titration: the CO2 was detected by potentiometric titration with an automatic titrator from Radiometer Copenhagen, model VIT90 Video Titrator and the H2S quantified by colorimetric titration with mercury acetate using dithizone for end point detection.
Results
Descriptive statistics
Soil CO2 concentration ranged between 0 and 20.7 vol.%, with the maximum values recorded at Faial Island (Table 2). Higher soil CO2 fluxes were also measured at Faial Island (339 g m−2 d−1). Maximum soil CO2 concentrations in Pico and São Jorge were lower (~13 vol.%; Figure 3).
Soil temperatures measured in the Pico and Faial fumaroles and in the surrounding areas are presented in Table 3 and maximum values were similar, i.e., 74.2 and 79.8°C, for Pico and Faial degassing areas, respectively. Higher temperatures were measured in the fumarolic/steam vents (Table 1).
Gas composition of the Pico fumarole, sampled in June 2016, is showed in Table 1. These emissions release essentially water vapor (~0.49 for the ratio dry gases/H2O), N2, CO2, O2, and Ar. The other elements, characteristic of hydrothermal fumaroles (CH4, He, H2, and H2S), were not detected.
CO2 degassing maps
Omnidirectional variograms of the soil CO2 concentration data were computed and modeled for each dataset (Figure 5). Degassing maps of the CO2 concentrations resulting from the sequential Gaussian simulation procedure for the sampled datasets are shown in Figures 6–8. The scarce data for the soil CO2 flux values do not allow the elaboration of an interpolated map, and for this reason only the distribution of the sampled points is showed (Figure 4).
Figure 6

E-type soil CO2 concentration map for Faial Island (cell size = 50 × 50 m; interpolation method: sGs). Numbers F1–F7 represent the DDS identified at Faial Island, which correspond to areas where CO2 is fed by hydrothermal sources (highlighted as orange and red—soil CO2 concentration >4 vol.%) (UTM(m)-WGS84, zone 26S). The blue arrow points to the location of the Capelinhos fumarolic field and red dotted lines represent degassing lineaments. Blue letters represent tectonic structures mentioned in the text: LB, Lomba de Baixo; LM, Lomba do Meio; F, Flamengos; E, Espalamaca; LG, Lomba Grande; CC, Chã da Cruz; R, Ribeirinha.
Figure 7

E-type soil CO2 concentration map for Pico Island (cell size = 50 × 50 m; interpolation method: sGs). Numbers P1–P4 represent the DDS identified at Pico Island (highlighted as orange and red colors) (UTM(m)-WGS84, zone 26S). The blue arrow represents the location of the Pico fumarolic field and red dotted line represents a degassing lineament. Blue letters represent tectonic structures mentioned in the text: LF, Lomba de Fogo-S. João; LC, Lagoa do Capitão; T, Topo.
Figure 8

E-type soil CO2 concentration map for São Jorge Island (cell size = 50 × 50 m; interpolation method: sGs). Numbers SJ1–SJ4 represent the DDS identified at São Jorge Island (highlighted as orange and red colors) (UTM(m)-WGS84, zone 26S). The red dotted line represents a degassing lineament. Blue letters represent tectonic structures mentioned in the text: PC, Pico do Carvão; RS, Ribeira Seca; USJ, Urze-S. João.
CO2 populations
Logarithmic probability plots of the soil CO2 flux/concentrations are presented in Figures 9, 10. Soil CO2 concentration shows bimodal populations for the data measured in the three islands, suggesting the presence of different CO2 sources (biogenic and volcanic-hydrothermal). In what concerns the soil CO2 flux data, polymodal distribution of the Faial Island data represents the overlapping of three log-normal populations (Figure 10).
Figure 9

Probability plots of soil CO2 concentrations measured at Faial (A), Pico (B), and São Jorge (C) islands. Green and red lines represent populations A and B, respectively.
Figure 10

Probability plots of soil CO2 fluxes measured at Faial (A) and Pico (B) islands. Green, orange and red lines represent populations A, B, and C, respectively.
Statistical parameters from the partitioned CO2 populations and the 90% confidence intervals of the mean for the different study sites are presented in Table 4. Populations named “A” are characterized by very low soil CO2 flux/concentration values. Population “C” is identified only for the case of Faial Island soil CO2 flux data.
Table 4
| Island | Variable | Population | CO2 source | Proportion of each population (%) | Mean CO2 | Mean CO2 90% confidence interval |
|---|---|---|---|---|---|---|
| Faial | CO2 flux (g m−2 d−1) | A | Mainly biogenic | 15 | 0.6 | – |
| B | Biogenic + hydrothermal | 81 | 22.2 | 20.44–24.51 | ||
| C | Hydrothermal | 4 | 537.5 | 159–421,776 | ||
| CO2 concentration (vol.%) | A | Mainly biogenic | 75 | 1.2 | 1.1–1.3 | |
| B | Biogenic + hydrothermal | 25 | 1.5 | 1.4–1.7 | ||
| Pico | CO2 flux (g m−2 d−1) | A | Mainly biogenic | 7 | 4.3 | 2.50–16.45 |
| B | Biogenic + hydrothermal | 93 | 22.7 | 21.36–24.40 | ||
| CO2 concentration (vol.%) | A | Mainly biogenic | 65 | 0.6 | 0.3–2.8 | |
| B | Biogenic + hydrothermal | 35 | 3.2 | 2.8–3.8 | ||
| São Jorge | CO2 concentration (vol.%) | A | Mainly biogenic | 65 | 0.7 | 0.5– 1.1 |
| B | Biogenic + hydrothermal | 35 | 1.3 | 1.2–1.5 |
Statistical parameters from the partitioned CO2 concentration/flux populations and the 90% confidence intervals of the mean.
Discussion
Visible secondary manifestations of volcanism were previously recognized in the islands of Faial, Pico, and São Jorge (Ferreira et al.,
The measured ranges of CO2 suggest different origins for the gases released from soils. Statistical methodologies applied to the recorded data showed different populations both for the CO2 concentrations and fluxes, suggesting biogenic and volcanic-hydrothermal sources. Discrimination of these sources is particularly relevant in the Azores archipelago due to the abundant vegetation that covers the volcanic soils. Lowest CO2 values (population “A,” Table 4) probably correspond to areas where biologic processes are responsible for the gas released from soils. Values very close to zero, or even zero, should represent areas where even the biologic production is scarce. Population “C,” only present in the soil CO2 flux dataset of Faial Island, is representative of CO2 fed by an endogenous source (volcanic-hydrothermal origin). Data associated with population “B” probably represent the mixture of biogenic and hydrothermal CO2 contributions. Most of the measured soil CO2 concentration data fall within populations “A” (65–75% of the data), corresponding to biogenic contribution for the emitted CO2 (Table 4). In what concerns the soil CO2 flux determinations, most of the data are interpreted as a mixture of biogenic and hydrothermal/volcanic contributions. The high percentage of CO2 flux data in population “B” is probably explained by the fact that the flux surveys were carried out in the areas where anomalous CO2 values were previously detected and mostly related with the presence of faults. The main objective of the 2011 campaigns was to select an area to install permanent soil CO2 flux stations in the islands of Pico and Faial.
The cut-off value defined for each dataset according to the criteria suggested by Sinclair (1974), i.e., the 95th percentile of a determined population, allows to establish thresholds to discriminate the biogenic contribution from a deep source (Table 5). The 95th percentile for population “B” is selected in this study to define the DDS (Chiodini et al.,
Table 5
| Variables | Statistics | Faial Island | Pico Island | São Jorge Island |
|---|---|---|---|---|
| Soil CO2 concentration (vol.%) | 95th percentile (population A) | 4.0 | 2.1 | 3.2 |
| 95th percentile (population B) (DDS threshold) | 4.0 | 4.9 | 3.5 | |
| Soil CO2 flux (g m−2 d−1) | 95th percentile (population A) | 12.6 | 12.6 | – |
| 95th percentile (population B) (DDS threshold) | 39.8 | 39.8 | – |
Values for the biogenic and DDS thresholds following the criteria of Sinclair (1974).
The thresholds proposed here are similar to data available in the literature. For instance, available data on CO2 production from a wide variety of ecosystems show CO2 fluxes ranging from ~0.5 to ~19 g m−2 d−1 (e.g., Raich and Schlesinger, 1992; Raich and Tufekcioglu, 2000) with maximum flux values attributed to grassland (~50 g m−2 d−1; e.g., Norman et al., 1992; Bajracharya et al.,
The CO2 degassing map for Faial Island shows an absence of CO2 emission in the youngest part of the island (Capelo Volcanic Complex), including the areas where the recent volcanic eruptions occurred (Figures 2, 6). Even in areas where thermal anomalous values were measured, which are associated to the only fumarolic (steam) field in Faial, in Capelinhos Volcano, no soil CO2 degassing (both concentration and flux) was detected. In the fumarolic field the maximum temperature measured was 91°C and the main gases released are the atmospheric constituents (N2, O2, and Ar), with minor amounts of CO2 and H2. These emissions are probably due to some remnant heat from the 1957–58 volcanic eruption that still heats some water and air in the system. Higher soil CO2 concentrations were measured in the floor of the summit caldera (Caldeira Volcano; DDS F1, Figure 6), where a phreatic eruption occurred during the Capelinhos 1957–58 volcanic eruption, however no thermal anomaly was identified in the area.
Two main fault systems trending WNW-ESE and NNW-SSE were recognized in the study islands by previous works (Madeira and Brum da Silveira, 2003; Madeira et al., 2015), and most of the DDS identified in this study occur associated with the previously mapped tectonic structures. A clear correlation between distribution of the CO2 degassing anomalies and one of the main trending pattern is not observed; instead, and similarly to the observed for volcanism that is tectonically controlled (Madeira et al., 2015 and references therein), the main degassing anomalies seem to be found out in the intersection of different fault systems.
Based on soil CO2 concentrations coupled with the soil CO2 flux values, seven DDS were then identified in the Faial island (Figure 6). Main CO2 anomalies are located in the eastern side of the island, essentially delimited by the Pedro Miguel graben structure that affects products from the Ribeirinha and Caldeira central volcanoes. The large DDS F3, which comprises more than half of the graben floor, lays between the important WNW-ESE tectonic structures that define the Pedro Miguel graben: between the Lomba Grande Fault (LG, Figure 2) that dips to the south and the north-dipping Espalamaca and Flamengos faults that define the southern part of the graben structure. The presence of deeply-derived CO2 in the Pedro Miguel graben suggests that the graben faults correspond to deep structures that link the surface to deep areas where CO2 is trapped and then channeled to the surface. DDS F5 probably results from the intersection of the graben faults with the mapped NNW-SSE structure identified in the area. In addition, this DDS contributes to define a degassing lineament (1′, Figure 6) that can be the prolongation of the Chã da Cruz Fault.
The graben faults do not constitute barriers to the gas traveling to the surface since diffuse degassing structures are not limited to the graben floor and, for instance, DDS F2 is identified in the upthrown block of the Flamengos Fault. Two degassing lineaments are also associated with this DDS, which argue to confirm the presence of the hidden N-S radial fracture (lineament 2′; Figure 6).
The DDS F6 is associated with Flamengos and Lomba de Baixo faults, which dip to the north. The main CO2 anomaly is however observed in the top of the scarp of Lomba de Baixo Fault and should result from the high permeability of the area that is crossed by several tectonic structures. Similarly, DDS F7 is also found out in the intersection of different trending tectonic structures (WNW-ESE, NNW-SSE, and NE-SW). As mentioned before, the CO2 anomalous areas occur, in general, where a network of faults/fractures are identified, and thus seem to result from the intersection of permeable structures that allow the gas to travel from depth and escape to the surface (see DDS F2, F5, F6, and F7 in Figure 6).
The highest soil CO2 flux and concentration values, representing a clear volcanic-hydrothermal origin, were measured at Praia do Almoxarife village (DDS F4) and were detected close to a building where in 2007 lethal indoor CO2 concentrations (>15 vol.%) were measured (Ferreira and Viveiros,
From a public health risk assessment perspective, the DDS F4, located at Praia do Almoxarife, highlights the need to perform a detailed CO2 flux map in that village to identify which buildings may be located over anomalous CO2 degassing sites, with the consequent impacts that it may have on the population. Even if the other DDS are not directly associated with inhabited areas, detailed CO2 degassing maps should also be performed to evaluate possible risks for future constructions.
A total of four DDS were identified at Pico Island. The main CO2 degassing anomalies are located in the summit area of the youngest volcano (DDS P2), Pico Volcano, where the only fumarolic emissions of the island are found (Figure 7). DDS P1 and P2 are associated with the Lomba de Fogo-S. João NNW-SSE trending fault that crosses the summit of Pico Volcano and fed the 1718 eruption. Similarly to the central volcano, the fumarolic field occurs in the intersection of the NNW-SSE, WNW-ESE, and NE-SW structures. Chemical analyses performed on Pico fumaroles in the summer of 2016 are similar to those previously obtained in 1994, showing the stability of its feeding system across two decades. The gases emitted are mainly composed of N2 and CO2, with minor amounts of O2 and Ar, whose composition suggests air contamination during the ascend of gas from the deep reservoir to the surface. This fact is probably correlated with the location of the fumaroles in the volcano's summit, allowing atmospheric air to be introduced in the system, interfering with the original gas composition in the hydrothermal reservoir.
Despite the DDS P1 and P2 the remaining anomalous CO2 degassing sites are also associated with tectonic structures previously identified by Madeira and Brum da Silveira (2003) in older volcanic areas. These important CO2 degassing anomalies are correlated with the Lagoa do Capitão and Topo areas (DDS P3 and DDS P4, respectively, Figure 7), where several faults intersect. DDS P3 is correlated with the important WNW-ESE trending Lagoa do Capitão Fault. This structure dips to the south, but CO2 anomalies are observed also in the upthrown block of the fault. This behavior may be associated with a higher permeability in the top of the scarp (due to fracture along the scarp) or with the presence of not mapped tectonic structures, such as the suggested lineament 1' that can be the elongation of the probable NNW-SSE fault defined by Madeira and Brum da Silveira (2003). In what concerns DDS P4, identified in the southwestern side of the WNW-ESE Topo trending fault, it probably results from the intersection of hidden NNW-SSE and NE-SW tectonic structures (Figure 7).
Similarly to Faial, faults/fractures intersection seems to favor higher permeability for the gas to travel from depth to the surface. Madeira and Brum da Silveira (2003) suggested several probable/hidden faults for the island trending WNW-ESE and NW-SE, and the present study shows that significant CO2 anomalies correlate well with the areas of intersection of hidden/probable NNW-SSE structures with the dominant WNW-ESE faults.
The areas associated with historical eruptions in Pico Island (1562–64, 1718 and 1720, Figure 2), similarly to what was observed for Faial (Capelo Peninsula), do not show CO2 emissions, what can be probably explained by the low or absent coverage of soils together with the high permeability of the recent lava flows that allow ingress of atmospheric air in the sub-superficial layers and dilution of the deep gases in their traveling to the surface. Nevertheless, it cannot be excluded that lava flows may act locally as barriers to the gas ascent to the surface when they are compact and low fractured. This behavior was also observed by Padrón et al. (2013) at Timanfaya Volcano (Lanzarote, Canary Islands), where the absence of CO2 emission along the eruptive fissures was justified by the capping nature of the lava flow fields of the 1730–36 eruptions. Even in the vent areas or in the eruptive fissures the CO2 degassing is low. Studies performed by Giammanco et al. (
From a hazard perspective, none of the main anomalous CO2 sites are located in areas with buildings; nevertheless, and considering the distances between sampling sites, a detailed survey should be carried out in the villages located close to the identified DDS.
The comparison of the two surveys carried out at Pico fumarolic emissions shows that the gas composition remained quite stable in the last two decades; this fact argues for the stability of the feeding systems, and consequently can be useful for any seismo-volcanic monitoring programme.
Four main DDS were identified in São Jorge Island, where no thermal anomalous zones are known. SJ3 is the only DDS located in the intersection of different WNW-ESE trending faults (Figure 8) and probably corresponds to an area of increased permeability that allows the gas to escape at the surface. The DDS SJ1 is aligned with a WNW-ESE trending probable fault located to the west and east of the degassing anomaly and mapped by Madeira and Brum da Silveira (2003). The alignment of the soil CO2 degassing anomalies (1′, Figure 8), once again, may contribute to define the location of hidden faults. The main WNW-ESE fault zones that dominate the younger western part of the island do not show any CO2 degassing, what can be explained by the superimposed lithological control, as mentioned also for Pico and Faial islands. In the eastern part of the island, and along the WNW-ESE north dipping Urze-São João Fault (Madeira et al., 2015) few anomalous soil CO2 concentrations are observed, even if not directly associated with the dipping direction, and were identified as DDS SJ4.
No tectonic structure has been identified in the area surrounding the SJ2 anomaly, thus further studies are needed to better understand the permeability in the area. These future works should include soil CO2 flux surveys that could contribute to better define the DDS here proposed and eventually identify other anomalous areas. This study does not show significant anomalous CO2 emissions close to the probable Ribeira Seca Fault, and thus cannot contribute to support its existence, that is still under debate.
Considering that CO2 may act as asphyxiant, a detailed survey at Ribeira Seca village is needed for land-use planning. Three of the anomalous zones (SJ2, SJ3, and SJ4) are located in the oldest volcanic system of the island (Topo Volcanic System) and, similarly to what is observed on Pico and Faial islands, no soil CO2 degassing was detected on the areas where historical eruptions (1580 and 1808) occurred (Figure 2), fact that should be explained by the same factors as mentioned above for Faial and Pico islands (absence of well-developed soils, permeability of the lava flows and existence of natural barriers for the gas migration and release at the surface).
Degassing phenomena during periods of quiescence may be a permanent risk in any volcanic area (Viveiros et al., 2009) and, even if no incidents were reported for the islands of Pico and São Jorge, a family was relocated in Faial Island in 2007 due to high indoor CO2 concentrations (>90 vol.%) (Ferreira and Viveiros,
Concluding remarks
Anomalous soil CO2 degassing areas were identified at Faial, Pico, and São Jorge islands and a total of 15 diffuse degassing structures (DDS) were defined based on the integration of soil CO2 flux and concentration measurements. Statistical methodologies applied to the measured values show distinct populations for the CO2 emitted at the surface, which probably represent different origins (biogenic and hydrothermal-volcanic). Thresholds were defined for each population highlighting the existence of deeply-derived CO2 that is transported to the surface through active faults. In fact, main WNW-ESE and NNW-SSE trending faults in the three study islands are usually well correlated with CO2 degassing and constitute the pathways for the migration of CO2 from deep reservoirs to the surface; few exceptions may be associated with the superimposed effect of the lithology. The main CO2 emissions seem to occur in the older volcanic systems, where the development of soils reduces possible air contamination and the gases migrate to the surface through the high permeability zones that correspond to the intersection of fractures/faults. The absence of soils and the permeability of the recent lava flows, which allow the dilution of the deep gases with atmospheric air, may explain the lack of CO2 emission in the younger volcanic systems. In addition, locally, some low permeable lava flows may act as barriers for the deep gas emissions.
For Pico and Faial islands two different methodologies were used in the present study to measure the soil CO2diffuse degassing. Even if no simultaneous measurements were performed, higher soil CO2 flux values were essentially coincident with areas with anomalous soil CO2 concentrations. The generally good spatial correlation observed between the soil CO2 flux and concentration values in different surveyed periods shows that the anomalous CO2 degassing areas remain stable along time, arguing for the use of these techniques as useful tools for volcanic monitoring programmes. Considering the need to establish a permanent gas monitoring network in these islands, several DDS are here suggested as potential sites to install a permanent soil CO2 flux station, taking into account the state of activity of the different volcanic systems. DDS F2 or F3 at Faial Island are two potential areas to install a permanent soil CO2 flux station since they correspond to the flank of the active Caldeira central volcano. At Pico Island both DDS P1 and P2 seem to be adequate to install a permanent equipment, even if detailed soil CO2 flux measurements should be performed to select the highest CO2 emission site.
Considering that (a) about 32% of the historical volcanic eruptions that affected the Azores archipelago occurred inland or offshore the islands of Pico, Faial, and São Jorge (Gaspar et al.,
Statements
Author contributions
FV revisited/processed the data and drafted the manuscript. MM, CF, FV, and CS carried out the soil diffuse degassing field surveys and contributed for the elaboration of the data. JG and TF provided important suggestions during the processing. All the authors have read and approved the final manuscript.
Acknowledgments
FV is supported by a Post-Doc Grant from Fundação para a Ciência e Tecnologia (FCT, Ref. SFRH/BPD/100032/2014). This study was also a part of the MARES Project (Ref. TDC/GEO-FIQ/1088/2014), funded by the FCT. The authors would like also to thank Vittorio Zanon for sampling the Pico fumaroles during June 2016 and Lucia Moreno for performing the chemical analyses.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AiuppaA.BurtonM.CaltabianoT.GiudiceG.GurrieriS.LiuzzoM.et al. (2010). Unusually large magmatic CO2 gas emissions prior to a basaltic paroxysm. Geophys. Res. Lett.37, L17303. 10.1029/2010GL043837
2
AllardP.CarbonnelleJ.DajlevicD.Le BronecJ.MorelP.RobeM. C.et al. (1991). Eruptive and diffuse emissions of CO2 from Mount Etna. Nature351, 387–391. 10.1038/351387a0
3
AsimowP. D.DixonJ. E.LangmuirC. H. (2004). A hydrous melting and fractionation model for mid-ocean ridge basalts: application to the Mid-Atlantic Ridge near the Azores. Geochem. Geophys. Geosys.5:Q01E16. 10.1029/2003GC000568
4
BajracharyaR. M.LalR.KimbleJ. M. (2000). Erosion effects on carbon dioxide concentration and carbon flux from an Ohio alfisol. Soil Sci. Soc. Am. J.64, 694–700. 10.2136/sssaj2000.642694x
5
BaubronJ. C.RigoA.ToutainJ. P. (2002). Soil gas profiles as a tool to characterise active tectonics areas: the Jaut Pass example (Pyrenees, France). Earth Planet. Sci. Lett.196, 69–81. 10.1016/S0012-821X(01)00596-9
6
BaxterP.BaubronJ. C.CoutinhoR. (1999). Health hazards and disaster potential of ground gas emissions at Furnas Volcano, São Miguel, Azores. J. Volcanol. Geotherm. Res.92, 95–106. 10.1016/S0377-0273(99)00070-0
7
BeaubienS. E.CiotoliG.LombardiS. (2003). Carbon dioxide and radon gas hazard in the Alban Hills area (central Italy). J. Volcanol. Geotherm. Res.123, 63–80. 10.1016/S0377-0273(03)00028-3
8
BlongR. J. (1984). Volcanic Hazards. A Sourcebook on the Effects of Eruptions. San Diego, CL: Academic Press.
9
BonattiE. (1990). Not so hot “hot spots” in the Oceanic mantle. Science250, 107–111. 10.1126/science.250.4977.107
10
CaliroS.ViveirosF.ChiodiniG.FerreiraT. (2015). Gas geochemistry of hydrothermal fluids of the S. Miguel and Terceira Islands, Azores. Geochim. Cosmochim. Acta168, 43–57, 10.1016/j.gca.2015.07.009
11
CannatM.BriaisA.DeplusC.EscartínJ.GeorgenJ.LinJ.et al. (1999). Mid-Atlantic Ridge-Azores hotspot interactions: along-axis migration of a hotspot-derived event of enhanced magmatism 10 to 4 Ma ago. Earth Planet. Sci. Lett.173, 257–269. 10.1016/s0012-821x(99)00234-4
12
CarapezzaM. L.GranieriD. (2004). CO2 soil flux at Vulcano (Italy): comparison between active and passive methods. Appl. Geochem.19, 73–88. 10.1016/s0883-2927(03)00111-2
13
CarapezzaM. L.InguaggiatoS.BruscaL.LongoM. (2004). Geochemical precursors of the activity of an open-conduit volcano: the Stromboli 2002–2003 eruptive events. Geophys. Res. Lett.31, L07620. 10.1029/2004gl019614
14
CardelliniC.ChiodiniG.FrondiniF. (2003). Application of stochastic simulation to CO2 flux from soil: mapping and quantification of gas release. J. Geophys. Res.108, 2425–2437. 10.1029/2002JB002165
15
ChiodiniG.CaliroS.CardelliniC.AvinoR.GranieriD.SchmidtA. (2008). Carbon isotopic composition of soil CO2 efflux, a powerful method to discriminate different sources feeding soil CO2 degassing in volcanic-hydrothermal areas. Earth Planet. Sci. Lett.274, 372–379. 10.1016/j.epsl.2008.07.051
16
ChiodiniG.CaliroS.CardelliniC.GranieriD.AvinoR.BaldiniA.et al. (2010a). Long-term variations of the Campi Flegrei, Italy, volcanic system as revealed by the monitoring of hydrothermal activity. J. Geophys. Res.115:B03205. 10.1029/2008JB006258
17
ChiodiniG.CioniR.GuidiM.RacoB.MariniL. (1998). Soil CO2 flux measurements in volcanic and geothermal areas. Appl. Geochem.13, 543–552. 10.1016/S0883-2927(97)00076-0
18
ChiodiniG.FrondiniF.CardelliniC.GranieriD.MariniL.VenturaG. (2001). CO2 degassing and energy release at Solfatara Volcano, Campi Flegrei, Italy. J. Geophys. Res.106, 16213–16221. 10.1029/2001JB000246
19
ChiodiniG.FrondiniF.KerrickD. M.RogieJ.ParelloF.PeruzziL.et al. (1999). Quantification of deep CO2 fluxes frm Central Italy. Examples of carbon balance for regional aquifers and soil diffuse degassing. Chem. Geol.159, 205–222. 10.1016/S0009-2541(99)00030-3
20
ChiodiniG.GranieriD.AvinoR.CaliroS.CostaA.MinopoliC.et al. (2010b). Non-volcanic CO2 Earth degassing: the case of Mefite di Ansanto (Southern Apennines), Italy. Geophys. Res. Lett.37, L11303. 10.1029/2010GL042858
21
ChovelonP. (1982). Evolution Volcanotectonique des Iles de Faial et de Pico. Thesis, University Paris-Sud, Orsay.
22
CostaA. C. G.HildenbrandA.MarquesF. O.SibrantA. L. R.Santos de CamposA. (2015). Catastrophic flank collapses and slumping in Pico Island during the last 130 kyr (Pico-Faial ridge, Azores Triple Junction). J. Volcanol. Geotherm. Res.302, 33–46. 10.1016/j.jvolgeores.2015.06.008
23
CruzJ. V.FreireP.CostaA. (2010). Mineral waters characterization in the Azores archipelago (Portugal). J. Volcanol. Geotherm. Res.190, 353–364. 10.1016/j.jvolgeores.2009.12.001
24
DavidM. (1977). Geostatistical Ore Reserve Estimation. Amsterdam: Elsevier Scientific Publishing Company.
25
DemandeJ.FabriolR.GérardA.IundtF. (1982). Prospection Géothermique des Iles de Faial e Pico (Açores). 82SGN003GTH. Bureau de Recherches Géologiques et Minières, Orléans.
26
de MoorJ. M.AiuppaA.AvardG.WehrmannH.DunbarN.MullerC.et al. (2016). Turmoil at turrialba volcano (Costa Rica): degassing and eruptive processes inferred from high-frequency gas monitoring. J. Geophys. Res.121, 5761–5775. 10.1002/2016JB013150
27
DeutschC. V.JournelA. G. (1998). GSLIB: Geostatistical Software Library and User's Guide. Applied Geostatistics Series. New York, NY; Oxford: Oxford Univerity Press.
28
DiasN. A.MatiasL.LourençoN.MadeiraJ.CarrilhoF.GasparJ. L. (2007). Crustal seismic velocity strcuture near Faial and Pico Islands (Azores), from local earthquake tomography. Tectonopphysics445, 301–317. 10.1016/j.tecto.2007.09.001
29
Di ChiaraA.SperanzaF.PorrecaM.PimentelA.D'Ajello CaraccioloF.PachecoJ. (2014). Constraining chronology and time-space evolution of Holocene volcanic activity on the Capelo Peninsula (Faial Island, Azores): the paleomagnetic contribution. Geol. Soc. Am. Bull. 126, 1164–1180. 10.1130/B30933.1
30
EvansB. M.StaudacherT. H. (2001). In situ measurements of gas discharges across fissures associated with lava flows at Réunion Island. J. Volcanol. Geotherm. Res.106, 255–263. 10.1016/S0377-0273(00)00242-0
31
FariaC. (2002). Estudo de Anomalias Geoquímicas Associadas a Processos de Desgaseificação Difusa na Ilha do Faial: Contribuição Para a Cartografia de Falhas Activas. Master Thesis, University of the Azores.
32
FéraudG.KaneokaI.AllègreC. J. (1980). K-Ar ages and stress pattern in the Azores: geodynamic implications. Earth Planet. Sci. Lett.46, 275–286. 10.1016/0012-821X(80)90013-8
33
FerreiraT. (1994). Contribuição Para o Estudo das Emanações Gasosas Associadas a Processos de Vulcanismo no Arquipélago dos Açores. Master Thesis, University of the Azores.
34
FerreiraT.GasparJ. L.ViveirosF.MarcosM.FariaC.SousaF. (2005). Monitoring of fumarole discharge and CO2 soil degassing in the Azores: contribution to volcanic surveillance and public health risk assessment. Ann. Geophys. 48, 787–796. 10.4401/ag-3234
35
FerreiraT.ViveirosF. (2007). Anomalia Geoquímica em Zona Residencial da Freguesia da Praia do Almoxarife (Ilha do Faial, Açores). Technical Report Ref. 01/CVARG/07
36
FischerT. P.ChiodiniG. (2015). Volcanic, magmatic and hydrothermal gases, in The Encyclopedia of Volcanoes, 2nd Edn., eds SigurdssonH.HoughtonB.McNuttS. R.RymerH.StixJ. (San Diego, CL: Elsevier Inc.), 779–797.
37
GasparJ. L.FerreiraT. (1995). Temperatura e composição química das emanações gasosas do Vulcão dos Capelinhos no decorrer da crise sísmica que se desenvolveu a NW da ilha do Faial em 1992/93. Açoriana8, 167–175.
38
GasparJ. L.QueirozG.FerreiraT.MedeirosA. R.GoulartC.MedeirosJ. (2015). Earthquakes and volcanic eruptions in the Azores region: geodynamic implications from major historical events and instrumental seismicity, in Volcanic Geology of S. Miguel Island (Azores archipelago), eds GasparJ. L.GuestJ. E.DuncanA. M.BarrigaF. J. A. S.ChesterD. K. (London: Geological Society, London, Memoirs), 33–49.
39
GiammancoS.GurrieriS.ValenzaM. (1999). Geochemical investigations applied to active fault detection in a volcanic area: the North-East Rift on Mt. Etna (Sicily, Italy). Geophys. Res. Lett.26, 2005–2008. 10.1029/1999gl900396
40
GiammancoS.GurrieriS.ValenzaM. (2006). Fault-controlled soil CO2 degassing and shallow magma bodies: summit and lower east rift of Kilauea volcano (Hawaii), 1997. Pure Appl. Geophys.163, 853–867. 10.1007/3-7643-7584-1_15
41
GiggenbachW. F. (1975). A simple method for the collection and analysis of volcanic gas samples. Bull.Volcanol.39, 132–145. 10.1007/BF02596953
42
GiggenbachW. F. (1996). Chemical composition of volcanic gases, in Monitoring and Mitigation of Volcano Hazards, eds ScarpaR.TillingR. I. (Berlin; Heidelberg: Springer-Verlag), 221–253.
43
GiggenbachW. F.GoguelR. L. (1989). Collection and Analysis of Geothermal Volcanic Water and Gas Discharges. Chemistry Division, DSIR, Report no. CD 2401.
44
GoovaertsP. (1999). Geostatistics in soil science: state-of-the-art and perspectives. Geoderma89, 1–45. 10.1016/S0016-7061(98)00078-0
45
GranieriD.ChiodiniG.MarzocchiW.AvinoR. (2003). Continuous monitoring of CO2 soil diffuse degassing at Phlegraean Fields (Italy): influence of environmental and volcanic parameters. Earth Planet. Sci. Lett.212, 167–179. 10.1016/S0012-821X(03)00232-2
46
GurrieriS.ValenzaM. (1988). Gas transport in natural porous mediums: a method for measuring CO2 flows from the ground in volcanic and geothermal areas. Rendiconti della Società Italiana di Mineralogia e Petrologia43, 1151–1158.
47
HansellA.OppenheimerC. (2004). Health hazards from volcanic gases: a systematic literature review. Arch. Environ. Health59, 628–639. 10.1080/00039890409602947
48
HernándezP. A.NotsuK.SalazarJ. M.MoriT.NataleG.OkadaH.et al. (2001). Carbon dioxide degassing by advective flow from Usu Volcano, Japan. Science292, 83–86. 10.1126/science.1058450
49
HildenbrandA.MadureiraP.MarquesF. O.CruzI.HenryB.SilvaP. (2008). Multi-stage evolution of a sub-aerial volcanic ridge over the last 1.3 Myr: S. Jorge Island, Azores Triple Junction. Earth Planet. Sci. Lett.273, 289–298. 10.1016/j.epsl.2008.06.041
50
HildenbrandA.MarquesF. O.CostaA. C. G.SibrandtA. L. R.SilvaP. F.HenryB.et al. (2012). Reconstructing the architectural evolution of volcanic islands from combined K/Ar, morphologic, tectonic, and magnetic data: the Faial Island example (Azores). J. Volcanol. Geotherm. Res. 241–242, 39–48. 10.1016/j.jvolgeores.2012.06.019
51
HipólitoA.MadeiraJ.CarmoR.GasparJ. L. (2013). Neotectonics of Graciosa Island (Azores): a contribution to seismic hazard assessment of a volcanic area in a complex geodynamic setting. Ann. Geophys.56:S0677. 10.4401/ag-6222
52
HutchisonW.MatherT. A.PyleD. M.BiggsJ.YirguG. (2015). Structural controls on fluid pathways in an active rift system: a case study of the Aluto volcanic complex. Geosphere11, 542–562. 10.1130/GES01119.1
53
IsaaksE. H.SrivastavaR. M. (1989). An Introduction to Applied Geostatistics.New York, NY; Oxford: Oxford University Press.
54
Jean-BaptisteP.AllardP.CoutinhoR.FerreiraT.FourréE.QueirozG.et al. (2009). Helium isotopes in hydrothermal volcanic fluids of the Azores archipelago. Earth Planet. Sci. Lett.281, 70–80. 10.1016/j.epsl.2009.02.009
55
JungN.-H.HanW. S.WatsonZ. T.GrahamJ. P.KimK.-Y. (2014). Fault-controlled CO2 leakage from natural reservoirs in the Colorado Plateau, East-Central Utah. Earth Planet. Sci. Lett.403, 358–367. 10.1016/j.epsl.2014.07.012
56
KanemasuE. T.PowersW. L.SijJ. W. (1974). Field chamber measurements of CO2 flux from soil surface. Soil Sci.118, 233–237. 10.1097/00010694-197410000-00001
57
LeeH.MuirheadJ. D.FischerT. P.EbingerC. J.KattenhornS. A.SharpZ. D.et al. (2016). Massive and prolonged deep carbon emissions associated with continental rifting. Nat. Geosci. 5, 145–149. 10.1038/ngeo2622
58
LiuzzoM.Di MuroA.GiudiceG.MichonL.FerrazziniV.GurrieriS. (2015). New evidence of CO2 soil degassing anomalies on Piton de la Fournaise volcano and the link with volcano tectonic structures. Geochem. Geophys. Geosyst. 16, 4388–4404. 10.1002/2015gc006032
59
LiuzzoM.GurrieriS.GiudiceG.GiuffridaG. (2013). Ten years of soil CO2 continuous monitoring on Mt. Etna: exploring the relationship between processes of soil degassing and volcanic activity. Geochem. Geophys. Geosyst.14, 2886–2899. !10.1002/ggge.20196
60
LourençoN.LuísJ.MirandaJ. M. (1997). Azores Triple Junction Bathymetry. Map edited by Science Faculty, Lisbon and Algarve Universities, 1 sheet.
61
LuoY.ZhouX. (2006). Soil Respiration and the Environment. San Diego, CA: Academic Press.
62
MachadoF.ParsonsW.RichardsA.MulfordJ. W. (1962). Capelinhos eruption of Fayal Volcano, Azores, 1957-1958. J. Geophys. Res., 67, 3519–3529. 10.1029/JZ067i009p03519
63
MadeiraJ. (1998). Estudos de Neotectónica nas Ilhas do Faial, Pico e S. Jorge: Uma Contribuição Para o Conhecimento Geodinâmico da Junção Tripla dos Açores. Ph.D. thesis, Faculdade de Ciências da Universidade de Lisboa, 481.
64
MadeiraJ.Brum da SilveiraA. (2003). Active tectonics and first paleoseismological results in Faial, Pico and S. Jorge islands (Azores, Portugal). Ann. Geophys.46, 733–761.
65
MadeiraJ.Brum da SilveiraA.HipólitoA.CarmoR. (2015). Active tectonics in the Central and Eastern Azores islands along the Eurasia-Nubia boundary: a review, in Volcanic Geology of S. Miguel Island (Azores archipelago), Vol. 44, eds GasparJ. L.GuestJ. E.DuncanA. M.BarrigaF. J. A. S.ChesterD. K. (London: Geological Society, London, Memoirs), 15–32.
66
MadeiraJ.RibeiroA. (1990). Geodynamic models for the Azores triple junction: a contribution from tectonics. Tectonophysics184, 405–415. 10.1016/0040-1951(90)90452-e
67
MarcosM. (2006). Estudo dos Processos de Desgaseificação Difusa nos Açores, com base na Análise de Variações Espaciais de CO2 e Temporais de 222Rn. Master Thesis, University of the Azores, 135.
68
MarquesF. O.CatalãoJ. C.DeMetsC.CostaA. C. G.HildenbrandA. (2013). GPS and tectonic evidence for a diffuse plate boundary at the Azores Triple Junction. Earth Planet. Sci. Lett.381, 177–187. 10.1016/j.epsl.2013.08.051
69
MatiasL.DiasN. A.MoraisI.ValesD.CarrilhoF.MadeiraJ.et al. (2007). The 9th of July 1998 Faial Island (Azores, North Atlantic) seismic sequence. J. Seismol.11, 275–298. 10.1007/s10950-007-9052-4
70
MendesV. B.MadeiraJ.Brum da SilveiraA.TrotaA.EloseguiP.PagareteJ. (2013). Present-day deformation in São Jorge Island, Azores, from episodic GPS measurements (2001-2011). Adv. Space Res.51, 1581–1592. 10.1016/j.asr.2012.10.019
71
MétrichN.ZanonV.CréonL.HildenbrandA.MoreiraM.MarquesF. O. (2014). Is the 'Azores Hotspot'a Wetspot? Insights from the geochemistry of fluid and melt inclusions in olivine of Pico Basalts. J. Petrol.55, 377–393. 10.1093/petrology/egt071
72
MirandaJ. M.LuísJ. F.LourençoN.FernandesR. M. S. (2015). The structure of the Azores Triple Junction: implications for São Miguel Island, in Volcanic Geology of S. Miguel Island (Azores archipelago), eds GasparJ. L.GuestJ. E.DuncanA. M.BarrigaF. J. A. S.ChesterD. K. (London: Geological Society, London, Memoirs), 5–13.
73
MoreiraM.DoucelanceR.KurzM. D.DupréB.AllègreC. J. (1999). Helium and lead isotope geochemistry of the Azores Archipelago. Earth Planet. Sci. Lett.169, 189–205. 10.1016/S0012-821X(99)00071-0
74
NakadaiT.YokozawaM.IkedaH.KoizumiH. (2002). Diurnal changes of carbon dioxide flux from bare soil in agricultural field in Japan. Appl. Soil Ecol.19, 161–171. 10.1016/S0929-1393(01)00180-9
75
NormanJ. M.GarciaR.VermaS. B. (1992). Soil surface CO2 fluxes and the carbon budget of a grassland. J. Geophys. Res.97, 18845–18853. 10.1029/92JD01348
76
NunesJ. C. (1999). A Actividade Vulcânica na Ilha do Pico do Plistocénico Superior ao Holocénico: Mecanismo Eruptivo e Hazard Vulcânico. Ph.D. thesis, University of the Azores, Ponta Delgada, 356.
77
PachecoJ. (2015). 14C Ages of the last eruptive stage of Faial Island Caldeira Volcano. Radiocarbon57, 1021–1027. 10.2458/azu_rc.57.17957
78
PachecoJ. M. R. (2001). Processos Associados ao Desenvolvimento de Erupções Vulcânicas Hidromagmátics Explosivas na ilha do Faial e Sua Interpretação Numa Perspectiva de Avaliação do Hazard e Minimização do Risco. Ph.D. thesis, Azores University, Ponta Delgada, 300.
79
PadrónE.PadillaG.HernándezP. A.PérezN. M.CalvoD.NolascoD.et al. (2013). Soil gas gechemistry in relation to eruptive fissures on Timanfaya volcano, Lanzarote Island (Canary Islands, Spain). J. Volcanol. Geotherm. Res.250, 91–99. 10.1016/j.jvolgeores.2012.10.013
80
PantaleoM.WalterT. R. (2014). The ring-shaped thermal field of Stefanos crater, Nisyros Island: a conceptual model. Solid Earth5, 183–198. 10.5194/se-5-183-2014
81
ParkinsonK. J. (1981). An improved method for measuring soil respiration in the field. J. Appl. Ecol.18, 221–228. 10.2307/2402491
82
PedoneM.ViveirosF.AiuppaA.GiudiceG.GrassaF.GaglianoA. L.et al. (2015). Total (fumarolic + diffuse soil) CO2 output from Furnas volcano. Earth Planets Space67, 174, 10.1186/s40623-015-0345-5
83
PeltierA.FinizolaA.DouilletG. A.BrothelandeE.GaraebitiE. (2012). Structure of an active volcano associated with a resurgent block inferred from thermal mapping: the Yasur-Yenkahe volcanic complex (Vanuatu). J. Volcanol. Geotherm. Res. 243–244, 59–68. 10.1016/j.jvolgeores.2012.06.022
84
QuartauR.MadeiraJ.MitchellN. C.TemperaF.SilvaP. F.BrandãoF. (2015). The insular shelves of the Faial-Pico Ridge (Azores archipelago): a morpholical record of its evolution. Geochem. Geophys. Geosyst.16, 1401–1420. 10.1002/2015GC005733
85
RaichJ. W.SchlesingerW. H. (1992). The global carbon dioxide flux in soil respiration and its relationship to vegetation and climate. Tellus44B, 81–99. 10.3402/tellusb.v44i2.15428
86
RaichJ. W.TufekciogluA. (2000). Vegetation and soil respiration: correlations and controls. Biogeochemistry48, 71–90. 10.1023/A:1006112000616
87
RinaldiA. P.VandemeulebrouckJ.TodescoM.ViveirosF. (2012). Effects of atmospheric conditions on surface diffuse degassing. J. Geophys. Res.117, B11201. 10.1029/2012JB009490
88
SchillingJ. G. (1975). Azores mantle blob: rare-earth evidence. Earth Planet. Sci. Lett.25, 103–115. 10.1016/0012-821X(75)90186-7
89
SchöpaA.PantaleoM.WalterT. R. (2011). Scale-dependent location of hydrothermal vents: stress field models and infrared field observations on the Fossa Cone, Vulcano Island, Italy. J. Volcanol. Geotherm. Res.203, 133–145. 10.1016/j.jvolgeores.2011.03.008
90
SearleR. (1980). Tectonic pattern of the Azores spreading centre and triple junction. Earth Planet. Sci. Lett.51, 415–434. 10.1016/0012-821X(80)90221-6
91
SerralheiroA.ForjazV. H.AlvesC. A. M.RodriguesB. (1989). Carta Vulcanológica dos Açores – ilha do Faial, scale 1:15000, Sheets 1, 2, 3 & 4. Ponta Delgada: Centro de Vulcanologia do INIC, Serviço Regional de Proteção Civil dos Açores e Universidade dos Açores.
92
SilvaC.FerreiraT.ViveirosF.AllardP. (2015). Soil radon (222Rn) monitoring at Furnas Volcano (São Miguel, Azores): applications and challenges. Eur. Phys. J. Spec. Top. 224, 659–686. 10.1140/epjst/e2015-02398-6
93
SilvaM. (2005). Caracterização da Sismicidae Histórica dos Açores com Base na Reinterpretação de Dados de Macrossísmica: Contribuição Para a Avaliação do Risco Sísmico nas Ilhas do Grupo Central. Master thesis, Universidade dos Açores, 146.
94
SinclairA. J. (1974). Selection of threshold values in geochemical data using probability graphs. J. Geochem. Explor.3, 129–149. 10.1016/0375-6742(74)90030-2
95
TrippaneraD.PorrecaM.RuchJ.PimentelA.AcocellaV.PachecoJ.et al. (2014). Relationships between tectonics and magmatism in a transtensive/transform setting: an example from Faial Island (Azores, Portugal). Geol. Soc. Am. Bull. 126, 164–181. 10.1130/B30758.1
96
ViveirosF. (2003). Contribuição Para o Estudo dos Processos de Desgaseificação Difusa nos Açores no mbito da Monitorização Sismovulcânica e da Avaliação do Risco: Discriminação de Factores que Influenciam a Variação do Fluxo do CO2. Master thesis, University of the Azores, 140.
97
ViveirosF.CardelliniC.FerreiraT.CaliroS.ChiodiniG.SilvaC. (2010). Soil CO2 emissions at Furnas volcano, São Miguel Island, Azores archipelago: volcano monitoring perspectives, geomorphologic studies, and land use planning application. J. Geophys. Res.115:B12208, 10.1029/2010JB007555
98
ViveirosF.FerreiraT.Cabral VieiraJ.SilvaC.GasparJ. L. (2008). Environmental influences on soil CO2 degassing at Furnas and Fogo volcanoes (São Miguel Island, Azores archipelago). J. Volcanol. Geotherm. Res.177, 883–893. 10.1016/j.jvolgeores.2008.07.005
99
ViveirosF.FerreiraT.SilvaC.GasparJ. L. (2009). Meteorological factors controlling soil gases and indoor CO2 concentration: a permanent risk in degassing areas. Sci. Tot. Environ.407, 1362–1372. 10.1016/j.scitotenv.2008.10.009
100
ViveirosF.GasparJ. L.FerreiraT.SilvaC. (2016a). Hazardous indoor CO2 concentrations in volcanic environments. Environ. Pollut.214, 776–786, 10.1016/j.envpol.2016.04.086
101
ViveirosF.GasparJ. L.FerreiraT.SilvaC.MarcosM.HipólitoA. (2015). Mapping of soil CO2 diffuse degassing at São Miguel Island and its public health implications, in Volcanic Geology of S. Miguel Island (Azores archipelago), eds GasparJ. L.GuestJ. E.DuncanA. M.BarrigaF. J. A. S.ChesterD. K. (London: Geological Society, London, Memoirs), 185–195.
102
ViveirosF.MorenoL.Carreiro-SilvaM.CoutoR.SilvaC.RangeP.et al. (2016b). Volcanic gas emissions offshore of São Miguel and Faial islands (Azores archipelago), in Actas das 4as Jornadas de Engenharia Hidrográfica, Instituto Hidrográfico (Lisbon), 320–323.
103
ViveirosF.VandemeulebrouckJ.RinaldiA. P.FerreiraT.SilvaC.CruzJ. V. (2014). Periodic behavior of soil CO2 emissions in diffuse degassing areas of the Azores archipelago: application to seismovolcanic monitoring. J. Geophys. Res.119, 7578–7597. 10.1002/2014JB011118
104
VogtP. R.JungW. Y. (2004). The Terceira rift as hyper-slow, hotspot dominated oblique spreading axis: a comparison with other slow-spreading plate boundaries. Earth Planet. Sci. Lett.218, 77–90. 10.1016/S0012-821X(03)00627-7
105
WeinsteinP.CookA. (2005). Volcanic emissions and health, in Essentials of Medical Geology – Impacts of the Natural Environment on Public Health, ed SelinusO. (Burlington, MA: Elsevier), 203–226.
106
WestonF. S. (1964). List of recorded volcanic eruptions in the Azores with brief reports. Bol. Mus. Lab. Min. Geol. Fac. Ciências Lisboa.10, 3–18.
107
ZanonV.FrezzottiM. L. (2013). Magma storage and ascent conditions beneath Pico and Faial islands (Azores Islands). A study on fluid inclusions. Geochem. Geophys. Geosys.14, 3494–3514. 10.1002/ggge.20221
108
ZanonV.KueppersU.PachecoJ.CruzI. (2013). Volcanism from fissure zones and the Caldeira central volcano at Faial Island, Azores archipelago: geochemical processes in multiple feeding systems. Geol. Mag.150, 536–555. 10.1017/s0016756812000702
Summary
Keywords
volcanic gases, diffuse degassing structures, CO2 degassing, tectonics, Azores archipelago, Faial Island, Pico Island, São Jorge Island
Citation
Viveiros F, Marcos M, Faria C, Gaspar JL, Ferreira T and Silva C (2017) Soil CO2 Degassing Path along Volcano-Tectonic Structures in the Pico-Faial-São Jorge Islands (Azores Archipelago, Portugal). Front. Earth Sci. 5:50. doi: 10.3389/feart.2017.00050
Received
13 November 2016
Accepted
01 June 2017
Published
14 June 2017
Volume
5 - 2017
Edited by
Benjamin Alexander Black, City College of New York (CUNY), United States
Reviewed by
José Madeira, Universidade de Lisboa, Portugal; J. Maarten De Moor, OVSICORI-UNA, Costa Rica
Updates

Check for updates
Copyright
© 2017 Viveiros, Marcos, Faria, Gaspar, Ferreira and Silva.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fátima Viveiros maria.fb.viveiros@azores.gov.pt
This article was submitted to Geochemistry, a section of the journal Frontiers in Earth Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.