Geology Reports⌕ Search

Geology topics

Research about Italy

Source-linked reports with geographic coverage including Italy.

At least 19 recordsLinked to original sources

Simulating past and future refugia for temperate trees in northern Italy

During the Quaternary, trees responded to the climatic changes of glacial–interglacial cycles with large-scale range shifts. Over cold glacials, temperate tree species contracted their ranges and survived in areas known as refugia. Several studies point to the Euganean Hills (Colli Euganei), in Veneto, northern Italy, as one of the northernmost European refugia of temperate tree species during the Last Glacial Maximum (LGM, ca 23 000–19 000 calibrated years BP). Using LandClim, a spatially explicit, dynamic forest landscape model, we demonstrate that climate conditions during the LGM likely allowed temperate tree species to persist in the Euganean Hills. The identified refugial locations lie at intermediate to high elevations and in sheltered valleys within the hilly complex. Therefore, the combined palaeoecological and modelling evidence suggests that today's temperate forests of the Euganean Hills have a full glacial legacy. Simulations under future climate conditions suggest a collapse of the sub-mediterranean and oro-mediterranean deciduous forests that are prevalent today and the expansion of thermo-mediterranean evergreen forests (with e.g. Quercus ilex , Q. suber , Olea europaea and Pinus sp.). Specifically, the extrazonal population of oro-mediterranean Fagus sylvatica , which is unique to the Po Plain and likely persisted locally through several glacial–interglacial cycles, is predicted to sharply decline and face local extinction, underscoring a conservation hazard.

Ecography↗

Temporal and spatial changes in seismic attenuation associated with inferred fluid migration in the 2016 central Apennines earthquake sequence

Prior work suggests that high‐frequency seismic attenuation acts as a highly sensitive proxy for crustal permeability and fluid mobility in fractured media. We test the hypothesis that the fault system responsible for the 2016–2017 Amatrice–Visso–Norcia–Capitignano sequence acted as an impermeable seal, compartmentalizing pressurized fluids until dynamic rupture triggered widespread fluid diffusion. By tracking across the sequence the spatiotemporal evolution of the S ‐wave anelastic attenuation parameter, we identify large, positive low‐frequency attenuation anomalies emerging within the hanging wall following the Amatrice mainshock and strictly preceding subsequent large ruptures. Conversely, we observe weaker, negative anomalies in the footwall, anticorrelated in time with those of the hanging wall, revealing a massive asymmetry in fluid redistribution and permeability evolution across the fault system. Furthermore, aftershock migration rates reveal distinct linear alignments in a distance‐reduced time space, allowing us to explicitly track and quantify episodes of lateral and upward fluid migration. These physically consistent patterns suggest that stress‐driven fluid diffusion directly weakens adjacent fault patches, dictating the spatiotemporal migration of seismicity. We conclude that near‐real‐time monitoring of seismic attenuation may help detect fluid redistribution in active fault systems and may provide useful information for time‐dependent seismic hazard assessment.

central Apennines↗

The morphological response to the 13 May 2022 major explosive event at Stromboli volcano, Italy, characterised by high-resolution UAS surveys

Since 2019, the frequency of major explosive eruptions at Stromboli volcano (Italy) has increased, heightening the exposure of population and scientists to the hazards posed by ejecta. Morphological changes can directly alter the hazard potential associated with these phenomena. Here, we present a quantitative morphological analysis of changes of the crater terrace area linked to the 13 May 2022 major explosive event. High resolution (2.5 cm pixel -1 ) aerial imagery was acquired by unoccupied aircraft systems 2 hours before and 19 hours after the event. The 13 May 2022 major explosive event consisted of a minimum of seven explosions from four vents located in the south-central crater area. The opportune timing of this campaign enabled the quantification of morphological changes at Stromboli related to a single major explosive event at high temporal and spatial resolution. A total of 12.7 × 10 3 m 3 was excavated and 5.5 × 10 3 m 3 deposited. Via the mapping and classification of bomb distributions we observe that angular blocks make up the largest fraction of ballistics >0.2 m, from which we infer a strong interaction with wall rock and/or fragmentation of solidified plugs in the shallow plumbing system. The morphological changes observed provide valuable constraints on how much material is displaced, and the shift in location and the number of active vents during major explosive events at Stromboli.

Stromboli volcano↗

Characterization and validation of tidally calibrated strains from the Alto Tiberina Near Fault Observatory Strainmeter Array (TABOO-NFO-STAR)

Six horizontal borehole tensor strainmeters (TSM1-6) installed from Fall 2021 to Spring 2022 comprise the Alto Tiberina Near Fault Observatory Strainmeter Array (STAR), providing an unprecedented opportunity to investigate seismic and aseismic deformation from hazardous high- and low-angle normal faults in Italy. Prior to use in tectonic applications, they require in-situ calibration and correction for non-tectonic signals. We tidally calibrate the instruments, characterize the calibration uncertainty, and test the results against environmental and earthquake signals originating from local to teleseismic distances. The STAR sites demonstrably deviate from assumptions common to the standard manufacturer's calibrations, including negative areal coupling at TSM3-6. While the tidally calibrated strains have ~3-56% uncertainty, the calibrated dynamic strains show interstation precision and accuracy to nanostrain levels, and static coseismic offsets in the array footprint are within uncertainty. TSM3 records a complex series of strains that may arise from dynamically triggered near-borehole fracture slip and fluid flow that does not appear to affect its sensitivity to lower strain rate deformation. Future calibration improvement may be afforded with longer stable timeseries, particularly for TSM4. Overall, our analyses demonstrate expanded geodetic capability for detecting deformation in the Alto Tiberina Near Fault Observatory.

Alto Tiberina Near Fault Observatory↗

Seismic response comparison of a historical masonry church subject to real and simulated ground motions

In recent years, advanced numerical models and high-performance computing have facilitated the utilization of ground motion time series in the assessment of the non-linear dynamic behavior of historic masonry structures. Since recorded accelerograms can be sparse for specific analysis conditions, stochastic ground motion simulations have become a viable alternative to overcome this limitation. This study simulates the recorded acceleration time series of the Central Italy 2016 earthquake event at the closest station to the town of Macerata using a site-based stochastic approach. The simulated motions are seismologically evaluated using a goodness-of-fit method in terms of various intensity measures. The simulated records, in conjunction with real records, are used to study the non-linear dynamic behavior of San Filippo Neri church located in Macerata. The church of San Filippo represents an important example of Baroque religious architecture in central Italy, which was damaged and closed off to the public after the 2016 earthquake events. The construction was investigated with a vast diagnostic campaign which included on-site testing and dynamic identification tests. The collected data is used to calibrate the dynamic response of a three-dimensional finite element model of the church. The model is finally used to compare the non-linear seismic responses under real and simulated ground motions with the site recorded damage. The results of structural responses demonstrate a strong agreement between the real and simulated records, providing evidence to support the validation of the site-based stochastic simulation.

Macerata↗

A benchmark dataset and workflow for landslide susceptibility zonation

Landslide susceptibility shows the spatial likelihood of landslide occurrence in a specific geographical area and is a relevant tool for mitigating the impact of landslides worldwide. As such, it is the subject of countless scientific studies. Many methods exist for generating a susceptibility map, mostly falling under the definition of statistical or machine learning. These models try to solve a classification problem: given a collection of spatial variables, and their combination associated with landslide presence or absence, a model should be trained, tested to reproduce the target outcome, and eventually applied to unseen data. Contrary to many fields of science that use machine learning for specific tasks, no reference data exist to assess the performance of a given method for landslide susceptibility. Here, we propose a benchmark dataset consisting of 7360 slope units encompassing an area of about 4,100 km 2 "> 4,100 km 2 in Central Italy. Using the dataset, we tried to answer two open questions in landslide research: (1) what effect does the human variability have in creating susceptibility models; (2) how can we develop a reproducible workflow for allowing meaningful model comparisons within the landslide susceptibility research community. With these questions in mind, we released a preliminary version of the dataset, along with a “call for collaboration,” aimed at collecting different calculations using the proposed data, and leaving the freedom of implementation to the respondents. Contributions were different in many respects, including classification methods, use of predictors, implementation of training/validation, and performance assessment. That feedback suggested refining the initial dataset, and constraining the implementation workflow. This resulted in a final benchmark dataset and landslide susceptibility maps obtained with many classification methods. Values of area under the receiver operating characteristic curve obtained with the final benchmark dataset were rather similar, as an effect of constraints on training, cross–validation, and use of data. Brier score results show larger variability, instead, ascribed to different model predictive abilities. Correlation plots show similarities between results of different methods applied by the same group, ascribed to a residual implementation dependence. We stress that the experiment did not intend to select the “best” method but only to establish a first benchmark dataset and workflow, that may be useful as a standard reference for calculations by other scholars. The experiment, to our knowledge, is the first of its kind for landslide susceptibility modeling. The data and workflow presented here comparatively assess the performance of independent methods for landslide susceptibility and we suggest the benchmark approach as a best practice for quantitative research in geosciences.

Earth-Science Reviews↗

Ore mineralization in the Mofete and San Vito geothermal fields, Campi Flegrei volcanic complex, Naples, Italy

The Mofete and San Vito geothermal fields, located west of Naples, Italy, are part of the Campi Flegrei volcanic complex. In the 1970s, exploratory wells were drilled to a depth of ~3000 m in an attempt to locate high-enthalpy fluids for potential power production. Drill core samples from Mofete wells (MF1, MF2, and MF5) and from San Vito wells (SV1 and SV3) contain authigenic ore mineralization. Pyrite, pyrrhotite, and galena are abundant. Less common are chalcopyrite, sphalerite, arsenopyrite, and scheelite; rare are millerite, violarite, native bismuth, tellurobismuthite, cassiterite, molybdenite, and acanthite. Mineral chemistry was determined by electron microprobe wavelength dispersive spectroscopy aided by a scanning electron microscope equipped with energy-dispersive spectroscopy. The mineral assemblage suggests a low sulfidation environment and the absence of pyrrhotite in the MF1 well and upper part of the SV1 well indicates variable sulfur activity. Both molybdenite and scheelite were identified in samples SV1–2860 and SV3–2353 and scheelite in the SV3 well is zoned with variable Mo 6+ content; low Mo 6+ zones show blue cathodoluminescence, whereas, zones with high Mo 6+ content are yellow to brown. Zoned scheelite and the occurrence of both Mo-bearing minerals attest to the variability of ƒO 2 and ƒS 2 in the geothermal fluid.

Naples↗

Pathology of tissue loss in three key gorgonian species in the Mediterranean Sea

The Mediterranean is known for its marine biodiversity, especially gorgonian forests. Unfortunately, these are experiencing rapid declines due to climate change, manifested by repeated marine heat waves resulting in mass mortality events since the early 1990 s. To better understand why gorgonians are declining, more systematic approaches to investigate the exact causes are needed, and pathology may aid in this goal. We described gross and microscopic pathology of tissue loss in three key gorgonian species in the Mediterranean region, Paramuricea clavata , Eunicella cavolini , and Leptogorgia sarmentosa , that were all experiencing various degrees of acute to subacute tissue loss characterized by exposed axial skeleton sometimes partly colonized by epibionts and thinning of adjacent tissues. The most significant variety of lesions was seen in P. clavata followed by L. sarmentosa and E. cavolini. For all species, dissociation of gastrodermal cells was the dominant microscopic lesion followed by necrosis of the gastrodermis. Ciliates invading gastrodermis and associated with necrosis of polyps were seen only in E. cavolini. Epidermal tissue loss was seen only in L. sarmentosa , while P. clavata was distinguished by a prominent inflammatory response and unidentified dark round structures within the tentacle epidermis and gastrodermis with no host response. Further work to understand the cause of death in gorgonians is needed, particularly to elucidate the role of ciliates and environmental co-factors or infectious agents not visible on light microscopy, as well as applications of additional tools such as cytology.

Journal of Invertebrate Pathology↗

Detection and quantification of preferential flow using artificial rainfall with multiple experimental approaches

Preferential flow in the unsaturated zone strongly influences important hydrologic processes, such as infiltration, contaminant transport, and aquifer recharge. Because it entails various combinations of physical processes arising from the interactions of water, air, and solid particles in a porous medium, preferential flow is highly complex. Major research is needed to improve the ability to understand, quantify, model, and predict preferential flow. Toward a solution, a combination of diverse experimental measurements at multiple scales, from laboratory scale to mesoscale, has been implemented to detect and quantify preferential paths in carbonate and karstic unsaturated zones. This involves integration of information from (1) core samples, by means of mercury intrusion porosimeter, evaporation, quasi-steady centrifuge and dewpoint potentiometer laboratory methods, to investigate the effect of pore-size distribution on hydraulic characteristics and the potential activation of preferential flow, (2) field plot experiments with artificial sprinkling, to visualize preferential pathways related to secondary porosity, through use of geophysical measurements, and (3) mesoscale evaluation of field data through episodic master recession modeling of episodic recharge. This study demonstrates that preferential flow processes operate from core scale to two different field scales and impact on the qualitative and quantitative groundwater status, by entailing fast flow with subsequent effects on recharge rate and contaminant mobilizing. The presented results represent a rare example of preferential flow detection and numerical modeling by reducing underestimation of the recharge and contamination risks.

Bari↗

UAS-Based Observations of Infrasound Directionality at Stromboli Volcano, Italy

Infrasound (low frequency sound waves) can be used to monitor and characterize volcanic eruptions. However, infrasound sensors are usually placed on the ground, thus providing a limited sampling of the acoustic radiation pattern that can bias source size estimates. We present observations of explosive eruptions from a novel uncrewed aircraft system (UAS)-based infrasound sensor platform that was strategically hovered near the active vents of Stromboli volcano, Italy. We captured eruption infrasound from short-duration explosions and jetting events. While potential vertical directionality was inconclusive for the short-duration explosion, we find that jetting events exhibit vertical sound directionality that was observed with a UAS close to vertical. This directionality would not have been observed using only traditional deployments of ground-based infrasound sensors, but is consistent with jet noise theory. This proof-of-concept study provides unique information that can improve our ability to characterize and quantify the directionality of volcanic eruptions and their associated hazards.

Stromboli Volcano↗

On the use of high-resolution and deep-learning seismic catalogs for short-term earthquake forecasts: Potential benefits and current limitations

Enhanced earthquake catalogs provide detailed images of evolving seismic sequences. Currently, these data sets take some time to be released but will soon become available in real time. Here, we explore whether and how enhanced seismic catalogs feeding into established short-term earthquake forecasting protocols may result in higher predictive skill. We consider three enhanced catalogs for the 2016–2017 Central Italy sequence, featuring a bulk completeness lower by at least two magnitude units compared to the real-time catalog and an improved hypocentral resolution. We use them to inform a set of physical Coulomb Rate-and-State (CRS) and statistical Epidemic-Type Aftershock Sequence (ETAS) models to forecast the space-time occurrence of M3+ events during the first 6 months of the sequence. We track model performance using standard likelihood-based metrics and compare their skill against the best-performing CRS and ETAS models among those developed with the real-time catalog. We find that while the incorporation of the triggering contributions from new small magnitude detections of the enhanced catalogs is beneficial for both types of forecasts, these models do not significantly outperform their respective near real-time benchmarks. To explore the reasons behind this result, we perform targeted sensitivity tests that show how (a) the typical spatial discretizations of forecast experiments ( 2 km) hamper the ability of models to capture highly localized secondary triggering patterns and (b) differences in earthquake parameters (i.e., magnitude and hypocenters) reported in different catalogs can affect forecast evaluation. These findings will contribute toward improving forecast model design and evaluation strategies for next-generation seismic catalogs.

Journal of Geophysical Research--Solid Earth↗

8,000 years of climate, vegetation, fire and land-use dynamics in the thermo-mediterranean vegetation belt of northern Sardinia (Italy)

Knowledge about the vegetation history of Sardinia, the second largest island of the Mediterranean, is scanty. Here, we present a new sedimentary record covering the past ~ 8,000 years from Lago di Baratz, north-west Sardinia. Vegetation and fire history are reconstructed by pollen, spores, macrofossils and charcoal analyses and environmental dynamics by high-resolution element geochemistry together with pigment analyses. During the period 8,100–7,500 cal BP , when seasonality was high and fire and erosion were frequent, Erica arborea and E. scoparia woodlands dominated the coastal landscape. Subsequently, between 7,500 and 5,500 cal BP , seasonality gradually declined and thermo-mediterranean woodlands with Pistacia and Quercus ilex partially replaced Erica communities under diminished incidence of fire. After 5,500 cal BP , evergreen oak forests expanded markedly, erosion declined and lake levels increased, likely in response to increasing (summer) moisture availability. Increased anthropogenic fire disturbance triggered shrubland expansions (e.g. Tamarix and Pistacia ) around 5,000–4,500 cal BP . Subsequently around 4,000–3,500 cal BP evergreen oak-olive forests expanded massively when fire activity declined and lake productivity and anoxia reached Holocene maxima. Land-use activities during the past 4,000 years (since the Bronze Age) gradually disrupted coastal forests, but relict stands persisted under rather stable environmental conditions until ca. 200 cal BP , when agricultural activities intensified and Pinus and Eucalyptus were planted to stabilize the sand dunes. Pervasive prehistoric land-use activities since at least the Bronze Age Nuraghi period included the cultivation of Prunus , Olea europaea and Juglans regia after 3,500–3,300 cal BP , and Quercus suber after 2,500 cal BP . We conclude that restoring less flammable native Q. ilex and O. europaea forest communities would markedly reduce fire risk and erodibility compared to recent forest plantations with flammable non-native trees (e.g. Pinus , Eucalyptus ) and xerophytic shrubland (e.g. Cistus , Erica ).

Sardinia↗

Incorporating the effects of complex soil layering and thickness local variability into distributed landslide susceptibility assessments

Incorporating the influence of soil layering and local variability into the parameterizations of physics-based numerical models for distributed landslide susceptibility assessments remains a challenge. Typical applications employ substantial simplifications including homogeneous soil units and soil-hydraulic properties assigned based only on average textural classifications; the potential impact of these assumptions is usually disregarded. We present a multi-scale approach for parameterizing the distributed Transient Rainfall Infiltration and Grid-Based Regional Slope-Stability (TRIGRS) model that accounts for site-specific spatial variations in both soil thickness and complex layering properties by defining homogeneous soil properties that vary spatially for each model grid cell. These effective properties allow TRIGRS to accurately simulate the timing and distribution of slope failures without any modification of the model structure. We implemented this approach for the carbonate ridge of Sarno Mountains (southern Italy) whose slopes are mantled by complex layered soils of pyroclastic origin. The urbanized foot slopes enveloping these mountains are among the most landslide-prone areas of Italy and have been subjected to repeated occurrences of damaging and deadly rainfall-induced flow-type shallow landslides. At this scope, a primary local-scale application of TRIGRS was calibrated on physics-based rainfall thresholds, previously determined by a coupled VS2D (version 1.3) hydrological modeling and slope stability analysis. Subsequently, by taking into account the spatial distribution of soil thickness and vertical heterogeneity of soil hydrological and mechanical properties, a distributed assessment of landslide hazard was carried out by means of TRIGRS. The combination of these approaches led to the spatial assessment of landslide hazard under different hypothetical rainfall intensities and antecedent hydrological conditions. This approach to parameterizing TRIGRS can be adapted to other spatially variable soil layering and thickness to improve hazard assessments.

Mount Vesuvius↗

Evidence of preferential flow activation in the vadose zone via geophysical monitoring

Preferential pathways allow rapid and non-uniform water movement in the subsurface due to strong heterogeneity of texture, composition, and hydraulic properties. Understanding the importance of preferential pathways is crucial, because they have strong impact on flow and transport hydrodynamics in the unsaturated zone. Particularly, improving knowledge of the water dynamics is essential for estimating travel time through soil to quantify hazards for groundwater, assess aquifer recharge rates, improve agricultural water management, and prevent surface stormflow and flooding hazards. Small scale field heterogeneities cannot be always captured by the limited number of point scale measurements collected. In order to overcome these limitations, noninvasive geophysical techniques have been widely used in the last decade to predict hydrodynamic processes, due to their capability to spatialize hydrogeophysical properties with high resolution. In the test site located in Bari, Southern Italy, the geophysical approach, based on electrical resistivity tomography (ERT) monitoring, has been implemented to detect preferential pathways triggered by an artificial rainfall event. ERT-derived soil moisture estimations were obtained in order to quantitatively predict the water storage (m 3 m −3 ), water velocity (ms −1 ), and spread (m 2 ) through preferential pathways by using spatial moments analysis.

Bari↗

Regional occurrence of aqueous tungsten and relations with antimony, arsenic and molybdenum concentrations (Sardinia, Italy)

T ungsten (W) is rarely found in natural waters, yet it can be introduced into the food chain and cause potentially toxic effects. Uptake of W by plants and vegetables, or trace presence of W in drinking water are possible vectors for ingestion of W by humans. The latter is recognized as a possible cause of lymphatic leukemia. Increased uses of W might result in a degradation of water resources, with attendant adverse effects on biota and human health. Therefore, this study was aimed at investigating regional occurrence and speciation of W in aquatic systems in Sardinia, Italy, factors affecting W mobility and possible relations with other oxyanion-forming trace elements such as Sb, As and Mo. Although our results are specifically from Sardinia, the implications are broader and should prompt future studies in other areas with known high W concentrations. A total of 350 sample sites are reported here, including surface waters, groundwaters, mine drainages, thermal waters and local seawater. The waters were analyzed for major and trace components, including W, Sb, As and Mo. The waters showed a variety of major chemical compositions and W concentrations. High concentrations of W were found in some mine waters and drainages from slag heaps, with W, Sb and As up to 140, 5000 and 800 μg L −1 , respectively. The highest concentrations of W occurred under slightly alkaline pH and oxygenated conditions, and were likely due to the dissolution of scheelite [CaWO 4 ] hosted in materials with which the water came into contact. High W concentrations also were observed in thermal waters, under alkaline pH and reducing conditions, and sometimes coincided with relatively high concentrations either of As or Mo. Previous studies of W geochemistry have focused on WO 4 2− as the major dissolved form of W. For this study, we have augmented the thermodynamic database in PHREEQC to include possible formation of many other W-bearing complexes gleaned from the literature. The results of the speciation calculations with the newly added complexation reactions shows that the neutral species CaWO 4 ° and MgWO 4 ° are particularly dominant in most W-bearing waters and lead to undersaturation with respect to scheelite and other W-bearing minerals. Assessing W contamination in water systems and establishing W limits in drinking water may prevent potential adverse effects of W on human and ecosystem health.

Sardinia↗

Ten years of volcanic activity at Mt Etna: High-resolution mapping and accurate quantification of the morphological changes by Pleiades and Lidar data

The topography of Mt. Etna, Italy, is subjected to continuous modifications depending on intensity and magnitude of eruptions that frequently occur at the volcano summit and flanks. In order to make high-resolution maps of morphological changes and accurately calculate the overall volume of the erupted products (e.g., lava flows, tephra fall out, scoriae cones) in ten years, we have compared the altimetry models of Mt. Etna derived from 2005 Airborne Laser Scanning data and 2015 Pleiades stereo satellite imagery. Both models cover a common area of 400 km 2 with spatial resolution of 2 m and comparable vertical accuracy (RMSE < 0.8 m). The results show that the area most affected by the erupted products is the mid-upper portion of the volcano with an altitude ranging from 1300 m to more than 3300 m a.s.l., value reached at the summit of the North East crater. In particular, this portion changes dramatically in the eastern sector due to the birth and growth of the New South-East Crater, the invasion of dozens of lava flows in the Valle del Bove, and the formation of the 2014 scoriae cones and lava field at the base of the North-East Crater. The total volume of products erupted in the investigated period results in 284.3±15.8 x 10 6 m 3 with a yearly average volume of 28.4 x 10 6 m 3 /y comparable with the previous decades. In addition, the products emitted by the 2014 sub-terminal eruption are mapped and quantified including, for the first time, the volume of the 2014 scoriae cones generated on the eastern flank of North-East Crater This study demonstrates how a rigorous comparison between digital elevation models derived from different remote sensing techniques produce high accurate mapping and quantifications of morphological changes applicable for worldwide active volcanoes. This allows to quantify volumes and areas of erupted products reducing the error estimations, a crucial point to provide precise data often used as key parameters for many volcanic hazard studies.

Mt. Etna, Sicily↗

Towards the understanding of hydrogeochemical seismic responses in karst aquifers: A retrospective meta-analysis focused on the Apennines (Italy)

Earthquakes are known to affect groundwater properties, yet the mechanisms causing chemical and physical aquifer changes are still unclear. The Apennines mountain belt in Italy presents a rich literature of case studies documenting hydrogeochemical response to seismicity, due to the high frequency of seismic events and the presence of different regional aquifers in the area. In this study, we synthesize published data from the last 30 years in the Apennine region in order to shed light on the main mechanisms causing earthquake induced water changes. The results suggest the geologic and hydrologic setting specific to a given spring play an important role in spring response, as well as the timing of the observed response. In contrast to setting, the main focal mechanisms of earthquake and the distance between epicenter and the analyzed springs seems to present a minor role in defining the response. The analysis of different response variables, moreover, indicates that an important driver of change is the degassing of CO 2 , especially in thermal springs, whereas a rapid increase in solute concentration due to permeability enhancement is observable in different cold and shallow springs. These findings also leave open the debate regarding whether earthquake precursors can be recognized beyond site-specific responses. Such responses can be understood more comprehensively through the establishment of a regional long-term monitoring system and continuous harmonization of data and sampling strategies, achievable in the Apennine region through the set-up of a monitoring network.

Apennines mountain belt↗

Combined seismic and geodetic analysis before, during and after the 2018 Mt. Etna eruption

In December 2018, Etna volcano experienced one of the largest episodes of unrest since the installation of geophysical monitoring networks in 1970. The unrest culminated in a short eruption with a small volume of lava erupted, a significant seismic crisis and deformation of the entire volcanic edifice of magnitude never recorded before at Mount Etna. Here we describe the evolution of the 2018 eruptive cycle from the analysis of seismic and geodetic data collected in the months preceding, during, and following the intrusion. We model the space-time evolution of high-rate deformation data starting from the active source previously identified from deformation data and the propagation of seismicity in a 3-D velocity model. The intrusion model suggests emplacement of two dikes: a smaller dike located beneath the eruptive fissure and a second, deeper dike between 1 and 5 km below sea level that opened ~2 m. The rise and eruption of magma from the shallower dike did not interrupt the pressurization of a long-lasting deeper reservoir (~6 km) that induced continuous inflation and intense deformation of the eastern flank. Shortly after the intrusion, on 26 December 2018, a M L 4.8 earthquake occurred near Pisano, destroying buildings and roads in two villages. We propose a time-dependent intrusion model that supports the hypothesis of the inflation inducing flank deformation and that this process has been active since September 2018.

Mt. Etna↗