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

Inherited crustal deformation along the East Gondwana margin revealed by seismic anisotropy tomography

Acknowledgments We thank Mallory Young for providing phase velocity measurements in mainland Australia and Tasmania. Robert Musgrave is thanked for making available his tilt-filtered magnetic intensity map. In the short term, data may be made available by contacting the authors (S.P. or N.R.). A new database of passive seismic data recorded in Australia is planned as part of a national geophysics data facility for easy access download. Details on the status of this database may be obtained from the authors (S.P., N.R., or A.M.R.). There are no restrictions on access for noncommercial use. Commercial users should seek written permission from the authors (S.P. or N.R.). Ross Cayley publishes with the permission of the Director of the Geological Survey of Victoria.Peer reviewedPublisher PD

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

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

Inherited crustal deformation along the East Gondwana margin revealed by seismic anisotropy tomography

The mechanisms of continental growth are a crucial part of plate tectonic theory, yet a clear understanding of the processes involved remains elusive. Here we determine seismic Rayleigh wave phase anisotropy variations in the crust beneath the southern Tasmanides of Australia, a Paleozoic accretionary margin. Our results reveal a complex, thick-skinned pervasive deformation that was driven by the tectonic interaction between the proto-Pacific Ocean and the ancient eastern margin of Gondwana. Stress-induced effects triggered by the collision and entrainment of a microcontinent into the active subduction zone are evident in the anisotropy signature. The paleofracturing trend of failed rifting between Australia and Antarctica is also recorded in the anisotropy pattern as well as a tightly curved feature in central Tasmania. The observed patterns of anisotropy correlate well with recent geodynamic and kinematic models of the Tasmanides and provide a platform from which the spatial extent of deformational domains can be refined

Inherited crustal deformation along the East Gondwana margin revealed by seismic anisotropy tomography

The mechanisms of continental growth are a crucial part of plate tectonic theory, yet a clear understanding of the processes involved remains elusive. Here we determine seismic Rayleigh wave phase anisotropy variations in the crust beneath the southern Tasmanides of Australia, a Paleozoic accretionary margin. Our results reveal a complex, thick-skinned pervasive deformation that was driven by the tectonic interaction between the proto-Pacific Ocean and the ancient eastern margin of Gondwana. Stress-induced effects triggered by the collision and entrainment of a microcontinent into the active subduction zone are evident in the anisotropy signature. The paleofracturing trend of failed rifting between Australia and Antarctica is also recorded in the anisotropy pattern as well as a tightly curved feature in central Tasmania. The observed patterns of anisotropy correlate well with recent geodynamic and kinematic models of the Tasmanides and provide a platform from which the spatial extent of deformational domains can be refined. ©2016. American Geophysical Union. All Rights Reserved

Crustal and Mantle Deformation Inherited From Obduction of the Semail Ophiolite (Oman) and Continental Collision (Zagros)

Abstract: A common deviation from typical subduction models occurs when thrust sheets of oceanic crust and upper‐mantle rocks are emplaced over more buoyant continental lithosphere. The archetypal example of ophiolite obduction is the Semail ophiolite in the United Arab Emirates (UAE)‐Oman orogenic belt, formed and obducted onto the Arabian continental margin during the Late Cretaceous. The Strait of Hormuz syntaxis, the northern extent of the UAE‐Oman mountains, marks the transition from ocean‐continent convergence in the Gulf of Oman to continental collision along the Zagros Mountains. Based on new seismic data from a focused recording network, we infer continental crustal and mantle deformation in the northeastern corner of the Arabian plate (including the southern Zagros and the UAE‐Oman mountains), using observations from anisotropic tomography and shear‐wave splitting (SWS) measurements. We recover a change of ∼90° (from approximately WNW to nearly NS) in the axis of fast‐anisotropic orientations in the crust from the Zagros to the UAE‐Oman mountain belt, consistent with the dominant strike of the orogenic belts. We also find evidence in our SWS parameters for localized fossil deformation in the lithospheric mantle underlying the UAE‐Oman mountain range, possibly related to stress‐induced tectonism triggered by north‐east oriented underthrusting of the proto‐Arabian continental margin beneath the overriding Semail ophiolite. Shear‐wave‐splitting anisotropy orientations along two transects across the northern Musandam peninsula, averaging 15° anticlockwise from the north, provide the first geophysical verification of previous geological evidence that suggests a NE polarity of the Late Cretaceous Oman subduction zone system

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

Transdimensional inversion of receiver functions and surface wave dispersion

International audienceWe 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