High precision locations of LP events on Mt. Etna: Reconstruction of the fluid-filled volume

During 1991–93 at Mount Etna, long-period (LP) events occurring in swarms characterized the evolution of the eruption. The presence of multiplets i.e. groups of events with similar waveform signatures, has been recognized within this activity. Traditional techniques for locating LP events do not allow obtaining reliable hypocenters, which have only succeeded in placing earthquakes in a roughly 1 km² area slightly east of the Mt. Etna Northeast Crater. Hypocenters have been relocated in two steps: the absolute location has been improved using Thurber’s code and a complex 3D velocity model; a highly precise relative location has been applied on multiplets to define the source geometry. 3D locations and high precision analysis suggest that during the 1991–93 eruption the resonator producing LP events was a part of the uppermost Northeast Crater conduit, measuring 210 meters in height and 45–50 meters in diameter.     

Compositionally zoned crystals and real-time degassing data reveal changes in magma transfer dynamics during the 2006 summit eruptive episodes of Mt. Etna

One of the major objectives of volcanology remains relating variations in surface monitoring signals to the magmatic processes at depth that cause these variations. We present a method that enables compositional and temporal information stored in zoning of minerals (olivine in this case) to be linked to observations of real-time degassing data. The integrated record may reveal details of the dynamics of gradual evolution of a plumbing system during eruption. We illustrate our approach using the 2006 summit eruptive episodes of Mt. Etna. We find that the history tracked by olivine crystals, and hence, most likely the magma pathways within the shallow plumbing system of Mt. Etna, differed considerably between the July and October eruptions. The compositional and temporal record preserved in the olivine zoning patterns reveal two mafic recharge events within months of each other (June and September 2006), and each of these magma supplies may have triggered the initiation of different eruptive cycles (July 14–24 and August 31–December 14). Correlation of these observations with gas monitoring data shows that the systematic rise of the CO₂/SO₂ gas values is associated with the gradual (pre- and syn-eruptive) supply of batches of gas-rich mafic magma into segments of Etna’s shallow plumbing system, where mixing with pre-existing and more evolved magma occurred.     

Detection and Analysis of Near-Surface Explosions on the Kola Peninsula

Seismic and infrasonic observations of signals from a sequence of near-surface explosions at a site on the Kola Peninsula have been analyzed. NORSAR’s automatic network processing of these events shows a significant scatter in the location estimates and, to improve the automatic classification of the events, we have performed full waveform cross-correlation on the data set. Although the signals from the different events share many characteristics, the waveforms do not exhibit a ripple-for-ripple correspondence and cross-correlation does not result in the classic delta-function indicative of repeating signals. Using recordings from the ARCES seismic array (250 km W of the events), we find that a correlation detector on a single channel or three-component station would not be able to detect subsequent events from this source without an unacceptable false alarm rate. However, performing the correlation on each channel of the full ARCES array, and stacking the resulting traces, generates a correlation detection statistic with a suppressed background level which is exceeded by many times its standard deviation on only very few occasions. Performing f-k analysis on the individual correlation coefficient traces, and rejecting detections indicating a non-zero slowness vector, results in a detection list with essentially no false alarms. Applying the algorithm to 8 years of continuous ARCES data identified over 350 events which we confidently assign to this sequence. The large event population provides additional confidence in relative travel-time estimates and this, together with the occurrence of many events between 2002 and 2004 when a temporary network was deployed in the region, reduces the variability in location estimates. The best seismic location estimate, incorporating phase information for many hundreds of events, is consistent with backazimuth measurements for infrasound arrivals at several stations at regional distances. At Lycksele, 800 km SW of the events, as well as at ARCES, infrasound is detected for most of the events in the summer and for few in the winter. At Apatity, some 230 km S of the estimated source location, infrasound is detected for most events. As a first step to providing a Ground Truth database for this useful source of infrasound, we provide the times of explosions for over 50 events spanning 1 year.     

Analysis of Signals from an Unique Ground-Truth Infrasound Source Observed at IMS Station IS26 in Southern Germany

Quantitative modeling of infrasound signals and development and verification of the corresponding atmospheric propagation models requires the use of well-calibrated sources. Numerous sources have been detected by the currently installed network of about 40 of the final 60 IMS infrasound stations. Besides non-nuclear explosions such as mining and quarry blasts and atmospheric phenomena like auroras, these sources include meteorites, volcanic eruptions and supersonic aircraft including re-entering spacecraft and rocket launches. All these sources of infrasound have one feature in common, in that their source parameters are not precisely known and the quantitative interpretation of the corresponding signals is therefore somewhat ambiguous. A source considered well-calibrated has been identified producing repeated infrasound signals at the IMS infrasound station IS26 in the Bavarian forest. The source results from propulsion tests of the ARIANE-5 rocket’s main engine at a testing facility near Heilbronn, southern Germany. The test facility is at a range of 320 km and a backazimuth of ~280° from IS26. Ground-truth information was obtained for nearly 100 tests conducted in a 5-year period. Review of the available data for IS26 revealed that at least 28 of these tests show signals above the background noise level. These signals are verified based on the consistency of various signal parameters, e.g., arrival times, durations, and estimates of propagation characteristics (backazimuth, apparent velocity). Signal levels observed are a factor of 2–8 above the noise and reach values of up to 250 mPa for peak amplitudes, and a factor of 2–3 less for RMS measurements. Furthermore, only tests conducted during the months from October to April produce observable signals, indicating a significant change in infrasound propagation conditions between summer and winter months.     

Volcanic Tremor at Mt. Etna,Italy, Preceding and Accompanying the Eruption of July – August, 2001

The July 17 – August 9, 2001 flank eruption of Mt. Etna was preceded and accompanied by remarkable changes in volcanic tremor. Based on the records of stations belonging to the permanent seismic network deployed on the volcano, we analyze amplitude and frequency content of the seismic signal. We find considerable changes in the volcanic tremor which mark the transition to different styles of eruptive activity, e.g., lava fountains, phreatomagmatic activity, Strombolian explosions. In particular, the frequency content of the signal decreases from 5 Hz to 3 Hz at our reference station ETF during episodes of lava fountains, and further decreases at about 2 Hz throughout phases of intense lava emission. The frequency content and the ratios of the signal amplitude allow us to distinguish three seismic sources, i.e., the peripheral dike which fed the eruption, the reservoir which fed the lava fountains, and the central conduit. Based on the analysis of the amplitude decay of the signal, we highlight the migration of the dike from a depth of ca. 5 km to about 1 km between July 10 and 12. After the onset of the effusive phase, the distribution of the amplitude decay at our stations can be interpreted as the overall result of sources located within the first half kilometer from the surface. Although on a qualitative basis, our findings shed some light on the complex feeding system of Mt. Etna, and integrate other volcanological and geophysical studies which tackle the problem of magma replenishment for the July–August, 2001 flank eruption. We conclude that volcanic tremor is fundamental in monitoring Mt. Etna, not only as a marker of the different sources which act within the volcano edifice, but also of the diverse styles of eruptive activity. An erratum to this article is available at .     

Review of Microgravity Observations at Mt. Etna: A Powerful Tool to Monitor and Study Active Volcanoes

Microgravity observations at Mt. Etna have been routinely performed as both discrete (since 1986) and continuous (since 1998) measurements. In addition to describing the methodology for acquiring and reducing gravity data from Mt. Etna, this paper provides a collection of case studies aimed at demonstrating the potential of microgravity to investigate the plumbing system of an active volcano and detect forerunners to paroxysmal volcanic events. For discrete gravity measurements, results from 1994–1996 and 2001 are reported. During the first period, the observed gravity changes are interpreted within the framework of the Strombolian activity which occurred from the summit craters. Gravity changes observed during the first nine months of 2001 are directly related to subsurface mass redistributions which preceded, accompanied and followed the July-August 2001 flank eruption of Mt. Etna. Two continuous gravity records are discussed: a 16-month (October 1998 to February 2000) sequence and a 48-hour (26–28 October, 2002) sequence, both from a station within a few kilometers of the volcano's summit. The 16-month record may be the longest continuous gravity sequence ever acquired at a station very close to the summit zone of an active volcano. By cross analyzing it with contemporaneous discrete observations along a summit profile of stations, both the geometry of a buried source and its time evolution can be investigated. The shorter continuous sequence encompasses the onset of an eruption from a location only 1.5 km from the gravity station. This gravity record is useful for establishing constraints on the characteristics of the intrusive mechanism leading to the eruption. In particular, the observed gravity anomaly indicates that the magma intrusion occurred “passively” within a fracture system opened by external forces.     

Locating the Tohoku-Oki 2011 tsunami source using acoustic–gravity waves

The giant Tohoku-Oki earthquake of 11 March 2011 in offshore Japan did not only generate tsunami waves in the ocean but also infrasound (or acoustic–gravity) waves in the atmosphere. We identified ultra-long-period signals (>500 s) in the recordings of infrasound stations in northeast Asia, the northwest Pacific, and Alaska. Their source was found close to the earthquake epicenter. Therefore, we conclude that in general, infrasound observations after a large offshore earthquake are evidence that the surface and the floor of the sea have been significantly vertically displaced by the earthquake and that a tsunami must be expected. Since infrasound is traveling faster than the tsunami, such information may be used for tsunami early warnings.     

Volcanic tremor at Mt. Etna: Inferences on magma dynamics during effusive and explosive activity

During 1999, the volcanic activity at Mt. Etna was both explosive and effusive at the summit craters: Strombolian activity, lava fountains and lava flows affected different areas of the volcano, involving three of the four summit craters. Results from analysis of the 1999 volcanic tremor features are shown at two different time scales. First, the long-term time variation of the features of the volcanic tremor (including spectral and polarization parameters), during the entire year, was compared with the evolution of the eruptive activity. This approach demonstrated the good agreement between tremor data and observed eruptive activity; the activation of different tremor sources was suggested. Then, a more refined analysis of the volcanic tremor, recorded during 14 lava fountain eruptions, was performed. In particular, a shift of the dominant frequencies towards lower values was noted which corresponds with increasing explosive activity. Similar behaviour in the frequency content has already been observed in other explosive eruptions at Mt. Etna as well as on other volcanoes. This behaviour has been explained in terms of either an increase in the tremor source dimension or a decrease in the sound speed in the magma within the conduit. These results confirm that the volcanic tremor is a powerful tool for better understanding the physical processes controlling explosive eruptions at Mt. Etna volcano.     

The role of syn-eruptive vesiculation on explosive basaltic activity at Mt. Etna,Italy

We investigated the dynamics of explosive activity at Mt. Etna between 31 August and 15 December 2006 by combining vesicle studies in the erupted products with measurements of the gas composition at the active, summit crater. The analysed scoria clasts present large, connected vesicles with complex shapes and smaller, isolated, spherical vesicles, the content of which increases in scoriae from the most explosive events. Gas geochemistry reports CO₂/SO₂ and SO₂/HCl ratios supporting a deep-derived gas phase for fire-fountain activity. By integrating results from scoria vesiculation and gas analysis we find that the highest energy episodes of Mt. Etna activity in 2006 were driven by a previously accumulated CO₂-rich gas phase but we highlight the lesser role of syn-eruptive vesicle nucleation driven by water exsolution during ascent. We conclude that syn-eruptive vesiculation is a common process in Etnean magmas that may promote a deeper conduit magma fragmentation and increase ash formation.     

High precision locations of multiplets on south-eastern flank of Mt. Etna (Italy): reconstruction of fault plane geometry

On 9 January 2001, a seismic swarm, located on the south-eastern flank of Mt. Etna and with nearly identical waveforms, caused some damage to Zafferana Etnea village, 3 km from the epicentral area.An analysis of the seismicity occurring in the last 8 years in this area has revealed other earthquakes with the same characteristics; some pre-empted and followed (up to a few months) the 2001 January swarm, others were recorded more than five years beforehand.Using similarity of waveforms, these earthquakes were classified into three families.The use of a multiplet-technique has allowed to obtain the spatial distribution of the events with higher precision (mean error of 10-20 m) with respect to traditional localization techniques.Mt. Etna earthquakes relocation clearly describes the geometry of the seismogenic tectonic structure; the hypocenters lie on a NE-SW oriented plane that is coincident with one of the focal planes obtained by first-arrival polarities. This alignment is also coherent with one of the main regional tectonic trends cutting the Mt. Etna area, and can be interpreted as a right-lateral strike seismic source on the south-eastern flank of Mt. Etna, distant from eruptive centres, which repeats from time to time and is able to produce strong energy releases.     

New Insights on Mt. Etna’s Crust and Relationship with the Regional Tectonic Framework from Joint Active and Passive P-Wave Seismic Tomography

In the Central Mediterranean region, the production of chemically diverse volcanic products (e.g., those from Mt. Etna and the Aeolian Islands archipelago) testifies to the complexity of the tectonic and geodynamic setting. Despite the large number of studies that have focused on this area, the relationships among volcanism, tectonics, magma ascent, and geodynamic processes remain poorly understood. We present a tomographic inversion of P-wave velocity using active and passive sources. Seismic signals were recorded using both temporary on-land and ocean bottom seismometers and data from a permanent local seismic network consisting of 267 seismic stations. Active seismic signals were generated using air gun shots mounted on the Spanish Oceanographic Vessel ‘Sarmiento de Gamboa’. Passive seismic sources were obtained from 452 local earthquakes recorded over a 4-month period. In total, 184,797 active P-phase and 11,802 passive P-phase first arrivals were inverted to provide three different velocity models. Our results include the first crustal seismic active tomography for the northern Sicily area, including the Peloritan–southern Calabria region and both the Mt. Etna and Aeolian volcanic environments. The tomographic images provide a detailed and complete regional seismotectonic framework and highlight a spatially heterogeneous tectonic regime, which is consistent with and extends the findings of previous models. One of our most significant results was a tomographic map extending to 14 km depth showing a discontinuity striking roughly NW–SE, extending from the Gulf of Patti to the Ionian Sea, south-east of Capo Taormina, corresponding to the Aeolian–Tindari–Letojanni fault system, a regional deformation belt. Moreover, for the first time, we observed a high-velocity anomaly located in the south-eastern sector of the Mt. Etna region, offshore of the Timpe area, which is compatible with the plumbing system of an ancient shield volcano located offshore of Mt. Etna.     

Upwards migration of seismic focii: A forerunner of the 1989 eruption of Mt Etna (Italy)

Data from a portable array of three-component digital stations, run at Mt Etna from 1988 to early 1990, highlight the seismic behaviour of the volcano before the 1989 eruption, one of the most significant in terms of energy of the last two decades. After a two-year period of weak and discotinuous seismicity, the depth of the seismically active volumes was observed to become shallower a few months before the volcanic event. The overall migration of the events, inferred by hypocentral locations and decreases of S-P time differences at two stations, agrees with other geophysical forerunners and allows further insights into the changes in the stress field leading to the eruption.     

Detection of Impact Points of Fragments of Spent Launch Vehicle Stages Using Infrasound Direction-Finding Methods

The paper describes the principles and techniques used to detect signals propagating in the atmosphere in the infrasonic frequency range. Such signals can be generated by different sources: ground and atmospheric explosions, as well as objects moving in the atmosphere at supersonic speed (aircraft, rockets, bolides, fragments of spent stages of launch vehicles). Portable infrasound monitoring stations are described, each of which includes three spaced infrasonic microphones. Each such station makes it possible to determine three basic parameters of the detected infrasound signal, which are subsequently used to solve the direction- finding problem: the time of arrival of an infrasonic wave, the azimuth to the source in the horizontal plane, and the wave approach angle from the source of infrasonic waves to the Earth’s surface in the vertical plane. An acoustic detector used to extract useful signals against a noise background is described. The detector is based on an algorithm similar to the STA/LTA detection algorithm known in seismology. Examples of the operation of an acoustic detector with data obtained during real measurements are given. Passive infrasound direction-finding technology is described. It is based on mathematical modeling of the of infrasonic wave propagation in the atmosphere, which are generated by objects moving along possible trajectories; comparison of theoretical signals with real ones recorded by monitoring stations; and determination of the realized trajectories. The paper gives examples of experimental verification of the effectiveness of passive infrasound direction-finding technology for determining the impact points of the first and second stages of launch vehicles. It is shown that infrasound direction-finding systems makes it possible to reduce the estimated search area for launch vehicle fragments that fall to the Earth, significantly decrease the time and costs for their search and utilization, and mitigate the negative environmental impact of the rocket and space industry.     

High precision tilt observation at Mt. Etna Volcano,Italy

In 2007–2008, we installed on Mt. Etna two deep tilt stations using high resolution, self-leveling instruments. These installations are a result of accurate instrument tests, site selection, drilling and sensor positioning that has allowed detecting variations related to the principal diurnal and semidiurnal tides for first time on Mt. Etna using tilt data.     

Monitoring Seismo-volcanic and Infrasonic Signals at Volcanoes: Mt. Etna Case Study

Volcanoes generate a broad range of seismo-volcanic and infrasonic signals, whose features and variations are often closely related to volcanic activity. The study of these signals is hence very useful in the monitoring and investigation of volcano dynamics. The analysis of seismo-volcanic and infrasonic signals requires specifically developed techniques due to their unique characteristics, which are generally quite distinct compared with tectonic and volcano-tectonic earthquakes. In this work, we describe analysis methods used to detect and locate seismo-volcanic and infrasonic signals at Mt. Etna. Volcanic tremor sources are located using a method based on spatial seismic amplitude distribution, assuming propagation in a homogeneous medium. The tremor source is found by calculating the goodness of the linear regression fit (R ²) of the log-linearized equation of the seismic amplitude decay with distance. The location method for long-period events is based on the joint computation of semblance and R ² values, and the location method of very long-period events is based on the application of radial semblance. Infrasonic events and tremor are located by semblance–brightness- and semblance-based methods, respectively. The techniques described here can also be applied to other volcanoes and do not require particular network geometries (such as arrays) but rather simple sparse networks. Using the source locations of all the considered signals, we were able to reconstruct the shallow plumbing system (above sea level) during 2011.     

Observations of volcanic tremor during the January–February 2005 eruption of Mt. Veniaminof,Alaska

Mt. Veniaminof, Alaska Peninsula, is a stratovolcano with a summit ice-filled caldera containing a small intracaldera cone and active vent. From January 2 to February 21, 2005, Mt. Veniaminof erupted. The eruption was characterized by numerous small ash emissions (VEI 0 to 1) and accompanied by low-frequency earthquake activity and volcanic tremor. We have performed spectral analyses of the seismic signals in order to characterize them and to constrain their source. Continuous tremor has durations of minutes to hours with dominant energy in the band 0.5–4.0 Hz, and spectra characterized by narrow peaks either irregularly (non-harmonic tremor) or regularly spaced (harmonic tremor). The spectra of non-harmonic tremor resemble those of low-frequency events recorded simultaneously with surface ash explosions, suggesting that the source mechanisms might be similar or related. We propose that non-harmonic tremor at Mt. Veniaminof results from the coalescence of gas bubbles while low-frequency events are related to the disruption of large gas pockets within the conduit. Harmonic tremor, characterized by regular and quasi-sinusoidal waveforms, has duration of hours. Spectra containing up to five harmonics suggest the presence of a resonating source volume that vibrates in a longitudinal acoustic mode. An interesting feature of harmonic tremor is that frequency is observed to change over time; spectral lines move towards higher or lower values while the harmonic nature of the spectra is maintained. Factors controlling the variable characteristics of harmonic tremor include changes in acoustic velocity at the source and variations of the effective size of the resonator.     

Entablature: fracture types and mechanisms

Stratigraphic and morphostructural analyses have been carried out in the Mt. Etna volcanic region (eastern Sicily) to investigate in detail the deformation events that have affected the sedimentary successions forming the substratum of the volcano. In the foredeep, Quaternary submarine sedimentation ended 600 ka ago when the whole area emerged in response to homogeneous regional uplift. The irregular distribution of a stratigraphic marker, recognized through the analysis of more than 250 borehole logs, suggests that local dynamics also affected the area. We identify both compressional tectonic dynamics and volcano-related tectonic activity, and discriminate among their associated deformations. In particular, we quantify the vertical deformation component of the compressional structures (thrusts and related folds) and recognize for the first time a vertical component of deformation whose pattern clearly indicates a doming process acting at Mt. Etna. The comparison between long-term and short-term rates suggests that the doming has acted consistently over space and time through the last 600 ka and provides clues to the source of uplift. This component, defined by a specific Quaternary sedimentary horizon, has been compared with vertical deformation obtained by analytical inversion of morphological substratum data, and localizes the source at a depth of ～16 km, at the mantle-crust transition. This uplift may be the consequence of hydration occurring in the altered ocean-like crust.     

Geochemistry of some volcanic rocks from south-eastern sicily: Rare earth and other trace element distribution

Some trace element data for volcanic rocks found at different levels, from Tertiary to Holocene, in south-eastern Sicily (Iblean Plateau and Mt. Etna) are presented and discussed in the present paper in order to better the information about the origin and relationships of the various rock types. Four groups of volcanic rocks have been recognized on the basis of their major element chemistry: 1) low-K tholeiites, 2) associated alkali basalts to nephelinites of the Iblean Plateau (Upper Pliocene to Lower Pleistocene), 3) the basal subalkaline lavas of Mt. Etna, and 4) the alkalic suite rocks that make up the bulk of the volcano. The distribution of Rb, Sr, Ni, Cr, Co, Cu, REE, Th and Sc suggests:

1. an origin of the Iblean magmas by a different degree of partial melting of a Rb-poor and possibly slightly hetereogeneous mantle;
2. quite distinct source compositions for the Etnean magmas, relative to those of the Iblean area, on the basis of their Rb and Sr contents;
3. an origin of the alkalic rocks of Mt. Etna from independently generated magma(s) rather than by crystal fractionation of the Etnean subalkaline magmas or of a magma having the geochemical features of the Iblean alkali basalts; evidence for this is given by the distribution features of the incompatible elements showing an origin for these rocks from compositionally different parent magmas and/or an evolution under widely variable environmental conditions;
4. the primary character for the chemical differences observed in some of the Etnean subalkaline rocks that can be accounted for by different physico-chemical conditions at their source rather than by crystal fractionation processes.

     

Building vulnerability and seismic risk analysis in the urban area of Mt. Etna volcano (Italy)

The tectonic system of the eastern flank of Mt. Etna volcano (Sicily, Italy) is the source of most of the strongest earthquakes occurring in the area over the last 205 years. A total of 12 events with epicentre intensities ≥VIII EMS have occurred at Mt. Etna, 10 of which were located on the eastern flank. This indicates a mean recurrence time of about 20 years. This area is highly urbanised, with many villages around the volcano at altitudes up to 700 m a.s.l. The southern and eastern flanks are particularly highly populated areas, with numerous villages very close to each other. The probabilistic seismic hazard due to local faults for Mt. Etna was calculated by adopting a site approach to seismic hazard assessment. Only the site histories of local volcano-tectonic earthquakes were considered, leaving out the effects due to strong regional earthquakes that occurred in north-eastern and south-eastern Sicily. The inventory used in this application refers to residential buildings. These data were extracted from the 1991 census of the Italian National Institute of Statistics, and are grouped according to the census sections. The seismic vulnerability of the elements at risk belonging to a given building typology is described by a vulnerability index, in accordance with a damage model based on macroseismic intensities. For the estimation of economic losses due to physical damage to buildings, an integrated impact indicator was used, which is equivalent to the lost building volume. The expected annualised economic earthquake losses were evaluated both in absolute and in relative terms, and were compared with the geographical distribution of seismic hazard and with similar evaluations of losses for other regions.