Sismicité du Massif Armoricain : relocalisations et interprétation tectonique

The Armorican Massif is an outcropping segment of the Hercynian belt of Western Europe. Nowadays it constitutes an intraplate deformation domain in a passive continental margin context. This deformation finds its expression in a moderate seismic activity characterized by low magnitude earthquakes. The geographical distribution of the seismicity, provided by French seismological organizations, exhibit mainly a diffuse behaviour, with large uncertainties on hypocentral parameters. A new analysis of all seismograms recorded between 1980 and 2004, combined with a stochastic inversion of the arrival times visually picked in this work, enabled me to relocate approximately 1500 earthquakes. A joint interpretation of the latter results and of geophysical data available to date, leads to a better understanding of the seismotectonic processes which trigger this intraplate seismicity.Le Massif Armoricain est un segment affleurant de la Chaîne Hercynienne d'Europe Occidentale et constitue actuellement un domaine de déformation intraplaque en contexte de marge continentale passive. Cette déformation se manifeste par une activité sismique modérée caractérisée par des séismes de faible magnitude. Les bulletins des organismes nationaux décrivent - dans cette région peu couverte par les réseaux sismologiques - une répartition diffuse de la sismicité, avec une grande imprécision sur la localisation des événements. Une nouvelle analyse de l'ensemble des sismogrammes, combinée à une inversion stochastique des temps d'arrivée pointés visuellement au cours de ce travail, a permis la relocalisation d'environ 1500 séismes survenus entre 1980 et 2004. L'interprétation jointe de ces résultats et des données géophysiques existantes conduit à une meilleure compréhension des processus sismotectoniques à l'origine de cette sismicité intraplaque

High-frequency ambient noise tomography of southeast Australia: New constraints on Tasmania's tectonic past

The island of Tasmania, at the southeast tip of Australia, is an ideal natural laboratory for ambient noise tomography, as the surrounding oceans provide an energetic and relatively even distribution of noise sources. We extract Rayleigh wave dispersion curves from the continuous records of 104 stations with ∼15 km separation. Unlike most passive experiments of this type, which observe very little coherent noise below a 5 s period, we clearly detect energy at periods as short as 1 s, thanks largely to the close proximity of oceanic microseisms on all sides. The main structural elements of the eastern and northern Tasmanian crust are revealed by inverting the dispersion curves (between 1 and 12 s period) for both group and phase velocity maps. Of particular significance is a pronounced band of low velocity, observed across all periods, that underlies the Tamar River Valley and continues south until dissipating in southeast Tasmania. Together with evidence from combined active source and teleseismic tomography and heat flow data, we interpret this region as a diffuse zone of strong deformation associated with the mid-Paleozoic accretion of oceanic crust along the eastern margin of Proterozoic Tasmania, which has important implications for the evolution of the Tasman Orogen of eastern Australia. In the northwest, a narrower low-velocity anomaly is seen in the vicinity of the Arthur Lineament, which may be attributed to local sediments and strong deformation and folding associated with the final phases of the Tyennan Orogeny

Transdimensional inversion of receiver functions and surface wave dispersion

We present a novel method for joint inversion of receiver functions and surface wave dispersion data, using a transdimensional Bayesian formulation. This class of algorithm treats the number of model parameters (e.g. number of layers) as an unknown in the problem. The dimension of the model space is variable and a Markov chain Monte Carlo (McMC) scheme is used to provide a parsimonious solution that fully quantifies the degree of knowledge one has about seismic structure (i.e constraints on the model, resolution, and trade-offs). The level of data noise (i.e. the covariance matrix of data errors) effectively controls the information recoverable from the data and here it naturally determines the complexity of the model (i.e. the number of model parameters). However, it is often difficult to quantify the data noise appropriately, particularly in the case of seismic waveform inversion where data errors are correlated. Here we address the issue of noise estimation using an extended Hierarchical Bayesian formulation, which allows both the variance and covariance of data noise to be treated as unknowns in the inversion. In this way it is possible to let the data infer the appropriate level of data fit. In the context of joint inversions, assessment of uncertainty for different data types becomes crucial in the evaluation of the misfit function. We show that the Hierarchical Bayes procedure is a powerful tool in this situation, because it is able to evaluate the level of information brought by different data types in the misfit, thus removing the arbitrary choice of weighting factors. After illustrating the method with synthetic tests, a real data application is shown where teleseismic receiver functions and ambient noise surface wave dispersion measurements from the WOMBAT array (South-East Australia) are jointly inverted to provide a probabilistic 1D model of shear-wave velocity beneath a given station

Thermal Nature of Mantle Upwellings Below the Ibero‐Western Maghreb Region Inferred From Teleseismic Tomography

Independent models of P wave and S wave velocity anomalies in the mantle derived from seismic tomography help to distinguish thermal signatures from those of partial melt, volatiles, and compositional variations. Here we use seismic data from SW Europe and NW Africa, spanning the region between the Pyrenees and the Canaries, in order to obtain a new S‐SKS relative arrival‐time tomographic model of the upper mantle below Iberia, Western Morocco, and the Canaries. Similar to previous P wave tomographic results, the S wave model provides evidence for (1) subvertical upper‐mantle low‐velocity structures below the Canaries, Atlas Ranges, and Gibraltar Arc, which are interpreted as mantle upwellings fed by a common lower‐mantle source below the Canaries; and (2) two low‐velocity anomalies below the eastern Rif and Betics that we interpret as the result of the interaction between quasi‐toroidal mantle flow induced by the Gibraltar slab and the mantle upwelling behind it. The analysis of teleseismic P wave and S wave arrival‐time residuals and the conversion of the low‐velocity anomalies to temperature variations suggest that the upwellings in the upper mantle below the Canaries, Atlas Ranges, and Gibraltar Arc system may be solely thermal in nature, with temperature excesses in the range ~100–350 °C. Our results also indicate that local partial melting can be present at lithospheric depths, especially below the Atlas Ranges. The locations of thermal mantle upwellings are in good agreement with those of thinned lithosphere, moderate to high heat‐flow measurements, and recent magmatic activity at the surface

PRISM3D: a 3-D reference seismic model for Iberia and adjacent areas

We present PRISM3D, a 3-D reference seismic model of P- and S-wave velocities for Iberia and adjacent areas. PRISM3D results from the combination of the most up-to-date earth models available for the region. It extends horizontally from 15°W to 5°E in longitude, 34°N to 46°N in latitude and vertically from 3.5 km above to 200 km below sea level, and is modelled on a regular grid with 10 and 0.5 km of grid node spacing in the horizontal and vertical directions, respectively. It was designed using models inferred from local and teleseismic body-wave tomography, earthquake and ambient noise surface wave tomography, receiver function analysis and active source experiments. It includes two interfaces, namely the topography/bathymetry and the Mohorovičić (Moho) discontinuity. The Moho was modelled from previously published receiver function analysis and deep seismic sounding results. To that end we used a probabilistic surface reconstruction algorithm that allowed to extract the mean of the Moho depth surface along with its associated standard deviation, which provides a depth uncertainty estimate. The Moho depth model is in good agreement with previously published models, although it presents slightly sharper gardients in orogenic areas such as the Pyrenees or the Betic-Rif system. Crustal and mantle P- and S-wave wave speed grids were built separately on each side of the Moho depth surface by weighted average of existing models, thus allowing to realistically render the speed gradients across that interface. The associated weighted standard deviation was also calculated, which provides an uncertainty estimation on the average wave speed values at any point of the grid. At shallow depths (<10 km), low P and S wave speeds and high VP/VS are observed in offshore basins, while the Iberian Massif, which covers a large part of western Iberia, appears characterized by a rather flat Moho, higher than average VP and VS and low VP/VS. Conversely, the Betic-Rif system seems to be associated with low VP and VS, combined with high VP/VS in comparison to the rest of the study area. The most prominent feature of the mantle is the well known high wave speed anomaly related to the Alboran slab imaged in various mantle tomography studies. The consistency of PRISM3D with previous work is verified by comparing it with two recent studies, with which it shows a good general agreement.The impact of the new 3-D model is illustrated through a simple synthetic experiment, which shows that the lateral variations of the wave speed can produce traveltime differences ranging from –1.5 and 1.5 s for P waves and from –2.5 and 2.5 s for S waves at local to regional distances. Such values are far larger than phase picking uncertainties and would likely affect earthquake hypocentral parameter estimations. The new 3-D model thus provides a basis for regional studies including earthquake source studies, Earth structure investigations and geodynamic modelling of Iberia and its surroundings

A common deep source for upper-mantle upwellings below the Ibero-western Maghreb region from teleseismic P-wave travel-time tomography

Upper-mantle upwellings are often invoked as the cause of Cenozoic volcanism in the Ibero-western Maghreb region. However, their nature, geometry and origin are unclear. This study takes advantage of dense seismic networks, which cover an area extending from the Pyrenees in the north to the Canaries in the south, to provide a new high-resolution P-wave velocity model of the upper-mantle and topmost lower-mantle structure. Our images show three subvertical upper-mantle upwellings below the Canaries, the Atlas Ranges and the Gibraltar Arc, which appear to be rooted beneath the upper-mantle transition zone (MTZ). Two other mantle upwellings beneath the eastern Rif and eastern Betics surround the Gibraltar subduction zone. We propose a new geodynamic model in which narrow upper-mantle upwellings below the Canaries, the Atlas Ranges and the Gibraltar Arc rise from a laterally-propagating layer of material below the MTZ, which in turn is fed by a common deep source below the Canaries. In the Gibraltar region, the deeply rooted upwelling interacts with the Gibraltar slab. Quasi-toroidal flow driven by slab rollback induces the hot mantle material to flow around the slab, creating the two low-velocity anomalies below the eastern Betics and eastern Rif. Our results suggest that the Central Atlantic plume is a likely source of hot mantle material for upper-mantle upwellings in the Ibero-western Maghreb region

Optimal resolution tomography with error tracking and the structure of the crust and upper mantle beneath Ireland and Britain

The classical Backus–Gilbert method seeks localized Earth-structure averages at the shortest length scales possible, given a data set, data errors, and a threshold for acceptable model errors. The resolving length at a point is the width of the local averaging kernel, and the optimal averaging kernel is the narrowest one such that the model error is below a specified level. This approach is well suited for seismic tomography, which maps 3-D Earth structure using large sets of seismic measurements. The continual measurement-error decreases and data-redundancy increases have reduced the impact of random errors on tomographic models. Systematic errors, however, are resistant to data redundancy and their effect on the model is difficult to predict. Here, we develop a method for finding the optimal resolving length at every point, implementing it for surface-wave tomography. As in the Backus–Gilbert method, every solution at a point results from an entire-system inversion, and the model error is reduced by increasing the model-parameter averaging. The key advantage of our method stems from its direct, empirical evaluation of the posterior model error at a point. We first measure inter- station phase velocities at simultaneously recording station pairs and compute phase-velocity maps at densely, logarithmically spaced periods. Numerous versions of the maps with varying smoothness are then computed, ranging from very rough to very smooth. Phase-velocity curves extracted from the maps at every point can be inverted for shear-velocity (V S ) profiles. As we show, errors in these phase-velocity curves increase nearly monotonically with the map roughness. We evaluate the error by isolating the roughness of the phase-velocity curve that cannot be explained by any Earth structure and determine the optimal resolving length at a point such that the error of the local phase-velocity curve is below a threshold. A 3-D V S model is then computed by the inversion of the composite phase-velocity maps with an optimal resolution at every point. The estimated optimal resolution shows smooth lateral variations, confirming the robustness of the procedure. Importantly, the optimal resolving length does not scale with the density of the data coverage: some of the best-sampled locations display relatively low lateral resolution, probably due to systematic errors in the data. We apply the method to image the lithosphere and underlying mantle beneath Ireland and Britain. Our very large data set was created using new data from Ireland Array, the Irish National Seismic Network, the UK Seismograph Network and other deployments. A total of 11 238 inter-station dispersion curves, spanning a very broad total period range (4–500 s), yield unprecedented data coverage of the area and provide fine regional resolution from the crust to the deep asthenosphere. The lateral resolution of the 3-D model is computed explicitly and varies from 39 km in central Ireland to over 800 km at the edges of the area, where the data coverage declines. Our tomography reveals pronounced, previously unknown variations in the lithospheric thickness beneath Ireland and Britain, with implications for their Caledonian assembly and for the mechanisms of the British Tertiary Igneous Province magmatism

SISMIKO:emergency network deployment and data sharing for the 2016 central Italy seismic sequence

At 01:36 UTC (03:36 local time) on August 24th 2016, an earthquake Mw 6.0 struck an extensive sector of the central Apennines (coordinates: latitude 42.70° N, longitude 13.23° E, 8.0 km depth). The earthquake caused about 300 casualties and severe damage to the historical buildings and economic activity in an area located near the borders of the Umbria, Lazio, Abruzzo and Marche regions. The Istituto Nazionale di Geofisica e Vulcanologia (INGV) located in few minutes the hypocenter near Accumoli, a small town in the province of Rieti. In the hours after the quake, dozens of events were recorded by the National Seismic Network (Rete Sismica Nazionale, RSN) of the INGV, many of which had a ML > 3.0. The density and coverage of the RSN in the epicentral area meant the epicenter and magnitude of the main event and subsequent shocks that followed it in the early hours of the seismic sequence were well constrained. However, in order to better constrain the localizations of the aftershock hypocenters, especially the depths, a denser seismic monitoring network was needed. Just after the mainshock, SISMIKO, the coordinating body of the emergency seismic network at INGV, was activated in order to install a temporary seismic network integrated with the existing permanent network in the epicentral area. From August the 24th to the 30th, SISMIKO deployed eighteen seismic stations, generally six components (equipped with both velocimeter and accelerometer), with thirteen of the seismic station transmitting in real-time to the INGV seismic monitoring room in Rome. The design and geometry of the temporary network was decided in consolation with other groups who were deploying seismic stations in the region, namely EMERSITO (a group studying site-effects), and the emergency Italian strong motion network (RAN) managed by the National Civil Protection Department (DPC). Further 25 BB temporary seismic stations were deployed by colleagues of the British Geological Survey (BGS) and the School of Geosciences, University of Edinburgh in collaboration with INGV. All data acquired from SISMIKO stations, are quickly available at the European Integrated Data Archive (EIDA). The data acquired by the SISMIKO stations were included in the preliminary analysis that was performed by the Bollettino Sismico Italiano (BSI), the Centro Nazionale Terremoti (CNT) staff working in Ancona, and the INGV-MI, described below

Le attività del gruppo operativo INGV "SISMIKO" durante la sequenza sismica "Amatrice 2016",

SISMIKO è un gruppo operativo dell’Istituto Nazionale di Geofisica e Vulcanologia (INGV) che coordina tutte le Reti Sismiche Mobili INGVPublishedLecce3T. Sorgente sismica4T. Sismicità dell'Italia8T. Sismologia in tempo reale1SR TERREMOTI - Sorveglianza Sismica e Allerta Tsunami2SR TERREMOTI - Gestione delle emergenze sismiche e da maremoto3SR TERREMOTI - Attività dei Centr