Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We investigated the response of the crystal structure, lattice parameters, and unit-cell volume of hydrous layered silicate phlogopite at conditions relevant to subduction zone settings. We have used first principles simulation based on density functional theory to calculate the equation of state and full elastic constant tensor. Based on the generalized gradient approximation, the full single crystal elastic constant tensor with monoclinic symmetry shows significant anisotropy with the compressional elastic constants: c11 = 181 GPa, c22 = 185 GPa, c33 = 62 GPa, the shear elastic constants c44 = 14 GPa, c55 = 20 GPa, c66 = 68 Ga, and c46 = −6 GPa; the off diagonal elastic constants c12 = 48 GPa, c13 = 12 GPa, c23 = 12 GPa, c15 = −16 GPa, c25 = −5 GPa and c35 = −1 GPa at zero pressure. The elastic anisotropy of phlogopite is larger than most of the layered hydrous phases relevant in the subduction zone conditions. The shear anisotropy, AVS for phlogopite is ∼77% at zero pressure condition and although it decreases upon compression it remains relatively high compared to other hydrous phases relevant in the subduction zone settings. We also note that the shear elastic constants for phlogopite are relatively low. Phlogopite also has a high isotropic bulk VP/VS ratio ∼2.0. However, the VP/VS ratio also exhibits significant anisotropy with values as low as 1.49. Thus, phlogopite bearing metasomatized mantle could readily explain unusual VP/VS ratio as observed from seismological studies from the mantle wedge regions of the subduction zone.
ISSN: 0031-9201

hal-01115670
https://hal.archives-ouvertes.fr/hal-01115670
DOI : 10.1016/j.pepi.2014.05.004

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Talc is a hydrous magnesium rich layered silicate that is widely disseminated in the Earth from the seafloor to over 100 km depth, in ultra-high pressure metamorphism of oceanic crust. In this paper we determine the single crystal elastic constants at pressures from 0 to 12 GPa of talc triclinic (C 1¯) and monoclinic (C2/c) polytypes using ab initio methods. We find that talc has an extraordinarily high elastic anisotropy at zero pressure that reduces with increasing pressure. The exceptional anisotropy is complemented by a negative Poisson's ratio for many directions in crystal space. Calculations show that talc is not only one of very few common minerals to exhibit auxetic behaviour, but the magnitude of this effect may be the largest reported so far for a mineral. The compression (Vp) and shear (Vs) wave velocity anisotropy is 80% and 85% for the triclinic polytype. At pressures where talc is known be stable in the Earth (up to 5 GPa) the Vp and Vs anisotropy is reduced to about 40% for both velocities, which is still a very high value. Vp is slow parallel to the c-axis and fast perpendicular to it. This remains unchanged with increasing pressure and is observed in both polytypes. The shear wave splitting (difference between fast and slow S-wave velocities) at low pressure has high values in the plane normal to the c-axis, with a maximum near the alow asterisk-axis in the triclinic and the b-axis in the monoclinic polytype. The c-axis is the direction of minimum splitting. The pattern of shear wave splitting does not change significantly with pressure. The volume fraction of talc varies between 11 and 41% for hydrated mantle rocks, but the lack of data on the crystallographic preferred orientation (CPO) precludes a detailed analysis of the impact of talc on seismic anisotropy in subduction zones. However, it is highly likely that CPO can easily develop in zones of deformation due to the platy habit of talc crystals. For random aggregates of talc, the isotropic Vp, Vs and Vp/Vs ratio have significantly lower values than those of antigorite and may explain low-velocity regions in the mantle wedge. Vp/Vs ratios are more complex in anisotropic media because there are fast and slow S-waves, resulting in Vp/Vs1 and Vp/Vs2 ratios for every propagation direction, making interpretation difficult in deformed polycrystalline talc with a CPO. Talc on the subduction plate boundary can strongly influence guided wave velocity as CPO would develop in this region of intense shearing. The very low coefficient of friction (< 0.1) of talc above 100 °C could also explain silent earthquakes at shallow depths (ca 30 km) along the subduction plate boundaries, frequently responsible for tsunami.
ISSN: 0012-821X

hal-00411481
https://hal.archives-ouvertes.fr/hal-00411481
DOI : 10.1016/j.epsl.2008.07.047

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
Using first principle simulations we calculated the elasticity of chlorite. At a density ρ~ 2.60 g cm−3, the elastic constant tensor reveals significant elastic anisotropy: VP ~27%, VS1 ~56%, and VS2 ~43%. The shear anisotropy is exceptionally large for chlorite and enhances upon compression. Upon compression, the shear elastic constant component C44 and C55 decreases, whereas C66 shear component stiffens. The softening in C44 and C55 is reflected in shear modulus, G, and the shear wave velocity, VS. Our results on elastic anisotropy at conditions relevant to the mantle wedge indicates that a 10-20 km layer of hydrated peridotite with serpentine and chlorite could account for the observed shear polarization anisotropy and associated large delay times of 1-2 s observed in some subduction zone settings. In addition, chlorite could also explain the low VP/VS ratios that have been observed in recent high-resolution seismological studies.
ISSN: 0094-8276

hal-01053707
https://hal.archives-ouvertes.fr/hal-01053707
https://hal.archives-ouvertes.fr/hal-01053707/document
https://hal.archives-ouvertes.fr/hal-01053707/file/Mookherjee_et_al-2014-Geophysical_Research_Letters.pdf
DOI : 10.1002/2014GL059334

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : A differential effective medium (DEM) model is used to predict elastic properties for a set of porous and anisotropic aggregates, comprised of mixtures of calcite and muscovite. The DEM takes into consideration an anisotropic background medium with triclinic or higher symmetry, in which inclusions of idealized ellipsoidal shape are added incrementally. In general, the calculated elastic properties of a solid that contains inclusions representing "dry" pores/cracks are strongly dependent on the orientation and aspect ratio of the inclusions. Aspect ratios of inclusions in the synthetic aggregates, which consist of air-filled pores, are estimated from anisotropy of magnetic susceptibility (AMS) of samples whose pore space has been impregnated with a colloidal ferrofluid. The AMS derived pore shape geometry is used as an input value for inclusions in the DEM. Modeling results are compared with laboratory determined elastic properties, measured with ultrasonic waves. Calculated shear wave velocities agree in general well with laboratory measured S wave velocities, whereas calculated P wave velocities are typically 0.5-1.1 km/s higher than measured values. Differences between calculated and measured P wave velocities are attributed mainly to incomplete and biased ferrofluid saturation of pores. Spherical pores are preferably filled during imbibition, in comparison to thin cracks, which leads to overprediction of the calculated P wave velocities. The amount of ferrofluid that fills the pore space is dependent on the ratio of calcite to muscovite and the load used for compaction during sample synthesis. The permeability decreases with increasing muscovite content and increasing compaction load. Incomplete saturation of samples with high-muscovite content is confirmed by X-ray microtomography density contrast imaging of dry and ferrofluid saturated specimens.
ISSN: 0148-0227

hal-00617717
https://hal.archives-ouvertes.fr/hal-00617717
DOI : 10.1029/2010JB007750

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
In this paper, we describe an outcrop to characterize the effect of fracture spacing and type on larger scale effective elasticity, which is measured for the first time in-situ with a Schmidt hammer. The outcrop is dominated by lime mudstones and belongs to the deformation zone of the St Clement fault, in southern France. Our results suggest that small spacing of faults, open fractures and styolites leads to lesser effective Young's modulus, whereas small sealed fracture spacing leads to greater effective Young's modulus. These relationships are compatible with theoretical models of effective elasticity. Using Amadei and Savage (1993) approach, we define a non-linear model that relates Schmidt hammer rebound to spacing by fracture type. A hemisphere with a radius of 40 to similar to 200 cm is the rheological volume characterized by the Schmidt hammer. Results of model inversion demonstrate that variations of Schmidt hammer rebound over the outcrop can be used to estimate fracture type and stiffness. Stiffness of sealed fractures is 2-3 orders of magnitude greater than the stiffness of faults, stylolites and open fractures. This result is consistent with an increase of the rate of interseismic stress build-up of major faults with sealing of fractures in their damage zone.
ISSN: 0191-8141

hal-00669954
https://hal.archives-ouvertes.fr/hal-00669954
DOI : 10.1016/j.jsg.2011.09.009

Éditeur(s) : HAL CCSD Elsevier
Résumé : The dense hydrous magnesium silicate called the D phase is the most likely candidate to recycle hydrogen into the lower mantle in subduction zones. As seismology represents the preferred method to detect this mineral in subduction zones, the single crystal elastic tensor has been calculated using first principle methods. A triple cell (View the MathML source) was used to account for an ordered ideal composition of the D phase to lower mantle pressures. The elastic tensor of the mineral at ambient conditions is predicted to be C11 = 387.7, C33 = 287.7, C44 = 100.4, C12 = 108.0, C13 = 51.1, C14 = − 14.6 in GPa with bulk modulus at zero pressure Ko = 163 GPa and a statistical error of about 4.0 GPa for all constants. The pressure derivatives of the elastic constants have been determined to pressure of 80 GPa. Our tensor correctly predicts the published experimental measurements of linear compressibility along the a- and c-axes, the logarithmic pressure derivative of the c/a axial ratio and bulk modulus at zero pressure. It contradicts the experimentally measured elastic tensor reported by (Liu, L.-G., Okamoto, K., Yang, Y.-J., Chen, C.-C., Lin, C.-C., 2004. Elasticity of single-crystal phase D (a dense hydrous magnesium silicate) by Brillouin spectroscopy. Solid State Commun. 132, 517–520), which fails to reproduce the experimental compressibility data. Using the elastic constants and density predicted to 85 GPa, the seismic anisotropy of P- and S-waves has been calculated as a function of pressure. The single crystal has high P-wave velocities and shear wave splitting in the basal plane below 20 GPa. At high pressures, the maximum P-wave velocities have threefold symmetry normal to a-axes in the basal plane, whereas the maximum shear wave splitting is parallel to a-axes with twofold symmetry. The D phase has high Vp and Vs anisotropy at ambient pressure, 17.6% for Vp and 19.9% for Vs. The Vp anisotropy decreases with increasing pressure to about 30 GPa, while it remains constant at 8% from 30 to 85 GPa. The Vs anisotropy decreases to 17% at 20 GPa and then increases to 22% at 85 GPa. The isotropic velocity ratios of S to P (Rs/p), bulk sound to S (Rphi/s) and density to S (Rρ/s) are calculated as a function of pressure. The Rs/p shows very little variation with pressure with low value of 0.6, whereas both Rphi/s and Rρ/s increase with pressure to peak values of 2.8 and 2.0 respectively at about 60 GPa. The presence of the D phase in subducted horizontal ‘stagnant' plates in the circum Pacific could contribute to a number of seismic observations, such as the observed velocity heterogeneity at depths between 500 and 1000 km and the fast S-waves travelling horizontally with horizontal polarization in this depth range. Further the lower density of the D phase compared with anhydrous minerals could influence buoyancy of the hydrated ‘stagnant' slabs. The transformation of the D phase to perovskite plus ferripericlase during thermal annealing may be responsible for horizontal slab seismogenesis.
ISSN: 0012-821X

hal-00406379
https://hal.archives-ouvertes.fr/hal-00406379
DOI : 10.1016/j.epsl.2007.04.053

Éditeur(s) : HAL CCSD Elsevier
Résumé : To contribute to the understanding of how opening-mode fractures (joints) form and open or close at depth in layered rocks, we present a 2D numerical study aiming to determine whether tensile stress can develop in pre-fractured elastic multilayers submitted to biaxial compressive strain conditions. First, we investigate the role of the elastic and geometrical properties of the layers on the development of tensile stress in models with five bonded layers and containing one open fracture in the central layer. Our results indicate that, in absence of elastic contrast (in Young's modulus) between the layers, no tensile stress develops in the models. However, when the fractured layer is stiffer than the two adjacent layers directly above and below, a lobe of horizontal tensile stress develops centered on the pre-existing fracture. The creation of this tensile stress is contingent upon the partial closing of the fracture. The levels of tensile stress and the thickness of the lobe of tensile stress increase logarithmically with an increase in the elastic contrast and are systematically larger for a larger Soft/Stiff ratio (ratio of the total thickness of the soft layers with the total thickness of the stiff layers). Second, we investigate the role of fracture interaction in the development of tensile stress in models containing a pair of open fractures. We observe that the levels of tensile stress in the region between the fractures are systematically higher than those observed in identical models containing a single fracture. This increase in tensile stress is very large for small elastic contrasts between the layers but diminishes when the elastic contrast increases. Furthermore, the spacing between the pre-existing fractures plays an important role in the stress distribution in the region between them. When the fracture spacing is equal to or lower than 1.15 times the height of the fractured layer for the experimental conditions chosen, the lobes of tensile stress centered on the fractures coalesce. This results in the formation of vast areas of tensile stress in models under remote compressive loading conditions. Such tensile areas are likely to allow the initiation and propagation of subsequent opening-mode fractures. The results obtained provide new insights into the formation of joints in layered rocks in compressive environments, with important consequences on fluid flow.
ISSN: 0040-1951

hal-00407576
https://hal.archives-ouvertes.fr/hal-00407576
DOI : 10.1016/j.tecto.2006.12.003

Éditeur(s) : HAL CCSD Geological Society of America
Résumé : International audience
We present the results of dynamic modeling of the western Mediterranean that accounts for observed global positioning system (GPS) surface deformation of the Alboran Sea and surrounding Rif and Betic Mountains as the result of the combined effects of relative motion of the Eurasian and Nubian plates, low strength in the Alboran Sea region and sub-lithospheric processes occurring beneath the External Rif domain. Assuming that the lithosphere behaves elastically over the short time period of the GPS observations, an elastic plate model is considered in our study, including an area of weak lithosphere (factor of 10) centered on the Alboran Sea and in which lateral boundary conditions consist of the Nubia-Eurasia oblique convergence. Sub-crustal processes are modeled by application of a horizontal traction on a small area (patch) at the base of the elastic plate. Our modeling studies demonstrate the need for sub-crustal or sub-lithospheric, southwestward-directed forcing to account for observed southwestward motion of the Rif and Betic domains. Based on the location, orientation, and small area of the traction patch, we hypothesize that forcing is associated with delamination and rollback of the subducted African continental lithospheric mantle beneath the External Rif zone, due to the pull of the oceanic part of the Western Mediterranean slab, a dynamic process that may be similar to that where the over-riding plate is driven toward the subduction zone during slab rollback.
ISSN: 0091-7613

hal-00524546
https://hal.archives-ouvertes.fr/hal-00524546
DOI : 10.1130/G30963.1

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Networks of localized conjugate compactive shear bands with large dihedral angles (corresponding to so-called shear-enhanced compaction bands) were generated in 3-layer plane strain finite-difference models. The central layer has elastic–plastic properties with a small positive internal friction coefficient, negative dilatancy factor, and the hardening modulus reducing during the deformation. The two other layers are elastic. The spatial organization of the networks in the central layer, resulting from the elastic–plastic instability of the 3-layer system, was shown to be dependent on the contrast of the elastic stiffness of the layers and on the friction between them. When the stiffness E of the elastic layers and/or the friction are small, the networks are rather symmetric (the bands of both orientation families are distributed homogeneously within the layer). An increase in E at high friction, results in the clustering of the bands of one orientation along a layer segment which is followed by another segment filled with the bands of the other orientation. Different network types are characterized by different band spacing which is smaller in the clustered networks. The intensity of inelastic deformation and hence of porosity reduction within the bands are also smaller in this case. The obtained band networks are very similar to those in nature where the same variability of the network organization is observed. Available geological and mechanical data confirm the conclusions from the modeling about the origin of this variability.
ISSN: 0040-1951

hal-01116201
https://hal.archives-ouvertes.fr/hal-01116201
DOI : 10.1016/j.tecto.2014.01.021

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The interseismic velocity field provided by geodetic methods is generally interpreted in the framework of a thick elastic lithosphere with a slipping fault at depth. Because lateral variations of lithospheric rheology play a key role in determining the geological strain distribution, I examine the idea that interseismic strain rate variations also occur in response to lateral variations in the elastic thickness of the lithosphere. Using a stress balance principle and some simplifying assumptions, I show using a 1D model that elastic thickness is inversely proportional to strain rate for the simple case of pure strike-slip faulting. Elastic thickness computed on three profiles crossing the San Andreas fault system (SAFS) suggests that the distribution of interseismic strain rate is compatible with a thick elastic lithosphere in the Great Basin-Sierra Nevada province and on the Pacific plate. Conversely, a thin plate with a shallow asthenosphere is needed on the SAFS to explain its high strain rate. A 2.5D Finite Element model of interseismic strain in the Carrizo Plain region in Central California shows how known vertical and horizontal variations of elastic properties refine 1D model predictions. In a case of a multiple fault system, I point out that the interseismic velocity is not causally tied to faults slip rate. Therefore, analysing the velocity field across the SAFS cannot reliably provide faults slip rate distribution as previously claimed. Rather, the apparent correlation between geologic slip rate and interseismic strain may only indicate that the elastic thickness plays a dominant role in controlling fault strength. Finally, I suggest that interseismic geodetic strain could be a new way to infer effective elastic plate thickness on the continents.
ISSN: 0012-821X

hal-00412155
https://hal.archives-ouvertes.fr/hal-00412155
DOI : 10.1016/j.epsl.2008.01.046

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Upon subduction, the oceanic crust transforms to blueschists and eclogites, with seismic properties that gradually become similar to those of the surrounding mantle. In order to evaluate the anisotropy of blueschists and glaucophane-bearing eclogites, the elastic constants of glaucophane single-crystal plates from the Sesia-Lanzo Zone (Aosta Valley, Western Alps) were measured using Brillouin spectroscopy at ambient conditions. The mean P- and S-wave velocities are 7.8 and 4.6 km s−1 respectively, and the anisotropy is high (38.1% (AVP) and 27.3% (AVS)). Glaucophane develops strong LPO, characterized by the [001]-axes concentrated sub-parallel to the lineation, and the {110} poles concentrated sub-perpendicular to the foliation in both blueschist and eclogite rocks. The measured LPO is in good agreement with viscoplastic selfconsistent numerical models. Seismic properties of glaucophane-bearing blueschist and eclogite are calculated by combining measured LPO and the single-crystal elastic moduli of glaucophane with the other main mineral constituents of the rock: mostly epidote for blueschist, and omphacite and garnet for eclogite. Blueschists present stronger anisotropies (AVP=16.1% and AVS=10.3%) than eclogites (AVP=2.9% and AVS=1.7%). The shear-wave splitting and resulting delay times for a 7-km thick layer of eclogite or blueschist are low for the eclogite (b0.03 s), but signi!cant for blueschist (0.16 s). Application to the subducted oceanic crust yields a decrease of velocity contrast with the surrounding mantle and of anisotropy at depth depending on the temperature of the slab. The details of the velocity variations in subducted oceanic crust are dif!cult to reconcile with the blueschist-eclogite transition as probed by exhumed metamorphic rocks and may require additional phases such as lawsonite or chemical variations such as oxidation state.
ISSN: 0040-1951

hal-00564518
https://hal.archives-ouvertes.fr/hal-00564518
DOI : 10.1016/j.tecto.2010.09.011

Éditeur(s) : HAL CCSD Nature Publishing Group
Résumé : International audience
The variation of elastic- wave velocities as a function of the direction of propagation through the Earth's interior is a widely documented phenomenon called seismic anisotropy. The geometry and amount of seismic anisotropy is generally estimated by measuring shearwave splitting, which consists of determining the polarization direction of the fast shear- wave component and the time delay between the fast and slow, orthogonally polarized, waves. In subduction zones, the teleseismic fast shear- wave component is oriented generally parallel to the strike of the trench(1), although a few exceptions have been reported (Cascadia(2) and restricted areas of South America(3,4)). The interpretation of shear- wave splitting above subduction zones has been controversial and none of the inferred models seems to be sufficiently complete to explain the entire range of anisotropic patterns registered worldwide(1). Here we show that the amount and the geometry of seismic anisotropies measured in the forearc regions of subduction zones strongly depend on the preferred orientation of hydrated faults in the subducting oceanic plate. The anisotropy originates from the crystallographic preferred orientation of highly anisotropic hydrous minerals (serpentine and talc) formed along steeply dipping faults and from the larger- scale vertical layering consisting of dry and hydrated crust - mantle sections whose spacing is several times smaller than teleseismic wavelengths. Fault orientations and estimated delay times are consistent with the observed shear- wave splitting patterns in most subduction zones.
ISSN: 0028-0836

hal-00412016
https://hal.archives-ouvertes.fr/hal-00412016
DOI : 10.1038/nature07376

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The seismic anisotropy of the deep Earth is reviewed as a profile from the upper mantle to the solid inner core at the centre of the Earth. The upper mantle is by far the most anisotropic region of the Earth, followed the D” layer above the core mantle boundary. In contrast it is shown that many minerals that are present in the upper mantle, transition zone, lower mantle, D” layer and solid inner core are elastically anisotropic to different degrees. The anisotropy of hydrous phases present in subduction zones is briefly introduced. The basic concepts of single crystal and polycrystalline elasticity and the extrapolation of elastic properties to high temperature and pressure are presented. The specific features of elastic wave propagation in an anisotropic medium of any arbitrary symmetry are illustrated using mantle and inner core phases. The roles of crystal preferred orientation, water and melt in producing seismic anisotropy in the upper mantle are discussed.
Treatese of Geophysics, vol.2

hal-00408321
https://hal.archives-ouvertes.fr/hal-00408321
DOI : 10.1016/B978-044452748-6.00045-6

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
Progress in seismic methodology and ambitious large-scale seismic projects are enabling high-resolution imaging of the continental crust. The ability to constrain interpretations of crustal seismic data is based on laboratory measurements on rock samples and calculations of seismic properties. Seismic velocity calculations and their directional dependence are based on the rock microfabric, which consists of mineral aggregate properties including crystallographic preferred orientation (CPO), grain shape and distribution, grain boundary distribution, and misorientation within grains. Single-mineral elastic constants and density are crucial for predicting seismic velocities, preferably at conditions that span the crust. However, high-temperature and high-pressure elastic constant data are not as common as those determined at standard temperature and pressure (STP; atmospheric conditions). Continental crust has a very diverse mineral composition; however, a select number of minerals appear to dominate seismic properties because of their high-volume fraction contribution. Calculations of microfabric-based seismic properties and anisotropy are performed with averaging methods that in their simplest form takes into account the CPO and modal mineral composition, and corresponding single crystal elastic constants. More complex methods can take into account other microstructural characteristics, including the grain shape and distribution of mineral grains and cracks and pores. Dynamic or active wave propagation schemes have recently been developed, which offer a complementary method to existing static averaging methods generally based on the use of the Christoffel equation. A challenge for the geophysics and rock physics communities is the separation of intrinsic factors affecting seismic anisotropy, due to properties of crystals within a rock and apparent sources due to extrinsic factors like cracks, fractures, and alteration. This is of particular importance when trying to deduce crustal composition and the state of deformation from seismic parameters.
ISSN: 8755-1209

Droits : info:eu-repo/semantics/OpenAccess
hal-01685568
https://hal.archives-ouvertes.fr/hal-01685568
https://hal.archives-ouvertes.fr/hal-01685568/document
https://hal.archives-ouvertes.fr/hal-01685568/file/Almqvist_et_al-2017-Reviewssur Geophysics.pdf
DOI : 10.1002/2016RG000552

Éditeur(s) : HAL CCSD Springer Verlag (Germany)
Résumé : International audience
The anisotropy of vibrational properties influences the acoustic behaviour of wooden pieces and their dependence on grain angle (GA). As most pieces of wood include some GA, either for technological reasons or due to grain deviations inside trunks, predicting its repercussions would be useful. This paper aims at evaluating the variability in the anisotropy of wood vibrational properties and analysing resulting trends as a function of orientation. GA dependence is described by a model based on transformation formulas applied to complex compliances, and literature data on anisotropic vibrational properties are reviewed. Ranges of variability, as well as representative sets of viscoelastic anisotropic parameters, are defined for mean hardwoods and softwoods and for contrasted wood types. GA-dependence calculations are in close agreement with published experimental results and allow comparing the sensitivity of different woods to GA. Calculated trends in damping coefficient (tanδ) and in specific modulus of elasticity (E′/ρ) allow reconstructing the general tanδ-E′/ρ statistical relationships previously reported. Trends for woods with different mechanical parameters merge into a single curve if anisotropic ratios (both elastic and of damping) are correlated between them, and with axial properties, as is indicated by the collected data. On the other hand, varying damping coefficient independently results in parallel curves, which coincide with observations on chemically modified woods, either “artificially”, or by natural extractives.
Wood Science and Technology

hal-01096868
https://hal-agroparistech.archives-ouvertes.fr/hal-01096868
DOI : 10.1007/S00226-010-0393-8

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
A set of synthetic aggregates, containing mixtures of calcite and muscovite, were compacted with varying uniaxial loads ranging from 20 MPa up to 400 MPa. Their elastic properties have been measured using compressional and shear waves in confining hydrostatic pressures up to 475 MPa. Measured seismic velocities are shown to depend on 1) the ratio of calcite to muscovite, 2) the uniaxial load used during sample manufacturing, and 3) the porosity amount prior to velocity measurements. The matrix framework may also affect the seismic velocities, but this effect is not easily quantifiable. In general, measured seismic velocities decrease with increasing muscovite content and porosity. Elastic properties have been calculated based on texture measurements, which were obtained by neutron diffraction. The calculated velocities are based on the calcite and muscovite single crystal elastic tensors, their orientation distribution functions, and their modal fractions. A large discrepancy is apparent between calculated and measured velocities, where the former always overestimate the actual velocities. Shear waves display less of a difference between calculated and measured values, than do compressional waves, indicating that waves propagating by shearing are less affected by sample porosity and matrix framework. A better agreement between measured and calculated seismic velocities arises when the initial porosity is taken into account, using a differential effective medium model. Seismic anisotropy is evident from both measured and calculated velocities, and is low for samples containing pure calcite, but becomes prominent as the muscovite concentration increases. The intensity of anisotropy further depends on the uniaxial load used during sample compaction and the initial porosity.
ISSN: 0148-0227

hal-00512666
https://hal.archives-ouvertes.fr/hal-00512666
DOI : 10.1029/2009JB006523

Éditeur(s) : HAL CCSD Springer-Verlag
Résumé : The seismology is the most effective method to explore the Structure of subduction zones to great depth. The distinguishing feature of the mantle in the subduction regions is the presence of hydrated phases, which transport water into the Earth's interior and release it with dramatic local consequences, triggering earthquakes and melting. The seismological detection of these hydrous phases and geo-dynamic interpretation of flow in the hydrous mantle depend on knowledge of the anisotropic elastic properties and the characteristics of the wave propagation in anisotropic media. In this paper we briefly recall the distinguishing features of anisotropic wave propagation and the observable parameters. We suggest the Vp/Vs ratio is a physically sound parameter than can be observed by seismology in anisotropic regions of the Earth, whereas the Poisson's ratio, which is often quoted, is not directly observable and does not correspond to the characteristics of wave propagation in an anisotropic or isotropic medium. We report for the first time the ratios of Vp/Vs1 and Vp/Vs2, where Vs1 is the fastest and Vs2 slowest S-wave velocities of an anisotropic media. We present the current knowledge of the anisotropic seismic properties of hydrous minerals in the upper mantle, transition zone, and lower mantle that are stable along low temperature geotherms associated with subduction, and identify which minerals are likely to influence seismological observations because they have very high volume fractions, or very high anisotropies, or both of these.;In the upper mantle antigorite and talc are exceptionally anisotropic (Vp 71% Vs 68% and Vp 65%, Vs 68%, respectively) and chlorite is also very anisotropic for S-waves (Vp 35%, Vs 76%). Comparatively less anisotropic are hornblende (Vp 27%, Vs 31%), used as proxy for tremolite in subduction zones, and clinohumite (21.8%, Vs 15.9%). Brucite at 4 GPa (Vp 26.5%, Vs 30.9%) and the dense hydrogen magnesium silicate (DHMS) phase A at 9GPa (Vp 9.3%, Vs 17.6%) are the only hydrous minerals stable in the upper mantle that have had their elastic properties measured at in situ mantle pressures. Except for the phase A, all these minerals are more anisotropic than olivine for at least one parameter (Vp, Vs, Vp/Vs1, or Vp/Vs2). In the transition zone the major phases hydrous wadsleyite and ringwoodite C have moderate to weak anisotropies (Vp 16.3%, Vs 16.5%, and Vp 1.9%, Vs 4.4%, respectively). The DHMS Superhydrous B has a moderate anisotropy (Vp 6.9%, Vs 11.6%). At greater depth the DHMS phase D Lit 24 GPa (Vp 10.8%, Vs 18.0%) is the only hydrous phase that can transport hydrogen from the transition zone into the lower mantle. The phase D is more anisotropic for S than P waves like many of the hydrous phases. From these data it is clear that hydrous phases are in general very anisotropic. However, pressure can play a strong role in reducing anisotropy. It is the case for brucite and talc, in which increasing pressure from ambient to 4 GPa reduces the anisotropy by about 50% for both P and S waves. In contrast, the anisotropy of the DHMS phase A does not change significantly with pressure.;Our picture of the seismic anisotropy of hydrated minerals remains incomplete; the elastic properties of many have not been measured even at ambient conditions (e.g., 10 angstrom phase and phase E) or not measured in their true elastic symmetry (e.g., clinochlore). The majority of hydrous minerals have not been measured at high pressure. and none related to hydrated mantle have been measured at elevated temperature. We have shown in the few cases where hydrous minerals have been measured as a function of pressure, that this variable has an important effect oil the velocity distribution and in most cases reduces the degree of anisotropy, hence we Would expect seismic anisotropy to play a key role in the determination of the shallow structure of subduction zones in the upper mantle.
Subduction zone geodynamics

hal-00420893
https://hal.archives-ouvertes.fr/hal-00420893

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
We determine the elasticity of FeSiO3 perovskite for various spin configurations using density-functional theory calculations. The elastic moduli and the bulk seismic wave velocities are weakly affected by the spin transition. However we show that the intrinsic differences in seismic anisotropy between the high-spin and low-spin phases of Fe-bearing perovskite coupled with lattice preferred orientation that can develop due and during the convection may lead to distinct seismic signatures between the top and the bottom of the lower mantle. These signatures should be detectable in observations and they need to be taken into account in tomographic studies of the Earth's lower mantle.
ISSN: 0094-8276

hal-00512662
https://hal.archives-ouvertes.fr/hal-00512662
DOI : 10.1029/2010GL043320

Éditeur(s) : HAL CCSD
Résumé : Different applications of linear elasticity in structural geology are presented in this thesis through the development of three types of numerical computer codes. The first one uses forward modeling to study displacement and perturbed stress fields around complexly faulted regions. We show that incorporating inequality constraints, such as static Coulomb friction, enables one to explain the angle of initiation of jogs in extensional relays. Adding heterogeneous material properties and optimizations, such as parallelization on multicore architectures and complexity reduction, admits more complex models. The second type deals with inverse modeling, also called parameter estimation. Linear slip inversion on faults with complex geometry, as well as paleo-stress inversion using a geomechanical approach, are developed. The last type of numerical computer code is dedicated to restoration of complexly folded and faulted structures. It is shown that this technique enables one to check balanced cross-sections, and also to retrieve fault chronology. Finally, we show that this code allows one to smooth noisy 3D interpreted faulted and folded horizons using geomechanics.
Différentes applications de l'élasticité linéaire en géologie structurale sont présentées dans cette thèse à travers le développement de trois types de codes numériques. Le premier utilise la modélisation directe pour étudier les déplacements et champs de contraintes autour de zones faillées complexes. On montre que l'ajout de contraintes inégalitaires, telles que la friction de Coulomb, permet d'expliquer l'angle d'initiation des dominos dans les relais extensifs. L'ajout de matériaux hétérogènes et d'optimisations, telles la parallélisation sur processeurs multi-coeurs ainsi que la réduction de complexité des modèles, permettent l'étude de modèles beaucoup plus complexes. Le second type de code numérique utilise la modélisation inverse, aussi appelée estimation de paramètres. L'inversion linéaire de déplacements sur les failles ainsi que la détermination de paléo-contraintes utilisant une approche géomécanique sont développées. Le dernier type de code numérique concerne la restauration de structures complexes plissées et faillées. Il est notamment montré qu'une telle méthode permet de vérifier l'équilibre de coupes géologiques, ainsi que de retrouver la chronologie des failles. Finalement, nous montrons que ce même code permet de lisser des horizons 3D faillés, plissés et bruités en utilisant la géomécanique.
https://tel.archives-ouvertes.fr/tel-00537899

tel-00537899
https://tel.archives-ouvertes.fr/tel-00537899
https://tel.archives-ouvertes.fr/tel-00537899/document
https://tel.archives-ouvertes.fr/tel-00537899/file/thesemaerten.pdf

Éditeur(s) : HAL CCSD E Schweizerbart Science Publishers
Résumé : International audience
We present an electron tomography study of dislocations in deformed olivine single crystals, at ca. 0.5 T m, along two distinct orientations. The easiest slip systems are [001](100) and [001]{110}. Disorientating a single crystal away from easy glide conditions leads to massive cross-slip, which generates three-dimensional dislocations, and thus contributes to hardening. Fast motion of curved non-screw dislocation in those planes leaves long straight screw dislocations which bear lattice friction and control plastic strain. We have identified several hardening mechanisms. Non-screw [001] dislocations interact elastically to form dipoles. Recovery mechanisms leading to dipole annihilation are observed, but they are slow at those temperatures and produce numerous sessile loops. These loops represent obstacles for gliding dislocations. Interactions between dislocations and sessile loops produce sessile segments (super jogs), which efficiently impede dislocation motions.
ISSN: 0935-1221

hal-01277821
https://hal.archives-ouvertes.fr/hal-01277821
DOI : 10.1127/ejm/2015/0027-2481