Search

International audience; Unravelling the evolution of supradetachment basins developed in the hanging wall of low‐angle detachment faults may be an invaluable tool in reconstructing the tectonic evolution of highly extended terrains. These basins may record major regional tectonic events related to the exhumation of metamorphic core complexes, and the reconstruction of their evolution helps to quantify the amount of extension accommodated by such processes. Here we present stratigraphic and structural field evidence and micropaleontological constraints to the Neogene‐to‐Quaternary evolution of the supradetachment Gediz Graben that developed on top of the exhuming Central Menderes Massif (SW Turkey). This basin displays three different structural styles during its evolution: (i) it initiated as a ramp basin following the activation of the Gediz Detachment in the Middle Miocene; (ii) evolved as a half graben during the Late Miocene following the activation of high‐angle brittle faults at its southern margin; and (iii) reached its final symmetric graben configuration in Late Pliocene (?) – Quaternary times following the activation of its northern margin. New micropaleontological data document a short‐lived upper Tortonian marine episode in the basin, and major along‐strike variations in exhumation are documented on its southern margin. Our reconstruction shows how sedimentary basins originally formed in the hanging‐wall of detachment faults may eventually end up in tectonic contact with the mylonitic footwall. Finally, we highlight the importance of magmatism in localizing the deformation in highly extended terrains and in controlling the evolution of supradetachment systems.

Magmatic conditions prior to explosive eruption are often investigated using geochemical signatures in glassy components of pyroclastic deposits and related to magmatic processes at depth. One important process is fractional crystallization, which causes systematic changes to the SiO 2 /Al 2 O 3 ratio of the residual melt that can be determined by observation of the mineralogy of fully crystallized lavas, by experimental petrology, and by magmatic modeling. However, for many alkaline-mafic pyroclastic deposits, the record of residual melt compositions is obscured by alteration, commonly affecting more than 50% of pyroclastic rock components including reactive glass and some susceptible minerals. In this study, melt signatures of SiO 2 /Al 2 O 3 represented heterogeneously by the scarce fresh glass and abundant, zeolitized proxy-glass in the alkaline deposits of a major, caldera-forming eruption were used in conjunction with a model system (Rhyolite-MELTS) to reconstruct residual melt compositions and characteristics that existed immediately prior to explosive eruption. Through the model, full major oxide compositions of residual melts and fractionally crystallizing minerals become accessible, with associated constraints on volatiles and physical characteristics (melt temperature, density, viscosity). The use of zeolitized proxy-glass signatures relies on established and deposit-specific evidence for "hydrologically closed" systems that suggests the SiO 2 /Al 2 O 3 ratio is closely retained through initial alteration reactions and therefore closely representative of SiO 2 /Al 2 O 3 in the precursor glass (erupted melt). The relationship is supported by a review of available, paired data (R 2 = 0.94). Therefore, magmatic system data for the abundant and pervasive fine ash fraction of pyroclastic deposits can be investigated using this method and can progress more deeply beyond the widely used simple affiliation to igneous rock classification. Model-predicted magmatic mineral compositions (clinopyroxene, spinel, and nepheline as demonstrated here) serve to validate a case study reconstruction by comparison with compositions reported from natural and experimental samples. This predictive capability of the novel procedure is demonstrated in the case of a major caldera-forming eruption, the 355 ka Villa Senni event of the quiescent Colli Albani volcano, Rome, Italy, and its pervasively zeolitized Tufo Lionato deposit (>50 km 3 ). The key finding is that a more-evolved residual melt fraction has been revealed, based on a reconstructed SiO 2 /Al 2 O 3 ratio of 2.05 relative to that of the parent magma at 2.68, with implications for a reappraisal of pre-eruptive conditions and eruption mechanisms, and potentially for similar patterns across the volcanic stratigraphy and for other alkaline volcanoes.

International audience; Trench motions represent the surface expression of the interaction between subducting plates and the underlyingmantle, but the inherent dynamics are not fully understood. One interesting case is the migration of the Izu-Bonin Marianas trench (IBM) that accommodates the subduction of the Pacific beneath the Philippine Sea Plate (PSP), which is in turn subducting beneath the Eurasian plate along the Ryukyu trench. The history of the IBM trench is dominated by fast, episodic retreat from 40 to 15 Ma. However, around 10–5 Ma, the IBM trenchreversed its motion from retreat to advance. The switch in trench motion occurred soon after the breakoff of thePSP slab along the Ryukyu trench and the onset of new subduction, and represents a fundamental change in thedynamics of the western Pacific subduction zones. Here, Following the modelling study of Čížková and Bina(2015), which suggested a link between IBM trench advance and Ryukyu subduction, we run 2-D numericalexperiments to test the influence of a newly formed Ryukyu slab on the established IBM subduction zone, we runtwo-dimensional numerical experiments to test the influence of this newly formed Ryukyu slab on the IBMsubduction zone. The results from our geodynamic model compare favourably with the reconstructed trend, indicating that the switch in trench motion along the IBM trench may indeed be related to the onset of a new subduction zone along the Ryukyu trench. Our analysis substantiates the idea that advancing trench motions in the western Pacific are due to the establishment of a double subduction system. Further analysis of such dynamics provides insights for the mechanisms controlling subducting plate and trench motions and mantle force transmission.

The Arabian Plate shows a strong asymmetry between its Shield and Platform, in terms of topography, seismic velocity and density structure of the upper mantle. This asymmetry also results in significant rheological differences between these blocks, as revealed by the effective elastic thickness (EET) estimates, obtained using a spectral gravity method. However, these estimates may be biased due to various factors: therefore, other approaches based on a direct rheological modeling of the lithospheric structure should be employed to verify these results. In this study, we use a recent model of the lithosphere, based on an integrative interpretation of the gravity field and seismic tomography, to correct an initial thermal model obtained from the inversion of seismic velocity, assuming a uniform composition. The results are used together with the most recent crustal model of the Arabian Plate to construct two alternative models of strength and EET of the lithosphere. The first model (Model I) assumes a constant value of 10− 15 s− 1 for the strain rates. In the second model (Model II), we used the strain rates obtained from a global mantle flow model. Model I confirms the asymmetry in the rigidity of the Shield and Platform. In contrast, Model II shows that the influence of the variable strain rates causes a significant increase in the strength and EET of the central and eastern part of the Shield and in contrast to previous studies, reveals that most of the Arabian Plate is a long-term stable tectonic feature, predominantly characterized by large EET values (≥ 70 km). Do you want to read the rest of this article?Request full-text Strength and elastic thickness variations in the Arabian Plate: A combination of temperature, composition and strain rates of the lithosphere. Available from: https://www.researchgate.net/publication/314832634_Strength_and_elastic_thickness_variations_in_the_Arabian_Plate_A_combination_of_temperature_composition_and_strain_rates_of_the_lithosphere [accessed Apr 18, 2017].

Litho-structural control on the spatial and temporal evolution of landslides is one of the major typical aspects on slopes constituted of structurally complex sequences. Mainly focused on instabilities of the earth flow type, a semi-quantitative analysis has been developed with the purpose of identifying and characterizing litho-structural control exerted by bedding on slopes and its effects on landsliding. In quantitative terms, a technique for azimuth data interpolation, Non-continuous Azimuth Distribution Methodological Approach (NADIA), is presented by means of a GIS software application. In addition, processed by NADIA, two indexes have been determined: (i) Δ, aimed at defining the relationship between the orientation of geological bedding planes and slope aspect, and (ii) C, which recognizes localized slope sectors in which the stony component of structurally complex formations is abundant and therefore operates an evolutive control of landslide masses. Furthermore, some Litho-Structural Models (LSMs) of slopes are proposed aiming at characterizing recurrent forms of structural control in the source, channel and deposition areas of gravitational movements. In order to elaborate evolutive models controlling landslide scenarios, LSMs were qualitatively related and compared with Δ and C quantitative indexes. The methodological procedure has been applied to a lithostructurally complex area of Southern Italy where data about azimuth measurements and landslide mapping were known. It was found that the proposed methodology enables the recognition of typical control conditions on landslides in relation to the LSMs. Different control patterns on landslide shape and on style and distribution of the activity resulted for each LSM. This provides the possibility for first-order identification to be made of the spatial evolution of landslide bodies.

The Iranian ophiolites are part of the vast orogenic suture zones that mark the Alpine-Himalayan convergence zone. Few petrological and geochronological data are available from these ophiolitic domains, hampering a full assessment of the timing and regimes of subduction zone metamorphism and orogenic construction in the region. This paper describes texture, geochemistry, and the pressure-temperature path of the Early Cretaceous mafic granulites that occur within the Tertiary Sabzevar ophiolitic suture zone of NE Iran. Whole rock geochemistry indicates that the Sabzevar granulites are likely derived from a MORB-type precursor. They are thus considered as remnants of a dismembered dynamo-thermal sole formed during subduction of a back-arc basin (proto-Sabzevar Ocean) formed in the upper-plate of the Neotethyan slab. The metamorphic history of the granulites suggests an anticlockwise pressure-temperature loop compatible with burial in a hot subduction zone, followed by cooling during exhumation. Transition from a nascent to a mature stage of oceanic subduction is the geodynamic scenario proposed to accomplish for the reconstructed thermobaric evolution. When framed with the regional scenario, results of this study point to diachronous and independent tectonic evolutions of the different ophiolitic domains of central Iran, for which a growing disparity in the timing of metamorphic equilibration and of pressure-temperature paths can be expected to emerge with further investigations.

The formation of mountain belts or rift zones is commonly attributed to interactions between plates along their boundaries, but the widely distributed deformation of Asia from Himalaya to the Japan Sea and other back‐arc basins is difficult to reconcile with this notion. Through comparison of the tectonic and kinematic records of the last 50 Ma with seismic tomography and anisotropy models, we show that the closure of the former Tethys Ocean and the extensional deformation of East Asia can be best explained if the asthenospheric mantle transporting India northward, forming the Himalaya and the Tibetan Plateau, reaches East Asia where it overrides the westward flowing Pacific mantle and contributes to subduction dynamics, distributing extensional deformation over a 3,000‐km wide region. This deep asthenospheric flow partly controls the compressional stresses transmitted through the continent‐continent collision, driving crustal thickening below the Himalayas and Tibet and the propagation of strike‐slip faults across Asian lithosphere further north and east, as well as with the lithospheric and crustal flow powered by slab retreat east of the collision zone below East and SE Asia. The main shortening direction in the deforming continent between the collision zone and the Pacific subduction zones may in this case be a proxy for the direction of flow in the asthenosphere underneath, which may become a useful tool for studying mantle flow in the distant past. Our model of the India‐Asia collision emphasizes the role of asthenospheric flow underneath continents and may offer alternative ways of understanding tectonic processes., Key Points Coeval India‐Asia collision and Pacific back‐arc extension explained by asthenospheric flowMantle flow correlates with long‐term kinematics from India to Pacific borderIndian mantle reaches Pacific border and Asian lithosphere spreads above ponding slabs

We consider a system made of n coplanar faults and study the conditions under which it may originate a sequence of n seismic events of similar magnitudes. The system is placed in an elastic half-space and is subject to a constant and uniform strain rate by tectonic motions. The state of the system is described by n variables that are the Coulomb stresses of the faults. Coulomb stresses increase steadily during the interseismic intervals and change suddenly in connection with seismic events. A permutation of the first n integers can be associated with each state of the system, expressing the order of Coulomb stresses according to their magnitudes. The permutation controls the evolution of the system. The state space can be divided into n! subsets corresponding to each of the n! permutations. The representative point of the system moves continuously in one of the subsets during the interseismic intervals and switches to a different subset in connection with seismic events. The order of events in a sequence can be also expressed as a permutation of the first n integers, so that the number of possible sequences is equal to n!. The order of events is determined by the initial stress state, the stress drops and the stress transfers associated with each event. If the sequence is made of n similar events and the differences between Coulomb stresses are always greater than the mutual stress transfers, the order is given by the initial permutation. However, the order implicit in the initial permutation is generally modified during the sequence, due to changes in the differences between Coulomb stresses and to different stress drops of the events. The model allows the retrieval of the stress states of a fault system from the observation of a seismic sequence. As an example, the model is applied to the 2012 Emilia (Italy) seismic sequence and enlightens the complex interplay between the fault dislocations that produced the observed order of events.

Seismological, tsunami and geodetic observations have shown that subduction zones are complex systems where the properties of earthquake rupture vary with depth as a result of different pre-stress and frictional conditions. A wealth of earthquakes of different sizes and different source features (e.g. rupture duration) can be generated in subduction zones, including tsunami earthquakes, some of which can produce extreme tsunamigenic events. Here, we offer a geological perspective principally accounting for depth-dependent frictional conditions, while adopting a simplified distribution of on-fault tectonic pre-stress. We combine a lithology-controlled, depth-dependent experimental friction law with 2D elastodynamic rupture simulations for a Tohoku-like subduction zone cross-section. Subduction zone fault rocks are dominantly incohesive and clay-rich near the surface, transitioning to cohesive and more crystalline at depth. By randomly shifting along fault dip the location of the high shear stress regions (“asperities”), moderate to great thrust earthquakes and tsunami earthquakes are produced that are quite consistent with seismological, geodetic, and tsunami observations. As an effect of depth-dependent friction in our model, slip is confined to the high stress asperity at depth; near the surface rupture is impeded by the rock-clay transition constraining slip to the clay-rich layer. However, when the high stress asperity is located in the clay-to-crystalline rock transition, great thrust earthquakes can be generated similar to the Mw 9 Tohoku (2011) earthquake.

Recent declines in ice shelf and sea ice extent experienced in polar regions highlight the importance of evaluating variations in local weather patterns in response to climate change. Airborne mineral particles (dust) transported through the atmosphere and deposited on ice sheets and glaciers in Antarctica and Greenland can provide a robust set of tools for resolving the evolution of climatic systems through time. Here we present the first high time resolution radiogenic isotope (strontium and neodymium) data for Holocene dust in a coastal East Antarctic ice core, accompanied by rare earth element composition, dust concentration, and particle size distribution during the last deglaciation. We aim to use these combined ice core data to determine dust provenance, with variations indicative of shifts in either dust production, sources, and/or transport pathways. We analyzed a series of 17 samples from the Taylor Dome (77°47′47″S, 158°43′26″E) ice core, 113–391 m in depth from 1.1–31.4 ka. Radiogenic isotopic and rare earth element compositions of dust during the last glacial period are in good agreement with previously measured East Antarctic ice core dust records. In contrast, the Holocene dust dataset displays a broad range in isotopic and rare earth element compositions, suggesting a shift from long-range transported dust to a more variable, local input that may be linked to the retreat of the Ross Ice Shelf during the last deglaciation. Observed changes in the dust cycle inferred from a coastal East Antarctic ice core can thus be used to infer an evolving local climate.

Nonstructural components (NSCs) should be subjected to a careful and rational seismic design, in order to reduce the economic loss and to avoid threats to the life safety, as well as what concerns structural elements. The design of NSCs is based on the evaluation of the maximum inertia force, which is related to the floor spectral accelerations. The question arises as to whether Eurocode 8 is able to predict actual floor response spectral accelerations occurring in structures designed according to Eurocode 8. A parametric study is conducted on five RC frame structures in order to evaluate the floor response spectra. The structures, designed according to Eurocode 8, are subjected to a set of earthquakes, compatible with the design response spectrum. Time-history analyses are performed both on elastic and inelastic models of the considered structures. Eurocode formulation for the evaluation of the seismic demand on NSCs does not well fit the numerical results. Some comments on the target spectrum provided by AC 156 for the seismic qualification of NSC are also included.

Stability of isolation rubber bearings is a topic widely studied and concerns (1) the critical load capacity in the undeformed configuration, under long-term load (gravity) effects; and (2) the stability condition under short-term vertical pressure (due to gravity plus seismic loads) at large lateral deformations. In this paper the problem of elastomeric bearing stability under large lateral displacements is addressed through FEM parametric analysis; rubber bearings typical of current design practice, characterized by S1 = 20 are considered, with the value of the secondary shape factor S2 being varied between 1.5 and 6.2 in order to assess the effect of slenderness on the mechanical behavior, failure mode, and interaction vertical pressure–shear deformation. The analysis results show that the sensitivity of the shear response to the applied vertical pressure is directly related to the value of S2 and that S2 has an overwhelming effect on the stability behavior and shear response of the bearing, as compared to the effect of the primary shape factor and of the rubber shear modulus. Finally, S2 is the parameter which governs the failure mode of the bearing in the seismic condition (vertical pressure–shear deformation). On the basis of these results, design implications are discussed.

Here we report on the first assessment of volatile fluxes from the hyperacid crater lake hosted within the summit crater of Copahue, a very active volcano on the Argentina-Chile border. Our observations were performed using a variety of in situ and remote sensing techniques during field campaigns in March 2013, when the crater hosted an active fumarole field, and in March 2014, when an acidic volcanic lake covered the fumarole field. In the latter campaign, we found that 566 to 1373 t d -1 of SO 2 were being emitted from the lake in a plume that appeared largely invisible. This, combined with our derived bulk plume composition, was converted into flux of other volcanic species (H 2 O ~ 10989 t d -1 , CO 2 ~ 638 t d -1 , HCl ~ 66 t d -1 , H 2 ~ 3.3 t d -1 , and HBr ~ 0.05 t d -1 ). These levels of degassing, comparable to those seen at many open-vent degassing arc volcanoes, were surprisingly high for a volcano hosting a crater lake. Copahue's unusual degassing regime was also confirmed by the chemical composition of the plume that, although issuing from a hot (65°C) lake, preserves a close-to-magmatic signature. EQ3/6 models of gas-water-rock interaction in the lake were able to match observed compositions and demonstrated that magmatic gases emitted to the atmosphere were virtually unaffected by scrubbing of soluble (S and Cl) species. Finally, the derived large H 2 O flux (10,988 t d -1 ) suggested a mechanism in which magmatic gas stripping drove enhanced lake water evaporation, a process likely common to many degassing volcanic lakes worldwide. Key Points First volatile flux record of crater lake gas emissions from Copahue volcano Magmatic gases breaching through the lake surface Magmatic gas stripping drive enhanced lake water evaporation Fil: Tamburello, G.. Università degli Studi di Palermo; Italia Fil: Caselli, Alberto Tomás. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Estudios Andinos "Don Pablo Groeber". Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Estudios Andinos "Don Pablo Groeber"; Argentina Fil: Tassi, F.. Università degli Studi di Firenze; Italia Fil: Vaselli, O.. Università degli Studi di Firenze; Italia Fil: Calabrese, S.. Università degli Studi di Palermo; Italia Fil: Rouwet, D.. Istituto Nazionale di Geofisica e Vulcanologia; Italia Fil: Capaccioni, B.. Universidad de Bologna; Italia Fil: Di Napoli, R.. Università degli Studi di Palermo; Italia Fil: Cardellini, C.. Università di Perugia; Italia Fil: Chiodini, G.. Istituto Nazionale di Geofisica e Vulcanologia; Italia Fil: Bitetto, M.. Università degli Studi di Palermo; Italia Fil: Brusca, L.. Istituto Nazionale Geofisica e Vulcanologia; Italia Fil: Bellomo, S.. Istituto Nazionale Geofisica e Vulcanologia; Italia Fil: Aiuppa, A.. Istituto Nazionale Geofisica e Vulcanologia; Italia

The Val d'Agri area provides the opportunity to analyse active structures in a seismic region for which a large amount of subsurface data is available. This area, which was struck in 1857 by one of the most destructive earthquakes in Italy (MW = 7.03), represents a unique natural laboratory to gain new insights into geometry, modes and rates of faulting controlling crustal deformation in an actively extending orogen. In this study, a crustal geological section through the southern Apennines is discretized into a finite element model (FEM). We present a 2D elastoplastic FEM that reproduces stress perturbations and strain field around the Val d'Agri active fault system. The influence of fault strand activity on interseismic crustal deformation is tested by a series of computer models, whose predictions are compared with the horizontal velocity components of continuous GPS sites in the region and with stress directions and geological data. The best fit with available geological and geophysical constraints is obtained with a 300 km long, 29 km deep model formed by a multilayer including three components having different rheological characteristics and including several shallow, locked fault segments, which branch into a freely slipping major basement fault at depth. Finite element modelling provides new insights into the controversial and widely debated active tectonic setting of the study area, pointing out the fundamental role played by a structural reactivation process involving inherited, long-lived, mature fault systems at depth. Our FEM, reconciling apparently contrasting geological and geophysical constraints from the study area, points to maximum stress build up and strain accumulation at a depth of 15 ± 5 km. Such a depth range is suggested as the most likely one for the nucleation of large events such as the 1857 Val d'Agri earthquake.

Over the years, the belief that the first mode of vibration governs the seismic response of shear-type frame structures has been widely accepted and proved to be effective for preliminary structural design. Indeed, most of the actual seismic design procedures are based on drift profiles which are typically an approximation of the shape of the fundamental mode of vibration. In this paper, an analytical study on the dynamic properties of multi-storey shear-type frames is carried out with the purpose of precisely identifying the contribution of the modes of vibration to the seismic response of such structures, both in terms of maximum inter-storey displacement profiles (which govern the beams and columns maximum actions) and maximum inter-storey velocity profiles (which govern the viscous dampers maximum forces, of fundamental importance for building structures equipped with additional viscous dampers). A new parameter, referred to as Seismic Modal Contribution Factor, which represents the contribution of the generic mode to the seismic response of the structure, is introduced. With respect to the well-known Modal Contribution Factor, grounded on the concept of modal static response, the Seismic Modal Contribution Factor explicitly takes into account also the dynamic nature of the response due to earthquake excitation. The Seismic Modal Contribution Factor could be a meaningful parameter to be implemented in a professional structural design software and used in conjunction with the common modal participating mass ratios to identify the number of modes to be included in the analyses.

The Eastern Cordillera of Colombia rose to maximum elevations of >5 km during the Cenozoic by inversion of a Mesozoic rift basin. Previous studies proposed that the exhumation of the Eastern Cordillera increased from ~6 Ma to the present due to the interplay between tectonic shortening and climate. In this study, we integrate new field observations, structural data, low-temperature thermochronology, thermobarometry, and vitrinite reflectance along a section through the Tablazo, Cocuy, and Llanos regions to estimate the amount of shortening and the exhumation history. Our results indicate that shortening started as early as the latest Maastrichtian-Paleocene in the Tablazo and Cocuy regions. Exhumation migrated eastward, starting in the Paleocene in the west and continuing in the Miocene in the east. The amount and rate of exhumation peaked in the Cocuy region with values of

We study the thermomechanical evolution of a stationary lava body cooling from the upper surface. The body progressively solidifies and develops a boundary layer made of an upper solid thermoelastic crust and a lower viscoelastic layer. The presented approach to the problem is one-dimensional and quasi-static. Under the assumption that the crust is a brittle solid with a given strength, we study the generation of tensile fractures. We consider a single fracture event and we assume that a fracture is produced when the thermal stress exceeds the strength over a certain depth interval, inferred from the observation of striae. The fracture propagates abruptly across the brittle layer, driven by the difference between the thermal stress and the tensile strength. The model can contribute to explain the formation of striae, surfaces caused by crack advance, in particular, of the first fracture. We find that the Young modulus controls the crack growth. The model provides an explanation for the presence of the rough bands on striae. For its higher values the fracture extends partially into the viscoelastic region; while for lower values the crack is confined to the brittle layer. The maximum opening width is typically in the order of a few thousandths of the vertical length of the crack. The resumption of crack process must await further cooling of the lava body.

The European Alps show major along-strike variations in the partitioning of crustal shortening between the upper and the lower plates, either explained by lateral variation in temperature and rheology, or by along-strike changes in crustal-scale lithology. Here, we present the first 3D thermomechanical geodynamic model to test the potential impact of along-strike lithological (and associated rheological) changes and lateral temperature variations on crustal shortening during continental collision in the Western and Central Alps. Numerical modeling results suggest a minor impact of lateral thermal variations and lithological changes in the lower plate for the Alpine orogeny. In contrast, lateral lithological changes in the upper plate structure, defined by the presence of the buried Ivrea Body in the Western Alps, should play a major role in determining the present-day architecture of the Alpine belt. Our models suggest that the Ivrea Body acts as a strong, rigid buttress constraining crustal shortening to the lower plate, promoting the formation of the Jura fold-and-thrust belt, and supporting the higher surface elevation observed in the Western Alps. On the other hand, in the Central Alps where the Ivrea Body is absent, much broader distribution of crustal shortening results in lower surface elevation. Numerical modeling results also suggest that crustal shortening in the Central Alps is focused within the upper plate during the final stages of continental collision, favored by the indentation of thickened upper crust of the lower plate, and by detachment of the Adriatic upper crust from the underlying lithosphere.

The Apennines is a well-studied orogeny formed by the accretion of continental slivers during the subduction of the Adriatic plate, but its deep structure is still a topic of controversy. Here we illuminated the deep structure of the Northern Apennines belt by combining results from the analysis of active seismic (CROP03) and receiver function data. The result from combining these two approaches provides a new robust view of the structure of the deep crust/upper mantle, from the back-arc region to the Adriatic subduction zone. Our analysis confirms the shallow Moho depth beneath the back-arc region and defines the top of the downgoing plate, showing that the two plates separate at depth about 40km closer to the trench than reported in previous reconstructions. This spatial relationship has profound implications for the geometry of the shallow subduction zone and of the mantle wedge, by the amount of crustal material consumed at trench.

Microbes play a fundamental role in the precipitation of silicate biominerals, thereby affecting the Si geochemical cycle. The fne mechanisms ruling biomineralization are not yet fully understood, and their microscopic structures can offer deep insight into their processes of formation, reactivity and stability. In this study, a Zn silicate biomineral, extracellularly produced by cyanobacterium Leptolingbya frigida, was investigated combining nuclear magnetic resonance (NMR), Zn K-edge X-ray absorption spectroscopy (XAS) and other complementary techniques. 29Si magic angle spinning and 29Si/1H cross polarization magic angle spinning analysis, Fourier transform infrared spectroscopy (FTIR) and XAS analysis revealed a poorly crystalline phase closely resembling hemimorphite [Zn4Si2O7(OH)2•H2O]. Zn K-edge extended X-ray absorption fne structure (EXAFS) provided further structural details, revealing that the Zn-O-Si interatomic distances were 7-8% shorter than the abiotic mineral. 13C NMR spectra analysis was conducted to investigate the composition of the Zn silicate biomineral organic matrix, and results revealed that C atoms occurred in several functional groups such as carbonyl carbons, C rings, O-aliphatic chains, N-aliphatic chains, and aliphatic chains. Under slightly alkaline conditions, bacterial cell walls exhibited fundamental control on the biomineralization process by binding Zn ions and forming Zn-O-Si bonds. In this way, L. frigida cell walls served as a reactive surface for the precipitation of this Zn sorosilicate, hindering the condensation of silicon dimers. Moreover, we found a 29Si NMR band at 85 ppm that could be attributed to a (C3H6O3)2Si complex. This complex could play a role in the control of silicon polymerization, with implications for Si biomineralization processes.

The 2017 Mw 7.3 Iran–Iraq earthquake occurred in a region where the pattern of major plate convergence is well constrained, but limited information is available on the seismogenic structures. Geological observations, interpretation of seismic reflection profiles, and well data are used in this paper to build a regional, balanced cross section that provides a comprehensive picture of the geometry and dimensional parameters of active faults in the hypocentral area. Our results indicate (i) the coexistence of thin- and thick-skinned thrusting, (ii) the reactivation of inherited structures, and (iii) the occurrence of weak units promoting heterogeneous deformation within the palaeo-Cenozoic sedimentary cover and partial decoupling from the underlying basement. According to our study, the main shock of the November 2017 seismic sequence is located within the basement, along the low-angle Mountain Front Fault. Aftershocks unzipped the up-dip portion of the same fault. This merges with a detachment level located at the base of the Paleozoic succession, to form a crustal-scale fault-bend anticline. Size and geometry of the Mountain Front Fault are consistent with a down-dip rupture width of 30 km, which is required for an Mw 7.3 earthquake.

In the original publication of the article, under Section 2.1, the sentence in the eighth paragraph that reads as “It has been demonstrated…very disturbed periods.” should read as “It has been demonstrated (Miro Amarante et al. 2007) that the ionosphere shows small scale variations that global or regional empirical models describe hardly, especially for very disturbed periods.” In the eleventh paragraph under Section 3.2, the sentence that reads as “In formulas (3)–(4),…respectively.” should read as “In formulas (3)–(4), Ulow2j;kand Uhigh2j;k(similarly for 2j - 1 terms) are the coefficient values calculated for IG12 or R12 equal to 0 (low solar activity) and equal to 100 (high solar activity), respectively.” Finally, under Section 3.3.2, the sentence in the third paragraph that reads as “In this case, the… or small values of h.” should read as “In this case, the experimental variogram has a parabolic behavior near the origin, that is, it is proportional to h2for small values of h.”.

Although earthquakes generally form clusters in both space and time, only mainshocks, usually the largest magnitude events within clusters, are considered by probabilistic seismic hazard analysis (PSHA; Cornell, 1968). On the other hand, aftershock PSHA (APSHA), based on the modified Omori law, allows the quantification of the aftershock threat (Yeo and Cornell, 2009). Classical PSHA often describes event occurrence via a homogeneous Poisson process, whereas APSHA describes occurrence of aftershocks via cluster-specific nonhomogeneous Poisson processes, the rate of which is a function of the mainshock magnitude. It is easy to recognize that clusters, each of which is made of the mainshock and the following aftershocks, occur at the same rate of mainshocks. This recently allowed the generalization of the hazard integral to account for aftershocks in PSHA (i.e., Iervolino et al., 2014), which resulted in the formulation of the so-called sequence-based PSHA (SPSHA). In the present study, SPSHA is applied to Italy countrywide, using the same source model (Stucchi et al., 2011) that lies at the base of the official PSHA used for structural design, to quantitatively assess the increase in seismic design actions for structures when accounting for the aftershocks.

This Special Issue hosts a selection of papers presented at the Third International Conference on Metrology for Archaeology and Cultural Heritage (Lecce, Italy, 23-25 October 2017), related to the new perspectives on geophysics for archaeology. In recent years, archaeological prospecting has seen major advances through a variety of remote sensing and computing technologies. Geophysical instrumentation continues to improve in sensitivity and acquisition speed, and new multi-sensor arrays, e.g., drawn by carts over land, now permit vast areas to be rapidly covered. On the other hand, the availability of highresolution remote sensing techniques provides multiscale and multi-temporal approaches to the study of ancient settlements and landscapes, and proves fundamental in the reconstruction of their development over centuries. Nowadays, research in landscape archaeology needs the integration of different high-resolution remote sensing techniques such as satellite (optical and radar data), aerial (photographic, infrared and lidar data from aircraft and unmanned aerial vehicles), but also land acquisitions (integration of different geophysical techniques, field walking, and differential GPS topographical surveys). All these investigations are based on a geoarchaeological approach, with several aims ranging from historical reconstruction to preventive archaeology and from the preservation of archaeological and monumental heritage to noninvasive diagnosis through micro-geophysical techniques.

Natural remanent magnetization and self-demagnetization on high-susceptibility bodies are two important factors affecting magnetic data inversion. We propose a framework for the inversion and interpretation of magnetic anomalies, in which significant remanent magnetization and self-demagnetization are present simultaneously. The framework is based on the assumptions that the external applied field and internal self-demagnetization field are uniform and the deflection of self-demagnetization in the total magnetization direction is negligible. First, the magnetization vector distributions are obtained from magnetic data by estimating the magnetization direction, then inverting for the magnetization intensity distribution, using the inferred magnetization direction as a constraint. Based on a priori information about the Koenigsberger ratio derived from petrophysical measurements, the direction and intensity of the remanent magnetization are obtained. The self-demagnetization factor is then computed using the finite volume method. Finally, the true-susceptibility distribution is achieved by correcting for the self-demagnetization effect. The method is first applied to synthetic magnetic data produced by a prism-shaped source model that has significant remanent magnetization and high susceptibility. In a case study of the Daye iron-ore deposit, Hubei province, China, the true susceptibility and remanent magnetization are reconstructed. The remanence direction information reveals that local geological activities such as synclines and faults lead to changes in the remanence directions at different local deposits.

The natural formation of rock slopes can result in the initiation of complex failure processes, including the deformation and fracturing of intact rock and the displacement of large individual blocks. The presence of non-persistent discontinuities may also promote stress concentrations that would result in stress induced fractures leading to the formation of continuous failure surfaces, thus creating the condition of kinematic freedom for previously finite and tapered blocks that may result in slope failure. This type of mechanism cannot be analyzed using a limit equilibrium approach and requires models capable of considering rock mass plastic yield, simulation of shear/tensile fracturing and ability to consider kinematic processes. In this context, hybrid finite-discrete element methods (FDEM) provide an attractive option for modelling rock slope failure mechanisms. In this paper a FDEM approach was applied to study a rock face located near the Passo della Morte (Carnian Alps, Italy), on the left flank of Tagliamento River valley. The rock face is approximately 130 m wide and 250 m high (650 to 900 m a.s.l.). The study includes 2D numerical analysis of a critical section for which the potential failure mechanism would be controlled by a series of faults bounding a densely stratified rock mass crossed by very persistent sub-vertical joints. Overall, the rock mass is divided into a large number of blocks of variable size, that are isolated by a dense network of discontinuities and random fractures. The importance of damage in the rock slope mass has been recognized by several researchers using remote sensing methods.

Ground‐motion prediction equations (GMPEs) are calibrated to predict the intensity of ground shaking at any given location, based on earthquake magnitude, source‐to‐site distance, local soil amplifications, and other parameters. GMPEs are generally assumed to be independent of time; however, evidence is increasing that large earthquakes modify the shallow soil conditions and those of the fault zone for months or years. These changes may affect the intensity of shaking and result in time‐dependent effects that can potentially be resolved by analyzing between‐event residuals (residuals between observed and predicted ground motion for individual earthquakes averaged over all stations). Here, we analyze a data set of about 65,000 recordings for about 1400 earthquakes in the moment magnitude range 2.5–6.5 that occurred in central Italy from 2008 to 2017 to capture the temporal variability of the ground shaking at high frequency. We first compute between‐event residuals for each earthquake in the Fourier domain with respect to a GMPE developed ad hoc for the analyzed data set. The between‐events show large changes after the occurrence of mainshocks such as the 2009 Mw 6.3 L'Aquila, the 2016 Mw 6.2 Amatrice, and Mw 6.5 Norcia earthquakes. Within the time span of a few months after the mainshocks, the between‐event contribution to the ground shaking varies by a factor 7. In particular, we find a large drop in the between‐events in the aftermath of the L'Aquila earthquake, followed by a slow positive trend that leads to a recovery interrupted by a new drop at the beginning of 2014. We also quantify the frequency‐dependent correlation between the Brune stress drop Δσ and the between‐events. We find that the temporal changes of Δσ resemble those of the between‐event residuals; in particular, during the period when the between‐events show the positive trend, the average logarithm of Δσ increases with an annual rate of 0.19 (i.e., the amplification factor for Δσ is 1.56 per year). Breakpoint analysis located a change in the linear trend coefficients of Δσ versus time in February 2014, although no large earthquakes occurred at that time. Finally, the temporal variability of Δσ mirrors the relative seismic‐velocity variations observed in previous studies for the same area and period, suggesting that both crack healing along the main fault system and healing of microcracks distributed at shallow depths throughout the surrounding region might be necessary to explain the wider observations of postearthquake recovery.

Kamarina, located in southern Sicily (Italy), was an important Greek colony since its foundation in the sixth century BC. Archaeological excavations, carried out since the twentieth century, uncovered only limited portions of the site so far. Despite the importance of the Greek colony, the presence of remarkable buildings that archaeologists expected to bring to light has not found fully correspondence in the archaeological excavations. Consequently, the integrated geophysical prospection carried out in the study area is aimed to support and address the future archaeological investigations. After the photographic and thermographic survey obtained by an unmanned aerial vehicle, we performed a systematic survey through ground magnetic and GPR methods over an area of 6200 m2. The acquisition procedures have been optimized in order to get the best results combining high resolution and elevated speed of acquisition. The results derived from the three geophysical techniques have been conveniently combined by means of a cluster analysis, allowing us to clearly identify a series of buried archaeological features. Because of their geometrical characteristics, often in good agreement with the spatial arrangement of the archaeological remains at the surface, these buried archaeological features can be interpreted as roads, walls, or buildings foundations in which the various construction phases of the city can be clearly recognized. The integrated approach has proven to be essential for a robust interpretation of the archaeogeophysical investigation.

"Super-deep" diamonds are thought to have a sub-lithospheric origin (i.e., below ~300 km depth) because some of the mineral phases entrapped within them as inclusions are considered to be the products of retrograde transformation from lower-mantle or transition-zone precursors. CaSiO3-walstromite, the most abundant Ca-bearing mineral inclusion found in super-deep diamonds, is believed to derive from CaSiO3-perovskite, which is stable only below ~600 km depth, although its real depth of origin is controversial. The remnant pressure (Pinc) retained by an inclusion, combined with the thermoelastic parameters of the mineral inclusion and the diamond host, allows calculation of the entrapment pressure of the diamond-inclusion pair. Raman spectroscopy, together with X ray diffraction, is the most commonly used method for measuring the Pinc without damaging the diamond host. In the present study we provide, for the first time, a calibration curve to determine the Pinc of a CaSiO3-walstromite inclusion by means of Raman spectroscopy without breaking the diamond. To do so, we performed high-pressure micro-Raman investigations on a CaSiO3-walstromite crystal under hydrostatic stress conditions within a diamond-Anvil cell. We additionally calculated the Raman spectrum of CaSiO3-walstromite by ab initio methods both under hydrostatic and non-hydrostatic stress conditions to avoid misinterpretation of the results caused by the possible presence of deviatoric stresses causing anomalous shift of CaSiO3-walstromite Raman peaks. Last, we applied single-inclusion elastic barometry to estimate the minimum entrapment pressure of a CaSiO3-walstromite inclusion trapped in a natural diamond, which is ~9 GPa (~260 km) at 1800 K. These results suggest that the diamond investigated is certainly sub-lithospheric and endorse the hypothesis that the presence of CaSiO3-walstromite is a strong indication of super-deep origin.

Pacaya volcano in Guatemala is one of the most active volcanoes of the Central American Volcanic Arc (CAVA). However, its magmatic gas signature and volatile output have received little attention to date. Here, we present novel volcanic gas information from in-situ (Multi-GAS) and remote (UV camera) plume observations in January 2016. We find in-plume H2O/SO2 and CO2/SO2 ratios of 2-20 and 0.6-10.5, and an end-member magmatic gas signature of 80.5 mol. % H2O, 10.4 mol. % CO2, and 9.0 mol. % SO2. The SO2 flux is evaluated at 885 ± 550 tons/d. This, combined with co-acquired volcanic plume composition, leads to H2O and CO2 fluxes of 2,230 ± 1,390 and 700 ± 440, and a total volatile flux of ∼3,800 tons/d. We use these results in tandem with previous SO2 flux budgets for Fuego and Santiaguito to estimate the total volcanic CO2 flux from Guatemala at ∼1,160 ± 600 tons/day. This calculation is based upon CO2/total S (St) ratios for Fuego (1.5 ± 0.75) and Santiaguito (1.4 ± 0.75) inferred from a gas (CO2/St ratio) versus trace-element (Ba/La ratio) CAVA relationship. The H2O-poor and low CO2/St ratio (∼1.0-1.5) signature of Pacaya gas suggests dominant mantle-wedge derivation of the emitted volatiles. This is consistent with 3He/4He ratios in olivine hosted fluid inclusions (FIs), which range between 8.4 and 9.0 Ra (being Ra the atmospheric 3He/4He ratio) at the upper limit of MORB range (8 ± 1 Ra). These values are the highest ever measured in CAVA and among the highest ever recorded in arc volcanoes worldwide, indicating negligible 4He contributions from the crust/slab.

The long-term carbon budget has major implications for Earth’s climate and biosphere, but the balance between carbon sequestration during subduction, and outgassing by volcanism is still poorly known. Although carbon-rich fluid inclusions and minerals are described in exhumed mantle rocks and xenoliths, compelling geophysical evidence of large-scale carbon storage in the upper mantle is still lacking. Here, we use a geophysical surface-wave seismic tomography model of the mantle wedge above the subducted European slab to document a prominent shear-wave low-velocity anomaly at depths greater than 180 km. We propose that this anomaly is generated by extraction of carbonate-rich melts from the asthenosphere, favoured by the breakdown of slab carbonates and hydrous minerals after cold subduction. The resulting transient network of carbon-rich melts is frozen in the mantle wedge without producing volcanism. Our results provide the first in-situ observational evidence of ongoing carbon sequestration in the upper mantle at a plate-tectonic scale. We infer that carbon sequestered during cold subduction may partly counterbalance carbon outgassed from ridges and oceanic islands. However, subducted carbon may be rapidly released during continental rifting, with global effects on long-term climate trends and the habitability of planet Earth.

Andic soils have unique morphological, physical, and chemical properties that induce both considerable soil fertility and great vulnerability to land degradation. Moreover, they are the most striking mineral soils in terms of large organic C storage and long C residence time. This is especially related to the presence of poorly crystalline clay minerals and metal–humus complexes. Recognition of andic soils is then very important.Here we attempt to show, through a combined analysis of 35 sampling points chosen in accordance to specific physical and vegetation rules, that some andic soils have an utmost ecological importance.More specifically, in Italian non-volcanic mountain ecosystems ( > 600 m a.s.l.) combining low slope (Our results are hoped to be a starting point for better understanding of the ecological importance of andic soils and also possibly to better consider pedological information in C balance calculations.

An integrated geophysical survey has been conducted at the Tarragona’s Cathedral (Catalonia, NE Spain) with the aim to confirm the potential occurrence of archaeological remains of the Roman Temple dedicated to the Emperor Augustus. Many hypotheses have been proposed about its possible location, the last ones regarding the inner part of the Cathedral, which is one of the most renowned temples of Spain (twelfth century) evolving from Romanesque to Gothic styles. A geophysical project including electrical resistivity tomography (ERT) and ground probing radar (GPR) was planned over 1 year considering the administrative and logistic difficulties of such a project inside a cathedral of religious veneration. Finally, both ERT and GPR have been conducted during a week of intensive overnight surveys that provided detailed information on subsurface existing structures. The ERT method has been applied using different techniques and arrays, ranging from standard Wenner–Schlumberger 2D sections to full 3D electrical imaging with the advanced Maximum Yield Grid array. Electrical resistivity data were recorded extensively, making available many thousands of apparent resistivity data to obtain a complete 3D image after a full inversion. In conclusion, some significant buried structures have been revealed providing conclusive information for archaeologists. GPR results provided additional information about shallowest structures. The geophysical results were clear enough to persuade religious authorities and archaeologists to conduct selected excavations in the most promising areas that confirmed the interpretation of geophysical data. In conclusion, the significant buried structures revealed by geophysical methods under the cathedral were confirmed by archaeological digging as the basement of the impressive Roman Temple that headed the Provincial Forum of Tarraco, seat of the Concilium of Hispania Citerior Province.

International audience; Slab rollback results in the development of low-angle normal faults (detachments) and metamorphic core complexes (MCCs) in back-arc domains. Although the mechanical consequences of slab dynamics on lithospheric and crustal behaviors have already been studied, thermal effects have not been investigated yet. This study shows that slab rollback produces lithospheric-scale thermal perturbations intrinsically associated with emplacement of amagmatic high-enthalpy geothermal systems. Using a multi-scale numerical modeling approach, with lithospheric-scale 3-D thermo-mechanical models of subduction, and 2-D models of fluid flow at the scale of detachments, we demonstrate that subduction-induced extensional tectonics controls the genesis and distribution of crustal-scale thermal domes from the base of the crust, and the location of high-energy geothermal systems. We find that when slab tearing occurs, Moho temperatures can temporarily increase by up to 250 °C due to significant shear heating in the flowing upper mantle. Associated thermal anomalies (with characteristic width and spacing of tens and hundreds of km, for crustal and lithospheric scales, respectively) then migrate systematically toward the retreating trench. These thermal domes weaken the crust, localize deformation and enhance the development of crustal-scale detachments. These thermo-mechanical instabilities mimic genesis of high-temperature MCCs with migmatitic cores in the back-arc domain, such as those of the Menderes (western Anatolia, Turkey) and Larderello (southern Tuscany) provinces in the Mediterranean realm, and those in the Basin and Range (western United States), where detachments control the bulk of the heat transport. At the scale of MCCs, the bulk fluid flow pattern is controlled by topography-driven flow while buoyancy-driven flow dominates within the permeable detachments, focusing reservoir location of high-energy geothermal systems at shallow depth beneath the detachments.

Dataset from the paper "Relation Between Near-Fault Ground Motion Impulsive Signals and Source Parameters" - Scala et al. 2018 (in review)

This article presents the main results of the validation phase of the PRENOLIN project. PRENOLIN is an international benchmark on 1D nonlinear (NL) site‐response analysis. This project involved 19 teams with 23 different codes tested. It was divided into two phases; with the first phase verifying the numerical solution of these codes on idealized soil profiles using simple signals and real seismic records. The second phase described in this article referred to code validation for the analysis of real instrumented sites. This validation phase was performed on two sites (KSRH10 and Sendai) of the Japanese strong‐motion networks KiK‐net and Port and Airport Research Institute (PARI), respectively, with a pair of accelerometers at surface and depth. Extensive additional site characterizations were performed at both sites involving in situ and laboratory measurements of the soil properties. At each site, sets of input motions were selected to represent different peak ground acceleration (PGA) and frequency content. It was found that the code‐to‐code variability given by the standard deviation of the computed surface‐response spectra is around 0.1 (in log10 scale) regardless of the site and input motions. This indicates a quite large influence of the numerical methods on site‐effect assessment and more generally on seismic hazard. Besides, it was observed that site‐specific measurements are of primary importance for defining the input data in site‐response analysis. The NL parameters obtained from the laboratory measurements should be compared with curves coming from the literature. Finally, the lessons learned from this exercise are synthesized, resulting also in a few recommendations for future benchmarking studies, and the use of 1D NL, total stress site‐response analysis.

In recent years, the use of radiographic inspection with cosmic-ray muons has spread into multiple research and industrial fields. This technique is based on the high-penetration power of cosmogenic muons. Specifically, it allows the resolution of internal density structures of large scale, geological objects through precise measurement of the muon absorption rate. So far in many previous works, this muon absorption rate has been considered to depend solely on the density of traversed material (under the assumption of a standard rock) but the variation in chemical composition has not been taken seriously into account. However, from our experience with muon tomography in Alpine environments we find that this assumption causes a substantial bias on the muon flux calculation, particularly where the target consists of high {Z2/A} (like basalts) or low {Z2/A} (e.g. dolomites) rocks and where the material thickness exceeds 300 m. In this paper we derive an energy loss equation for different minerals and we additionally derive a related equation for mineral assemblages that can be used for any rock type on which mineralogical data is available. Thus, for muon tomography experiments in which high/low {Z2/A} rock thicknesses can be expected, it is advisable to plan an accompanying geological field campaign to determine a realistic rock model.

Understanding the thermal behavior of pyroclastic density currents (PDCs) is crucial for forecasting impact scenarios for exposed populations as it affects their lethality and destructiveness. Here we report the emplacement temperatures of PDC deposits produced during the paroxysmal explosive eruption of Merapi (Central Java) on 5 November 2010 based on the reflectance of entombed charcoal fragments. This event was anomalously explosive for Merapi, and destroyed the summit dome that had been rapidly growing, with partial collapses and associated PDCs, since October 26. Results show mean reflectance values mainly between 0.17 and 0.41. These new data provide a minimum temperature of the flow of 240–320 °C, consistent with previous estimations determined from independent field, engineering, and medical observations published in the literature for this eruption. A few charcoal fragments recorded higher values, suggestive of temperatures up to 450 °C, and we suggest that this is due to the thermal disequilibrium of the deposits, with larger block-size clasts being much hotter than the surrounding ash matrix. Charring temperatures show no major differences between proximal and distal PDC deposits and are significantly lower than those that may be associated with a fast growing dome dominated by dense and vitric non-vesicular rocks. We therefore infer that the decrease in temperature from that at fragmentation (>900 °C) occurred in the very initial part of the current

A multi-2D imaging of the Solfatara Crater inside the Campi Flegrei Caldera, was obtained by the joint interpretation of geophysical evidences and the new active seismic dataset acquired during the RICEN experiment (EU project MEDSUV) in 2014. We used a total of 17,894 first P-wave arrival times manually picked on pre-processed waveforms, recorded along two 1D profiles criss-crossing the inner Solfatara crater, and performed a tomographic inversion based on a multi-scale strategy and a Bayesian estimation of velocity parameters. The resulting tomographic images provide evidence for a low velocity (500–1500 m/s) water saturated deeper layer at West near the outcropping evidence of the Fangaia, contrasted by a high velocity (2000–3200 m/s) layer correlated with a consolidated tephra deposit. The transition velocity range (1500–2000 m/s) layer suggests a possible presence of a gas-rich, accumulation volume. Thanks to the mutual P-wave velocity model, we infer a detailed image for the gas migration path to the Earth surface. The gasses coming from the deep hydrothermal plume accumulate in the central and most depressed area of the Solfatara being trapped by the meteoric water saturated layer. Therefore, the gasses are transmitted through the buried fault toward the east part of the crater, where the ring faults facilitate the release as confirmed by the fumaroles. Starting from the eastern surface evidence of the gas releasing in the Bocca Grande and Bocca Nuova fumaroles, and the presence of the central deeper plume we suggest a fault situated in the central part of the crater which seems to represent the main buried conduit among them plays a key role.

In the last two decades, several studies addressed the revaluation and homogenization of the Italian instrumental seismic catalog, but all of them refer to the time interval from 1981, that is, the starting year of the Catalogo Strumentale dei Terremoti Italiani (CSTI). At the time, the CSTI was conceived as the continuation of the catalog of the Progetto Finalizzato Geodinamica (PFG) but, over time, the PFG catalog was almost totally forgotten, and presently it is even difficult to obtain because it is not provided by any website. In this work, we integrate a genuine copy of PFG, with additional locations from the bulletins of the Istituto Nazionale di Geofisica (ING, now known as INGV) and of the International Seismological Centre (ISC) and with local magnitudes from two couples of Wood–Anderson (WA) seismometers operational in Italy in the 1970s and 1980s, mostly derived from a careful scrutiny of paper bulletins of the Osservatorio Geofisico Sperimentale (OGS) and of the ING. We restrict our analysis to the time interval from 1960 to 1980 because, based on various evidence, we can infer that within such period most instrumental magnitudes reported by the PFG catalog are reasonably coherent with the Richter’s definition. Magnitudes provided by WA stations and other data sources are calibrated with respect to Mwby general orthogonal regressions. The final catalog from 1960 to 1980 contains 8536 earthquakes, of which we compute a true or proxy Mwmagnitude with related uncertainty for 6407. The analysis of the frequency–magnitude distribution indicates completeness for about Mw≥4:0. This work extends the time coverage of the Italian instrumental catalog to about 55 yrs before the present, allowing the statistical study of some important seismic periods that occurred, for example, in 1962 (Irpinia), 1968 (Belice Valley), 1976 (Friuli), 1979 (Umbria), and 1980 (Irpinia).

The realization of the international terrestrial reference frame (ITRF) is currently based on the data provided by four space geodetic techniques. The accuracy of the different technique-dependent materializations of the frame physical parameters (origin and scale) varies according to the nature of the relevant observables and to the impact of technique-specific errors. A reliable computation of the ITRF requires combining the different inputs, so that the strengths of each technique can compensate for the weaknesses of the others. This combination, however, can only be performed providing some additional information which allows tying together the independent technique networks. At present, the links used for that purpose are topometric surveys (local/terrestrial ties) available at ITRF sites hosting instruments of different techniques. In principle, a possible alternative could be offered by spacecrafts accommodating the positioning payloads of multiple geodetic techniques realizing their co-location in orbit (space ties). In this paper, the GNSS–SLR space ties on-board GPS and GLONASS satellites are thoroughly examined in the framework of global reference frame computations. The investigation focuses on the quality of the realized physical frame parameters. According to the achieved results, the space ties on-board GNSS satellites cannot, at present, substitute terrestrial ties in the computation of the ITRF. The study is completed by a series of synthetic simulations investigating the impact that substantial improvements in the volume and quality of SLR observations to GNSS satellites would have on the precision of the GNSS frame parameters.

International audience; Volcanogenic halogens, in particular bromine, potentially play an important role in the ozone depletion of the atmosphere. Understanding bromine behaviour in magmas is therefore crucial to properly evaluate the contribution of volcanic eruptions to atmospheric chemistry and their environmental impact. To date, bromine partitioning between silicate melts and the gas phase is very poorly constrained, with the only relevant experimental studies limited to investigation of synthetic melt with silicic compositions. In this study, fluid/melt partitioning experiments were performed using natural silicate glasses with mafic, intermediate and silicic compositions. For each composition, experiments were run with various Br contents in the initial fluid (H2O–NaBr), at T–P conditions representative of shallow magmatic reservoirs in volcanic arc contexts (100–200 MPa, 900–1200 °C). The resulting fluid/melt partition coefficients (DBrf/m) are: 5.0 ± 0.3 at 1200 °C–100 MPa for the basalt, 9.1 ± 0.6 at 1060 °C–200 MPa for the andesite and 20.2 ± 1.2 at 900 °C–200 MPa for the rhyodacite. Our experiments show that DBrf/m increases with increasing SiO2 content of the melt (as for chlorine) and suggest that it is also sensitive to melt temperature (increase of DBrf/m with decreasing temperature). We develop a simple model to predict the S–Cl–Br degassing behaviour in mafic systems, which accounts for the variability of S–Cl–Br compositions of volcanic gases from Etna and other mafic systems, and shows that coexisting magmatic gas and melt evolve from S-rich to Cl–Br enriched (relative to S) upon increasing degree of degassing. We also report first Br contents for melt inclusions from Etna, Stromboli, Merapi and Santorini eruptions and calculate the mass of bromine available in the magma reservoir prior to the eruptions under consideration. The discrepancy that we highlight between the mass of Br in the co-existing melt and fluid prior to the Merapi 2010 eruption (433 and 73 tons, respectively) and the lack of observed BrO (from space) hints at the need to investigate further Br speciation in ‘ash-rich’ volcanic plumes. Overall, our results suggest that the Br yield into the atmosphere of cold and silicic magmas will be much larger than that from hotter and more mafic magmas.

A new and unusual eruptive event occurred on 29 February 2016 within the summit crater of the Mount Nyiragongo volcano. Based on field campaigns performed between July 2015 and September 2017, and building on a previously published buoyancy‐driven bidirectional magma flow model explaining the progression of Mount Nyiragongo lava lake level, we provide the first quantitative estimations of volumes of erupted lava outpouring from the new spatter cone. Besides matching field data of the lava lake level covering the period December 2002 to September 2017, numerical solutions of the model reveal that the most probable dike path is one originating from the shallow magma reservoir, and not from the lava lake or branching from its feeding conduit. According to these simulations, the reservoir and erupted lava volumes are respectively estimated to the order of 10 km3 and 20 M m3. Magma overpressure at the level of the shallow reservoir is estimated in the range 12 to 16 MPa, high enough to potentially initiate new erupting events. A new and unusual eruptive event occurred on 29 February 2016 within the summit crater of the Mount Nyiragongo volcano. Based on field campaigns performed between July 2015 and September 2017, and building on a previously published buoyancy-driven bidirectional magma flow model explaining the progression of Mount Nyiragongo lava lake level, we provide the first quantitative estimations of volumes of erupted lava outpouring from the new spatter cone. Besides matching field data of the lava lake level covering the period December 2002 to September 2017, numerical solutions of the model reveal that the most probable dike path is one originating from the shallow magma reservoir, and not from the lava lake or branching from its feeding conduit. According to these simulations, the reservoir and erupted lava volumes are respectively estimated to the order of 10 km3 and 20 M m3. Magma overpressure at the level of the shallow reservoir is estimated in the range 12 to 16 MPa, high enough to potentially initiate new erupting events.

Highlights • Active tectonics in Sicilian Channel • Analogue models of transcurrent structures • 3D fault pattern reconstruction Seismo-stratigraphic and structural analysis of a large number of multichannel seismic reflection profiles acquired in the northern part of the Sicilian Channel allowed a 3-D reconstruction of a regional NS-trending transfer zone which displays a transcurrent tectonic regime, and that is of broad relevance for its seismotectonic and geodynamic implications. It is constituted of two major transcurrent faults delimiting a 30-km-wide, mostly undeformed basin. The western fault (Capo Granitola) does not show clear evidence of present-day tectonic activity, and toward the south it is connected with the volcanic area of the Graham Bank. The eastern fault (Sciacca) is structurally more complex, showing active deformation at the sea-floor, particularly evident along the Nerita Bank. The Sciacca Fault is constituted of a master and splay faults compatible with a right-lateral kinematics. Sciacca Fault is superimposed on an inherited weakness zone (a Mesozoic carbonate ramp), which borders to the east a 2.5-km-thick Plio-Quaternary basin, and that was reactivated during the Pliocene. A set of scaled claybox analogue models was carried out in order to better understand the tectonic processes that led to the structural setting displayed by seismic data. Tectonic structures and uplift/subsidence patterns generated by the models are compatible with the 3-D model obtained from seismic reflection profiles. The best fit between the tectonic setting deriving from the interpretation of seismic profiles and the analogue models was obtained considering a right-lateral movement for the Sciacca Fault. Nevertheless, the stress field in the study area derived from GPS measurements does not support the present-day modelled right-lateral kinematics along the Sciacca Fault. Moreover, seismic events along this fault show focal mechanisms with a left-lateral component. We ascribe the slip change along the Sciacca Fault, from a right-lateral transcurrent regime to the present-day left-lateral kinematics to a change of principal horizontal stress direction starting from Late Pliocene.

Nakamura (Q Rep Railway Tech Res Inst 30:25–33, 1989) popularized the application of the horizontal-to-vertical spectral ratio (HVSR) analysis of microtremor (seismic noise or ambient vibration) recordings to estimate the predominant frequency and amplification factor of earthquake shaking. During the following quarter century, popularity in the microtremor HVSR (MHVSR) method grew; studies have verified the stability of a site’s MHVSR response over time and validated the MHVSR response with that of earthquake HVSR response. Today, MHVSR analysis is a popular reconnaissance tool used worldwide for seismic microzonation and earthquake site characterization in numerous regions, specifically, in the mapping of site period or fundamental frequency and inverted for shear-wave velocity depth profiles, respectively. However, the ubiquity of MHVSR analysis is predominantly a consequence of its ease in application rather than our full understanding of its theory. We present the state of the art in MHVSR analyses in terms of the development of its theoretical basis, current state of practice, and we comment on its future for applications in earthquake site characterization.