Volcanic evolution of the island of Tenerife (Canary Islands) in the light of new K-Ar data

New age determinations from Tenerife, together with those previously published (93 in all), provide a fairly comprehensive picture of the volcanic evolution of the island. The oldest volcanic series, with ages starting in the late Miocene, are formed mainly by basalts with some trachytes and phonolites which appear in Anaga, Teno and Roque del Conde massifs. In Anaga (NE), three volcanic cycles occurred: one older than 6.5 Ma, a second one between 6.5 and 4.5 Ma, with a possible gap between 5.4 and 4.8 Ma, and a late cycle around 3.6 Ma. In Teno (NW), after some undated units, the activity took place between 6.7 and 4.5 Ma, with two main series separated by a possible pause between 6.2 and 5.6 Ma. In the zone of Roque del Conde (S), the ages are scattered between 11.6 and 3.5 Ma. Between 3.3 and 1.9 Ma, the whole island underwent a period of volcanic quiescence and erosion.The large Cañadas volcano, made up of basalts, trachytes and phonolites, was built essentially between 1.9 and 0.2 Ma. To the NE of this central volcano, linking it with Anaga, is a chain of basaltic emission centers, with a peak of activity around 0.8 Ma. The Cañadas Caldera had several collapse phases, associated with large ignimbrite emissions. There were, at least, an older phase more than 1 Ma old, on the western part of the volcano, and a younger one, less than 0.6 Ma old, in the eastern side. The two large “valleys” of Guimar and la Orotava were formed by large landslides less than 0.8 Ma ago, and probably before 0.6 Ma ago. The present Cañadas caldera was formed by another landslide, less than 0.2 Ma ago. This caldera was later filled by the huge Teide volcano, which has been active even in historic times. During the same period a series of small volcanoes erupted at scattered locations throughout the island.The average eruptive rate in Tenerife was 0.3 km³/ka, with relatively small variations for the different eruptive periods. This island and La Gomera represent a model of growth by discontinuous pulses of volcanic activity, separated by gaps often coinciding with episodes of destruction of the edifices and sometimes extended for several million years. The neighbouring Gran Canaria, on the other hand, had an initial, rapid “shield-building phase” during which more than 90% of the island was built, and a series of smaller pulses at a much later period.A comparison between these three central islands indicates that the previously postulated westward displacement in time of a gap in the volcanic activity is valid only as a first approximation. Several gaps are present on each island, overlapping in time and not clearly supporting either of the models proposed to explain the evolution of the Canaries.     

Ar–Ar ages and geochemistry of the basaltic shield stage of Tenerife, Canary Islands, Spain

We report the first ³⁹Ar–⁴⁰Ar ages from the three early basic shield-like massifs of Tenerife, Canary islands, and couple these with detailed major and trace element chemistry to constrain the nature and timing of the mantle melting processes. The massifs have chemically different sources, and independent evolutionary histories. The Teno and Roque del Conde massifs appear chemically to represent the products of single mantle melting cycles, with progressive decrease in mean melt fraction and increase in mean melting depth in younger rocks. The Teno massif (NW) was erupted in a short time period around 6.0–6.4 Ma, while at least the lower half of the Roque del Conde massif (SW) is older than 11 Ma. In contrast, the Anaga massif (NE) is polygenetic, with ³⁹Ar–⁴⁰Ar ages ranging from 8.0–4.2 Ma, and no simple stratigraphic chemical progression. These ages run counter to published suggestions of progressive younging of Canary shield stages to the southwest. Basic rocks in all three massifs are the result of much deeper melting and smaller melt fractions than equivalent units in Gran Canaria, but nevertheless the melting column must have extended significantly into the spinel facies, requiring substantial disruption of the local lithosphere. The age and melting relationships broadly support the mantle blob model for Canary magmatism proposed by Hoernle and Schmincke (Hoernle, K., Schminke, H.-U., 1993. The role of partial melting in the 15-Ma geochemical evolution of Gran Canaria: a blob model for the Canary hotspot. J. Petrol. 34, 599–626). In all three massifs, extensive fractional crystallisation has taken place at crustal levels so that mean MgO contents are only some 6–7%. The fractionation sequence is olivine–clinopyroxene–magnetite in basaltic compositions, with the involvement of plagioclase, amphibole and apatite only to generate the infrequent more evolved hawaiites to benmoreites. Despite the abundance of basanitic magmas in the Tenerife older massifs, these follow a differentiation trend towards weakly undersaturated benmoreite rather than to phonolite. This probably reflects early crystallisation of magnetite, perhaps resulting from somewhat high oxygen fugacity. The chemical evidence for replenished magma chambers in Tenerife described by Neumann et al. (Neumann, E.R., Wulff-Oedersen, E., Simonsen, S.L., Pearson, N.J., Martí, J., Mitjavila, J., 1999. Evidence for fractional crystallisation of periodically refilled magma chambers in Tenerife, Canary Islands. J. Petrol. 40, 1089–1123) is a consequence of treating as a single cogenetic suite the products of several magmatic systems that differ in parental melt fraction.     

Implications for the early shield-stage evolution of Tenerife from K/Ar ages and magnetic stratigraphy

The combined use of field geology, radioisotopic dating and magnetic stratigraphy applied to the old shield volcanoes of Tenerife provides a reliable time framework for the early, shield-stage evolution of the island. The greater part of this new set of ages, obtained from sequences of lava flows is in agreement with the astronomical polarity time scale. This approach illustrates that previous K-Ar data collected without a comprehensive stratigraphy should be viewed with caution, and in some cases discarded altogether. The shield volcanoes of Tenerife encompass a relatively small number of magnetozones, an observation consistent with the relatively short periods of growth shown by the new ages (1-2 my). The island was constructed by the aggregation of three successive shields: the Roque del Conde (Central shield), between about 11.9 and 8.9 Ma, and the Teno (6.2-5.6 Ma) and Anaga (4.9-3.9 Ma) volcanoes. This new oldest subaerial age of Tenerife fits with the others obtained in the Canaries in a clear west to east monotonous age progression, one of the main restrictions for hotspot-related island chains.     

Dykes, faults and palaeostresses in the Teno and Anaga massifs of Tenerife (Canary Islands)

A structural field study was made of 578 sheet intrusions (mostly dykes) and 153 (mostly normal) faults dissecting the Anaga and Teno massifs, where a complex volcanic succession of Tertiary age (the ‘Old Basaltic Series’) representing the shield-building stage of Tenerife (Canary Islands) crops out. Many of the intrusions, mostly sub-vertical mafic dykes, are emplaced by multiple magma injections, with cumulative thicknesses mostly less than 2 m. Dyke tips are exposed and preserved for 12% of the dykes. Three differently oriented sets of dykes exist in the Anaga massif (NNW–SSE, NNE–SSW, E–W), whereas there is only one main set in Teno, trending NNW–SSE. Dyke swarms and other structural features having similar orientations also exist in other Canary Islands. A minimum value of the horizontal component of extension induced by dykes is computed using a step of 5° of azimuth, accounting also for the dip of dykes. The cumulative crustal dilation is at least 300 m (4%) in Anaga and 270 m (6%) in Teno; the maximum extension peaks at N75° in Anaga and N60° in Teno, indicating a general prevailing extension in direction ENE–WSW. Most of the measured faults are normal and strike NNW–SSE. Computation of palaeostresses from inversion of fault-slip data sets suggests the existence of a polyphase brittle deformation due to an extensional stress field with the minimum compressive principal axes trending NE–SW and WNW–ESE.     

Evolution of the Cañadas edifice and its implications for the origin of the Cañadas Caldera (Tenerife, Canary Islands)

The volcano-stratigraphic and geochronologic data presented in this work show that the Tenerife central zone has been occupied during the last 3 Ma by shield or central composite volcanoes which reached more than 3000 m in height. The last volcanic system, the presently active Teide-Pico Viejo Complex began to form approximately 150 ka ago. The first Cañadas Edifice (CE) volcanic activity took place between about 3.5 Ma and 2.7 Ma. The CE-I is formed mainly by basalts, trachybasalts and trachytes. The remains of this phase outcrop in the Cañadas Wall (CW) sectors of La Angostura (3.5–3.0 Ma and 3.0–2.7 Ma), Boca de Tauce (3.0 Ma), and in the bottom of some external radial ravines (3.5 Ma). The position of its main emission center was located in the central part of the CC. The volcano could have reached 3000 m in height. This edifice underwent a partial destruction by failure and flank collapse, forming debris-avalanches during the 2.6–2.3 Ma period. The debris-avalanche deposits can be seen in the most distal zones in the N flank of the CE-I (Tigaiga Breccia). A new volcanic phase, whose deposits overlie the remains of CE-I and the former debris-avalanche deposits, constituted a new volcanic edifice, the CE-II. The dyke directions analysis and the morphological reconstruction suggest that the CE-II center was situated somewhat westward of the CE-I, reaching some 3200 m in height. The CE-II formations are well exposed on the CW, especially at the El Cedro (2.3–2.00 Ma) sector. They are also frequent in the S flank of the edifice (2.25–1.89 Ma) in Tejina (2.5–1.87 Ma) as well as in the Tigaiga massif to the N (2.23 Ma). During the last periods of activity of CE-II, important explosive eruptions took place forming ignimbrites, pyroclastic flows, and fall deposits of trachytic composition. Their ages vary between 1.5 and 1.6 Ma (Adeje ignimbrites, to the W). In the CW, the Upper Ucanca phonolitic Unit (1.4 Ma) could be the last main episode of the CE-II. Afterwards, the Cañadas III phase began. It is well represented in the CW sectors of Tigaiga (1.1 Ma–0.27 Ma), Las Pilas (1.03 Ma–0.78 Ma), Diego Hernández (0.54 Ma–0.17 Ma) and Guajara (1.1 Ma–0.7 Ma). The materials of this edifice are also found in the SE flank. These materials are trachybasaltic lava-flows and abundant phonolitic lava and pyroclastic flows (0.6 Ma–0.5 Ma) associated with abundant plinian falls. The CE-III was essentially built between 0.9 and 0.2 Ma, a period when the volcanic activity was also intense in the ‘Dorsal Edifice' situated in the easterly wing of Tenerife. The so called ‘valleys' of La Orotava and Güimar, transversals to the ridge axis, also formed during this period. In the central part of Tenerife, the CE-III completed its evolution with an explosive deposit resting on the top of the CE, for which ages from 0.173 to 0.13 Ma have been obtained. The CC age must be younger due to the fact that the present caldera scarp cuts these deposits. On the controversial origin of the CC (central vertical collapse vs. repeated flank failure and lateral collapse of mature volcanic edifices), the data discussed in this paper favor the second hypothesis. Clearly several debris-avalanche type events exist in the history of the volcano but most of the deposits are now under the sea. The caldera wall should represent the proximal scarps of the large slides whose intermediate scarps are covered by the more recent Teide-Pico Viejo volcanoes.     

Rift zone reorganization through flank instability in ocean island volcanoes: an example from Tenerife, Canary Islands

The relationship between rift zones and flank instability in ocean island volcanoes is often inferred but rarely documented. Our field data, aerial image analysis, and ⁴⁰Ar/³⁹Ar chronology from Anaga basaltic shield volcano on Tenerife, Canary Islands, support a rift zone—flank instability relationship. A single rift zone dominated the early stage of the Anaga edifice (~6–4.5 Ma). Destabilization of the northern sector led to partial seaward collapse at about ~4.5 Ma, resulting in a giant landslide. The remnant highly fractured northern flank is part of the destabilized sector. A curved rift zone developed within and around this unstable sector between 4.5 and 3.5 Ma. Induced by the dilatation of the curved rift, a further rift-arm developed to the south, generating a three-armed rift system. This evolutionary sequence is supported by elastic dislocation models that illustrate how a curved rift zone accelerates flank instability on one side of a rift, and facilitates dike intrusions on the opposite side. Our study demonstrates a feedback relationship between flank instability and intrusive development, a scenario probably common in ocean island volcanoes. We therefore propose that ocean island rift zones represent geologically unsteady structures that migrate and reorganize in response to volcano flank instability.Editorial responsibility: T. DruittThis revised version was published online in February 2005 with typographical corrections and a changed wording.     

The felsic dikes of La Gomera (Canary Islands): identification of cone sheet and radial dike swarms

On the northern part of La Gomera there exists a great abundance of trachytic–phonolitic dikes showing a broad diversity in dip and strike. Several methods have been applied in order to separate these dikes in different sets, localise the area from where they derive, and reconstruct the geometry of the swarms. The oldest dikes correspond to a radial swarm dated at 8 Ma. The felsic activity migrated then southwestwards and a second radial swarm and a cone sheet complex were developed between 7.5 and 6.4 Ma ago. The cone sheet complex is 10 km in diameter and shared its centre with that of the second radial structure. The cone sheets exhibit an outward decrease of dip angle whilst every individual sheet maintains a constant inclination. This geometry reflects the existence of an ancient single dome-shaped shallow magma chamber situated some 1650 m below present sea level. The eastern radial swarm represents a felsic episode that could mark the ending of the Lower Old Basalts, the earlier subaerial activity of La Gomera. The two other dike swarms represent a younger episode coeval with the Upper Old Basalts.     

Subsurface geometry and structural evolution of La Gomera island based on gravity data

We hereby present a new Bouguer gravity map of the La Gomera island (Canarian Archipelago), which is analysed and interpreted by means of a 3-dimensional inversion, in order to contribute to the knowledge of the structural setting of this volcanic island and its evolutionary history. A land gravity data set covering the whole island of La Gomera is used in combination with offshore measurements to achieve a better determination of the gravity field in areas near the coasts.The study of this map let us to shed some light on the hypothesis established about the volcanism of this island. Moreover, it shows the information that is hidden from a geological surface exploration, modelling deep sections of the crust in La Gomera, which have been unknown until now.A first interpretation of the Bouguer gravity anomaly is achieved from 1) the residual gravity map calculated by removing a regional component and 2) the total horizontal gradient of the gravity. These residual and derivative maps allow us to identify the horizontal location and borders of the shallowest gravity sources. This provides a useful tool to study the structures associated to the latest periods of the volcanism in the area. Moreover, the information so obtained supports the hypothesis about the migration of volcanic activity towards the south of the island.Subsequently, an inversion process is carried out looking for the 3D-modelisation of the sources of the observed gravity field, which provides a comprehensive view of the structures in volcanic environments. The inversion technique used is based on a genetic algorithm (GA) applied upon a prismatic partition of the subsoil volume, and adopting a priori values of density contrast (positive and negative). The main advantage of this method is that let us to model deep and shallow bodies which exhibit very different geometries and density contrasts. So, results indicate that this inversion strategy can be very effective for characterization of volcanic structures, improving the information from previous geologic and volcanologic studies. The inversion model obtained shows correlation between several sources of the gravity field and the volcanic units associated with the growth of La Gomera Island. The main gravity source of this model is associated with the oldest unit, called the Basal Complex. This unit corresponds to the first submarine growth stage and it is modelled as the most important and deepest high density structure. According to previous geological studies, the following edifice (Old Edifice) was also submarine in its initial phases, later being represented by a wide basaltic shield volcano. The original location and morphology of this Old Edifice is deduced from the distribution of positive density contrasts that appears in the model. Moreover, other gravity field sources are identified and associated to several feeding systems of this stage of the volcanism in La Gomera. The shallowest sections of the model let us recognise the distribution of light material inside the Vallehermoso caldera, surrounded by high density structures. This gives us some insight into the internal structure and morphology of the caldera, pointing to a vertical collapse origin followed by erosion and other destructive processes. Finally, other conclusions are obtained from the correlation found between the sources of the gravity field and the migration of the volcanic activity towards the southern area of the island.     

Deformation structures associated with the Tazo landslide (La Gomera, Canary Islands)

Deformation structures below the basal plane of gravitational slides can provide useful information about the state of stress undergone by rocks prior to the sliding process and about the triggering forces acting at each particular sliding event. In the present work we conducted a structural analysis of the rocks below the surface of the gravitational slide of Tazo (La Gomera, Canary Islands) and determined the epigenetic processes involved in the filling of the amphitheatre. We also inferred the possible triggering phenomena related to the Tazo landslide. The rocks located below the surface of the gravitational slide of Tazo -i.e., the basaltic lava flows, sills and dikes of the Lower Old Edifice and the submarine volcanic rocks, gabbros, pyroxenites and dikes of the Basal Complex of La Gomera- are strongly deformed close to this sliding surface. The lava flows and dikes of the Lower Old Edifice are folded, with fault breccias and gouges, and locally foliated, defining the sliding surface. The dikes of the Basal Complex are also folded, and the gabbros and pyroxenites are affected by a large number of small faults. In the Basal Complex, the sliding surface is defined by a foliated granular gouge. In the damage zone, the Basal Complex rocks show an incipient fracture cleavage. The sliding amphitheatre has been filled by the debris avalanche or cohesive debris flow generated within the slide, as well as by later debris flows, hyperconcentrated flows, sheet flows, and by interspersed lava flows from the Upper Old Edifice. We suggest here that the collapse of the north-western flank of the Lower Old Edifice at Tazo could in part have been triggered by continuous magma injection, associated with the emplacement of dikes in a rift zone with an ENE-WSW direction, enhanced by the mechanical weakness of the Basal Complex unit, which was affected by hydrothermal metamorphism under greenschist facies conditions and by the displacement along the Montaña de Alcalá and Guillama normal faults, which are deeply entrenched in the altered rocks of the Basal Complex.     

Cladistic analysis applied to the classification of volcanoes

Cladistics is a systematic method of classification that groups entities on the basis of sharing similar characteristics in the most parsimonious manner. Here cladistics is applied to the classification of volcanoes using a dataset of 59 Quaternary volcanoes and 129 volcanic edifices of the Tohoku region, Northeast Japan. Volcano and edifice characteristics recorded in the database include attributes of volcano size, chemical composition, dominant eruptive products, volcano morphology, dominant landforms, volcano age and eruptive history. Without characteristics related to time the volcanic edifices divide into two groups, with characters related to volcano size, dominant composition and edifice morphology being the most diagnostic. Analysis including time based characteristics yields four groups with a good correlation between these groups and the two groups from the analysis without time for 108 out of 129 volcanic edifices. Thus when characters are slightly changed the volcanoes still form similar groupings. Analysis of the volcanoes both with and without time yields three groups based on compositional, eruptive products and morphological characters. Spatial clusters of volcanic centres have been recognised in the Tohoku region by Tamura et al. (Earth Planet Sci Lett 197:105–106, 2002). The groups identified by cladistic analysis are distributed unevenly between the clusters, indicating a tendency for individual clusters to form similar kinds of volcanoes with distinctive but coherent styles of volcanism. Uneven distribution of volcano types between clusters can be explained by variations in dominant magma compositions through time, which are reflected in eruption products and volcanic landforms. Cladistic analysis can be a useful tool for elucidating dynamic igneous processes that could be applied to other regions and globally. Our exploratory study indicates that cladistics has promise as a method for classifying volcanoes and potentially elucidating dynamic and evolutionary volcanic processes. Cladistics may also have utility in hazards assessment where spatial distributions and robust definitions of a volcano are important, as in locating sensitive facilities such as nuclear reactors and repositories.     

Late Pleistocene flank collapse of Zempoala volcano (Central Mexico) and the role of fault reactivation

Zempoala is an extinct Pleistocene (∼ 0.7–0.8 Ma) stratovolcano that together with La Corona volcano (∼ 0.9 Ma) forms the southern end of the Sierra de las Cruces volcanic range, Central Mexico. The volcano consists of andesitic and dacitic lava flows and domes, as well as pyroclastic and epiclastic sequences, and has had a complex history with several flank collapses. One of these collapses occurred during the late Pleistocene on the S–SE flank of the volcano and produced the Zempoala debris avalanche deposit. This collapse could have been triggered by the reactivation of two normal fault systems (E–W and NE–SW), although magmatic activity cannot be absolutely excluded. The debris avalanche traveled 60 km to the south, covers an area of 600 km² and has a total volume of 6 km³, with a calculated Heim coefficient (H/L) of 0.03. Based on the textural characteristics of the deposit we recognized three zones: proximal, axial, and lateral distal zone. The proximal zone consists of debris avalanche blocks that develop a hummocky topography; the axial zone corresponds with the main debris avalanche deposit made of large clasts set in a sandy matrix, which transformed to a debris flow in the lateral distal portion. The deposit is heterolithologic in composition, with dacitic and andesitic fragments from the old edifice that decrease in volume as bulking of exotic clasts from the substratum increase. Several cities (Cuernavaca, Jojutla de Juárez, Alpuyeca) with associated industrial, agricultural, and tourism activities have been built on the deposit, which pose in evidence the possible impact in case of a new event with such characteristics, since the area is still tectonically active.     

The characteristics of magma reservoir failure beneath a volcanic edifice

Eruptions fed from subsurface reservoirs commonly construct volcanic edifices at the surface, and the growth of an edifice will in turn modify the subsurface stress state that dictates the conditions under which subsequent rupture of the inflating reservoir can occur. We re-examine this problem using axisymmetric finite element models of ellipsoidal reservoirs beneath conical edifices, explicitly incorporating factors (e.g., full gravitational loading conditions, an elastic edifice instead of a surface load, reservoir pressures sufficient to induce tensile rupture) that compromise previous solutions to illustrate why variations in rupture behavior can occur. Relative to half-space model results, the presence of an edifice generally rotates rupture toward the crest of a spherical reservoir, with increasing flank slope (for an edifice of constant volume) and larger edifices (or greater reservoir scaled depths) normally serving to enhance this trend. When non-spherical reservoirs are considered, the presence of an edifice amplifies previously identified half-space failure characteristics, shifting rupture to the crest more rapidly for prolate reservoirs while forcing rupture closer to the midpoint of oblate reservoirs. Rupture is always observed to occur in the σ_t orientation, and depending on where initial failure occurs rupture favors the initial emplacement of either lateral sills, circumferential intrusions or vertically ascending dikes. Ultimately, integration of our numerical model results with other information, for instance the sequence of intrusion/eruption events observed at a given volcano, can provide useful new insight into how a volcano's subsurface magma plumbing system evolved. We demonstrate this process through application of our model to Summer Coon, a well-studied stratocone on Earth, and Ilithyia Mons, a large conical shield volcano on Venus.     

New radiometric and petrological constraints on the evolution of the Pichincha volcanic complex (Ecuador)

Fieldwork, radiometric (⁴⁰Ar/³⁹Ar and ¹⁴C) ages and whole-rock geochemistry allow a reconstruction of eruptive stages at the active, mainly dacitic, Pichincha Volcanic Complex (PVC), whose eruptions have repeatedly threatened Quito, most recently from 1999 to 2001. After the emplacement of basal lavas dated at ∼1100 to 900 ka, the eruptive activity of the old Rucu Pichincha volcano lasted from ∼850 ka to ∼150 ka before present (BP) and resulted in a 15 × 20 km-wide edifice, which comprises three main building stages: (1) A lower stratocone (Lower Rucu, ∼160 km³ in volume) developed from ∼850 to 600 ka; (2) This edifice was capped by a steeper-sided and less voluminous cone (the Upper Rucu, 40–50 km³), the history of which started 450–430 ka ago and ended around 250 ka with a sector collapse; (3) A smaller (8–10 km³) but more explosive edifice grew in the avalanche amphitheatre and ended Rucu Pichincha's history about 150 ka ago. The Guagua Pichincha volcano (GGP) was developed from 60 ka on the western flank of Rucu with four growth stages separated by major catastrophic events. (1) From ∼60 to 47 ka, a basal effusive stratocone developed, terminating with a large ash-and-pumice flow event. (2) This basal volcano was followed by a long-lasting dome building stage and related explosive episodes, the latter occurring between 28–30 and 22–23 ka. These first two stages formed the main GGP (∼30 km³), a large part of which was removed by a major collapse 11 ka BP. (3) Sustained explosive activity and viscous lava extrusions gave rise to a new edifice, Toaza (4–5 km³ in volume), which in turn collapsed around 4 ka BP. (4) The ensuing amphitheatre was partly filled by the ∼1-km³ Cristal dome, which is the historically active centre of the Pichincha complex. The average output rate for the whole PVC is 0.29 km³/ka. Nevertheless, the chronostratigraphic resolution we obtained for Lower Rucu Pichincha and for the two main edifices of Guagua Pichincha (main GGP and Toaza), leads to eruptive rates of 0.60–0.65 km³/ka during these construction stages. These output rates are compared to those of other mainly dacitic volcanoes from continental arcs. Our study also supports an overall SiO₂ and large-ion lithophile elements enrichment as the PVC develops. In particular, distinctive geochemical signatures indicate the involvement of a new magma batch at the transition between Rucu and Guagua. At the GGP, the same phenomenon occurs at each major collapse event marking the onset of the ensuing magmatic stage. Since the 11-ka-BP collapse event, this magmatic behaviour has led to increasingly explosive activity. Four explosive cycles of between 100 and 200 years long have taken place at the Cristal dome in the past 3.7 ka, and repose intervals between these cycles have tended to decrease with time. As a consequence, we suggest that the 1999–2001 eruptive period may have initiated a new eruptive cycle that might pose a future hazard to Quito (∼2 million inhabitants).     

The geological evolution of Merapi volcano, Central Java, Indonesia

Merapi is an almost persistently active basalt to basaltic andesite volcanic complex in Central Java (Indonesia) and often referred to as the type volcano for small-volume pyroclastic flows generated by gravitational lava dome failures (Merapi-type nuées ardentes). Stratigraphic field data, published and new radiocarbon ages in conjunction with a new set of ⁴⁰K–⁴⁰Ar and ⁴⁰Ar–³⁹Ar ages, and whole-rock geochemical data allow a reassessment of the geological and geochemical evolution of the volcanic complex. An adapted version of the published geological map of Merapi [(Wirakusumah et al. 1989), Peta Geologi Gunungapi Merapi, Jawa Tengah (Geologic map of Merapi volcano, Central Java), 1:50,000] is presented, in which eight main volcano stratigraphic units are distinguished, linked to three main evolutionary stages of the volcanic complex—Proto-Merapi, Old Merapi and New Merapi. Construction of the Merapi volcanic complex began after 170?ka. The two earliest (Proto-Merapi) volcanic edifices, Gunung Bibi (109?±?60?ka), a small basaltic andesite volcanic structure on Merapi’s north-east flank, and Gunung Turgo and Gunung Plawangan (138?±?3?ka; 135?±?3?ka), two basaltic hills in the southern sector of the volcano, predate the Merapi cone sensu stricto. Old Merapi started to grow at ~30?ka, building a stratovolcano of basaltic andesite lavas and intercalated pyroclastic rocks. This older Merapi edifice was destroyed by one or, possibly, several flank failures, the latest of which occurred after 4.8?±?1.5?ka and marks the end of the Old Merapi stage. The construction of the recent Merapi cone (New Merapi) began afterwards. Mostly basaltic andesite pyroclastic and epiclastic deposits of both Old and New Merapi (<11,792?±?90 ¹⁴C years BP) cover the lower flanks of the edifice. A shift from medium-K to high-K character of the eruptive products occurred at ~1,900 ¹⁴C years BP, with all younger products having high-K affinity. The radiocarbon record points towards an almost continuous activity of Merapi since this time, with periods of high eruption frequency interrupted by shorter intervals of apparently lower eruption rates, which is reflected in the geochemical composition of the eruptive products. The Holocene stratigraphic record reveals that fountain collapse pyroclastic flows are a common phenomenon at Merapi. The distribution and run-out distances of these flows have frequently exceeded those of the classic Merapi-type nuées ardentes of the recent activity. Widespread pumiceous fallout deposits testify the occurrence of moderate to large (subplinian) eruptions (VEI 3–4) during the mid to late Holocene. VEI 4 eruptions, as identified in the stratigraphic record, are an order of magnitude larger than any recorded historical eruption of Merapi, except for the 1872?AD and, possibly, the October–November 2010 events. Both types of eruptive and volcanic phenomena require careful consideration in long-term hazard assessment at Merapi.     

Major and trace element geochemistry of neutral and acidic thermal springs at El Chichón volcano,Mexico: Implications for monitoring of the volcanic activity

Four groups of thermal springs with temperatures from 50 to 80 °C are located on the S–SW–W slopes of El Chichón volcano, a composite dome-tephra edifice, which exploded in 1982 with a 1 km wide, 160 m deep crater left. Very dynamic thermal activity inside the crater (variations in chemistry and migration of pools and fumaroles, drastic changes in the crater lake volume and chemistry) contrasts with the stable behavior of the flank hot springs during the time of observations (1974–2005). All known groups of hot springs are located on the contact of the basement and volcanic edifice, and only on the W–SW–S slopes of the volcano at almost same elevations 600–650 m asl and less than 3 km of direct distance from the crater. Three groups of near-neutral (pH ≈ 6) springs at SW–S slopes have the total thermal water outflow rate higher than 300 l/s and are similar in composition. The fourth and farthest group on the western slope discharges acidic (pH ≈ 2) saline (10 g/kg of Cl) water with a much lower outflow rate (< 10 l/s).     

Growth and evolution of an emergent tuff cone: Considerations from structural geology, geomorphology and facies analysis of São Roque volcano, São Miguel (Azores)

The syn-eruptive and post-eruptive history of São Roque tuff cone, its geological setting and volcanological features were studied in detail to understand the role played by the different factors that contributed to the morphological evolution of this relatively simple and small volcanic edifice.In addition, attention was also focused on the series of natural changes that affected the tuff cone during the course of the years and that finally led to its structural disassembly. A novel model is proposed to explain this process.The São Roque volcanic centre, located on the island of São Miguel (Azores), consists of two well-consolidated bodies and numerous small islets that formed more than 4700 years ago during the hydromagmatic activity that took place along an intruding dyke, whose NNW–SSE trend is in agreement with the regional tectonic pattern. The eruptive vents probably migrated progressively from SSE to NNW, forming small edifices through the rapid accumulation of sediments during alternating phases of “dry” and “wet” magmatic emissions. Syn-eruptive partial collapses greatly modified the original morphological structure of these edifices, probably allowing sea water to continuously flow into the vents. The complex interaction of these factors controlled the depth of magma fragmentation, producing different types of deposits, in which the ash-lapilli ratio varies considerably. The high-water saturation degree of these deposits caused syn-eruptive and post-eruptive remobilization which resulted in collapses and some small-scale landslides.Post-eruptive, WNW–ESE trending transtensional and extensional tectonic activities operated during the initial dissection of the cone, generating instability. Furthermore, the rapid accumulation of “wet” tephra, and its following consolidation, caused selective collapses that favoured the fragmentation of the deposit and caused the formation of numerous islets separated by radially-arranged channels. Collapses also involved the lava units emplaced in more recent times around the tuff cone, which show that brittle deformation has been significant in the area for a prolonged period.     

Relationships between volcano gravitational spreading and magma intrusion

Volcano spreading, with its characteristic sector grabens, is caused by outward flow of weak substrata due to gravitational loading. This process is now known to affect many present-day edifices. A volcano intrusive complex can form an important component of an edifice and may induce deformation while it develops. Such intrusions are clearly observed in ancient eroded volcanoes, like the Scottish Palaeocene centres, or in geophysical studies such as in La Réunion, or inferred from large calderas, such as in Hawaii, the Canaries or Galapagos volcanoes. Volcano gravitational spreading and intrusive complex emplacement may act simultaneously within an edifice. We explore the coupling and interactions between these two processes. We use scaled analogue models, where an intrusive complex made of Golden syrup is emplaced within a granular model volcano based on a substratum of a ductile silicone layer overlain by a brittle granular layer. We model specifically the large intrusive complex growth and do not model small-scale and short-lived events, such as dyke intrusion, that develop above the intrusive complex. The models show that the intrusive complex develops in continual competition between upward bulging and lateral gravity spreading. The brittle substratum strongly controls the deformation style, the intrusion shape and also controls the balance between intrusive complex spreading and ductile layer-related gravitational spreading. In the models, intrusive complex emplacement and spreading produce similar structures to those formed during volcano gravitational spreading alone (i.e. grabens, folds, en échelon fractures). Therefore, simple analysis of fault geometry and fault kinetic indicators is not sufficient to distinguish gravitational from intrusive complex spreading, except when the intrusive complex is eccentric from the volcano centre. However, the displacement fields obtained for (1) a solely gravitational spreading volcano and for (2) a gravitational spreading volcano with a growing and spreading intrusive complex are very different. Consequently, deformation fields (like those obtained from geodetic monitoring) can give a strong indication of the presence of a spreading intrusive complex. We compare the models with field observations and geophysical evidence on active volcanoes such as La Réunion Island (Indian Ocean), Ometepe Island (Nicaragua) and eroded volcanic remnants such as Ardnamurchan (Scotland) and suggest that a combination between gravitational and intrusive complex spreading has been active.     

Volcano collapse promoted by progressive strength reduction: new data from Mount St. Helens

Rock shear strength plays a fundamental role in volcano flank collapse, yet pertinent data from modern collapse surfaces are rare. Using samples collected from the inferred failure surface of the massive 1980 collapse of Mount St. Helens (MSH), we determined rock shear strength via laboratory tests designed to mimic conditions in the pre-collapse edifice. We observed that the 1980 failure shear surfaces formed primarily in pervasively shattered older dome rocks; failure was not localized in sloping volcanic strata or in weak, hydrothermally altered rocks. Our test results show that rock shear strength under large confining stresses is reduced ∼20% as a result of large quasi-static shear strain, as preceded the 1980 collapse of MSH. Using quasi-3D slope-stability modeling, we demonstrate that this mechanical weakening could have provoked edifice collapse, even in the absence of transiently elevated pore-fluid pressures or earthquake ground shaking. Progressive strength reduction could promote collapses at other volcanic edifices.     

A new scenario for the last magmatic eruption of La Soufrière of Guadeloupe (Lesser Antilles) in 1530 A.D. Evidence from stratigraphy radiocarbon dating and magmatic evolution of erupted products

La Soufrière of Guadeloupe is a dangerous volcano characterized over the last decade by moderate seismic and fumarolic unrest. In the last 15,000 years it has experienced phreatic and magmatic eruptions and unusually numerous flank collapse events sometimes associated with a magmatic eruption. We propose a new age of 1530 A.D. and a new eruptive scenario for the last magmatic eruption on the basis of a novel statistical analysis of radiocarbon age dates, and new field and geochemical data. This eruption is the only magmatic eruption likely to have occurred in Guadeloupe during the last 1400 years. The eruption mainly involved an andesitic magma which, in the first phase of the eruption, partially mixed with a slightly more differentiated magma stored in a small and shallow magma chamber. Ascent of magma to the surface generated a partial collapse of the hydrothermally altered edifice that increased the magma discharge and led to a sub-plinian phase with scoria fallout and column-collapse pyroclastic flows followed by near-vent pyroclastic scoria fountains. The eruption ended with growth of a lava dome. Our revised interpretation of the last magmatic eruption of La Soufrière constitutes the most likely key to a future magmatic eruption scenario for this volcano which displays strong evidence of unrest since 1992.