Guilhem Barruol

A quantitative evaluation of the contribution of crustal rocks to the shear-wave splitting of teleseismic SKS waves

Abstract The seismic (elastic) properties of the crust have been modelled using hypothetical polycrystals with typical crustal compositions and commonly measured petrofabrics. This modelling allows us to better evaluate and quantify the crustal contribution to splitting of teleseismic SKS-waves. Mafic rocks such as anorthosite or pyroxene-bearing gabbros have complex S-waves properties, i.e. with small shear-wave anisotropies (

Seismic anisotropy in the eastern United States: Deep structure of a complex continental plate

We have analyzed shear wave splitting recorded by portable and permanent broadband and long-period stations located in the eastern United States. Teleseismic shear waves (SKS, SKKS, and PKS) were used to retrieve the splitting parameters: the orientation of the fast wave polarization plane ϕ and the delay time δt. In total, 120 seismic events were processed, allowing for more than 600 splitting measurements. Within the Appalachians, stations located in the western (external) part are characterized by δt≈1s and ϕ trending N50°– 70°E in the south and central regions and N30°–40°E in the north, closely following the trend of the orogenic belt in these areas. The transition region between north and central is characterized by δt≈1–1.3 s and by E-W trending ϕ that are at a high angle to the regional geologic trend. Measurements at two stations located in the eastern (internal) part of the belt indicate very weak anisotropy. The large-scale pattern of anisotropy is not consistent with that predicted for simple asthenospheric flow beneath the plate. Splitting along the southern and eastern margins of the continent is consistent with that expected for Grenvillian deformation, an alternative model of asthenospheric flow around the cratonic keel cannot be ruled out. Within the cratonic core, the correlation between δt and lithospheric thickness suggests a lithospheric anisotropy. Smaller-length-scale variations also argue for a significant contribution of lithospheric structures. The fabric responsible for shear wave splitting may have formed during tectonic episodes that affected the eastern United States, i.e., the Grenville and Appalachian orogenies and the subsequent rifting of the North Atlantic Ocean. Our observations in the western Appalachians suggest that the anisotropy may be preserved since the Grenvillian orogeny. The absence of detectable splitting in the two stations in the eastern Appalachians is attributed to the igneous intrusions related to the Atlantic rifting. The measurements in the transition between the northern and central southern Appalachians, constitute an intriguing anomaly, whose E-W ϕ have little obvious relation to the regional surface geology. We suggest two possible causes: (1) the local dominance of asthenospheric flow, motivated by the proximity of a pervasive low-velocity anomaly and (2) lithospheric deformation in a transcontinental strike-slip fault zone active during the Appalachian collision.

Seismic anisotropy and shear-wave splitting in lower-crustal and upper-mantle rocks from the Ivrea Zone—experimental and calculated data

To quantify the relationships between anisotropy, S-wave splitting and tectonics, we determined the seismic properties of lower-crustal and upper-mantle rocks outcropping in the Ivrea Zone (Northern Italy). We obtained P- and S-wave seismic velocities by laboratory direct velocity measurements and/or by calculations based on the modal compositions of the rocks, the lattice preferred orientations (LPOs), and the single crystal stiffness coefficients. Measured P- and S-wave velocities (6.0-7.5 km s-1 and 3.6-4.2 km s-1) are typical of the lower crust. The P-wave anisotropy is in the range 0-10%. Shear-wave birefringence is in the range 0.0-0.6 km s-1, with typical values between 0.0 and 0.2 km s-1. In many cases, the birefringence is clearly related to fabric elements (foliation, lineation). Mafic rocks such as anorthosite or pyroxene-bearing gabbros exhibit low P-wave anisotropies ( < 5%) and low shear-wave birefringences (less than 0.1 km s-1). In contrast, the seismic properties of felsic rocks such as biotite-bearing gneisses and mafic rocks such as amphibolites display high Vp anisotropy (10%) and strong birefringence (0.3 km s-1). Biotite and amphibole preferred orientations clearly control seismic anisotropy and particularly shear-wave splitting. In these rocks, maximum splitting is observed in directions parallel to the foliation with the fast split shear wave polarized parallel to the foliation plane. To have an overview of the seismic properties of this lower-crustal section at a broader scale, we calculated from our data the anisotropic seismic properties of several hypothetical samples that are perhaps more representative of the regional anisotropy than each sample individually. For instance, the average lower-crustal sample displays an anisotropy of 5.5% for P waves and a birefringence around 0.14 km s-1 for S waves propagating parallel to the foliation. We observe little splitting for waves propagating at high angle to the foliation.

Upper mantle anisotropy beneath the Geoscope stations

Seismic anisotropy has been widely studied this last decade, particularly by measuring splitting of vertically propagating core shear waves. The main interest in this technique is to characterize upper mantle flow beneath seismic stations. On the other hand, the major restriction in this method is that a single station gives a single anisotropy measurement. Alternative methods have been developed in order to avoid this restriction. An accurate determination of upper mantle seismic anisotropy beneath a seismic station may allow one, by doing anisotropy correction, to characterize remote or deeper anisotropy. The Geoscope network is ideal for this purpose because it is composed of a large set (about 26) of high-quality, broadband seismometers globally distributed and because some of these stations have run for more than 10 years and most of them for more than 5 years. We selected about 100 events at each site, generally of magnitude (mb) > 6.0, and we performed systematic measurements of the splitting parameters (fast polarization direction ϕ and delay time δt) on SKS, SKKS, and PKS phases. Splitting on oceanic islands has been difficult to observe owing to the low quality of the signal but also perhaps owing to complex upper mantle structures beneath the stations. Station KIP (Kipapa, Hawaii) in the Pacific is the only oceanic Geoscope station with a clear anisotropy. We determined well-constrained splitting parameters for 10 of the 17 continental stations that may be explained by a single anisotropic layer. The poor correlation between fast polarization directions and the absolute plate motion together with the apparent incoherence between the plate velocities and the observed delay times suggest that a simple drag-induced asthenospheric flow alone fails to explain most of the observations. For some stations located on or near major lithospheric structures (TAM, Tamanrasset, Algeria, for instance), we observe a good correlation between fast polarization directions and regional structures. At station SCZ (Santa Cruz, California), we found clear variations of the splitting parameters as a function of the event backazimuth, compatible with two layers of anisotropy. Three stations (CAN (Canberra), HYB (Hyderabad, India) and SSB (Saint Sauveur Badole, France)) seem to be devoid of detectable anisotropy.

Rheological heterogeneity, mechanical anisotropy and deformation of the continental lithosphere

Abstract This paper aims to present an overview on the influence of rheological heterogeneity and mechanical anisotropy on the deformation of continents. After briefly recapping the concept of rheological stratification of the lithosphere, we discuss two specific issues: (1) as supported by a growing body of geophysical and geological observations, crust/mantle mechanical coupling is usually efficient, especially beneath major transcurrent faults which probably crosscut the lithosphere and root within the sublithospheric mantle; and (2) in most geodynamic environments, mechanical properties of the mantle govern the tectonic behaviour of the lithosphere. Lateral rheological heterogeneity of the continental lithosphere may result from various sources, with variations in geothermal gradient being the principal one. The oldest domains of continents, the cratonic nuclei, are characterized by a relatively cold, thick, and consequently stiff lithosphere. On the other hand, rifting may also modify the thermal structure of the lithosphere. Depending on the relative stretching of the crust and upper mantle, a stiff or a weak heterogeneity may develop. Observations from rift domains suggest that rifting usually results in a larger thinning of the lithospheric mantle than of the crust, and therefore tends to generate a weak heterogeneity. Numerical models show that during continental collision, the presence of both stiff and weak rheological heterogeneities significantly influences the large-scale deformation of the continental lithosphere. They especially favour the development of lithospheric-scale strike-slip faults, which allow strain to be transferred between the heterogeneities. An heterogeneous strain partition occurs: cratons largely escape deformation, and strain tends to localize within or at the boundary of the rift basins provided compressional deformation starts before the thermal heterogeneity induced by rifting are compensated. Seismic and electrical conductivity anisotropies consistently point towards the existence of a coherent fabric in the lithospheric mantle beneath continental domains. Analysis of naturally deformed peridotites, experimental deformations and numerical simulations suggest that this fabric is developed during orogenic events and subsequently frozen in the lithospheric mantle. Because the mechanical properties of single-crystal olivine are anisotropic, i.e. dependent on the orientation of the applied forces relative to the dominant slip systems, a pervasive fabric frozen in the mantle may induce a significant mechanical anisotropy of the whole lithospheric mantle. It is suggested that this mechanical anisotropy is the source of the so-called tectonic inheritance, i.e. the systematic reactivation of ancient tectonic directions; it may especially explain preferential rift propagation and continental break-up along pre-existing orogenic belts. Thus, the deformation of continents during orogenic events results from a trade-off between tectonic forces applied at plate boundaries, plate geometry, and the intrinsic properties (rheological heterogeneity and mechanical anisotropy) of the continental plates.

EBSD-measured lattice-preferred orientations and seismic properties of eclogites

We investigated the deformation mechanisms and the seismic properties of 10 eclogite samples from different localities (Alps, Norway, Mali and eastern China) through the analysis of their microstructures and lattice-preferred orientations (LPO). These samples are representative of various types and intensity of deformation under eclogitic metamorphic conditions. Omphacite and garnet LPO were determined from electron backscatter diffraction (EBSD) technique. Garnet appears to be almost randomly oriented whereas omphacite develops strong LPO, characterized by the [001]-axes concentrated sub-parallel to the lineation, and the (010)-poles concentrated sub-perpendicular to the foliation. In order to analyze the deformation mechanisms that produced such omphacite LPO, we compare our observations to LPO simulated by viscoplastic self-consistent numerical models. A good fit to the measured LPO is obtained for models in which the dominant slip systems are 1/2h110i{11 ¯ 0}, [001] {110} and [001] (100). Dominant activation of these slip systems is in agreement with TEM studies of naturally deformed omphacite. Seismic properties of eclogite are calculated by combining the measured LPO and the single crystal elastic constants of omphacite and garnet. Although eclogite seismic anisotropies are very weak (less than 3% for both P-and S-wave), they are generally characterized by a maximum P-wave velocity sub-parallel to the lineation and by a minimum velocity approximately normal to foliation. The mean P-and S-wave velocities are high (respectively, 8.6 and 4.9 km/s). The S-wave anisotropy pattern displays complex relationships with the structural frame but the fast polarization plane generally tends to be parallel to the foliation. Calculated reflection coefficients show that an eclogite/crust interface is generally a good reflector (Rc > 0.1), whereas an eclogite body embedded in the upper mantle would be hardly detectable.

The Kaapvaal craton seismic anisotropy: Petrophysical analyses of upper mantle kimberlite nodules

A dense network of seismic stations has been deployed on the Kaapvaal craton (South Africa) to investigate the upper mantle seismic structures. In order to bring independent petrophysical constraints, we analyze a direct sampling of the cratonic upper mantle and determine the seismic properties of 48 mantle nodules brought up to the Earths surface by kimberlite eruptions. Seismic properties of these nodules are calculated from the olivine and pyroxene crystal preferred orientations and the single crystal elastic constants. Despite variations in the nodules compositions, microstructures and crystallographic preferred orientations, seismic anisotropy is rather homogeneous throughout the craton. Mean S-wave anisotropy is weak (2.64 %), which is compatible with the small measured SKS wave splitting (mean delay time of 0.62 s).

Upper mantle deformation and seismic anisotropy in continental rifts

Abstract We investigate the seismic anisotropy signature of the continental rifting process. Several sources of anisotropy are considered: the lithospheric deformation, the asthenospheric flow, and the occurrence of oriented meltpockets in the asthenospheric mantle. Our results show that contrasted anisotropy patterns should be associated with the various conceptual models of rifting. Thus seismic anisotropy measurements may allow one to discriminate between these models. Anisotropy measurements in the Rio Grande, East-African and Rhine rifts suggest that these rifts formed by a transtensional deformation of the lithospheric mantle rather than by homogeneous extension of the lithosphere. Alignment of melt-lenses in the asthenospheric wedge may also account for a significant part of the seismic anisotropy recorded in the internal domains of these rifts.

Lithospheric anisotropy beneath the Pyrenees from shear wave splitting

We investigate upper mantle anisotropy beneath the Pyrenean range along three N-S profiles across the mountain belt. The results of a first profile that operated in 1993 in the central part of the belt have been presented elsewhere. We present the results of two other profiles that ran in 1995-1996 and 1996-1997 in the eastern and western part of the belt, respectively and propose an interpretation of the whole results. Teleseismic shear waves (SKS, SKKS, and PKS) are used to determine splitting parameters: the fast polarization direction φ and the delay time δt. Teleseismic shear wave splitting in the eastern Pyrenees displays homogeneous φ values trending N100°E and δt values in the range 1.1 to 1.5 s. A station located in the southern Massif Central, 100 km north of the range, is characterized by different splitting parameters (φ = N70°E, δt = 0.7 s). In the western part of the belt, anisotropy parameters are similar across the whole belt (φ = N110°E and δt = 1.3 to 1.5 s). Most of the measured delay times, including those obtained in the central part of the range, are above the global average of the SKS splitting (around 1 s). At the belt scale, φ is generally poorly correlated with recent estimations of the absolute plate motion, which predicts a fast direction ranging between N50°E and N80°E. Instead, the orientation of φ (N100°E) is parallel to the trend of the Pyrenean belt but also to Hercynian preexisting structures. This parallelism supports an anisotropy primarily related to frozen or active lithospheric structures. We show that a signature related to the Pyrenean orogeny is likely for the stations located in the internal domains of the belt. By contrast, the anisotropy measured at the stations located on the external parts of the belt could reflect a pre-Pyrenean (Hercynian) deformation. We suggest that a late Hercynian strike-slip deformation is responsible for this frozen upper mantle anisotropy and that the Pyrenean tectonic fabric developped parallel to this preexisting fabric. Finally, no particularly strong splitting is related to the North Pyrenean Fault, commonly believed to represent the plate boundary between Iberia and Eurasia.

A Tertiary asthenospheric flow beneath the southern French Massif Central indicated by upper mantle seismic anisotropy and related to the west Mediterranean extension

Abstract Upper mantle flow beneath the French Massif Central is investigated using teleseismic shear wave splitting induced by seismic anisotropy. About 25 three-component stations (short period, intermediate and broadband) were installed during the period 1998–1999 in the southern Massif Central, from the Clermont Ferrand volcanic area to the Mediterranean Sea. Teleseismic shear waves (SKS, SKKS and PKS) were used to determine the splitting parameters: the fast polarization direction and the delay time. Delay times ranging between 0.7 and 1.5 s have been observed at most of the sites. The azimuths of the fast split shear waves trend homogeneously NW–SE in the southern Massif Central suggesting a homogeneous mantle flow beneath this area. The observed NW–SE direction differs from the N100°E Pyrenean anisotropy further south. It does not appear to be correlated to Hercynian structures nor to the present-day motion of the plate but is well correlated to the Tertiary extension direction. We propose that the opening of the western Mediterranean induced by the rotation of the Corsica–Sardinia lithospheric block and the roll-back to the SE of the Tethys slab may have generated a large asthenospheric mantle flow beneath the southern Massif Central and a deflection of the up going plume centered beneath the northern Massif Central toward the SE.