Shear-wave Velocity Structure around Teide Volcano: Results Using Microtremors with the SPAC Method and Implications for Interpretation of Geodetic Results

Deformation analysis and simulation of volcanic edifices require the construction of models of elastic properties of those structures. In this paper we present an analysis of microtremor measurements recorded during the performance tests of two temporary seismic arrays installed in the eastern portion of the Teide caldera in 1994. We take advantage of recent developments of the SPAC method and use spatial cross-correlation computations to estimate phase velocity dispersion of Rayleigh waves at the location of the arrays. We show that the extension of the standard SPAC method is valid in the case of our data, justifying its use and supporting the generalization of the SPAC method to single station pairs. The phase velocity dispersion curve obtained was inverted to recover the shear-wave profile at the site of the arrays. Our results indicate that the subsoil structure of the caldera is laterally homogeneous at the scale of a few km about the location of the arrays. We obtained about 315 m of volcanic sediments overlying rocks with a shear-wave velocity of 2 km/s. These results are robust and are a starting point to further modelling of deformation, permanent or transient, at this volcanic edifice, which can be useful in the interpretation of different observed fields. In fact, the computation of deformations and gravity changes due to possible volcanic intrusions in two models; one considering the volcanic sediments and the other without considering them, provided different results in the near field.     

Magma storage conditions of the last eruption of Teide volcano (Canary Islands, Spain)

Phase equilibrium experiments were performed to determine the pre-eruptive conditions of the phonolitic magma responsible for the last eruption (about 1,150 yr B.P.) of Teide volcano. The Lavas Negras phonolite contains 30 to 40 wt% of phenocrysts, mainly anorthoclase, diopside, and magnetite. We have investigated pressures from 100 to 250 MPa, temperatures from 750 to 925°C, water contents from 1.3 to 10 wt%, at an oxygen fugacity (fO₂) of 1 log unit above the Ni-NiO solid buffer. Comparison of the natural and experimental phase proportions and compositions indicates that the phonolite was stored at 900 ± 20°C, 150 ± 50 MPa, 3 ± 0.5 wt% dissolved H₂O in the melt. The fO₂ was probably close to the fayalite-magnetite-quartz solid buffer judging from results of other experimental studies. These conditions constrain the magma storage depth at about 5 ± 1 km below current summit of Teide volcano. Given that the island has not suffered any major structural or topographic changes since the Lavas Negras eruption, any remaining magma from this event should still be stored at such depth and probably with a similar thermal and rheological state.     

Array analyses of volcanic earthquakes and tremor recorded at Las Cañadas caldera (Tenerife Island, Spain) during the 2004 seismic activation of Teide volcano

We analyze data from three seismic antennas deployed in Las Cañadas caldera (Tenerife) during May–July 2004. The period selected for the analysis (May 12–31, 2004) constitutes one of the most active seismic episodes reported in the area, except for the precursory seismicity accompanying historical eruptions. Most seismic signals recorded by the antennas were volcano-tectonic (VT) earthquakes. They usually exhibited low magnitudes, although some of them were large enough to be felt at nearby villages. A few long-period (LP) events, generally associated with the presence of volcanic fluids in the medium, were also detected. Furthermore, we detected the appearance of a continuous tremor that started on May 18 and lasted for several weeks, at least until the end of the recording period. It is the first time that volcanic tremor has been reported at Teide volcano. This tremor was a small-amplitude, narrow-band signal with central frequency in the range 1–6 Hz. It was detected at the three antennas located in Las Cañadas caldera. We applied the zero-lag cross-correlation (ZLCC) method to estimate the propagation parameters (back-azimuth and apparent slowness) of the recorded signals. For VT earthquakes, we also determined the S–P times and source locations. Our results indicate that at the beginning of the analyzed period most earthquakes clustered in a deep volume below the northwest flank of Teide volcano. The similarity of the propagation parameters obtained for LP events and these early VT earthquakes suggests that LP events might also originate within the source volume of the VT cluster. During the last two weeks of May, VT earthquakes were generally shallower, and spread all over Las Cañadas caldera. Finally, the analysis of the tremor wavefield points to the presence of multiple, low-energy sources acting simultaneously. We propose a model to explain the pattern of seismicity observed at Teide volcano. The process started in early April with a deep magma injection under the northwest flank of Teide volcano, related to a basaltic magma chamber inferred by geological and geophysical studies. The stress changes associated with the injection produced the deep VT cluster. In turn, the occurrence of earthquakes permitted an enhanced supply of fresh magmatic gases toward the surface. This gas flow induced the generation of LP events. The gases permeated the volcanic edifice, producing lubrication of pre-existing fractures and thus favoring the occurrence of VT earthquakes. On May 18, the flow front reached the shallow aquifer located under Las Cañadas caldera. The induced instability constituted the driving mechanism of the observed tremor.     

High CO2 Levels in Boreholes at El Teide Volcano Complex (Tenerife, Canary Islands): Implications for Volcanic Activity Monitoring

Emissions of CO₂ have been known for more than a hundred years as fumarolic activity at the terminal crater of El Teide volcano and as diffuse emissions at numerous water prospection drillings in the volcanic island of Tenerife. Large concentrations of CO₂ (>10% in volume) have been found inside galleries, long horizontal tunnels excavated for water mining. However, CO₂ concentrations of only 2900 ppm have been observed at the surface of the central region of the island (Las Cañadas del Teide caldera). In this work we analysed CO₂ concentrations in the subsurface of Las Cañadas caldera, in an attempt to study the vertical distribution of carbon dioxide and, in particular, the low emissions at the surface. This has been done through a series of 17 vertical profiles in two deep boreholes excavated in the Caldera. We found high levels of CO₂, varying in time from 13 vol% up to 40 vol% in different profiles directly above the water table, while no significant concentrations were detected above the thermal inversion that takes places in both boreholes at approximately 100 m from the water table. Water analyses also showed high dissolved CO₂ levels in equilibrium with the air, and an average 13C value in DIC of +4.7 (PDB), apparently induced by fast CO₂ degassing in the bicarbonated water.     

Methodology for the computation of volcanic susceptibility: An example for mafic and felsic eruptions on Tenerife (Canary Islands)

A new method to calculate volcanic susceptibility, i.e. the spatial probability of vent opening, is presented. Determination of volcanic susceptibility should constitute the first step in the elaboration of volcanic hazard maps of active volcanic fields. Our method considers different criteria as possible indicators for the location of future vents, based on the assumption that these locations should correspond to the surface expressions of the most likely pathways for magma ascent. Thus, two groups of criteria have been considered depending on the time scale (short or long term) of our approach. The first one accounts for long-term hazard assessment and corresponds to structural criteria that provide direct information on the internal structure of the volcanic field, including its past and present stress field, location of structural lineations (fractures and dikes), and location of past eruptions. The second group of criteria concerns to the computation of susceptibility for short term analyses (from days to a few months) during unrest episodes, and includes those structural and dynamical aspects that can be inferred from volcano monitoring. Thus, a specific layer of information is obtained for each of the criteria used. The specific weight of each criterion on the overall analysis depends on its relative significance to indicate pathways for magma ascent, on the quality of data and on their degree of confidence. The combination of the different data layers allows to create a map of the spatial probability of future eruptions based on objective criteria, thus constituting the first step to obtain the corresponding volcanic hazards map. The method has been used to calculate long-term volcanic susceptibility on Tenerife (Canary Islands), and the results obtained are also presented.     

Time Evolution of Deformation Using Time Series of Differential Interferograms: Application to La Palma Island (Canary Islands)

Differential interferometry is a very powerful tool for detecting changes in the Earth’s crust where coherence conditions are good, but is difficult to employ in some volcanic areas due to dense vegetation. We apply two differential InSAR methods using the time series associated with the interferograms to perform a phase analysis on a data set for La Palma island (Canary Islands) from the ERS-1 and ERS-2 European Space Agency (ESA) satellites for the time period 1992 to 2000. Both methods involve choosing a master image from the database and creating a series of interferograms with respect to this image. The “Coherent Pixel Time Series” (CPTS) technique chooses pixels with good average coherence, aligns the unwrapped interferograms with a stable area and then performs an inversion to calculate the linear velocity to quantify the deformation. The Coherent Target Modeling (CTM) method calculates the temporal coherence of each pixel to identify stable targets and then determines the best velocity for each pixel by using a linear fit that maximizes the temporal coherence. Using these two methods we have been able to detect deformation on La Palma Island that has been previously undetectable by conventional InSAR methods. There is a roughly circular region on the Southern part of the island that is actively deforming at ~ −4 to −8 mm/yr. This region is located near the Teneguia valcano, the host of the last known eruption on La Palma in 1971. A thorough investigation of the possible sources for this deformation revealed that it was most likely created by a subsurface thermal source.     

Reply to ‘Comment on: Rainfall,fog and throughfall dynamics in a subtropical ridge top cloud forest,National Park of Garajonay (La Gomera,Canary Islands,Spain)’ by A. Ritter and C.M. Regalado

This Reply aims to clarify and address many of the errors introduced to the discussion of our original paper (García‐Santos and Bruijnzeel, 2011; referred to henceforth as GSB2011) in the comment by Ritter and Regalado (2012) (henceforth referred to as RRR2012) as explained more fully under replies nos. 1–6, 8–12 and 14–18 below. Two (largely inconsequential) shortcomings in our original paper as detected by RRR2012 are acknowledged (see replies nos. 7 and 13 below). It is concluded that the majority of the comments offered by RRR2012 are unfounded as well as unnecessary. In summary, the two errors in GSB2011 spotted by RRR2012 (see points 7 and 13 below) could easily have been addressed in an erratum. Copyright © 2013 John Wiley & Sons, Ltd.     

Temporal variations in magma composition at Merapi Volcano (Central Java, Indonesia): magmatic cycles during the past 2000 years of explosive activity

Merapi Volcano (Central Java, Indonesia) has been frequently active during Middle to Late Holocene time producing basalts and basaltic andesites of medium-K composition in earlier stages of activity and high-K magmas from 1900 ¹⁴C yr BP to the present. Radiocarbon dating of pyroclastic deposits indicates an almost continuous activity with periods of high eruption rates alternating with shorter time spans of distinctly reduced eruptive frequency since the first appearance of high-K volcanic rocks. Geochemical data of 28 well-dated, prehistoric pyroclastic flows of the Merapi high-K series indicate systematic cyclic variations. These medium-term compositional variations result from a complex interplay of several magmatic processes, which ultimately control the periodicity and frequency of eruptions at Merapi. Low eruption rates and the absence of new influxes of primitive magma from depth allow the generation of basaltic andesite magma (56–57 wt% SiO₂) in a small-volume magma reservoir through fractional crystallisation from parental mafic magma (52–53 wt% SiO₂) in periods of low eruptive frequency. Magmas of intermediate composition erupted during these stages provide evidence for periodic withdrawal of magma from a steadily fractionating magma chamber. Subsequent periods are characterised by high eruption rates that coincide with shifts of whole-rock compositions from basaltic andesite to basalt. This compositional variation is interpreted to originate from influxes of primitive magma into a continuously active magma chamber, triggering the eruption of evolved magma after periods of low eruptive frequency. Batches of primitive magma eventually mix with residual magma in the magmatic reservoir to decrease whole-rock SiO₂ contents. Supply of primitive magma at Merapi appears to be sufficiently frequent that andesites or more differentiated rock types were not generated during the past 2000 years of activity. Cyclic variations also occurred during the recent eruptive period since AD 1883. The most recent eruptive episode of Merapi is characterised by essentially uniform magma compositions that may imply the existence of a continuously active magma reservoir, maintained in a quasi-steady state by magma recharge. The whole-rock compositions at the upper limit of the total SiO₂ range of the Merapi suite could also indicate the beginning of another period of high eruption rates and shifts towards more mafic compositions.     

A study of melt inclusions at Vulcano (Aeolian Islands, Italy): insights on the primitive magmas and on the volcanic feeding system

 This work presents the results of a microthermometric and EPMA-SIMS study of melt inclusions in phenocrysts of rocks of the shoshonitic eruptive complex of Vulcano (Aeolian Islands, Italy). Different primitive magmas related to two different evolutionary series, an older one (50–25 ka) and a younger one (15 ka to 1890 A.D.), were identified as melt inclusions in olivine Fo_88–91 crystals. Both are characterized by high Ca/Al ratio and present very similar Rb/Sr, B/Be and patterns of trace elements, with Nb and Ti anomalies typical of a subduction zone. The two basalts present the same temperature of crystallization (1180±20  °C) and similar volatile abundances. The H₂O, S and Cl contents are relatively high, whereas magmatic CO₂ concentrations are very low, probably due to CO₂ loss before low-pressure crystallization and entrapment of melt inclusions. The mineral chemistry of the basaltic assemblages and the high Ca/Al ratio of melt inclusions indicate an origin from a depleted, metasomatized clinopyroxene-rich peridotitic mantle. The younger primitive melt is characterized with respect to the older one by higher K₂O and incompatible element abundances, by lower Zr/Nb and La/Nb, and by higher Ba/Rb and LREE enrichment. A different degree of partial melting of the same source can explain the chemical differences between the two magmas. However, some anomalies in Sr, Rb and K contents suggest either a slightly different source for the two magmas or differing extents of crustal contamination. Low-pressure degassing and cooling of the basaltic magmas produce shoshonitic liquids. The melt inclusions indicate evolutionary paths via fractional crystallization, leading to trachytic compositions during the older activity and to rhyolitic compositions during the recent one. The bulk-rock compositions record a more complex history than do the melt inclusions, due to the syneruptive mixing processes commonly affecting the magmas erupted at Vulcano. The composition and temperature data on melt inclusions suggest that in the older period of activity several shallow magmatic reservoirs existed; in the younger one a relatively homogeneous feeding system is active. The shallow magmatic reservoir feeding the recent eruptive activity probably has a vertical configuration, with basaltic magma in the deeper zones and differentiated magmas in shallower, low-volume, dike-like reservoirs. Received: 11 March 1998 / Accepted: 14 July 1998     

Comment on ‘García‐Santos G,Bruijnzeel LA. 2011. Rainfall,fog and throughfall dynamics in a subtropical ridge top cloud forest,National Park of Garajonay (La Gomera,Canary Islands,Spain). Hydrological Processes 25: 411–417’

A previous study aiming to characterize the water dynamics of a cloud forest in the Garajonay National Park (La Gomera) from measurements taken in a plot located in the upper part of a selected watershed within the park is here commented. Reported magnitudes of hydrological variables and conclusions based on them are in disagreement with those of numerous studies carried out previously at the same site. Large data dispersion and inapplicability of some of the hypothesis assumed are shown to invalidate most of the results. Copyright © 2011 John Wiley & Sons, Ltd.