Geodetic techniques applied to the study of the Etna volcano area (Italy)

Volcanic behaviour of Mt. Etna is due to the complex interaction between both the local and the regional stress field involving the eastern Sicily. Eruptions could trigger (be triggered?) during crust extension and/or compression, which are strictly linked with the dynamics of the lower mantle. In this study, very long baseline interferometry (VLBI) space geodesy technique has been used in order to study Etna volcano’s activity by means of the crustal deformations between Noto and Matera (located on the African and the Eurasian Plates, respectively). By analysing VLBI data, we obtained the behaviour of the baseline which crosses the Etnean area, from 1990 December to 2003 March, representing the time variations of the distance between the two geodetic stations; the linear trend of the baseline shows a general increasing, pointing out an extension of the crust between them. A detailed analysis of the Noto-Matera baseline allows the identification of three parts of the VLBI curve in the considered period. In the first part of the curve (from 20/12/90 to 09/02/94), VLBI data are rather poor and therefore no reliable consideration about correlation between crust movements and volcanic and seismicity activity has been made. In the second part of the curve (from 09/02/94 to 04/09/00), VLBI data are more frequent and show slightly fluctuations in the distance. Increasing in the extension and compression were observed in the central and in the final part of this period. In the third period (from 04/09/00 to 25/03/03), VLBI data are very sparse even if the time series was quite long; therefore, to fill gaps in the information, we analysed global positioning system (GPS) data. GPS technique performs continuous observations, and we were able to highlight both extensions and compressions in detail. Comparisons between the trend of Noto-Matera baseline length variations, volcanic activity and seismicity in the Etna area show the complexity of the development over time and space of the phenomenology determined by a deep cause which can be traced, in our opinion, to the interaction between the asthenospheric mantle, deep crust and surface crust. Therefore, we state that crustal distension and compression are determined by the lower pulsating mantle.     

Exploratory seismic site response surveys in a complex geologic area: a case study from Mt. Etna volcano (southern Italy)

A preliminary study targeting to evaluate the local seismic response was performed in the eastern flank of Mt. Etna (southern Italy) using ambient noise measurements. The obtained spectral ratios were subdivided through cluster analysis into different classes of fundamental frequency permitting to draw an iso-frequency contour map. The analysis set into evidence the extreme heterogeneity of lava sequences, which makes difficult to identify a single seismic bedrock formation. Another important outcome, concerning the local seismic effects in terms of frequency and azimuth, is the important role played by the fracture fields associated with the main structural systems of the area. The existence of two zones with strong directional effects striking WNW–ESE and NW–SE, nearly orthogonal to the orientation of the main fracture fields, corroborate such hypothesis.     

Exploratory seismic site response surveys in a complex geologic area: a case study from Mt. Etna volcano (southern Italy)

A preliminary study targeting to evaluate the local seismic response was performed in the eastern flank of Mt. Etna (southern Italy) using ambient noise measurements. The obtained spectral ratios were subdivided through cluster analysis into different classes of fundamental frequency permitting to draw an iso-frequency contour map. The analysis set into evidence the extreme heterogeneity of lava sequences, which makes difficult to identify a single seismic bedrock formation. Another important outcome, concerning the local seismic effects in terms of frequency and azimuth, is the important role played by the fracture fields associated with the main structural systems of the area. The existence of two zones with strong directional effects striking WNW–ESE and NW–SE, nearly orthogonal to the orientation of the main fracture fields, corroborate such hypothesis.

     

Late Pleistocene to Holocene record of changing uplift rates in southern Calabria and northeastern Sicily (southern Italy, Central Mediterranean Sea)

A combination of published and new radiometric dates on uplifted Holocene fossil beaches from northeastern Sicily and southern Calabria (southern Italy) is compared with the altitude of the inner margin of the Last Interglacial (LIg) (Late Pleistocene, 124 ka) and older marine terraces in order to gain a regional-scale outline of uplift rates and their temporal changes in a region which is one of the fastest uplifting sectors of the Central Mediterranean Sea. Late Holocene radiocarbon dates from Ioppolo (southern Calabria) and Ganzirri (northeast Sicily), two newly discovered sites are here presented for the first time. The Holocene uplift rates are highest at St. Alessio and Taormina in eastern Sicily (2.4 mm/y) and at Scilla in southwestern Calabria (2.1 mm/y), two sites located across the Messina Straits and which separate the island of Sicily from mainland Italy. Uplift rates decrease towards the south and north from this centre of uplift. Late Holocene uplift rates show an apparent increase of between 64 and 124% when compared with the longer-term uplift rates calculated from the LIg highstand terraces. Furthermore, we discovered that the locations of fastest Late Pleistocene and Late Holocene uplift rates spatially coincide. To what extent the Holocene increase in uplift rates results from incomplete elastic strain release along the major extensional faults which frame the seismotectonic of the area, or indicate a true change in regional tectonic processes, is not resolved. Nonetheless, the heterogeneity of uplift, with a well-defined centre that crosses the Messina Straits, and its persistence at different time-scales indicates a tight connection between wider regional processes and fault-related displacement in controlling crustal instability in this area.     

Atmospheric sources and sinks of volcanogenic elements in a basaltic volcano (Etna, Italy)

This study reports on the first quantitative assessment of the geochemical cycling of volcanogenic elements, from their atmospheric release to their deposition back to the ground. Etna’s emissions and atmospheric depositions were characterised for more than 2 years, providing data on major and trace element abundance in both volcanic aerosols and bulk depositions. Volcanic aerosols were collected from 2004 to 2007, at the summit vents by conventional filtration techniques. Precipitation was collected, from 2006 to 2007, in five rain gauges, at various altitudes around the summit craters. Analytical results for volcanic aerosols showed that the dominant anions were S, Cl, and F, and that the most abundant metals were K, Ca, Mg, Al, Fe, and Ti (1.5-50 μg m⁻³). Minor and trace element concentrations ranged from about 0.001 to 1 μg m⁻³. From such analysis, we derived an aerosol mass flux ranging from 3000 to 8000 t a⁻¹. Most analysed elements had higher concentrations close to the emission vent, confirming the prevailing volcanic contribution to bulk deposition. Calculated deposition rates were integrated over the whole Etna area, to provide a first estimate of the total deposition fluxes for several major and trace elements. These calculated deposition fluxes ranged from 20 to 80 t a⁻¹ (Al, Fe, Si) to 0.01-0.1 t a⁻¹ (Bi, Cs, Sc, Th, Tl, and U). Comparison between volcanic emissions and atmospheric deposition showed that the amount of trace elements scavenged from the plume in the surrounding of the volcano ranged from 0.1% to 1% for volatile elements such as As, Bi, Cd, Cs, Cu, Tl, and from 1% to 5% for refractory elements such as Al, Ba, Co, Fe, Ti, Th, U, and V. Consequently, more than 90% of volcanogenic trace elements were dispersed further away, and may cause a regional scale impact. Such a large difference between deposition and emission fluxes at Mt. Etna pointed to relatively high stability and long residence time of aerosols in the plume.     

The Upper Pleistocene to Holocene sediments on the Mediterranean island of Lampedusa (Italy)

The island of Lampedusa lies on the northern edge of the African continental shelf, but during some Quaternary marine lowstands it was joined to the African continent. The study and dating of the aeolian, alluvial, detrital sediments, calcareous crusts and speleothems have established that the climatic–environmental variations recorded on the island can be related chronologically to those known for northern Libya, Tunisia and the Italian peninsula. During the Last Glacial Maximum, phases of Saharan dust accumulation on Lampedusa occurred, and were coeval with dust accumulation in crater lakes and on high mountains in central‐southern Italy, and with phases of glacial advance in the Apennines and in the Alps. During the late Holocene, accumulation of Saharan dust on Lampedusa occurred but there was little accumulation of dust on the northern side of the Mediterranean Sea. With the new data from Lampedusa, it is possible to envisage two different scenarios of atmospheric circulation relating to the Last Glacial Maximum and to the late Holocene. During the Last Glacial Maximum, southerly atmospheric circulation brought rainfall to the southern slopes of the Alps and to the Apennines. During the late Holocene, a prevalent westerly atmospheric circulation became established in the northern Mediterranean. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling the long‐term deformation of the sedimentary substrate of Mt. Etna volcano (Italy)

This study investigates in detail the deformation events that have affected the sedimentary successions forming the substrate of Mt. Etna volcano (Italy). Based on the geometric reconstruction of a buried sedimentary marker, we have been able to identify and quantify the effects of three different mechanisms of deformation that have affected the area in the last 600 ka. Numerical results from Finite Element Method (FEM) applied to model viscoelastic deformation suggest the occurrence of a crustal doming process originating at the mantle‐crust transition (~16 km). We propose that the source of deformation is related to the diapiric uprise of hydrothermal material originating in altered ocean‐like crust and its emplacement at a shallower level in the crust. This process has great relevance in the volcanic system and should be considered for the full assessment of its origin and evolution.     

Repeating volcano-tectonic earthquakes at Mt. Etna volcano (Sicily,Italy) during 1999–2009

Repeating volcano-tectonic (VT) earthquakes, taking place at Mt. Etna during 1999–2009, were detected and analyzed to investigate their behavior. We found 735 families amounting to 2479 VT earthquakes, representing ~ 38% of all the analyzed VT earthquakes. The number of VT earthquakes making up the families ranges from 2 to 23. Over 70% of the families comprise 2 or 3 VT earthquakes and only 20 families by more than 10 events. The occurrence lifetime is also highly variable ranging from some minutes to ten years. In particular, more than half of the families have a lifetime shorter than 0.5 day and only ~ 10% longer than 1 year. On the basis of these results, most of the detected families were considered “burst-type”, i.e., show swarm-like occurrence, and hence their origin cannot be explained by a temporally constant tectonic loading. Indeed, since the analyzed earthquakes take place in a volcanic area, the rocks are affected not only by tectonic stresses related to the fairly steady regional stress field but also by local stresses, caused by the volcano, such as magma batch intrusions/movements and gravitational loading. We focused on the five groups of families characterized by the longest repeatability over time, namely high number of events and long lifetime, located in the north-eastern, eastern and southern flanks of the volcano. Unlike the first four groups, which similarly to most of the detected families show swarm-like VT occurrences, group “v”, located in the north-eastern sector, exhibits a more “tectonic” behavior with the events making up such a group spread over almost the entire analyzed period. It is clear how both occurrence and slip rates do not remain constant but vary over time, and such changes are time-related to the occurrence of the 2002–2003 eruption. Finally, by FPFIT algorithm a good agreement between directions identified by nodal planes and the earthquake epicentral distribution was generally found.     

Seismological constraints on the 2018 Mt. Etna (Italy) flank eruption and implications for the flank dynamics of the volcano

3D earthquake locations, focal mechanisms and stress tensor distribution in a 16‐month interval covering the 2018 Mt. Etna flank eruption, enabled us to investigate the relationship between magma intrusion and structural response of the volcano and shed light on the dynamic processes affecting the instability of Mt. Etna. The magma intrusion likely caused tension in the flanks of the volcano, leading to significant ground deformation and redistribution of stress on the neighbouring faults at the edge of Mt. Etna's unstable sector, encouraging the ESE sliding of the eastern flank of the volcano. Accordingly, FPSs of the post‐eruptive events show strike slip faulting mechanisms, under a stress regime characterized by a maximum compressive σ₁, NE‐SW oriented. In this perspective, any flank eruption could temporarily enhance the sliding process of both the southern and eastern flanks of the volcano.     

Natural and anthropogenic factors affecting groundwater quality of an active volcano (Mt. Etna,Italy)

New geochemical data on dissolved major and minor constituents in 276 groundwater samples from Etna aquifers reveal the main processes responsible for their geochemical evolution and mineralisation. This topic is of particular interest in the light of the progressive depletion of water resources and groundwater quality in the area. Multivariate statistical analysis reveal 3 sources of solutes: (a) the leaching of the host basalt, driven by the dissolution of magma-derived CO₂; (b) mixing processes with saline brines rising from the sedimentary basement below Etna; (c) contamination from agricultural and urban wastewaters. The last process, highlighted by increased concentrations of SO₄, NO₃, Ca, F and PO₄, is more pronounced on the lower slopes of the volcanic edifice, associated with areas of high population and intensive agriculture. However, this study demonstrates that natural processes (a) and (b) are also very effective in producing highly mineralised waters, which in turn results in many constituents (B, V, Mg) exceeding maximum admissible concentrations for drinking water.     

Rapid late Pleistocene/Holocene uplift and coastal evolution of the southern Arabian (Persian) Gulf

The coastline along the southern Arabian Gulf between Al Jubail, Kingdom of Saudi Arabia, and Dubai, UAE, appears to have risen at least 125 m in the last 18,000 years. Dating and topographic surveying of paleo-dunes (43–53 ka), paleo-marine terraces (17–30 ka), and paleo-marine shorelines (3.3–5.5 ka) document a rapid, > 1 mm/a subsidence, followed by a 6 mm/a uplift that is decreasing with time. The mechanism causing this movement remains elusive but may be related to the translation of the coastal area through the backbasin to forebulge hinge line movement of the Arabian plate or, alternatively, by movement of the underlying Infracambrian-age Hormuz salt in response to sea-level changes associated with continental glaciation. Independent of the mechanism, rapid and episodic uplift may impact the design of engineering projects such as nuclear power plants, airports, and artificial islands as well as the interpretation of sedimentation and archeology of the area.     

Late-Quaternary paleoenvironmental evolution of Lesina lagoon (southern Italy) from subsurface data

Integrated sedimentological and micropaleontological (foraminifers and ostracods) analyses of two 55 m long borehole cores (S3 and S4) drilled in the subsurface of Lesina lagoon (Gargano promontory—Italy) has yielded a facies distribution characteristic of alluvial, coastal and shallow-marine sediments. Stratigraphic correlation between the two cores, based on strong similarity in facies distribution and AMS radiocarbon dates, indicates a Late Pleistocene to Holocene age of the sedimentary succession.Two main depositional sequences were deposited during the last 60-ky. These sequences display poor preservation of lowstand deposits and record two major transgressive pulses and subsequent sea-level highstands. The older sequence, unconformably overlying a pedogenized alluvial unit, consists of paralic and marine units (dated by AMS radiocarbon at about 45–50,000 years BP) that represent the landward migration of a barrier-lagoon system. These units are separated by a ravinement surface (RS1). Above these tansgressive deposits, highstand deposition is characterised by progradation of the coastal sediments.The younger sequence, overlying an unconformity of tectonic origin, is a 10 m-thick sedimentary body, consisting of fluvial channel sediments overlain by transgressive–regressive deposits of Holocene age. A ravinement surface (RS2), truncating the transgressive (lagoonal and back-barrier) deposits in core S4, indicates shoreface retreat and landward migration of the barrier/lagoon system. The overlying beach, lagoon and alluvial deposits are the result of mid-Holocene highstand sedimentation and coastal progradation.     

VS crustal models of the Roccamonfina volcano and relationships with Neapolitan volcanoes (southern Italy)

Shear wave velocities of the crust and upper mantle are defined beneath the Roccamonfina volcano and surrounding Apennines (southern Italy) from the simultaneous nonlinear inversion of the local group velocity dispersion data, obtained from seismic events recorded in 1988–2004 at Roccamonfina station of the INGV-RSNC network, and regional dispersion data obtained in previous studies. The main features of the representative V_S models are a carbonatic basement and a low velocity zone at 6–10?km of depth. The sedimentary succession is ~5?km thick below the Roccamonfina volcano and lays above a high V_S (3.8?km/s) ascribable to solidified magma body, while it is ~10?km thick below the surrounding Apennines. A low velocity layer with an average thickness of 10?km is detected below the Roccamonfina volcano which can be associated with the presence of partial melting and interpreted as magmatic reservoir. Such low velocity layer, also found below the surrounding Apennines but with a reduced thickness of 2–3?km, extends to the Campanian Plain and to the Neapolitan volcanic area, from Campi Flegrei to Somma-Vesuvius.     

Carbonatites in a subduction system: The Pleistocene alvikites from Mt. Vulture (southern Italy)

We report here, for the first time, on the new finding of extrusive calciocarbonatite (alvikite) rocks from the Pleistocene Mt. Vulture volcano (southern Italy). These volcanic rocks, which represent an outstanding occurrence in the wider scenario of the Italian potassic magmatism, form lavas, pyroclastic deposits, and feeder dikes exposed on the northern slope of the volcano. The petrography, mineralogy and whole-rock chemistry attest the genuine carbonatitic nature of these rocks, that are characterized by high to very high contents of Sr, Ba, U, LREE, Nb, P, F, Th, high Nb/Ta and LREE/HREE ratios, and low contents of Ti, Zr, K, Rb, Na and Cs. The O–C isotope compositions are close to the “primary igneous carbonatite” field and, thus, are compatible with an ultimate mantle origin for these rocks. The Sr–Nd–Pb–B isotope compositions, measured both in the alvikites and in the silicate volcanic rocks, indicate a close genetic relationship between the alvikites and the associated melilitite/nephelinite rocks. Furthermore, these latter products are geochemically distinct from the main foiditic-phonolitic association of Mt. Vulture. We propose a petrogenetic/geodynamic interpretation which has important implications for understanding the relationships between carbonatites and orogenic activity. In particular, we propose that the studied alvikites are generated through liquid unmixing at crustal levels, starting from nephelinitic or melilititic parent liquids. These latter were produced in a hybrid mantle resulting from the interaction through a vertical slab window, between a metasomatized mantle wedge, moving eastward from the Tyrrhenian/Campanian region, and the local Adriatic mantle. The occurrence of carbonatite rocks at Mt. Vulture, that lies on the leading edge of the Southern Apennines accretionary prism, is taken as an evidence for the carbonatation of the mantle sources of this volcano. We speculate that mantle carbonatation is related to the introduction of sedimentary carbon from the Adriatic lithosphere during Tertiary subduction.     

The 18.7 ka phreatomagmatic flank eruption on Etna (Italy): relationship between eruptive activity and sedimentary basement setting

This paper documents a phreatomagmatic flank eruption that occurred 18 700 ± 100 a BP , on the lower north-eastern slope of Etna during the Ellittico volcano activity, which produced fall and surge deposits. This type of eruption is connected to a sedimentary basement ridge at Etna. The interaction between the rising magma and the shallow groundwater hosted in the volcanic pile overlying the impermeable sediments resulted in phreatomagmatic instead of strombolian activity. Three eruptive phases are distinguished based on field and analytical data: (i) an explosive phreatomagmatic opening, (ii) a main phase producing coarse lithic-rich fallout and a strombolian deposit, and (iii) the final pulsating surge-forming phase. The discovery of this phreatomagmatic flank eruption, which occurred at lower altitude, raises important issues for previous hazard assessments at Etna.     

Focal parameters of seismic sources during the 1981 and 1983 eruption at Mt. Etna volcano (Sicily,Italy)

In this study 50 seismic events, preceding and accompanying the eruptions occurring in 1981 and 1983, have been considered. Seismic moments, fault radii, stress drops and seismic energies have been calculated using Brune’s model (J Geophys Res 75:4997–5009, 1970; J Geophys Res 76:5002, 1971); site, anelastic attenuation along the propagation path, geometrical spreading and interaction with the free surface effects are taken into account. For each event we have also estimated the equivalent Wood–Anderson magnitude (MWAeq) (Scherbaum and Stoll in Bull Seism Soc Am 73:1321–1343, 1983); relations among all these source parameters have been determined. Furthermore, the hypothesis of self-similarity (Aki in J Geophys Res 72:1217–1231, 1967) is not verified for events with seismic moments <10¹² N-m: in fact the relationship between log-stress drop and log-moment is linear up to a moment of 10¹² N-m (events of 1981 eruption), while for higher moments (events of 1983 eruption) the slope of the regression line is not significantly different from zero. We suppose that such a behaviour is related to a heterogeneous medium with barriers on the faults. Finally, the main conclusion is that eruptions of 1981 and 1983 differ from one another both in eruptive and seismic aspects; analysis of seismic energies indicates an increase in Mt. Etna’s activity, confirmed by studies performed on the following lateral eruption of 1991–1993 (Patanè et al. in Bull Volcanol 47:941–952, 1995), occurring on the same structural trend.     

Hydrogeochemistry of Roccamonfina volcano (Southern Italy)

This is the first hydro-geochemical investigation carried out on the Roccamonfina Volcanic Complex groundwaters. The chemistry of Roccamonfina waters is defined by water–rock and water–rock–gas interactions. In fact, interactions between rocks of the first eruptive high-K formations and circulating groundwaters are recognized by high K concentrations. On the other hand, inverse concentration of calcium versus alkali metals is related to two different rock interactions occurring in different areas of the volcano: (a) within the caldera where groundwaters flow within latite and pyroclastic formations releasing calcium, and (b) similarly at the base of the volcano where groundwaters flowing from surrounding carbonates got strongly enriched in Ca. These geochemical processes are also associated with K (SE of caldera) and Mg/Ca (in sites located at the NE base of the volcano) decrease. Completely different dynamics occurs at Riardo groundwaters (SE). Here waters are the result of a mix between the Roccamonfina deep aquifer and the carbonate aquifer of the Riardo plain. Rich-CO₂ emissions make these waters strongly mineralized. Minor elements show a similar geochemical behavior of major ions and are crucial defining interactions processes. The evolution of Roccamonfina groundwaters is also evident along the simultaneous enrichment of Ba, Sr, and Ca. Ba increase is the result of deep local carbonate alteration enhanced by CO₂ emissions and, the lower Sr/Ca ratio, from 10 to 2 (ppb/ppm), is also due to the same process. In the light of our results the Roccamonfina aquifer can be schematically divided into two main reservoirs: (a) a superficial aquifer which basically follows the volcanic structure morphology and tectonics and (b) a deeper reservoir, originating within the oldest Roccamonfina volcano ultra potassic lavas and then flowing into the carbonate aquifers of the neighboring plain. Eventually, the chemistry of the Roccamonfina aquifer does not show any specific and visible pollution, contrary to what happens in the volcano surrounding plains. In fact, only 14% of the samples we collected (206) show a NO₃ content >30 mg/l. These sites are all located at the base of the volcano, near the plain.     

New stratigraphical data on the Middle-Late Jurassic biosiliceous sediments from the Sicanian basin, Western Sicily (Italy)

The reported data present the stratigraphy of several sections across a Middle-Late Jurassic Radiolaritic Unit, well exposed in different thrust sheets pertaining to the Maghrebian chain of Southwestern Sicily. The aim was to define the chronostratigraphical distribution of the Jurassic biosiliceous sedimentation in the Sicanian palaeogeographical zone, a deep water basin belonging to the Southern Tethys continental margin. The radiolarian biostratigraphy indicates that the switching from carbonate to siliceous sedimentation in the Sicanian Basin is referable to the Bajocian, as shown by the section of Campofiorito, near Corleone. The biostratigraphical dataset allows the correlation between the onset of biosiliceous sedimentation and the fall of biodiversity in the Sicanian basin with the carbonate productivity crisis, indicated by the highest eutrophication that affected Western Tethys during Middle Jurassic times. Editorial handling: J.-P. Billon-Bruyat & M. Chiari (Guest)