New geomorphic evidence for anticlinal growth driven by blind-thrust faulting along the northern Marche coastal belt (central Italy)

The Northern Marche coastal belt is characterised by a series of NW-SE trending, NE verging folds forming the easternmost edge of the Apennines thrust front. Several geomorphic features suggest that the folds are still growing and hence that the thrust front is active. The occurrence of several historical and instrumental earthquakes (e.g. 1672, 1690, 1786, 1875, 1916, 1930, 1972, all having M_e 5.2) suggests that the thrust faults are also seismogenic.We performed a geomorphological analysis to identify and characterise the faults driving the active folds. Our approach assumes that anomalous drainage patterns and deformed Middle-Late Pleistocene alluvial and coastal terraces are indicators of the vertical component of tectonic strain. We identified, mapped and correlated with sea-level fluctuations a sequence of alluvial and coastal terraces. Longitudinal profiles of six rivers (Conca, Foglia, Metauro, Cesano, Misa, and Esino) show that terraces (1) consistently converge downstream, suggesting that they result from regional uplift that dies out near the coast, and (2) some are slightly warped where they cross anticline axes. We interpreted as coastal terraces several land-surface remnants arranged parallel to the present coastline. Lower remnants clearly top off gently landward-tilted coastal deposits. Reconstructed coastal terraces also seem to be tectonically warped.Our results help characterise the geometry and segmentation of a system that generated the largest earthquakes of the region and suggest the loci of potential seismic gaps. We conclude that the earthquake potential of the densely populated northern Marche coastal belt may be substantially higher than currently estimated.     

The Seismotectonics of the Po Plain (Northern Italy): Tectonic Diversity in a Blind Faulting Domain

We present a systematic and updated overview of a seismotectonic model for the Po Plain (northern Italy). This flat and apparently quiet tectonic domain is, in fact, rather active as it comprises the shortened foreland and foredeep of both the Southern Alps and the Northern Apennines. Assessing its seismic hazard is crucial due to the concentration of population, industrial activities, and critical infrastructures, but it is also complicated because (a) the region is geologically very diverse, and (b) nearly all potential seismogenic faults are buried beneath a thick blanket of Pliocene–Pleistocene sediments, and thus can be investigated only indirectly. Identifying and parameterizing the potential seismogenic faults of the Po Plain requires proper consideration of their depth, geometry, kinematics, earthquake potential and location with respect to the two confronting orogens. To this end, we subdivided them into four main, homogeneous groups. Over the past 15 years we developed new strategies for coping with this diversity, resorting to different data and modeling approaches as required by each individual fault group. The most significant faults occur beneath the thrust fronts of the Ferrara-Romagna and Emilia arcs, which correspond to the most advanced and buried portions of the Northern Apennines and were the locus of the destructive May 2012 earthquake sequence. The largest known Po Plain earthquake, however, occurred on an elusive reactivated fault cutting the Alpine foreland south of Verona. Significant earthquakes are expected to be generated also by a set of transverse structures segmenting the thrust system, and by the deeper ramps of the Apennines thrusts. The new dataset is intended to be included in the next version of the Database of Individual Seismogenic Sources (DISS; http://diss.rm.ingv.it/diss/, version 3.2.0, developed and maintained by INGV) to improve completeness of potential sources for seismic hazard assessment.     

Sharp spatially constrained inversion with applications to transient electromagnetic data

Time‐domain electromagnetic data are conveniently inverted by using smoothly varying 1D models with fixed vertical discretization. The vertical smoothness of the obtained models stems from the application of Occam‐type regularization constraints, which are meant to address the ill‐posedness of the problem. An important side effect of such regularization, however, is that horizontal layer boundaries can no longer be accurately reproduced as the model is required to be smooth. This issue can be overcome by inverting for fewer layers with variable thicknesses; nevertheless, to decide on a particular and constant number of layers for the parameterization of a large survey inversion can be equally problematic. Here, we present a focusing regularization technique to obtain the best of both methodologies. The new focusing approach allows for accurate reconstruction of resistivity distributions using a fixed vertical discretization while preserving the capability to reproduce horizontal boundaries. The formulation is flexible and can be coupled with traditional lateral/spatial smoothness constraints in order to resolve interfaces in stratified soils with no additional hypothesis about the number of layers. The method relies on minimizing the number of layers of non‐vanishing resistivity gradient, instead of minimizing the norm of the model variation itself. This approach ensures that the results are consistent with the measured data while favouring, at the same time, the retrieval of horizontal abrupt changes. In addition, the focusing regularization can also be applied in the horizontal direction in order to promote the reconstruction of lateral boundaries such as faults. We present the theoretical framework of our regularization methodology and illustrate its capabilities by means of both synthetic and field data sets. We further demonstrate how the concept has been integrated in our existing spatially constrained inversion formalism and show its application to large‐scale time‐domain electromagnetic data inversions.     

A seismic source zone model for the seismic hazard assessment of the Italian territory

We designed a new seismic source model for Italy to be used as an input for country-wide probabilistic seismic hazard assessment (PSHA) in the frame of the compilation of a new national reference map.

We started off by reviewing existing models available for Italy and for other European countries, then discussed the main open issues in the current practice of seismogenic zoning.

The new model, termed ZS9, is largely based on data collected in the past 10 years, including historical earthquakes and instrumental seismicity, active faults and their seismogenic potential, and seismotectonic evidence from recent earthquakes. This information allowed us to propose new interpretations for poorly understood areas where the new data are in conflict with assumptions made in designing the previous and widely used model ZS4.

ZS9 is made out of 36 zones where earthquakes with M_w > = 5 are expected. It also assumes that earthquakes with M_w up to 5 may occur anywhere outside the seismogenic zones, although the associated probability is rather low. Special care was taken to ensure that each zone sampled a large enough number of earthquakes so that we could compute reliable earthquake production rates.

Although it was drawn following criteria that are standard practice in PSHA, ZS9 is also innovative in that every zone is characterised also by its mean seismogenic depth (the depth of the crustal volume that will presumably release future earthquakes) and predominant focal mechanism (their most likely rupture mechanism). These properties were determined using instrumental data, and only in a limited number of cases we resorted to geologic constraints and expert judgment to cope with lack of data or conflicting indications. These attributes allow ZS9 to be used with more accurate regionalized depth-dependent attenuation relations, and are ultimately expected to increase significantly the reliability of seismic hazard estimates.     

Experimental solidification of anhydrous latitic and trachytic melts at different cooling rates: The role of nucleation kinetics

Two sets of cooling experiments were run at atmospheric conditions for two anhydrous starting latitic and trachytic melts: 1) five cooling rates (25, 12.5, 3, 0.5, and 0.125 °C/min) between 1300° and 800 °C, and 2) a 0.5 °C/min cooling rate from 1300 °C with quench temperatures at 1200°, 1100°, 1000° and 900 °C. Trachytic run-products are invariably glassy. Nucleation is also suppressed in the latitic run-products at the three highest cooling rates. Conversely, in the 0.5 and 0.125 °C/min runs, latites have a crystal content of  90 vol.%. The phases are: plagioclase, clinopyroxene, glass and iron-bearing oxide (in order of abundance). The variable quench temperatures, investigated by coupling experiments with Pt wire and Pt capsule sample containers in set 2, again did not produce crystallization of trachyte, whereas latitic samples are characterized by  10 vol.% of oxides, pyroxenes and plagioclase (in order of appearance), at temperature < 1000 °C. Effects of (preferential) heterogeneous nucleation on sample holders, of superheating degree, and chemical species loss during cooling are absent for both melt compositions. The difference of solidification paths between these two silicate melts can be ascribed only to their small chemical differences. In comparison with calculated equilibrium conditions all the experimental latitic and trachytic run-products revealed strong kinetic effects, interpretable in the light of the nucleation theory. The glass-forming ability (GFA) of trachyte is higher, whereas their critical cooling rate (Rc) is lower (< 0.125 °C/min), in comparison to latitic melts (Rc > 0.5 °C/min). The experimental results carried out in this study can be applied to lava flows and domes; trachytic lavas are able to flow for longer period with respect to latitic ones in a metastable condition. Glass-rich terrestrial lavas, i.e. obsidians, can be the result of sluggish nucleation kinetics due to the relative high polymerisation of evolved silicate melts.     

Numerical simulation of the tsunami generated by a past catastrophic landslide on the volcanic island of Ischia, Italy

The island of Ischia, Gulf of Naples, Italy, like many other volcanic islands is affected by mass failures, that are mainly related to secondary volcanic processes such as slope steepening and seismic shaking. The block resurgence of its main relief, Mount Epomeo, has been recognised to contribute cyclically to mass instability and cause landslides, that occasionally may reach the sea and start tsunamis. In this work we explore the consequences of the Ischia Debris Avalanche (IDA), a flank collapse that occurred in historical times, and involved the whole Mount Epomeo edifice including its submarine portion, and that may have caused gigantic sea waves affecting all the coasts of Ischia and of the Gulf of Naples. The IDA and the generated tsunami have been taken as the worst-case scenario for the occurrence of a new tsunami in the area. They have been simulated through numerical codes developed and maintained by the University of Bologna. The simulation shows that the IDA-induced tsunami attacks severely all the coasts of the Gulf of Naples with the highest waves obtained for the island of Ischia, the island of Capri and the peninsula of Sorrento. The propagation pattern of the IDA tsunami can be used to get hints on the impact that such an event may have had on early populations habiting Gulf of Naples, but also to get clues on the area that could be most severely hit by a tsunami generated by a smaller-scale landslide that may occur in the same source zone.     

A note on maar eruption energetics: current models and their application

Hydromagmatic eruptions convert thermal into mechanical energy via the expansion of ground- and/or surface-water. Several models address the energetics of these eruptions based on different physical-volcanological approaches. Here we test different models with two case studies in the Colli Albani Volcanic District (central Italy): the monogenetic Prata Porci and the polygenetic Albano maars. Test results are mutually consistent, and show cumulative mechanical energy releases on the order of 10¹⁵–10¹⁷ J for single maars. The fraction of thermal energy turned into mechanical ranges from 0.45 (as calculated from the theoretical maximum mechanical energy), through 0.1 (calculated from country rock fragmentation, crater formation and ballistic ejection), to 0.03 (derived from magma fragmentation by thermohydraulic explosions). It appears that the energy released during the most intense hydromagmatic events may account for a dominant fraction of the total mechanical energy released during the whole maar eruptive histories. Finally, we consider the role of magmatic explosive activity intervening during maar eruptions in causing departures from predictions of the models evaluated.     

Backward automatic calibration for three-dimensional landslide models

Back-analysis is broadly used for approaching geotechnical problems when monitoring data are available and information about the soils properties is of poor quality.For landslide stability assessment back-analysis calibration is usually carried out by time consuming trial-and-error procedure.This paper presents a new automatic Decision Support System that supports the selection of the soil parameters for three-dimensional models of landslides based on monitoring data.The method considering a pool of possible solutions,generated through permutation of soil parameters,selects the best ten configurations that are more congruent with the measured displacements.This reduces the operator biases while on the other hand allows the operator to control each step of the computation.The final selection of the preferred solution among the ten best-fitting solutions is carried out by an operator.The operator control is necessary as he may include in the final decision process all the qualitative elements that cannot be included in a qualitative analysis but nevertheless characterize a landslide dynamic as a whole epistemological subject,for example on the base of geomorphological evidence.A landslide located in Northeast Italy has been selected as example for showing the system potentiality.The proposed method is straightforward,scalable and robust and could be useful for researchers and practitioners.