Seismic response of L��Aquila downtown from comparison between 2D synthetics spectral ratios of SH, P-SV and Rayleigh waves and observations of the 2009 earthquake sequence

This is the first part of a study on the seismic response of the L’Aquila city using 2D simulation and experimental data. We have studied two velocity-depth models with the aim of outlining the behavior of a velocity reversal in the top layer, which is associated with the stiff Brecce de L’Aquila unit (BrA). In this setting, the SMTH model is topped by a layer with about 2:1 impedance contrast with the underlying layer while the NORV model has no velocity reversal. We have simulated the propagation of SH and P-SV wavefields in the range 0–10 Hz for incidence 0°–90°. Earthquake spectral ratios of the horizontal and vertical components at six sites in L’Aquila downtown are compared to corresponding synthetics spectral ratios. The vertical component of P-SV synthetics enables us to investigate a remarkable amplification effect seen in the vertical component of the recorded strong motion. Sites AQ04 and AQ05 are best matched by synthetics from the NORV model while FAQ5 and AQ06 have a better match with synthetics spectral ratios from the SMTH model. All simulations show this behavior systematically, with horizontal and near-horizontal incident waves predicting the overall pattern of matches more clearly than vertical and near-vertical incidence. The model inferences are in agreement with new geological data reporting lateral passages in the top layer from the stiff BrA to softer sediments. Matches are good in terms of frequency of the first amplification peak and of spectral amplitude: the horizontal components have spectral ratio peaks predominantly at 0.5 Hz in the simulations and at 0.7 Hz in the data, both with amplitudes of 4, while the vertical component spectral ratios reach values of 6 at frequencies of about 1 Hz in both data and simulations. The vertical component spectral ratios are very well matched using Rayleigh waves with incidence at 90°. The NORV model without the velocity reversal predicts spectral ratio peaks for the horizontal components at frequencies up to 6 Hz. The reversal of velocity acts as a low-pass frequency filter on the horizontal components reducing the amplification effect of the sediment filled valley.     

Tier-based approaches for landslide susceptibility assessment in Europe

In the framework of the European Soil Thematic Strategy and the associated proposal of a Framework Directive on the protection and sustainable use of soil, landslides were recognised as a soil threat requiring specific strategies for priority area identification, spatial hazard assessment and management. This contribution outlines the general specifications for nested, Tier-based geographical landslide zonings at small spatial scales to identify priority areas susceptible to landslides (Tier 1) and to perform quantitative susceptibility evaluations within these (Tier 2). A heuristic, synoptic-scale Tier 1 assessment exploiting a reduced set of geoenvironmental factors derived from common pan-European data sources is proposed for the European Union and adjacent countries. Evaluation of the susceptibility estimate with national-level landslide inventory data suggests that a zonation of Europe according to, e.g. morphology and climate, and performing separate susceptibility assessments per zone could give more reliable results. To improve the Tier 1 assessment, a geomorphological terrain zoning and landslide typology differentiation are then applied for France. A multivariate landslide susceptibility assessment using additional information on landslide conditioning and triggering factors, together with a historical catalogue of landslides, is proposed for Tier 2 analysis. An approach is tested for priority areas in Italy using small administrative mapping units, allowing for relating socioeconomic census data with landslide susceptibility, which is mandatory for decision making regarding the adoption of landslide prevention and mitigation measures. The paper concludes with recommendations on further work to harmonise European landslide susceptibility assessments in the context of the European Soil Thematic Strategy.     

The influence of variable topography on the depositional behaviour of pyroclastic density currents: The examples of the Upper Pollara eruption (Salina Island,southern Italy)

Stratigraphic reconstruction of the Upper Pollara eruption has allowed for the inference of eruptive mechanisms and the distillation of a sedimentological model for pyroclastic density currents (PDCs) moving across variable topography. The pre-eruptive topography in the study area was characterised by a tuff ring-like morphology, with both inward and outward dipping slopes. Highly viscous, moderately porphyritic, dacitic to rhyolitic magmas fed the eruption, which was characterised by a Vulcanian eruptive style. The stratigraphic succession was divided into five eruption units (EUs), which result from different phases of the eruption separated by stases. Sustained columns occurred only during EU1, while PDC generation dominates EU2–5. Lithofacies analysis of the PDC deposits indicates the prevalence of massive coarse-grained deposits on the inner slopes of the Pollara crater, which are interpreted as the deposits of a flow-boundary zone dominated by granular flow or fluid escape regimes. Dune-bedded, massive to stratified lithofacies dominate the outer slopes of the Pollara crater, and are interpreted as the deposits of PDCs with flow-boundary zones in which traction played a major role. Thin, massive PDC deposits are exposed on the sub-horizontal Malfa terrace, and are interpreted as representative of flow-boundary zones dominated by a granular flow regime. The occurrence of stacked deposits indicates that most of the PDCs were characterised by unsteady pulsatory behaviour, with development of trains of pulses during their transport. The downcurrent lithofacies transitions observed for the Upper Pollara deposits have finally been compared with other similar lithofacies associations which have been described for short-lived PDCs at tuff rings, in order to discuss the influence of pre-eruptive topography on lithofacies association.     

Pairs of pulsars (PSR) and the origin of PSR

The results of research on pairs of PSR which were possibly components of disrupted binaries, seem to advance further arguments in favour of the idea, suggested by some authors, that most of PSR were formed in binary systems.     

Sources of Mw 5+ earthquakes in northeastern Italy and western Slovenia: An updated view based on geological and seismological evidence

We present an overview of the seismogenic sources of northeastern Italy and western Slovenia, included in the last version of the Database of Individual Seismogenic Sources (DISS 3.0.2) and a new definition of the geometry of the Montello Source that will be included in the next release of the database. The seismogenic sources included in DISS are active faults capable of generating M_w > 5.5 earthquakes. We describe the method and the data used for their identification and characterization, discuss some implications for the seismic hazard and underline controversial points and open issues.In the Veneto–Friuli area (NE Italy), destructive earthquakes up to M_w 6.6 are generated by thrust faulting along N-dipping structures of the Eastern Southalpine Chain. Thrusting along the mountain front responds to about 2 mm/a of regional convergence, and it is associated with growing anticlines, tilted and uplifted Quaternary palaeolandsurfaces and forced drainage anomalies. In western Slovenia, dextral strike–slip faulting along the NW–SE trending structures of the Idrija Fault System dominates the seismic release. Activity and style of faulting are defined by recent earthquakes (e.g. the M_s 5.7, 1998 Bovec–Krn Mt. and the M_w 5.2, 2004 Kobarid earthquakes), while the related recent morphotectonic imprint is still a debated matter.We reinterpreted a large set of tectonic data and developed a segmentation model for the outermost Eastern Southalpine Chain thrust front. We also proposed the association of the four major shocks of the 1976 Friuli earthquake sequence with individual segments of three major thrust fronts. Although several sub-parallel active strike–slip strands exist in western Slovenia, we were able to positively identify only two segments of the Idrija Fault System. A comparison of the regional GPS velocity with long-term geological slip-rates of the seismogenic sources included in DISS shows that from a quarter to half of the deformation is absorbed along the external alignment of thrust faults in Veneto and western Friuli. The partitioning of the deformation in western Slovenia among the different strike–slip strands could not be quantified.     

Pressure dependence of structural parameters of paragonite

The compressibility and structure of a 2M₁ paragonite with composition [Na_0.88K_0.10Ca_0.01Ba_0.01] [Al_1.97Ti_0.007Fe_0.01Mn_0.002Mg_0.006]Si_3.01Al_0.99O₁₀OH₂ were determined at pressures between 1 bar and 41 kbar, by single crystal X-ray diffraction using a Merrill-Bassett diamond anvil cell. Compressibility turned out to be largely anisotropic, linear compressibility coefficients parallel to the unit cell edges being β_a=3.5(1)·10^?4, β_b=3.6(1)·10^?4, β_c=8.3(3)·10^?4 kbar^?1 (β_a:β_b·β_c=1:1028:2.371). The isothermal bulk modulus, calculated as the reciprocal of the mean compressibility of the cell volume, was 650(20) kbar. The main features of the deformation mechanism resulting from structural refinements at pressures of 0.5, 25.4, 40.5 kbar were: –?variation in sheet thickness, showing that compression of the c parameter was mainly due to the interlayer thickness reduction from 3.07 Å at 0.5 kbar to 2.81 Å at 40.5 kbar; –?the compressibility of octahedra was greater than that of tetrahedra, the dimensional misfit between tetrahedral and octahedral sheets increased with P, so that tetrahedral rotation angel α increased from 15° at 0.5 kbar to 21.6° at 40.5 kbar; –?the basal surface corrugation (Δz) of the tetrahedral layer, due to the different dimensions of M1 and M2 octahedra and to the octahedral distortion, decreased with P (Δz=0.19 and 0.12 Å at 0.5 and 40.5 kbar respectively). Comparison of the new data on paragonite with those of a K-muscovite and a Na-rich muscovite (Comodi and Zanazzi 1995) revealed a clear trend toward decreasing of compressibility when Na substitutes for K atoms in the interlayer sites.     

GIS techniques and statistical models in evaluating landslide hazard

Geographical Information Systems (GIS) and numerical cartography may greatly facilitate the development and use of statistical models for the assessment of regional landslide hazard. From a small drainage basin located in Central Italy, relevant geological and geomorphological factors were collected and processed by applying GIS technology. In particular, modules were used which both generate high-fidelity digital terrain models and automatically partition the terrain into main slope-units. The resulting information was then analysed by discriminant analysis which enabled landslide hazard and risk to be evaluated in each slope-unit. Although not lacking drawbacks, the method proved to be a feasible and cost-effective approach to landslide susceptibility assessment and mapping.     

An Innovative Approach to Meteorite Analysis by Laser‐Induced Breakdown Spectroscopy

An innovative approach of double pulse laser‐induced breakdown spectroscopy (DP‐LIBS) coupled with optical microscopy was applied to the characterisation and quantitative analysis of the Agoudal iron meteorite in bulk sample and in petrographic thin section. Qualitative analysis identified the elements Ca, Co, Fe, Ga, Li and Ni in the thin section and the whole meteorite. Two different methods, calibration‐free LIBS and one‐point calibration LIBS, were used as complementary methodologies for quantitative LIBS analysis. The elemental composition data obtained by LIBS were in good agreement with the compositional analyses obtained by traditional methods generally applied for the analysis of meteorites, such as ICP‐MS and EDS‐SEM. Besides the recognised advantages of LIBS over traditional techniques, including versatility, minimal destructivity, lack of waste production, low operating costs, rapidity of analysis, availability of transportable or portable systems, etc., additional advantages of this technique in the analysis of meteorites are precision and accuracy, sensitivity to low atomic number elements such as Li and the capacity to detect and quantify Co contents that cannot be obtained by EDS‐SEM.     

Magnesiochloritoid: Compressibility and high pressure structure refinement

The response of magnesiochloritoid to pressure has been studied by single crystal X-ray diffraction in a diamond anvil cell, using crystals with composition Mg_1.3Fe_0.7Al₄Si₂O₁₀(OH)₄. The unit cell parameters decrease from a = 9.434 (3), b = 5.452 (2), c = 18.136 (5) Å, β = 101.42° (2) (1 bar pressure) to a = 9.370 (7), b = 5.419 (5), c = 17.88 (1) Å, β = 101.5° (1) (42 kbar pressure), following a slightly anisotropic compression pattern (linear compressibilities parallel to unit cell edges: β _a = 1.85, β _b = 1.74, β_c = 3.05 × 10^?4 kbar^?1) with a bulk modulus of 1480 kbar. Perpendicular to c, the most compressible direction, the crystal structure (space group C2/c) consists of two kinds of alternating octahedral layers connected via isolated SiO₄ tetrahedra. With increasing pressure the slightly wavy layer [Mg_1.3Fe_0.7AlO₂(OH)₄] tends to flatten. Furthermore, the octahedra in this layer, with all cations underbonded, are more compressible than the octahedra in the (A1₃O₈) layer with slightly overbonded aluminum. Comparison between high-pressure and high-temperature data yields the following equations: $$\begin{gathered} a_{P,T} = 9.434{\text{ }}{\AA} - 174 \cdot 10^{ - 5} {\text{ }}{\AA}{\text{kb}}^{{\text{ - 1}}} \cdot P \hfill \\ {\text{ }} + 9 \cdot 10^{ - 5} {\text{ }}{\AA}^\circ C^{ - 1} \cdot (T - 25^\circ C) \hfill \\ b_{P,T} = 5.452{\text{ }}{\AA} - 95 \cdot 10^{ - 5} {\text{ }}{\AA}{\text{kb}}^{{\text{ - 1}}} \cdot P \hfill \\ {\text{ }} + 5 \cdot 65 \cdot 10^{ - 5} {\text{ }}{\AA}^\circ C^{ - 1} \cdot (T - 25^\circ C) \hfill \\ c_{P,T} = 18.136{\text{ }}{\AA} - 549 \cdot 10^{ - 5} {\text{ }}{\AA}{\text{kb}}^{{\text{ - 1}}} \cdot P \hfill \\ {\text{ }} + 16 \cdot 2^{ - 5} {\text{ }}{\AA}^\circ C^{ - 1} \cdot (T - 25^\circ C) \hfill \\ \end{gathered} $$ with P in kbar and T in °C. These equations indicate that the unit cell and bond geometry of magnesiochloritoid at formation conditions do not differ greatly from those at the outcrop conditions, e.g. the calculated unitcell volume is 917.3 ų at P = 16 kbar and T=500 °C, whereas the observed volume at room conditions is 914.4 ų. In addition, they show that the specific gravity increases from formation at depth to outcrop at surface conditions.