Stratigraphy, structure, and volcanic evolution of the Pico Teide–Pico Viejo formation, Tenerife, Canary Islands

The structure and volcanic stratigraphy of the Pico Teide–Pico Viejo (PT–PV) formation, deriving from the basanite–phonolite stratovolcanoes PT and PV, and numerous flank vent systems, are documented in detail based on new field and photogeologic mapping, geomorphologic analysis, borehole data, and petrological and geochemical findings. Results provide insight into the structure and evolution of the PT–PV magma system, and the long-term, cyclic evolution of Tenerife's post-shield volcanic complex. The PT–PV formation comprises products of central volcanism, mainly emplaced into the Las Cañadas caldera (LCC), and contemporaneous products from adjacent rifts. PT–PV central volcanic products become more differentiated up-section with felsic lavas dominating the recent output of the system. This is attributed to the evolution of a shallow magma reservoir beneath PT that was emplaced early in the PT–PV cycle on the intra-caldera segment of Tenerife's post-shield rift system. The rift axis has been the focus of PT–PV intrusive and eruptive activity, and has controlled the location of the stratocones. The current geometry of the rifts reflects a major structural reorganisation defining the start of the PT–PV cycle at 0.18 Ma, namely the truncation of the north side of the LCC/LCE by the giant Icod landslide. The internal stratigraphy of the PT–PV formation suggests that PT developed early, with PV developing as a satellite vent. Activity has since alternated between PT and PV due to episodes of vent blockage or chamber sealing. These processes have allowed significant volumes of phonolitic magmas to develop and accumulate within the PT chamber, which have vented through radial dike systems during tumescence episodes and from the rift system, which has permitted lateral magma transport. The PT–PV magma system is a potentially hazardous source of future, felsic eruptive activity on Tenerife.     

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.     

Estimating building exposure and impact to volcanic hazards in Icod de los Vinos,Tenerife (Canary Islands)

Principal and subsidiary building structure characteristics and their distribution have been inventoried in Icod, Tenerife (Canary Islands) and used to evaluate the vulnerability of individual buildings to three volcanic hazards: tephra fallout, volcanogenic earthquakes and pyroclastic flows. The procedures described in this paper represent a methodological framework for a comprehensive survey of all the buildings at risk in the area around the Teide volcano in Tenerife. Such a methodology would need to be implemented for the completion of a comprehensive risk assessment for the populations under threat of explosive eruptions in this area. The information presented in the paper is a sample of the necessary data required for the impact estimation and risk assessment exercises that would need to be carried out by emergency managers, local authorities and those responsible for recovery and repair in the event of a volcanic eruption. The data shows there are micro variations in building stock characteristics that would influence the likely impact of an eruption in the area. As an example of the use of this methodology for vulnerability assessment, we have applied a deterministic simulation model of a volcanic eruption from Teide volcano and its associated ash fallout which, when combined with the vulnerability data collected, allows us to obtain the vulnerability map of the studied area. This map is obtained by performing spatial analysis with a Geographical Information System (GIS). This vulnerability analysis is included in the framework of an automatic information system specifically developed for hazard assessment and risk management on Tenerife, but which can be also applied to other volcanic areas. The work presented is part of the EU-funded EXPLORIS project (Explosive Eruption Risk and Decision Support for EU Populations Threatened by Volcanoes, EVR1-2001-00047).     

Internal structure of Tenerife (Canary Islands) based on gravity, aeromagnetic and volcanological data

Gravity and magnetic methods have been applied to the Tenerife Island, to provide new information about its internal structure. For this study, 365 gravity stations covering the central part of the island have been selected. The anomalous density maps at different depths were obtained by means of an inversion global adjustment, on fixed density contrast, to describe the three-dimensional (3D) geometry of the anomalous bodies. On the other hand, several analysis techniques, such as reduction to the pole, spectral analysis, low-pass filtering, terrain correction and forward modelling, were applied to process the high-resolution data obtained in an aeromagnetic survey, completed with marine and terrestrial data.The joint analysis of gravity and magnetic anomalies has shown tectonic and volcanic features that define some fundamental aspects of the structural framework and volcanic evolution of the island. A strong gravity anomaly produced by a large and deep source has been associated with an uplifted block of the basement beneath the southern part of Tenerife. The sources of the observed gravity highs from 8 km b.s.l. may be associated with the growth of the submarine shield stage that was clearly controlled by regional tectonic.The long-wavelength magnetic anomalies reveal highly magnetic sources, interpreted as gabbro-ultramafic cumulates associated with the root zone of a large dyke swarm. This intrusive body could be topped by the emplacement zone of magma chambers that correlate with a magnetic horizon at 5.7±0.8 km depth. Rooted in this highly magnetic zone, two dike–like structures can be associated with the magmatic feeding system of large recent basaltic volcanoes. A shallow magnetic horizon (1.4 km a.s.l.) can be correlated with the bottom phonolites of the Las Cañadas Edifice.In the central part of the island the coincidence of some gravity and magnetic lows is consistent with the presence of low-density and low-magnetic materials, that infill a collapsed caldera system. The structures close to the surface are characterised by low-density areas connected with the recent volcanism, in particular the minimum over the Teide volcano. Hydrothermal alteration is assumed to be the cause of a short-wavelength magnetic low over the Teide volcano.     

Physical-chemical conditions of the Teide volcanic system (Tenerife, Canary Islands)

The Teide volcano (3717 m) is the central structure of the island of Tenerife and at present its morphology is that of a stratovolcano which has grown on a large caldera with a collapse 17 km in diameter, which was generated some 0.6 million years ago.The different studies that have been carried out seem to indicate that, in a oversimplified model, there is an intermediate magma chamber with an approximate volume of 30 km³ and located 2–3 km below the actual base of the caldera, i.e., almost at sea level, with a temperature of 430 ± 50°C, and a pressure of 400 ± 100 bar.The summit fumarole emissions are 85°C and are formed mainly of CO₂ with small amounts of sulphur species, H₂, CH₄ and He. The water vapor (68–82%) emitted with the gases comes from the vaporization of a perched aquifer in the upper cone, as shown by the isotopic analyses.     

Zonation of the main volcanic hazards (lava flows and ash fall) in Tenerife, Canary Islands. A proposal for a surveillance network

The island of Tenerife is volcanically complex, and its eruptive history predominantly reflects the processes and products of two different eruptive styles: (1) non-explosive effusions of basaltic lavas from fissure vents mostly aligned along two ridges; and (2) less frequent but explosive salic eruptions from central vents associated with the Las Cañadas volcanic edifice and associated summit caldera. We have taken into account this fundamental distinction to develop a volcanic-hazards zonation (for lava flows and ash fall only) that includes: definition of the principal hazards; identification of the areas that have higher probability of containing emission centres; and numerical modelling of the vulnerable areas to be affected by volcanic hazards. Not only does the volcanic-hazards zonation map provide emergency-management officials with an updated assessment of the volcanic hazards, but it also represents a starting point for the preparation of a volcanic risk map for Tenerife. Finally, the hazards-zonation map also furnishes the basis for the design of a proposed volcano surveillance network.     

Underground Temperature Measurements as a Tool for Volcanic Activity Monitoring in the Island of Tenerife, Canary Islands

The spatial distribution of groundwater temperatures in the volcanic island of Tenerife, Canary Islands, has been inferred through measurements of water temperatures collected in the vast network of wells and subhorizontal tunnels, locally called “galleries,” which constitutes the main water supply of the island. The spatial coverage of the network of galleries allows us to reach from depth almost any geological feature of the island. The complex spatial distribution of temperatures in the interior of Tenerife is the result of the complex geological evolution of the island. Groundwater temperatures are greatly affected by groundwater flow and are considerably warmer in those galleries located in areas where water circulation is reduced due to the low permeability of materials and/or to the low infiltration rate of cooling meteoric water. In this sense, groundwater temperature should be characterized in quiescent conditions (background level), in order to facilitate monitoring changes in heat flow, such as those induced by ascending gases expected with an increase in volcanic activity.     

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.     

Probabilistic seismic hazard assessment for Romania and sensitivity analysis: A case of joint consideration of intermediate-depth (Vrancea) and shallow (crustal) seismicity

The earthquake risk on Romania is one of the highest in Europe, and seismic hazard for almost half of the territory of Romania is determined by the Vrancea seismic region, which is situated beneath the southern Carpathian Arc. The region is characterized by a high rate of occurrence of large earthquakes in a narrow focal volume at depth from 70 to 160 km. Besides the Vrancea area, several zones of shallow seismicity located within and outside the Romanian territory are considered as seismically dangerous. We present the results of probabilistic seismic hazard analysis, which implemented the “logic tree” approach, and which considered both the intermediate-depth and the shallow seismicity. Various available models of seismicity and ground-motion attenuation were used as the alternative variants. Seismic hazard in terms of macroseismic intensities, peak ground acceleration, and response spectra was evaluated for various return periods. Sensitivity study was performed to analyze the impact of variation of input parameters on the hazard results. The uncertainty on hazard estimates may be reduced by better understanding of parameters of the Vrancea source zone and the zones of crustal seismicity. Reduction of uncertainty associated with the ground-motion models is also very important issue for Romania.     

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.     

The last 5000 years of activity at Sete Cidades volcano (São Miguel Island,Azores): Implications for hazard assessment

Sete Cidades is a central volcano with a summit caldera at the western end of São Miguel Island, Azores. Its stratigraphy comprises two main geological groups: the Inferior Group, the units of which date from more than 200 000 years ago through to 36 000 years before present, consisting of thick lava flows and subaerial volcaniclastic deposits that built the base of the central volcano; and the Superior Group which comprises all the activity from the last 36 000 years, including pumice and scoria fallout and PDC deposits with minor lava flows. The volcanostratigraphy is divided into six main formations — Risco, Ajuda, Bretanha, Lombas, Santa Bárbara and Lagoas, each defined by different activity phases in the volcano's evolution.     

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.     

A simplified approach for vulnerability assessment in moderate-to-low seismic hazard regions: application to Grenoble (France)

Due to the moderate seismic risks in France, the building vulnerability assessment methods developed for high seismic risk countries could not easily be used here because of their cost and the low-risk perception among the public and officials. A light vulnerability assessment method is proposed and tested in Grenoble (France), based on classes and scores provided in the GNDT method but simplified in terms of visual screening and number of structural parameters used. Compared to the RiskUE method, the damage obtained by our approach shows that 90% of buildings have residuals smaller than 0.2, i.e. one grade of the EMS98 damage scale. A large scale survey is devised and conducted among the inhabitants of Grenoble in order to collect the main structural parameters. By comparing the results from the survey to the historical urbanization of Grenoble and to expert surveys performed in two urban districts, the information useful for the light method of vulnerability assessment can be rapidly collected by non-experts reducing substantially the estimate cost. The average damage is then computed using the GNDT formula considering the probable intensities which could be observed in Grenoble (VII and VIII). The average damage reaches 0.4 in the oldest part of Grenoble mainly made of masonry buildings and 0.2 in reinforced concrete suburbs where reinforced concrete predominates. The results are a relative vulnerability assessment that provides useful initial information for the urban zones of Grenoble where the vulnerability is higher. This method can be used to classify the seismic vulnerability in wide seismic-prone regions to a fair degree of accuracy and at low cost.     

Seismic hazard assessment for the Itoiz dam site (Western Pyrenees,Spain)

This work summarises the seismic hazard analysis performed for the complete characterisation of strong ground-motion at the site of the Itoiz dam (Western Pyrenees, Spain). The hazard analysis includes the compilation of a composite catalogue from French and Spanish agencies, the definition of an original hybrid seismogenic source model (including zones and major faults) and the selection of ground motion prediction equations (GMPEs). Hazard results are provided as hazard curves and acceleration response spectra on rock for the 1000- and 5000-year return periods, which correspond respectively to the operating basis earthquake (OBE) and safety evaluation earthquake (SEE). The impact of truncating GMPEs at a number of standard deviations (epsilon) has been found not critical here for the return periods targeted. Subsequently, an analysis of the contribution of each source to total hazard and a hazard disaggregation analysis are performed in order to establish the earthquake-source parameters for both the OBE and SEE scenarios consistently with the seismotectonics of the region. The European Strong Motion database is then searched and a selection of records is proposed for each of the scenarios. Our results suggest that seismic hazard in the region is underestimated by the official Spanish seismic hazard map included in the current version of the code (NCSE-02), which is the reference document for the definition of seismic actions for dam projects in the whole Pyrenees.     

Reconstruction and analysis of sub-plinian tephra dispersal during the 1530 A.D. Soufrière (Guadeloupe) eruption: Implications for scenario definition and hazards assessment

The last magmatic eruption of Soufrière of Guadeloupe dated at 1530 A.D. (Soufrière eruption) is characterized by an onset with a partial flank-collapse and emplacement of a debris-avalanche that was followed by a sub-plinian VEI 2–3 explosive short-lived eruption (Phase-1) with a column that reached a height between 9 and 12 km producing about 3.9 × 10⁶ m³ DRE (16.3 × 10⁶ m³ bulk) of juvenile products. The column recurrently collapsed generating scoriaceous pyroclastic flows in radiating valleys up to a distance of 5–6 km with a maximum interpolated bulk deposit volume of 11.7 × 10⁶ m³ (5 × 10⁶ m³ DRE). We have used HAZMAP, a numerical simple first-order model of tephra dispersal [Macedonio, G., Costa, A., Longo, A., 2005. A computer model for volcanic ash fallout and assessment of subsequent hazard. Comput. Geosci. 31, 837–845] to reconstruct to a first approximation the potential dispersal of tephra and associated tephra mass loadings generated by the sub-plinian Phase 1 of the 1530 A.D. eruption. We have tested our model on a deterministic average dry season wind profile that best-fits the available data as well as on a set of randomly selected wind profiles over a 5 year interval that allows the elaboration of probabilistic maps for the exceedance of specific tephra mass load thresholds. Results show that in the hypothesis of a future 1530 A.D. scenario, populated areas to a distance of 3–4 km west–southwest of the vent could be subjected to a static load pressure between 2 and 10 kPa in case of wet tephra, susceptible to cause variable degrees of roof damage. Our results provide volcanological input parameters for scenario and event-tree definition, for assessing volcanic risks and evaluating their impact in case of a future sub-plinian eruption which could affect up to 70 000 people in southern Basse-Terre island and the region. They also provide a framework to aid decision-making concerning land management and development. A sub-plinian eruption is the most likely magmatic scenario in case of a future eruption of this volcano which has shown, since 1992, increasing signs of low-energy seismic, thermal, and acid degassing unrest without significant deformation.     

A desktop image processing and photogrammetric method for rapid volcanic hazard mapping: application to air-photo interpretation of Mount Pelée,Martinique

 A desktop image processing and photogrammetric method was developed for digitizing black-and-white aerial photographs. The technique was applied to airborne optical images of Mt. Pelée, Martinique, a historically active volcano in the tropical Lesser Antilles island arc, to evaluate its utility for rapid geologic mapping and hazard assessment in vegetated areas. The digital approach provides several advantages over traditional air-photo interpretation by allowing for change detection in time-series images, morphologic characterization, development of digital elevation models from stereopairs, and geo-referencing with other digital data sets. A digital mosaic of Mt. Pelée was created from air photos acquired in 1951, which covered the region affected by the 1902 eruption. Severe mismatches occurred along edges of adjacent photographs prior to correction, which precluded quantitative morphologic analysis of the volcanic edifice. Geometric corrections and histogram equalization of digitized air photos allowed creation of a continuous mosaic. Comparison of the mosaic and a map based on differences in gray scale and texture to a volcanostratigraphic map revealed that not only the various deposits produced during the 1902 event were easily differentiated, but that older eruptive products were identified, suggesting that this approach may be used for rapid hazard evaluation of historically active tropical volcanoes. Received: 22 January 1996 / Accepted: 26 July 1996     

A methodology for the evaluation of long-term volcanic risk from pyroclastic flows in Campi Flegrei (Italy)

The volcanological history of Campi Flegrei suggests that the most frequent eruptions are characterized by the emplacement of pyroclastic flow and surge deposits erupted from different vents scattered over a 150-km² caldera. The evaluation of volcanic risk in volcanic fields is complex because of the lack of a central vent. To approach this problem, we subdivided the entire area of Campi Flegrei into a regular grid and evaluated the relative spatial probability of opening of vents based on geological, geophysical and geochemical data. We evaluated the volcanic risk caused by pyroclastic flows based on the formula proposed by UNESCO (1972), R=H×V×Va, where H is the hazard, V is the vulnerability and Va is the value of the elements at risk. The product H×V was obtained by performing simulations of type eruptions centered in each cell of the grid. The simulation is based on the energy cone scheme proposed by Sheridan and Malin [J. Volcanol. Geotherm. Res. 17 (1983) 187–202], hypothesizing a column collapse height of 100 m for eruptions of VEI=3 and 300 m for eruptions of VEI=4 with a slope angle of 6°. Each simulation has been given the relative probability value associated with the corresponding cell. We made use of the GIS software ArcView 3.2 to evaluate the intersection between the energy cone and the topography. The superposition of the areas invaded by pyroclastic flows (124 simulations for VEI=3 and 37 for VEI=4) was used to obtain the relative hazard map of the area. The relative volcanic risk map is obtained by superimposing the urbanization maps.     

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     

The earthquake of 6 September 2002 and the seismic history of Palermo (Northern Sicily,Italy): Implications for the seismic hazard assessment of the city

On September 6, 2002, a M_L = 5.6 earthquake, occurring some tens of kilometres offshore from the Northern Sicilian coast (Southern Tyrrhenian Sea), slightly damaged the city of Palermo and surroundings (degree 6 in the European Macroseismic Scale 1998). The macroseismic investigation of the shock and a detailed study of effects of the main earthquakes which affected Palermo in the past have been performed in order to evaluate the seismic response of the city. Moreover, the comparison of the recent event, which is instrumentally constrained, with historical earthquakes allows us to infer new insights on the seismogenic sources of the area, that seem located offshore in the Tyrrhenian sea.In the last 500 years, Palermo has never been completely destroyed but has suffered effects estimated between intensities 6 and 8 EMS-98 many times (1693, 1726, 1751, 1823, 1940, 1968, 2002). The damage scenarios of the analysed events have shown that damage distribution is strongly conditioned by soil response in the different parts of the city and by a high building vulnerability, mainly in the historical centre and in the south-eastern zone of the modern city. As a matter of fact, Palermo has always suffered greater effects than those reported for other nearby localities. The hazard assessment obtained using observed site intensities has shown that the probability of occurrence for intensity 8 (the strongest intensity observed in Palermo) exceeds 99% for 550 years, while the estimated mean return period is 152 ± 40 years. These results, in connection with building vulnerability due to the urban expansion before the introduction of seismic code, suggest that the city is exposed to a relatively high seismic risk.This paper has not been submitted elsewhere in identical or similar form, nor will it be during the first 3 months after its submission to Journal of Seismology.