3-D seismic velocities calculated from lattice-preferred orientation and reflectivity of a lower crustal section: examples of the Val Sesia section (Ivrea zone, northern Italy)

International audienceTo quantify the seismic properties of lower crustal rocks and to better constrain the origin of the lower crustal seismic reflectivity, we determined the complete 3-D seismic properties of a lower crustal section. Eight representative samples of the main lithologic and structural units outcropping in the Val Sesia (Ivrea zone) were studied in detail. The seismic velocities were calculated using the single crystal stiffness coefficients and the lattice preferred orientation (LPO) of each mineral in all samples. The 21 stiffness coefficients characterizing the elastic behaviour of each rock are determined. Mafic and ultramafic rocks such as pyroxenite and pyroxene-bearing gabbros display complex shear wave properties. These rocks are weakly birefringent (maximum 0.1 kms-1) and it is difficult to find consistent relationships between the seismic properties and the rock structure. On the other hand, seismic properties of deformed felsic rocks are essentially controlled by mica. They display strong S-wave birefringence (0.3 km s-1) and relatively high Vp anisotropy (7.6 per cent). Amphibole also strongly influences the rock birefringence patterns. For both kind of rocks, the foliation is highly birefringent and the fast polarized shear wave is systematically oriented parallel to the foliation. We show that the number of mineral phases in the rock strongly controls the anisotropy. The seismic anisotropy has a complex role in the P-wave reflectivity. Compared to the isotropic case, anisotropy enhances the reflection coefficient for about 60 per cent of the possible lithological interfaces. For 40 per cent of the interfaces, the reflection coefficient is much lower when one considers the medium as anisotropic

Mantle structural geology from seismic anisotropy

International audienceSeismic anisotropy is a ubiquitous feature of the subcontinental mantle. This can be inferred both from direct seismic observations of shear wave splitting from teleseismic shear waves, as well as the petrofabric analyses of mantle nodules from kimberlite pipes. The anisotropy is principally due to the strain-induced lattice preferred orientation (LPO) of olivine. The combined use of these mantle samples, deformation experiments on olivine, and numerical modeling of LPO, provides a critical framework for making inferences about mantle deformation from observed seismic anisotropy. In most cases there is a close correspondence between mantle deformation derived from seismic observations of anisotropy, and crustal deformation, from the Archean to the present. This implies that the mantle plays a major, if not dominant role in continental deformation. No clear evidence is found for a continental asthenospheric decoupling zone, suggesting that continents are probably coupled to general mantle circulation

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

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

Effects of crystallographic anisotropy on fracture development and acoustic emission in quartz

Transgranular microcracking is fundamental for the initiation and propagation of all fractures in rocks. The geometry of these microcracks is primarily controlled by the interaction of the imposed stress field with the mineral elastic properties. However, the effects of anisotropic elastic properties of minerals on brittle fracture are not well understood. This study examines the effects of elastic anisotropy of quartz on the geometry of brittle fracture and related acoustic emissions (AE) developed during indentation experiments on single crystals at ambient pressure and temperature. A Hertzian cone crack developed during blunt indentation of a single crystal of flawless Brazilian quartz parallel to the c axis shows geometric deviation away from predictions based on the isotropic case, consistent with trigonal symmetry. The visible cone crack penetration depth varies from 3 to 5 mm and apical angle from 53 to 40. Electron backscatter diffraction (EBSD) mapping of the crack tip shows that fracturing initiates along a ~40 μm wide process zone, comprising damage along overlapping en echelon high-index crystallographic planes, shown by discrete bands of reduced electron backscatter pattern (EBSP) quality (band contrast).Coalescence of these surfaces results in a stepped fracture morphology. Monitoring of AE during indentation reveals that the elastic anisotropy of quartz has a significant effect on AE location and focal mechanisms. Ninety-four AE events were recorded during indentation and show an increasing frequency with increasing load. They correspond to the development of subsidiary concentric cracks peripheral to the main cone crack. The strong and complex anisotropy in seismic velocity (~28% Vp, ~43% Vs with trigonal symmetry) resulted in inaccurate and high uncertainty in AE locations using Geiger location routine with an isotropic velocity model. This problem was overcome by using a relative (master event) location algorithm that only requires a priori knowledge of the velocity structure within the source volume. The AE location results correlate reasonably well to the extent of the observed cone crack. Decomposition of AE source mechanisms of the Geiger relocated events shows dominantly end-member behavior between tensile and compressive vector dipole events, with some double-couple-dominated events and no purely tensile or compressive events. The same events located by the master event algorithm yield greater percentage of vector dipole components and no double-couple events, indicating that AE source mechanism solutions can depend on AE location accuracy, and therefore, relocation routine that is utilized. Calculations show that the crystallographic anisotropy of quartz causes apparent deviation of the moment tensors away from double-couple and pure tensile/compressive sources consistent with the observations. Preliminary modeling of calcite anisotropy shows a response distinct from quartz, indicating that the effects of anisotropy on interpreting AE are complex and require detailed further study

A Physics‐Based Rock Friction Constitutive Law: Steady State Friction

Experiments measuring friction over a wide range of sliding velocities find that the value of the friction coefficient varies widely: friction is high and behaves according to the rate and state constitutive law during slow sliding, yet markedly weakens as the sliding velocity approaches seismic slip speeds. We introduce a physics‐based theory to explain this behavior. Using conventional microphysics of creep, we calculate the velocity and temperature dependence of contact stresses during sliding, including the thermal effects of shear heating. Contacts are assumed to reach a coupled thermal and mechanical steady state, and friction is calculated for steady sliding. Results from theory provide good quantitative agreement with reported experimental results for quartz and granite friction over 11 orders of magnitude in velocity. The new model elucidates the physics of friction and predicts the connection between friction laws to independently determined material parameters. It predicts four frictional regimes as function of slip rate: at slow velocity friction is either velocity strengthening or weakening, depending on material parameters, and follows the rate and state friction law. Differences between surface and volume activation energies are the main control on velocity dependence. At intermediate velocity, for some material parameters, a distinct velocity strengthening regime emerges. At fast sliding, shear heating produces thermal softening of friction. At the fastest sliding, melting causes further weakening. This theory, with its four frictional regimes, fits well previously published experimental results under low temperature and normal stress

Diffusion, deformation and mineral properties of the Earth's interior. Introduction

International audienceThis special issue in honour of the late Olivier Jaoul brings together contributions from experimental and theoretical mineralogy, the area of most of Olivier's own work, focused for the most part on the Earth's mantle. The volume is divided into three sections; diffusion, deformation and seismic properties. New technological developments have often had a strong impact on research in this area, as again shown by Olivier's constant search to improve the resolution of the short diffusion profiles typical of major cations in mantle silicates. In the section on diffusion the impact of improved computer performance for modelling diffusion and point defect configurations on one hand, and high pressure experimental techniques on the other hand, highlight the continued development in this area. The section on deformation deals with wide range of subjects, including the deformation of olivine and quartz and the relation between deformation and diffusion, both subjects close to Olivier's personal interests. This illustrates the new possibilities to undertake deformation experiments at mantle pressures and temperatures, as well as improvement in high-fidelity viscoelastic measurements at high temperature. In the last section on seismic properties we discover theoretical developments on predicting properties near the Earth's surface and a new experimental application, which allows the determination of sound velocity at very high pressures with implications for the Earth's inner core