1. Introduction
Many volcanic systems are characterized by the presence of fissure systems or rifts, radiating from the summit of a central volcano, following local weakness patterns that are linked to the gravitational stress (e.g. Borgia, Ferrari & Pasquaré, Reference Borgia, Ferrari and Pasquarè1992; Borgia, Reference Borgia1994; Walter, Klügel & Münn, Reference Walter, Klügel and Münn2006; Münn, Walter & Klügel, Reference Münn, Walter and Klügel2006). Magma eruptions from fissure zones or from summit vents in active volcanoes, such as Etna (Italy) and Kilauea (Hawaii, USA), respond to the pressure dynamics of multiple reservoirs (Ryan, Koyanagi & Fiske, Reference Ryan, Koyanagi and Fiske1981; Acocella & Neri, Reference Acocella and Neri2003; Corsaro & Pompilio, Reference Corsaro and Pompilio2004). Magmas feeding these eruptions are genetically linked to each other by evolution processes occurring in these reservoirs (Garcia et al. Reference Garcia, Pietruszka, Rhodes and Swanson2000; Corsaro et al. Reference Corsaro, Métrich, Allard, Andronico, Miraglia and Fourmentraux2009).
In other cases, such as the islands of the Sicily Channel (Italy), some Icelandic volcanoes and the East African Rift, magmas erupted from fissure zones and associated central volcanoes show variable degrees of evolution and are believed to have independent evolution histories and origins (Civile et al. Reference Civile, Lodolo, Tortorici, Lanzafame and Brancolini2008; Peccerillo et al. Reference Peccerillo, Barberio, Yirgu, Ayalew, Barbieri and Wu2003, Reference Peccerillo, Donati, Santo, Orlando, Yirgu and Ayalew2007; Lacasse et al. Reference Lacasse, Sigurdsson, Carey, Jóhannesson, Thomas and Rogers2007). The concurrent emission from different vents of magmas with contrasting compositions, during the same eruption, is interpreted as proof of the simultaneous tapping of distinct magma reservoirs (e.g. Água de Pau Volcano, ad1563 eruption, Azores; Storey, Reference Storey1981).
The Island of Faial is located in an area characterized by extensional tectonics (Searle, Reference Searle1980; Vogt & Jung, Reference Vogt and Jung2004). Here, there are two young fissure zones located E and W of a central volcano, respectively (Métrich, Bizouard & Varet, Reference Métrich, Bizouard and Varet1981). While the composition of the magmas emitted from the central volcano ranges from basalt to trachyte, the lavas erupted from the fissure zones are basalts to hawaiites. These magmatic suites are geographically well distinguishable from each other and, furthermore, there is no evidence of interaction between them. However, during the 1957–58 eruption, a fracture system, caused by dyke intrusion, propagated from the fissure zone towards the central volcano and ultimately triggered a phreatic explosion inside its caldera (Machado et al. Reference Machado, Parsons, Richards and Mulford1962).
In this paper, we studied the mineral chemistry as well as major and trace element compositions of volcanic rocks from the two fissure zones and the central volcano. The main objectives are (i) to define the geochemical characteristics of mafic magmas and the composition of mantle sources beneath Faial as compared to other Azorean volcanoes, (ii) to define the role of differentiation processes during magma ascent and (iii) to provide a model of the plumbing system and the dynamics of magma ascent at Faial.
2. Geodynamic setting and geological history
The nine islands of the Azores archipelago are the emergent portions of large volcanic edifices, which rise from the Azores plateau in the North Atlantic Ocean (Fig. 1b). Magmatism in this area is commonly explained either by the presence of a mantle plume head (e.g. Beier et al. Reference Beier, Haase, Abouchami, Krienitz and Hauff2008 and references therein) or by the melting of a volatile-enriched mantle (Schilling, Bergeron & Evans, Reference Schilling, Bergeron and Evans1980; Bonatti, Reference Bonatti1990). The archipelago is also subject to intense but low-magnitude seismicity, linked to the interaction among three major tectonic structures: the East Azores Fracture Zone; the transtensional fault system, known as the Terceira Rift; and the Mid-Atlantic Ridge (Searle, Reference Searle1980; Vogt & Jung, Reference Vogt and Jung2004).
Faial is the westernmost island of the central part of the archipelago (Fig. 1c), which rises from a 10–15 Ma crust formed by Mid-Atlantic Ridge eruptions (Müller et al. Reference Müller, Sdrolias, Gaina and Roest2008). There are two main tectonic lineaments, with directions WNW–ESE and NNW–SSE (Fig. 1a). The first is the expression of extensional tectonics with normal dextral faults, whereas the latter is a conjugate system constituting a series of transtensional faults on the northern flank of the Caldeira central volcano (Madeira & Silveira, Reference Madeira and Silveira2003). Recent stratigraphic studies and new K–Ar dates revealed the sequence of the geological events that characterized the history of the island (Hildenbrand et al. Reference Hildenbrand, Marques, Costa, Sibrant, Silva, Henry, Miranda and Madureira2012).
Eruptions took place both from the central volcano and from the fissure zones during the last 1 Ma, in overlapping periods. The oldest rocks, dated back to about 850 ka, crop out in the Pedro Miguel graben (ESE sector) and refer to an older stratovolcano, mostly buried under younger deposits originated from eruptions of the Caldeira Volcano. Nearby, the products of the Ribeirinha Volcano (active from 400 to 350 ka ago) are represented by strongly weathered lava flows dissected by faults. About 120 ka ago, the Caldeira Volcano, in the central part of the island, started its activity that continued until recent times. Two different lithostratigraphic units form this volcano (Pacheco, Reference Pacheco2001 and references therein). The Grupo Inferior is formed by lava flows cropping out in the northern and southern coastal cliffs of the island. The Grupo Superior is mostly constituted by pumice fallout deposits produced between about 16 and 1.2 ka ago. During these eruptions a small summit caldera was formed. Eruptions from the Horta Fissure Zone took place in the SE part of the island presumably during the last 11 ka, leading to the formation of a WNW–ESE-oriented basaltic ridge and a lava platform. Eruptions from the Capelo Fissure Zone occurred west of the central volcano in recent times (< 6 ka?; Dias et al. Reference Dias, Matias, Lourenço, Madeira, Carrilho and Gaspar2007), forming a chain of cinder cones and associated lava flows aligned along the WNW–ESE tectonic trend. Two historic eruptions in this area occurred in 1672–73 and 1957–58.
3. Sampling and analytical methods
Eighty-seven samples of lavas and pyroclastic rocks, collected from the four volcanic systems of Faial (Fig. 1a), were chosen for petrographic investigation, and to perform bulk rock and mineral analyses. The largest number of samples comes both from the Capelo (31) and Horta (14) fissure zones and from the Caldeira Volcano (38). Only three samples of the Ribeirinha Volcano were used in this study, owing to limited outcrops and their high degree of weathering.
Electron microprobe analyses were performed making use of a JEOL JXA 8200 Superprobe, equipped with five wavelength-dispersive spectrometers, an energy-dispersive X-ray spectrometer and cathodoluminescence detector, located in the Earth Sciences Department of the ‘Ardito Desio’ University of Milan (Italy). Generally, a spot size of 1 μm with a beam current of 15 nA was utilized, with the exception of biotites and alkali feldspars (5 μm, with a beam current of 5 nA). Glasses were analysed with a current of 2–4 nA and a ~ 10 μm wide defocused beam to minimize sodium loss. Counting times for major elements were 30 s on the peak and 10 s on each background. Natural and synthetic minerals and glasses used as standards were calibrated within 2% at 2σ standard deviation. Raw data were corrected by applying a Phi-Rho-Z quantitative analysis program. The typical detection limit for each element is 0.01%. The discrimination between phenocrysts and microphenocrysts was based on the relative crystal size (microphenocryst size is ~150–300 μm, while phenocrysts size is ≥ 300 μm).
Whole-rock compositions were obtained at the Activation Laboratories Ltd, Ontario (Canada), after pulverization with an agate ball mill. Rock powders were mixed with a flux of lithium metaborate and lithium tetraborate and fused inside an induction furnace. The melt was then poured into a solution of 5% nitric acid containing cadmium as an internal standard and continuously stirred until complete dissolution was achieved (~ 30 minutes). The samples were contemporaneously analysed by means of a Perkin Elmer 9000 inductively coupled plasma mass spectrometer (ICP-MS) and an Agilent 735 inductively coupled plasma atomic emission spectrometer (ICP-AES). FeO was determined by titration with permanganate solution. Three blanks and five controls were analysed per group of samples. Duplicates were fused and analysed every 15 samples. The instruments were recalibrated after 40 samples. Analytical precision (2σ) for trace elements was generally better than 1% for most major elements and 8% for most minor and trace elements. Fourteen international rock standards (online Supplementary Data Table S1 available at http://journals.cambridge.org/geo) were used to calibrate the two methods. Whole-rock and mineral compositions are reported in the online Supplementary Data Tables S1 and S2, respectively, at http://journals.cambridge.org/geo.
4. Classification, petrography and mineral chemistry
The composition of sampled rocks ranges from basalt to trachyte (Fig. 2a). Basalt and hawaiite are the most abundant rock types, whereas evolved rocks (mugearite to trachyte) are scarce and mostly owing to eruptions of the Caldeira Volcano in recent times. Trachytic-benmoreitic fall deposits constitute thick covers mantling most of the island topography and severely limit sampling of older rocks.
The investigated samples have a mild Na-alkaline affinity. The samples from the Capelo Fissure Zone are mafic rocks with variable and sometimes very high MgO contents. CIPW norms, calculated using FeO/Fe2O3 ratios according to Irvine & Baragar (Reference Irvine and Baragar1971), indicate that these rocks are all undersaturated in silica with moderate contents in normative nepheline. Horta rocks are also mafic, nepheline- or hypersthene-normative and show less variable MgO contents than the rocks from Capelo. The rocks of the Caldeira Volcano are more evolved (MgO ≤ 4.4 wt%), ne- to hy-normative and become quartz-normative in trachytes (qz ≤ 3.6%). The few Ribeirinha samples are hawaiites and are hy-normative (hy = 1.6–3.2%) (Fig. 2b). A brief petrographic description of the most representative samples for all volcanic systems is reported in Table 1.
Mineral abundances are indicated in subscripts and in percentage. Acronyms used: ph – phenocryst, micro – microphenocryst, gms – groundmass, cpx – clinopyroxene, ol – olivine, pl – plagioclase, af – alkali feldspar, amp – amphibole, ae – aegirine, ap – apatite, ilm – ilmenite, Ti-mag – titanomagnetite, ox – oxides (unspecified), bt – biotite, gl – glass.
4.a. Ribeirinha Volcano
All rocks sampled from this volcano are lavas showing a porphyritic texture both at the meso- and micro-scale. Total phenocryst content is about 30–40 vol.%, constituted by plagioclase, olivine and rare clinopyroxene, with Fe-Ti oxides as accessory phases. The groundmass is intersertal to intergranular.
Plagioclase phenocrysts are dominant, typically tabular, euhedral or subhedral; less commonly they are euhedral and acicular. The largest crystals are antecrysts (Davidson et al. Reference Davidson, Morgan, Charlier, Harlou and Hora2007) and show disequilibrium features (rounded shape, embayments, anhedral melt inclusions at rims and strong rim zoning). Microphenocrysts in the matrix are euhedral and their composition shows a limited variability (An58–63, Ab36–41, Or1–2). In addition, the groundmass plagioclase is compositionally homogeneous (An48–50, Ab48–50, Or2) (Fig. 3a). Many olivine phenocrysts and microphenocrysts (9–16 vol.%) show hopper or more rarely swallow-tail and hourglass morphologies (Donaldson, Reference Donaldson1976). Euhedral crystals show 10 to 30 μm thick iddingsitic films at their rims and along cracks. Their composition ranges from Fo72 to Fo80, while the groundmass composition is constant (Fo72). Augite clinopyroxenes (Wo43–45, En41–43, Fs13–16) are found only in the groundmass. Highly resorbed phenocrysts and anhedral microphenocrysts are rarely observed. Fe-Ti oxides (from Usp43 to Usp49) are ubiquitous in all samples as euhedral microphenocrysts in the matrix and as inclusions in olivines. Subhedral and/or anhedral ilmenite is also present in some samples. Apatite is in the groundmass.
4.b. Caldeira Volcano
Sampled lavas range from aphyric to porphyritic (total phenocryst content 0–40%) hypocrystalline, and contain feldspar, olivine, clinopyroxene and amphibole in decreasing amounts, set in an intergranular, pilotaxitic to doleritic groundmass. Most of the mafic rocks show variably porphyritic textures with plagioclase (up to 39%) predominating over other phases. Some other mafic rocks are almost aphyric with moderate amounts of microphenocrysts set in a micro- to cryptocrystalline groundmass. Pumices are hypocrystalline, aphyric to scarcely porphyritic (total crystal content ≤ 18%), with phenocrysts of plagioclase predominating over anorthoclase and sanidine, clinopyroxene (both augite and aegirine), amphibole and biotite. Pumices are typically banded, both at a meso- and micro-scale.
Plagioclase phenocrysts are euhedral, tabular, twinned and show oscillatory to normal zoning and embayments. Phenocryst and microphenocryst compositions from various samples with different degrees of evolution show a continuous trend from bytownite to andesine (An83–30, Ab17–66, Or0.4–4) in Figure 3a. In some pumices, plagioclase phenocrysts show resorption, indicating disequilibrium with the surrounding melt. Alkali feldspar (anorthoclase and Na-sanidine) is also present as euhedral, tabular and zoned phenocrysts or microphenocrysts in pumices and in some of the lavas (modal abundance < 10 vol.%). Many crystals in the pumices are antecrysts, showing evident resorption. Pyroxene appears as small anhedral and partly resorbed crystals with embayments and reaction rims. Its core compositions range from diopside and diopside-augite in lavas to augite in pumices (Wo42–48, En36–43, Fs12–19) (Fig. 3b). Aegirine-augite crystals are found in a few pumice samples. Euhedral magnetites (Usp48–83) and subhedral ilmenites (Fig. 3c) are commonly found as microphenocrysts. Fluorapatite (Cl < 0.3 wt%, F > 2 wt%), kaersutite amphibole and biotite are present in many evolved samples as both phenocrysts and microlites in the matrix. Quartz is present in many trachytic pumices.
4.c. Horta and Capelo fissure zones
Most of the rocks erupted from the Horta and Capelo fissure zones range from sub-aphyric (Fig. 4a) to porphyritic (Fig. 4b). A group of mafic rocks from the Capelo Fissure Zone shows porphyritic textures with up to 30–50% of large (up to 1 cm) phenocrysts of olivine, clinopyroxene and rare plagioclase, set in a microcrystalline to cryptocrystalline groundmass. Some large crystals show resorbed rims and embayments, indicating disequilibrium with the surrounding melt, and providing evidence for an antecrystic origin (Fig. 4c). These samples will be referred to as ankaramites. Other mafic samples are moderately porphyritic with euhedral clinopyroxene, plagioclase and olivine microphenocrysts set in a micro- to cryptocrystalline groundmass made of the same phases plus Fe-Ti oxides and some glass. Glomeroporphyritic aggregates of clinopyroxene + plagioclase ± olivine ± titanomagnetite are sometimes present.
Bytownite plagioclase is common in all samples as euhedral phenocrysts showing normal or oscillatory zoning and albite twinning. Labradorite is found in the groundmass as small acicular crystals. Olivine is common and occurs as rounded antecrysts in the ankaramites and as phenocrysts and groundmass phase in all mafic samples. Phenocrysts are generally euhedral/subhedral with polyhedral or granular morphology and hopper and swallow-tail morphologies in only minor amounts. Phenocryst composition (Fo72–86) is more Fo-rich than the groundmass (Fo63–76). Clinopyroxene occurs as large resorbed crystals with small embayments and normal to oscillatory zoning, mainly in the strongly porphyritic ankaramites, and as phenocrysts in most of the analysed samples. Composition is diopside to augite (Wo43–49, En37–45, Fs10–17), sometimes with aegirine-augite rims. Clinopyroxene megacrysts in the most MgO-rich ankaramite of the Capelo Fissure Zone have lower contents of Ti and Ca and higher contents of Mg and AlVI than other clinopyroxenes. The variations of the AlIV/AlVI ratio and Ti in clinopyroxenes of the central volcano and the Capelo Fissure Zone are negatively correlated with Si and form distinct trends (Fig. 5). According to Wass (Reference Wass1979), these negative trends are related to the variation of the crystallization pressure. In contrast, data from the Horta Fissure Zone are scattered and mostly show high and poorly variable AlIV/AlVI and a more pronounced variation in Si, as compared to the other volcanic systems. Fe-Ti oxides (1–9 vol.%) are common in all samples as euhedral crystals (Usp41–83) and as grains in the groundmass. Ilmenite forms subhedral acicular crystals and is also present in the groundmass together with titanomagnetite (Fig. 3c). Euhedral and small-size Cr-spinels occur embedded in many olivine antecrysts and phenocrysts.
5. Geochemistry
Variations in major elements versus MgO indicate an overall continuous variation from mafic to felsic rocks (Fig. 6), with a significant scattering of some elements for the intermediate samples (MgO = 2–4 wt%). A clear kink in the variation trend of some major oxides is observed at MgO values of 4 wt%.
The rocks from the Capelo and Horta fissure zones show large variations in MgO (between 5.2 and 12.5 wt%), with highly porphyritic ankaramites displaying MgO > 10 wt%. There is a small increase in CaO and a decrease in most other major elements with increasing MgO, all defining smooth trends. FeOtotal is almost constant. For MgO contents of 6–8 wt%, Capelo basalts show slightly higher values of Na2O and K2O (respective average 3.6 and 1.6 wt%) than those of Horta (respective average 3.4 and 1.5 wt%). Caldeira rocks show lower MgO contents as well as different variation trends to Horta and Capelo. Silica and TiO2 show a small compositional gap between mugearites and benmoreites, which might depend on the lack of exposure of intermediate rocks in the field, owing to the large post-16 ka pyroclastic cover.
The few Ribeirinha Volcano lava compositions show limited variations in MgO content (5.7–6.3 wt%). They show slightly higher P2O5, Na2O, SiO2 and lower CaO contents than samples of the younger fissure zones at the same MgO content.
Compatible elements such as Ni, Sc, Cr and Co show positive correlation with MgO, whereas large ion lithophile elements (LILE: Rb, Ba, K, Cs) and high-field-strength elements (HFSE: Ta, Nb, Zr, Hf, Th, U) have opposite trends (Fig. 7). Barium and Sr increase with decreasing MgO content in mafic and intermediate samples, and decrease in the trachytes of the Caldeira Volcano. Two highly porphyritic benmoreites have an anomalously high Ba content and fall outside the trend defined by the bulk of the samples. A group of trachytes appears strongly enriched in incompatible elements (Rb, U, Th and rare earth elements) and depleted in Sr, as compared to other trachytes that have similar major element compositions. At the same MgO content (8–6 wt%), Capelo basalts show slightly higher Sr, Rb and Ba and lower heavy rare earth element (HREE) contents as compared to those from Horta.
Mantle-normalized multi-element patterns (Fig. 8) from basalts to mugearites from all the volcanic systems are very similar to one another and show a smooth upward convex shape with small anomalies in K, Hf and Pb (when detected), which are more remarkable in samples from Caldeira. Patterns for the evolved rocks show negative spikes in Ba, Sr and Ti. REEs are moderately fractionated for both light and heavy elements in all the mafic rocks (not shown). Benmoreites to trachytes display poorly fractionated HREE patterns and sometimes negative Eu anomalies.
6. Thermobarometry calculations
Thermobarometric data, reported in Table 2, were calculated to obtain information on the physical conditions of magma storage and crystallization. Regarding the temperature calculation, we used, whenever possible, the composition of groundmass glass and of the coexisting mineral phases, either plagioclase or clinopyroxene (Putirka, Reference Putirka, Putirka and Tepley2008). However, despite the large mineral database, these conditions were rarely met. For this reason, we used alternative algorithms that consider either the equilibrium among clinopyroxene–olivine–spinel–ilmenite (Andersen, Lindsley & Davidson, Reference Andersen, Lindsley and Davidson1993), or solely focus on the chemical composition of the clinopyroxenes (Putirka, Reference Putirka, Putirka and Tepley2008). The presence of amphibole in the evolved magmas allowed the use of a recently developed thermometer based upon the chemical composition of amphibole (Ridolfi & Renzulli, Reference Ridolfi and Renzulli2012). Although it was calibrated for calcic amphibole, it was also proven quite reliable for alkaline amphibole.
Data are calculated either based upon the chemical composition of clinopyroxene or upon the clinopyroxene-liquid equilibrium (*), or upon plagioclase-liquid equilibrium (†). In a single case, the condition of equilibration of various coexisting phases was calculated by the use of QUILF software, while temperatures for trachytes were calculated using the thermometer based upon the chemical composition of amphibole. References are: (a) Andersen & Lindsley (Reference Andersen and Lindsley1985); (b) Andersen, Lindsley & Davidson (Reference Andersen, Lindsley and Davidson1993); (c) Putirka (Reference Putirka, Putirka and Tepley2008); (d) Ridolfi & Renzulli (Reference Ridolfi and Renzulli2012).
However, since the algorithm involving only the composition of the clinopyroxene is quite imprecise, and owing to the diffuse presence of antecrysts, we extrapolated the temperature path by comparing all available temperature estimations with the MgO content in the bulk rocks. With the exception of a few outliers (antecrystic clinopyroxenes), a good correlation was found (R2 = 0.98), suggesting the existence of a single liquid line of descent, common to all the volcanic series. For this reason, we considered these temperatures fully representative of the conditions of crystallization of the Faial magmas. The use of these temperatures in the clinopyroxene-liquid barometer provided reliable calculations of pressure values (Putirka, Reference Putirka, Putirka and Tepley2008).
7. Discussion
7.a. Nature of the magma source and characteristics of primitive melts
The geochemistry of the source regions beneath the Azores archipelago has been widely discussed during the last two decades to reveal the existence of considerable differences in the Sr–Nd–Pb isotopic signature. In particular, four isotopically distinct end-members have been recognized (Millet et al. Reference Millet, Doucelance, Baker and Schiano2009; Beier, Haase & Turner, Reference Beier, Haase and Turner2012). These are represented by the high 206Pb/204Pb, high 87Sr/86Sr magmas from São Miguel Island, the high 206Pb/204Pb, low 87Sr/86Sr magmas from Graciosa, Terceira and São Jorge islands, and by the low 206Pb/204Pb, low 87Sr/86Sr magmas from Faial Island and the western part of Pico Island. A fourth component, present in most of the islands, can be found at the intersection of the compositional field of São Miguel with that of Terceira, and is thought to consist of upper mantle material from the nearby Mid-Atlantic Ridge (Turner et al. Reference Turner, Hawkesworth, Rogers and King1997; Moreira et al. Reference Moreira, Doucelance, Dupré and Allègre1999; Millet et al. Reference Millet, Doucelance, Baker and Schiano2009).
The most common explanation for this large variety of end-members refers to a mixing process at the regional scale between enriched and depleted mantle types (França et al. Reference França, Tassinari, Cruz, Aparicio, Araña and Rodrigues2006; Millet et al. Reference Millet, Doucelance, Baker and Schiano2009). Enriched mantle types (EM-I, EM-II and HIMU end-members) show an important ‘crustal signature’, which is thought to be the result of either melting of a deep-seated eclogitic crust, related to subduction processes older than 2.5 billion years (Turner et al. Reference Turner, Tonarini, Bindeman, Leeman and Schaefer2007), or of contamination/assimilation with recycled basaltic material (Widom & Shirley, Reference Widom and Shirey1996; Beier, Stracke & Haase, Reference Beier, Stracke and Haase2007; Elliott et al. Reference Elliott, Blichert-Toft, Heumann, Koetsier and Forjaz2007). The depleted mid-ocean ridge basalt (MORB) mantle (or DMM) is the source for MORB, and uncontaminated by enriched components.
While the isotopic compositions of the easternmost islands, São Miguel and Terceira islands, indicate interaction between EM-II and HIMU end-members (Beier, Haase & Hansteen, Reference Beier, Haase and Hansteen2006; Elliott et al. Reference Elliott, Blichert-Toft, Heumann, Koetsier and Forjaz2007), lavas erupted on the islands of the central part of the archipelago, São Jorge and especially Pico, which is only a few kilometres from Faial, indicate possible mixing among DMM and EM-I end-members. The Sr–Nd–Pb fields of the basalts of Faial Island overlap with the least 206Pb/204Pb radiogenic compositions of Pico Island; however, Faial basalts show more radiogenic Sr isotopic compositions. The EM-I type signature of some Faial–Pico magmas is associated with very low 187Os/188Os ratios that are typical of unradiogenic reservoirs and suggest the presence of material from an ancient subducted oceanic plate (Schaefer et al. Reference Schaefer, Turner, Parkinson, Rogers and Hawkesworth2002; Turner et al. Reference Turner, Tonarini, Bindeman, Leeman and Schaefer2007). On the other hand, the depleted signatures of other samples from Faial match with the short distance of this island from the Mid-Atlantic Ridge (~ 100 km), suggesting that these magmas suffered extensive mixing with MORB melts (Millet et al. Reference Millet, Doucelance, Baker and Schiano2009).
To unravel the melting processes and the composition of the source rocks for the Faial magmas it is necessary, however, to investigate the geochemical characteristics of the primitive (i.e. poorly evolved) mafic rocks, which are commonly found at both the Capelo and Horta fissure zones. These rocks are aphyric-to-sub-aphyric at the hand specimen scale and reveal a poorly porphyritic texture under the microscope, with small amounts of euhedral microphenocrysts (~ 5 to 15 vol.%). The compositional range of these rocks spans from saturated (hy-normative) to slightly undersaturated (ne-normative) in silica (Fig. 2b) and, in addition, they have lower contents of MgO (~ 9–7 wt%), Ni (~ 80–150 ppm) and Cr (~ 100–300 ppm), as well as Mg no. (~ 55–65), than mantle-equilibrated melts. They can be assumed to be derivative liquids of primary melts that underwent moderate degrees of fractional crystallization of mafic phases, especially olivine with accessory spinel. These processes significantly modified compatible element abundances such as Ni and Cr, but caused a minor enrichment in incompatible elements, while leaving incompatible element ratios unaffected. Therefore, incompatible trace element ratios and, to a lesser extent, their abundances can be used in these samples to constrain the composition of primary magmas and mantle melting processes.
The variable degrees of silica saturation and LILE abundances of the Faial basalts can be explained by the assumption that they are related to each other by high-pressure fractional crystallization processes (O'Hara, Reference O'hara1968; Green & Ringwood, Reference Green and Ringwood1976; Gupta, Green & Taylor, Reference Gupta, Green and Taylor1987). However, similar abundances of MgO, Ni and Cr (Fig. 7) argue against this hypothesis, and the variable pristine composition of the primary melts is considered to be a better explanation.
These basalts show low variability of some elements, especially HFSEs (Fig. 7) and some element ratios (Ba/Nb, Rb/Ta, Th/Ta and LILE/LILE; Fig. 9). Only Rb, Ba and K, to a lesser extent, show greater variations; Ba/Nb and Rb/Ta ratios increase with normative nepheline (Fig. 10). REEs are fractionated for both light REEs (LREEs) and HREEs, but there is variable La/Yb at almost constant TbN/YbN (~ 1.7–2.0) (Fig. 9f).
Since garnet is one of the few discriminating mantle phases as regards HREEs (i.e. DGt Tb = 1.0–1.18; DGt Yb = 5.0–6.25) (Adam & Green, Reference Adam and Green2006), the observed HREE fractionation of the Faial basalts can be explained only by the assumption that some garnet was left in the residue during mantle melting. This explanation supports the hypothesis that the magma's parental rock was a garnet peridotite (i.e. P > 2.8 GPa; Robinson & Wood, Reference Robinson and Wood1998). As confirmation, Beier, Haase & Turner (Reference Beier, Haase and Turner2012) calculated that melting beneath Faial occurred at P > 3 GPa.
According to geochemical modelling by Bourdon, Turner & Ribe (Reference Bourdon, Turner and Ribe2005), the degree of partial melting in the mantle source for the Azores archipelago can be estimated from the variations of La/Yb versus Tb/Yb. Based on the geochemistry of the most primitive Faial basalts (MgO = 7–10%), the degree of partial melting for this island would be around 3–5% for the Horta Fissure Zone and 2–3% for the Capelo Fissure Zone (Fig. 11), with a slightly higher amount of garnet left during the formation of the Capelo magmas. It should be noted, however, that such a model assumes the existence of a common homogeneous source for all the magmas, while this hypothesis may not be always true, both at a regional scale and even at the level of a single island.
Variable degrees of partial melting of a homogeneous source should also produce a parallel increase in HFSEs (e.g. Ta, Nb), which are as greatly incompatible with mantle mineralogy as many LILEs and LREEs. This is not observed at Faial, as discussed earlier. Therefore, geochemical evidence argues against variable degrees of partial melting and rather supports the hypothesis that melting originated from a variably LILE-enriched heterogeneous source in the mantle. In effect, this could be connected to various degrees of metasomatic modifications, which heterogeneously enriched the mantle source in LILEs such as Rb, Ba, Na and K, but not in HFSEs. LILEs could have been stored in various hydrous phases (e.g. phlogopite, amphibole; McKenzie & O'Nions, Reference McKenzie and O'Nions1991; Zack, Foley & Jenner, Reference Zack, Foley and Jenner1997), whose breakdown enriched the source in these elements (Melzer & Wunder, Reference Melzer and Wunder2001).
Mafic rocks from Faial show similar compositions for some elements, but also display important geochemical differences with respect to other Azorean basalts, at a similar range of MgO (7–10 wt%) and SiO2 (White, Tapia & Schilling, Reference White, Tapia and Schilling1979; Turner et al. Reference Turner, Hawkesworth, Rogers and King1997; Widom et al. Reference Widom, Carlson, Gill and Schmincke1997; Claude-Ivanaj, Joron & Allégre, Reference Claude-Ivanaj, Joron and Allègre2001; Beier, Haase & Hansteen, Reference Beier, Haase and Hansteen2006; Beier, Stracke & Hansteen, Reference Beier, Stracke and Haase2007; Elliott et al. Reference Elliott, Blichert-Toft, Heumann, Koetsier and Forjaz2007; Beier et al. Reference Beier, Haase, Abouchami, Krienitz and Hauff2008; Hildenbrand et al. Reference Hildenbrand, Madureira, Ornelas Marques, Cruz, Henry and Silva2008; Prytulak & Elliott, Reference Prytulak and Elliott2009; Millet et al. Reference Millet, Doucelance, Baker and Schiano2009; Madureira et al. Reference Madureira, Mata, Mattielli, Queiroz and Silva2011; Beier, Haase & Turner, Reference Beier, Haase and Turner2012). All mafic rocks from the Azores have comparable abundances of HFSEs and HFSE/HFSE ratios. In contrast, they are slightly undersaturated in silica in most islands with moderate amounts of normative nepheline. Unlike at Faial, hy-normative basalts are scarce or absent. LILE abundances are slightly higher at Faial, with some LILE/HFSE displaying notably higher values than at other islands. Overall, the available trace element data indicate some differences in the enrichment mechanisms in the upper mantle beneath Faial, with more intensive degrees of LILE enrichment than in the other islands.
Compared to the islands of Capo Verde, Madeira, the Canary archipelago and St Helena, the LILE/HFSE ratios of Faial (Rb/Ta, Ba/Nb) are generally higher, whereas they are more similar to those of Tristan da Cunha and slightly lower than those of Gough. In contrast, La/Yb and Tb/Yb ratios are generally low (data from the ‘georoc database’ and references therein: http://georoc.mpch-mainz.gwdg.de/georoc/). However, a detailed discussion on such a complex pattern of element variations is beyond the scope of this paper. It can only be commented that this complexity likely reflects variable degrees of mantle metasomatism and of partial melting in the generation of mafic magmas in these islands. High LILE/HFSE values at Faial highlight important LILE enrichments that suggest a different style of metasomatic modification from that of other Atlantic islands. Such a preferential enrichment is not typical of intraplate environments, and, if confirmed by further detailed studies on other mafic rocks at a regional scale, it would have important implications regarding the mechanism of mantle enrichment beneath the Azores. It would also constitute a constraint on the solution to the dilemma of whether volcanism represents a deep mantle plume or a volatile-enriched mantle (Bonatti, Reference Bonatti1990; McKenzie & O'Nions, Reference McKenzie and O'Nions1995; Montelli et al. Reference Montelli, Nolet, Dahlen, Masters, Engdahl and Hung2004; Asimow, Dixon & Langmuir, Reference Asimow, Dixon and Langmuir2004).
There are various metasomatic agents invoked to explain the geochemical characteristics of mantle worldwide. Among them, the only one that can cause the variations observed between the two fissure zones (i.e. selective variations in Ba, Rb, Sr, K, LREEs and degree of silica saturation/undersaturation) is an aqueous fluid phase, coexisting or not with silicate melts (e.g. Scambelluri et al. Reference Scambelluri, Herman, Morten and Rampone2006). This phase, however, should play a role only at a local scale, as not all the magmas erupted at the Horta Fissure Zone, as an example, show the same characteristics. This particular geochemical signature should be acquired at the moment of melt formation, at the spinel/garnet transition. For this reason, the fluids released by serpentinized oceanic crust at the depth of the Moho can be ruled out. At the same time, it is improbable that these fluids could be provided by the 3 Ga material lying at depth. Lithospheric mantle beneath the Azores is mainly constituted by ultrarefractory lithologies (Simon et al. Reference Simon, Neumann, Bonadiman, Coltorti, Delpech, Gregoire and Widom2008); however, the local presence of hydrous lithologies cannot be fully excluded, and thus their dehydration (phlogopite or amphibole breakdown) could release the necessary fluid (e.g. Coltorti & Gregoire, Reference Coltorti, Gregoire, Coltorti and Gregoire2008).
7.b. Magma evolution processes: the role of fractional crystallization
The most primitive magmas of Faial experienced variable degrees of evolution during their ascent to the surface. The comprehension of the processes of magma evolution has important implications for understanding the origin of derivative magmas and for modelling the plumbing system of the volcano. In general, only magmas from the Capelo and Horta fissure zones experienced moderate degrees of fractional crystallization and evolved to hawaiites, whereas the evolution process of the lavas erupted from the Caldeira Volcano continued and completed with the formation of trachytes.
Magma evolution processes were tested both by trace element modelling and mass balance calculations, based on major elements and mineral chemistry (Stormer & Nicholls, Reference Stormer and Nicholls1978), aiming at determining (1) the amount and type of fractionation processes, (2) the nature of crystallizing mineral phases in the different volcanic centres and (3) the conditions of magma storage and fractionation.
Starting from a parent magma with the composition of FYS301, one of the most primitive and sub-aphyric rocks from the Horta Fissure Zone, the formation of the ne-normative alkali basalts of the fissure zones must have required about 20% removal of Fe-Mg phases plus some plagioclase. A higher amount of fractional crystallization (37%) of the same phases, plus oxides has been calculated for the hy-normative basalt-to-hawaiite suite (see Table 3 for details). The generation of the hawaiites of the Capelo Fissure Zone required the removal of an amount of mineral phases similar to those of Horta (24.5%), starting from the least evolved basaltic sample (FYS357) of the series. The most primitive rocks of the Caldeira Volcano are both hy- and ne-normative and show low abundances of MgO and Ni and Cr trace elements (Figs 2, 7). Mass balance calculations suggest that fractional crystallization of mafic phases, titanomagnetite and plagioclase accounts for the formation of the hy-normative alkali basalts of the Caldeira Volcano. An even higher degree of fractionation of these phases, along with amphibole, some biotite, apatite and alkali feldspars, favoured the generation of benmoreite and trachyte magmas (see Table 3 for details). Overall, about 80% total fractional crystallization is necessary to form trachytes from basaltic parental magmas. The formation of the undersaturated series requires significantly lower amounts (~ 60%) of fractional crystallization to obtain trachytic compositions. It must be noted, however, that the modelled undersaturated trachytes (FYS14) contain less silica and higher MgO than the modelled oversaturated ones (FYS17) (online Supplementary Data Table S1 at http://journals.cambridge.org/geo).
Samples from the Grupo Inferior, characterized by MgO contents in the range 3–4.4 wt%, show anomalously high Al2O3, CaO and Sr contents. While Al2O3 is negatively correlated with MgO, the contents of the other two elements show a positive correlation. These variations are linked to the accumulation of plagioclase in these rocks. Other samples of the Grupo Inferior, especially the most evolved ones, show high Ba contents and a small but significant positive Eu anomaly in REE patterns. These two parameters are positively correlated with each other (not shown). In particular, two benmoreites show very high Ba content (exceeding 1130 ppm), which is not linked to any other chemical parameter. These rocks contain a few alkali feldspar phenocrysts with disequilibrium features, indicating that the high Ba content is related to the accumulation of alkali feldspar.
The variation of trace elements is more sensitive to evolution processes and can provide further constraints on the model. Figure 12 shows the variation diagrams of Sr v. Ba and La/Yb v. Zr. Barium-rich (cumulated) samples have been excluded from the Ba–Sr plot, in which three different trends may be distinguished (Fig. 12a). A positive trend (trend 1) is defined by the bulk of mafic rocks, when fractionation is dominated by mafic phases, and both Ba and Sr are incompatible. It is notable, however, that some original variations of the primitive magmas also superimpose over this trend. A negative trend is observed for intermediate compositions when plagioclase becomes an important fractionating phase and Sr, but not Ba, becomes compatible (trend 2). Finally, both Ba and Sr decrease towards trachytic composition when both alkali feldspar and plagioclase become dominant fractionating phases (trend 3). This is in agreement with the results of mass balance calculations.
Various trends are also present in the plot where the ratios of elements with a different degree of incompatibility, such as La and Yb, are plotted against an incompatible element such as Zr (Fig. 12b). Poorly evolved compositions show limited variations of Zr, as compared to the values of the La/Yb ratio. As discussed in the previous section, this is an effect of the partial melting of sources with different degrees of metasomatism (i.e. variable enrichments in LREEs and LILEs), followed by little degrees of fractional crystallization. The other two trends regarding the evolved rocks indicate separation of different phenocrysts during fractional crystallization. A positive trend in Zr v. La/Yb, which is mostly formed by intermediate rocks of the Caldeira Volcano, requires the fractional crystallization of mainly clinopyroxene. The other trend, which shows poorly variable La/Yb at variable Zr, can be modelled by about 80% fractional crystallization, dominated by separation of feldspars.
7.c. Modelling of the storage and feeding system
The petrological and geochemical modelling of the magma evolution discussed above may aid in the comprehension of the dynamics governing magma ascent and storage at Faial Island and provide important indications about the structure of the plumbing system, both along the fissure zones and at the central volcano. The main steps regarding magma generation and evolution at Faial consist of (i) low degrees of partial melting of a variably metasomatized upper mantle in the garnet facies (> 85–90 km deep; Robinson & Wood, Reference Robinson and Wood1998); (ii) moderate degrees of fractional crystallization beneath the areas of the Horta and Capelo fissure zones, with the formation of hawaiites; and (iii) high degrees of fractional crystallization beneath the Caldeira Volcano, with the formation of benmoreites and trachytes.
Geobarometric data for basalts and hawaiites from the Horta and Capelo fissure zones (Table 2), suggest high to intermediate pressures of crystallization, corresponding to the uppermost mantle (~ 15 to 28 km). This is shallower than the zone of magma origin of basalts, located at a depth ranging on average from 3.1 to 3.4 GPa (Beier, Haase & Turner, Reference Beier, Haase and Turner2012), in the garnet peridotite facies (~ 100 to 110 km). However, primitive crystals formed at pressures greater than 1 GPa. A polybaric crystallization is also suggested by AlIV/AlVI–Si–Ti variation for clinopyroxenes only for the Capelo Fissure Zone and the central volcano (Fig. 5). These pressures are higher than the stability field of plagioclase, which explains the absence of this phase as a phenocryst in many rocks, including ankaramites, as well as the incompatible behaviour of Sr during fractionation of these magmas. The poorly variable and generally high AlIV/AlVI values in lavas from the Horta Fissure Zone suggest that clinopyroxene crystallization occurred at a shallower depth than for the other sites.
Hawaiites rose from storage areas located at the boundary between the crust and mantle and erupted through fractures that followed the regional trend. The crystallization of plagioclase and magnetite-ilmenite microlites in these rocks probably occurred during water degassing in shallow feeder dykes prior to eruption, producing the characteristic intersertal texture of these lavas. The low degree of evolution of the lavas erupted at the Horta and Capelo fissure zones indicates a lack or minor degree of magma ponding and associated crystallization within the crust.
As regards the Caldeira Volcano, hy-normative basalts, similar in composition to those of the Horta Fissure Zone, equilibrated at 760 MPa (~ 26 km) beneath the present edifice, until the probable generation of a Sr-rich mugearite–benmoreite series at 159 MPa (~ 9 km). A Sr-poor benmoreite–trachyte series formed after the ~ 70% fractionation of mafic phases, biotite, apatite and feldspars in a reservoir possibly located at a depth of 100–132 MPa (~ 3.5–5 km). Possibly, this ponding level indicates that the conditions of neutral buoyancy for these magmas were matched, owing to a sudden jump in density of the country rocks. This jump may be related to the shift from hydrothermalized MORB-type basalts, ~ 15 Ma (Müller et al. Reference Müller, Sdrolias, Gaina and Roest2008), to the younger (~ 3–10 Ma; Cannat et al. Reference Cannat, Briais, Deplus, Escartín, Georgen, Lin, Mercouriev, Meyzen, Muller, Pouliquen, Rabain and Da Silva1999; Escartín et al. Reference Escartín, Cannat, Pouliquen, Rabain and Lin2001; Maia, Goslin & Gente, Reference Maia, Goslin and Gente2007) and denser basaltic rocks of the Azores Plateau. A schematic petrology-based model of the plumbing system of Faial is outlined in Figure 13.
8. Conclusions
The Island of Faial was built by mostly lavas and subordinate pyroclastic rocks ranging in composition from basalt to trachyte. Eruptions of mafic magmas occurred along the fissure zones of Capelo and Horta. A suite of rocks of basaltic to trachytic composition built up the Caldeira Volcano starting from about 120 ka ago. Emission of trachytes started about 16 ka ago, when the two fissure zones also began their activity.
Mafic magmas range from silica-saturated to moderately silica-undersaturated. They show a positive correlation between normative nepheline and many LILE abundances and LILE/HFSE ratios. Geochemical evidence excludes the hypothesis that these correlations are the result of fractional crystallization and rather suggests that mantle-equilibrated parental magmas were generated in a heterogeneous source that had undergone variable degrees of metasomatism. Such a process generated variable enrichments in LILEs, but not in other incompatible elements such as HFSE (Nb, Ta, etc).
Mafic magmas underwent fractional crystallization beneath the fissure zones and the Caldeira Volcano. Moderate degrees of polybaric fractionation in deep magma reservoirs beneath the Capelo Fissure Zone generated the basalt-to-hawaiite magmas, which directly fed the eruptions at this zone. In the case of the Horta Fissure Zone, fractional crystallization occurred at a shallower depth to generate magmas with a similar degree of evolution as the Capelo Fissure Zone. Polybaric fractional crystallization at shallower depths beneath Caldeira led to the formation of basalt-to-benmoreite and trachyte series, with silica increasing over time.
The space and time distribution of volcanic activity in Faial, coupled with the variation of petrological and geochemical features of the erupted products, can be summarized in a conceptual model of the internal structure of the island, as already attempted for several other active volcanic systems (e.g. Vulcano Island, Italy; Zanon, Frezzotti & Peccerillo, Reference Zanon, Frezzotti and Peccerillo2003; Peccerillo et al. Reference Peccerillo, Frezzotti, De Astis and Ventura2006). According to the proposed model, the slightly undersaturated to saturated parental magmas evolved by fractional crystallization in the uppermost mantle, producing a hawaiite suite. These magmas reached the surface without significant stagnation and erupted along the two fissure zones. On the contrary, more complex evolution processes took place beneath the Caldeira Volcano, producing a suite of basalt-to-trachyte melts. Geochemical and petrological data indicate a polybaric fractionation of these magmas, and suggest the possible existence of reservoirs at different depths within the crust. The main emission of evolved magmas during the latest stages of Caldeira activity suggests that the eruptions tapped the shallower magma chambers where evolution had generated trachytic melts. The lack of ponding areas for basaltic-hawaiitic magmas in the crust may depend on the presence of a density filter at the base of the Moho. The presence of ultramafic cumulitic levels in the lithospheric mantle at the base of the crust could act as a density filter, stopping ascending magmas and favouring their accumulation (underplating) and the subsequent fractionation. This phenomenon has already been reported for the volcanoes of the Canary Islands and Madeira (Dañobeitia & Canales, Reference Dañobeitia and Canales2000; Schwarz, Klügel & Wohlgemuth-Ueberwasser, Reference Schwarz, Klügel and Wohlgemuth-Ueberwasser2004; Klügel, Hansteen & Galipp, Reference Klügel, Hansteen and Galipp2005) and therefore would not constitute a new finding in the framework of the oceanic islands of this area of the Atlantic Ocean.
The variable behaviour of magmas during their ascent beneath the central volcano and the rift zones could have been the result of different local conditions of the tectonic stress. Fissure volcanism developed along extensional WNW- and ESE-trending fissures, whereas the central volcano formed at the intersection between those faults and the NNW–SSE transtensional systems. In such a tectonic setting, the formation of magma reservoirs was favoured and trachytic magmas were formed.
Acknowledgements
This work has been funded by the Direcção Regional da Ciência e Tecnologia (DRCT – medida M2.1.2/004/2007 and programme RAA/DRCT). A. Risplendente and S. Poli of the Earth Sciences Department of the ‘Ardito Desio’ University of Milan (Italy) are gratefully acknowledged for their assistance during microprobe analyses. Thanks are also extended to R. Coutinho, who provided help during fieldwork. Discussions with C. Beier and criticisms from two anonymous reviewers proved very useful and contributed to the improvement of the manuscript.