The Deep Seismicity of the Tyrrhenian Sea*

The deep seismicity of the Tyrrhenian Sea is analysed using data from a new instrumental catalogue of the seismicity of the Italian area. We use algorithms for the determination of absolute and relative hypocentral locations and for the evaluation of the geometry and coherence of the state of stress within the subducting slab. The structure of the Benioff zone, although simpler than previously indicated, reveals anomalous traits both in the seismicity distribution and in the stress geometry, confirming that standard subduction models cannot be applied in the Tyrrhenian region. The velocity anomaly and the location of few isolated events indicate that the subducted slabs extend to the north along the Apenninic chain approximately to the latitude of the Irpinia region, in Central Italy (? 42°N).     

Seismic potential of Southern Italy

To improve estimates of the long-term average seismic potential of the slowly straining South Central Mediterranean plate boundary zone, we integrate constraints on tectonic style and deformation rates from geodetic and geologic data with the traditional constraints from seismicity catalogs. We express seismic potential (long-term average earthquake recurrence rates as a function of magnitude) in the form of truncated Gutenberg–Richter distributions for seven seismotectonic source zones. Seismic coupling seems to be large or even complete in most zones. An exception is the southern Tyrrhenian thrust zone, where most of the African–European convergence is accommodated. Here aseismic deformation is estimated to range from at least 25% along the western part to almost 100% aseismic slip around the Aeolian Islands. Even so, seismic potential of this zone has previously been significantly underestimated, due to the low levels of recorded past seismicity. By contrast, the series of 19 M6–7 earthquakes that hit Calabria in the 18th and 19th century released tectonic strain rates accumulated over time spans up to several times the catalog duration, and seismic potential is revised downward. The southern Tyrrhenian thrust zone and the extensional Calabrian faults, as well as the northeastern Sicilian transtensional zone between them (which includes the Messina Straits, where a destructive M7 event occurred in 1908), all have a similar seismic potential with minimum recurrence times of M ≥ 6.5 of 150–220 years. This potential is lower than that of the Southern Apennines (M ≥ 6.5 recurring every 60 to 140 years), but higher than that of southeastern Sicily (minimum M ≥ 6.5 recurrence times of 400 years). The high seismicity levels recorded in southeastern Sicily indicate some clustering and are most compatible with a tectonic scenario where the Ionian deforms internally, and motions at the Calabrian Trench are small. The estimated seismic potential for the Calabrian Trench and Central and Western Sicily are the lowest (minimum M ≥ 6.5 recurrence times of 550–800 years). Most zones are probably capable of generating earthquakes up to magnitudes 7–7.5, with the exception of Central and Western Sicily where maximum events sizes most likely do not exceed 7.     

Evaluating Earthquake-Triggered Landslide Hazard at the Basin Scale Through Gis in the Upper Sele River Valley

To evaluate techniques for assessing earthquake-triggeredlandslide hazard in the Southern Apennines (Italy), a GIS-based analysis was used to modelseismically induced slope deformations. Geological, geotechnical, geomorphological and seismologicaldata were integrated into a standard earthquake slope stability model. The model assessed thelandslide potential that existed during the 1980 Irpinian earthquake in the Upper Sele river Valley.The standard Newmark displacement analysis, widely used for predicting the location of shallowunstable slopes, does not take into account errors and/or uncertainties in the input parameters.Therefore, a probabilistic Newmark displacement analysis technique has been used. Probabilistictechniques allow, e.g., an estimation of the probability that a slope will exceed a certain criticalvalue of Newmark displacement. In our probabilistic method, a Monte-Carlo based simulation modelis used in conjunction with a GIS. The random variability of geotechnical data is modelled by probabilitydensity functions (pdfs), while for the seismic input three different regression laws wereconsidered. Input probability distributions are sampled and the resulting values input into empiricalrelations for estimating Newmark displacement. The outcome is a map in which to each siteis related a spatial probability distribution for the expected displacement in response to seismic loading.Results of the experiments show a high grade of uncertainty in the application of the Newmarkanalysis both for the deterministic and probabilistic approach in a complex geological setting suchas the high Sele valley, quite common in the Southern Apennines. They show a strong dependence onthe reliability of the spatial data used in input, so that, when the model is used at basin scale,results are strongly influenced by local environmental condition (e.g., topography, lithology, groundwatercondition) and decrease the model performance.     

Mars in depth

P Lognonné, T Spohn and D Giardini describe the state of knowledge of Mars' interior and make the case for deep seismology to address the unanswered questions.     

Introduction to ground surface self-potential tomography

A new approach to self-potential (SP) data interpretation for the recognition of a buried causative SP source system is presented. The general model considered is characterized by the presence of primary electric sources or sinks, located within any complex resistivity structure with a flat air-earth boundary. First, using physical considerations of the nature of the electric potential generated by any arbitrary distribution of primary source charges and the related secondary induced charges over the buried resistivity discontinuity planes, a general formula is derived for the potential and the electric field component along any fixed direction on the ground surface. The total effect is written as a sum of elementary contributions, all of the same simple mathematical form. It is then demonstrated that the total electric power associated with the standing natural electric field component can be written in the space domain as a sum of cross-correlation integrals between the observed component of the total electric field and the component of the field due to each single constitutive elementary charge. By means of the cross-correlation bounding inequality, the concept of a scanning function is introduced as the key to the new interpretation procedure. In the space domain, the scanning function is the unit strength electric field component generated by an elementary positive charge. Next, the concept of charge occurrence probability is introduced as a suitable function for the tomographic imaging of the charge distribution geometry underground. This function is defined as the cross-correlation product of the total observed electric field component and the scanning function, divided by the square root of the product of the respective variances. Using this physical scheme, the tomographic procedure is described. It consists of scanning the section, through any SP survey profile, by the unit strength elementary charge, which is given a regular grid of space coordinates within the section, at each point of which the charge occurrence probability function is calculated. The complete set of calculated grid values can be used to draw contour lines in order to single out the zones of highest probability of concentrations of polarized, primary and secondary electric charges. An extension to the wavenumber domain and to three-dimensional tomography is also presented and discussed. A few simple synthetic examples are given to demonstrate the resolution power of the new SP inversion procedure.     

Exhumation patterns along shallow low‐angle normal faults: an example from the Altotiberina active fault system (Northern Apennines,Italy)

A multi‐method approach (palaeothermal and thermochronological analyses; thermal modelling) is applied to reconstruct the exhumation history of the Altotiberina Fault (ATF), a representative example of crustal‐scale active low‐angle normal faulting in the Northern Apennines (Italy). Thermal maturity and thermochronological data yield similar burial histories but different exhumation patterns for the sedimentary successions in the hangingwall and the footwall of the ATF. Since 3.8 Ma, the ATF footwall has exhumed at rates of 0.90 mm a^?1. Exhumation led to bending and deactivation of the ATF uppermost portion as a result of tectonic unloading and isostatic adjustment, followed by migration of extension and the development of a set of domino‐like, east‐dipping normal faults, rooting on the buried portion of the ATF. ATF activity and isostatic rebound exhumed Triassic rock units from depths of about 4 km. We suggest that isostatic instability is accommodated at shallow crustal levels, in a similar way to what is observed on larger structures at mid‐low crustal levels.     

Considering large earthquake clustering in seismic risk analysis

Earthquakes cluster in space and time resulting in nonlinear damage effects. We compute earthquake interactions using the Coulomb stress transfer theory and dynamic vulnerability from the concept of ductility capacity reduction. We combine both processes in the generic multi-risk framework where risk scenarios are simulated using a variant of the Markov chain Monte Carlo method. We apply the proposed approach to the thrust fault system of northern Italy, considering earthquakes with characteristic magnitudes in the range ~[6, 6.5], different levels of tectonic loading \(\dot{\tau }\) = {10⁻⁴, 10⁻³, 10⁻²} bar/year and a generic stock of fictitious low-rise buildings with different ductility capacities μ _Δ = {2, 4, 6}. We describe the process’ stochasticity by non-stationary Poisson earthquake probabilities and by binomial damage state probabilities. We find that earthquake clustering yields a tail fattening of the seismic risk curve, the effect of which is amplified by damage-dependent fragility due to clustering. The impact of clustering alone is in average more important than dynamic vulnerability, the spatial extent of the former phenomenon being greater than of the latter one.

     

Exploiting historical rainfall and landslide data in a spatial database for the derivation of critical rainfall thresholds

Critical rainfall thresholds for landslides are powerful tools for preventing landslide hazard. The thresholds are commonly estimated empirically starting from rainfall events that triggered landslides in the past. The creation of the appropriate rainfall–landslide database is one of the main efforts in this approach. In fact, an accurate agreement between the landslide and rainfall information, in terms of location and timing, is essential in order to correctly estimate the rainfall–landslide relationships. A further issue is taking into account the average moisture conditions prior the triggering event, which reasonably may be crucial in determining the sufficient amount of precipitation. In this context, the aim of this paper is exploiting historical landslide and rainfall data in a spatial database for the derivation of critical rainfall thresholds for landslide occurrence in Sicily, southern Italy. The hourly rainfall events that caused landslides occurred in the twentieth century were specifically identified and reconstructed. A procedure was proposed to automatically convert rain guages charts recorded on paper tape into digital format and then to provide the cumulative rainfall hyetograph in digital format. This procedure is based on a segmentation followed by signal recognition techniques which allow to digitalize and to recognize the hyetograph automatically. The role of rainfall prior to the landslide events was taken into account by including in the analysis the rainfall occurred 5, 15 and 30 days before each landslide. Finally, cumulated rainfall duration thresholds for different exceedance probability levels were determined. The obtained thresholds resulted in agreement with the regional curves proposed by other authors for the same area; antecedent rainfall turned out to be particularly important in triggering landslides.     

Earthquakes in Switzerland and surrounding regions during 2006

This report of the Swiss Seismological Service summarizes the seismic activity in Switzerland and surrounding regions during 2006. During this period, 572 earthquakes and 91 quarry blasts were detected and located in the region under consideration. Of these earthquakes, two occurred in conjunction with the construction of the new Gotthard railway tunnel and 165 were induced artificially by the stimulation of a proposed geothermal reservoir beneath the city of Basel. With 20 events with M _L ≥ 2.5, five of which were artificially induced, the seismic activity in the year 2006 was far below the average over the previous 31 years. Nevertheless, six events were felt by the public, most prominently the strongest of the induced Basel events (M _L 3.4), which caused some non-structural building damage. Noteworthy are also the two earthquakes near Cortaillod (M _L 3.2), on the shore of Lake Neuchatel, and in Val Mora (M _L 3.5), between the Engadin and Val Müstair, as well as the 42 aftershocks of the M _L 4.9 Vallorcine earthquake, between Martigny and Chamonix, of September 2005. Editorial handling: Stefan Bucher