Éditeur(s) : HAL CCSD American Geophysical Union (AGU)
Résumé : International audience
Observations from deep boreholes at several locations worldwide, laboratory measurements of frictional strength on quartzo-feldspathic materials, and earthquake focal mechanisms indicate that crustal faults are strong (apparent friction μ ≥ 0.6). However, friction experiments on phyllosilicate-rich rocks and some geophysical data have demonstrated that some major faults are considerably weaker. This weakness is commonly considered to be characteristic of mature faults in which rocks are altered by prolonged deformation and fluid-rock interaction (i.e., San Andreas, Zuccale, and Nankai Faults). In contrast, in this study we document fault weakening occurring along a marly shear zone in its infancy (<30 m displacement). Geochemical mass balance calculation and microstructural data show that a massive calcite departure (up to 50 vol %) from the fault rocks facilitated the concentration and reorganization of weak phyllosilicate minerals along the shear surfaces. Friction experiments carried out on intact foliated samples of host marls and fault rocks demonstrated that this structural reorganization lead to a significant fault weakening and that the incipient structure has strength and slip behavior comparable to that of the major weak faults previously documented. These results indicate that some faults, especially those nucleating in lithologies rich of both clays and high-solubility minerals (such as calcite), might experience rapid mineralogical and structural alteration and become weak even in the early stages of their activity.
ISSN: 0278-7407

hal-01211481
https://hal.archives-ouvertes.fr/hal-01211481
DOI : 10.1002/2014TC003716

Éditeur(s) : HAL CCSD
Résumé : This work provides new constraints on the development and on the distribution of the deformation in the upper mantle and particularly beneath transform plate boundaries. USArray experiment and the remarkable increase of the dataset in California for the past ten years allowed us to scrutinize the lateral variations of the anisotropy in the vicinity of the San Andreas Fault zone. We have confirmed and increased the detection of two layers of anisotropy beneath this plate boundary. The first layer, located in the lithosphere, is related to the deformation induced at the fault, and the other one, located in the asthenosphere, is coherent with the anisotropy observed far from it, its origin is however less clear. We show that the deformation zone associated both to the San Andreas, Calaveras and Hayward Faults, is likely 40 km wide at 70 km depth. We then performed numerical thermomechanical modeling (ADELI) of the displacement of a transform plate boundary associated with the computation of the development of crystallographic fabrics using a viscoplastic selfconsistent approach (VPSC). We analyzed the distribution of the deformation in the model ant looked after the possible interactions at depth between deformation caused at surface by the strike-slip dynamic of the fault and the shearing at the base of the lithosphere caused by the horizontal displacement of the plates. Elastic properties derived from the crystallographic fabrics modeled, show that such interactions exist and induce, beneath the fault zone, a progressive rotation of the crystallographic fabrics with depth. Seismological signature of these smooth rotations is however not relevant with the presence of two anisotropic layers as proposed beneath California. We thus consider that a decoupling zone exists between the lithosphere and the asthenosphere beneath the California to account for the sharp separation between a lithospheric and an asthenospheric deformation. We furthermore estimate that anisotropy observed far form the San Andreas Fault in California cannot be explained only by the drag of the asthenosphere by the North America lithosphere as proposed in our article. Indeed, we can only expect few tenths of second of splitting delay from the anisotropic properties derived from the numerical modeling of a plate moving in the same direction and in the same velocity than the North American lithosphere only for 10 Ma of displacement. As delays observed in California rather reach 1.5 s, anisotropy in this region thus requires the existence of an active asthenospheric flow to be explained.
Le travail réalisé pendant cette thèse a permis d'apporter de nouvelles contraintes sur le développement et la distribution de la déformation dans le manteau supérieur et plus particulièrement au niveau des grandes limites de plaques décrochantes. Grâce à l'apport de l'expérience USArray et d'une dizaine d'années d'enregistrements sismologiques supplémentaires, nous avons pu étudier, de manière précise, les variations d'anisotropie dans le voisinage de la Faille de San Andreas. Nous avons confirmé et étendu l'observation de deux couches anisotropes sous cette limite de plaque. On y observe une première couche localisée dans la lithosphère marquant la déformation induite à la limite de plaque, et une autre, asthénosphérique, cohérente avec l'anisotropie observée loin de la faille et d'origine plus discutée. Nous avons montré que la zone de déformation associée aux failles de San Andreas, Calaveras et d'Hayward a, vraisemblablement, une largeur d'au moins 40 kilomètres en base de lithosphère, sous chacune de ces failles. Nous avons ensuite procédé à la modélisation thermomécanique (ADELI) de la migration d'une limite de plaques décrochante couplée à une modélisation du développement de fabriques cristallographiques par une approche viscoplastique auto-cohérente (VPSC). Ceci nous a permis d'y observer le développement de la déformation et les conséquences des possibles interactions entre la déformation décrochante en surface et le cisaillement en base de lithosphère dû au déplacement horizontal des plaques. Les propriétés élastiques déduites des fabriques cristallographiques modélisées montrent que de telles interactions existent et provoquent, sous la limite de plaques, une rotation des orientations cristallographiques avec la profondeur. Le signal associé à ces rotations progressives n'est toutefois pas cohérent avec la présence de deux couches d'anisotropie comme proposée sous la faille de San Andreas. Nous pensons par conséquent qu'il existe, sous la Californie, une zone de découplage entre la lithosphère et l'asthénosphère, permettant d'individualiser une déformation lithosphérique d'une déformation asthénosphérique. Nous estimons, en outre, que l'anisotropie observée dans l'asthénosphère sous la Californie ne peut être expliquée seulement par le cisaillement induit par le déplacement de la lithosphère Nord Amérique. En effet, les propriétés anisotropes obtenues par modélisation à partir d'une plaque se déplaçant dans une direction et une vitesse proche de celle de la plaque Amérique du Nord montrent qu'on ne peut espérer guère plus que quelques dixièmes de seconde de délai au bout de 10 Ma de déplacement. Les déphasages mesurés en Californie étant de l'ordre de 1,5 s, il est donc nécessaire d'invoquer la présence d'écoulements mantelliques actifs sous cette région.
https://tel.archives-ouvertes.fr/tel-00767124

tel-00767124
https://tel.archives-ouvertes.fr/tel-00767124
https://tel.archives-ouvertes.fr/tel-00767124/document
https://tel.archives-ouvertes.fr/tel-00767124/file/theseBonnin.pdf

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
While vertical motion induced by long-term geological loads is often used to estimate the flexural rigidity of the lithosphere, we intend to evaluate the shear rigidity of the lithosphere using horizontal motion. Our approach considers that the rigidity of the lithosphere may be defined as its resistance to horizontal tectonic lateral forces. In this case, a spatial distribution of an effective shear rigidity can be estimated from the analysis of the interseismic velocity fields. We consider the Western United States zone where weakly strained areas (e.g., the Sierra Nevada) are connected with areas of large strain rate (e.g. San Andreas Fault system). By inverting interseismic strain distribution measured by geodetic methods, we infer the effective shear rigidity of the lithosphere. The forward problem is defined using the equations of linear elasticity. The inversion relies on the minimization of the sum of a quadratic measure of the differences between measured and modelled velocity fields. The functional also includes regularization terms for the parameters of the model. The gradient of the functional with respect to the minimization parameters is computed using an adjoint formulation. This permits the treatment of large dimensional minimization problems. Finally, a measure of the uncertainty of our inversion is illustrated through the covariance matrix of the parameters at the optimum. The optimization chart is validated on two synthetic velocity distributions. Then, the effective shear rigidity variations of the Western United States are estimated using the CMM3 interseismic velocities. The inversion displays low effective rigidities along the San Andreas Fault system, the Mojave Desert and in the Eastern California Shear Zone, while rigid areas are found in the Sierra Nevada and in the South Basin and Range. Finally, we discuss the differences between our strain rate and rigidity maps with previously published results for the Western United States.
ISSN: 0040-1951

Droits : http://hal.archives-ouvertes.fr/licences/copyright/
hal-01505024
https://hal.archives-ouvertes.fr/hal-01505024
https://hal.archives-ouvertes.fr/hal-01505024/document
https://hal.archives-ouvertes.fr/hal-01505024/file/Paper_tectophysics.pdf
DOI : 10.1016/j.tecto.2017.03.015

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
In order to constrain the vertical and lateral extent of deformation and the interactions between lithosphere and asthenosphere in a context of a transpressional plate boundary, we performed teleseismic shear wave splitting measurements for 65 permanent and temporary broadband stations in central California. We present evidence for the presence of two anisotropic domains: (1) one with clear E–W trending fast directions and delay times in the range 1.5 to 2.0 s and (2) the other closely associated with the San Andreas Fault system with large azimuthal variations of the splitting parameters that can be modeled by two anisotropic layers. The upper of the two layers provides fast directions close to the strike of the main Californian faults and averaged delay times of 0.7 s; the lower layers show E–W directions and delay times in the range 1.5 to 2.5 s and thus can be compared to what is observed in stations that require a single layer. We propose the E–W trending anisotropic layer to be a 150 to 200 km thick asthenospheric layer explained by the shearing associated with the absolute plate motion of the North American lithosphere. The shallower anisotropic layer ought to be related to the dynamics of the San Andreas Fault system and thus characterized by a vertical foliation with lineation parallel to the strike of the faults localized in the lithosphere. We also propose that the anisotropic layer associated with each fault of the San Andreas Fault system is about 40 km wide at the base of the lithosphere.
ISSN: 0148-0227

hal-01236161
http://hal.univ-reunion.fr/hal-01236161
http://hal.univ-reunion.fr/hal-01236161/document
http://hal.univ-reunion.fr/hal-01236161/file/Bonnin_SKS_California_JGR2010.pdf
DOI : 10.1029/2009JB006438

Éditeur(s) : HAL CCSD Elsevier
Résumé : We present field data from segmented normal faults having particular displacement and overlapping geometries that may be related to down-dip branching of the fault segments. Based on a 3-D numerical modeling study of computed displacement and stress fields around different geometries of down-dip branched normal fault segments, we show that the bends of fault surfaces that coalesce at depth exert a significant influence on (I) the displacement distribution on the fault segment branches and (2) the quasi-static stress fields around the relay zones. The asymmetry in the displacement profiles and low fault interaction at relay zones modeled are consistent with the fault segment geometries observed in field data. As an application, we model the geometry and interaction of the Vallo di Diano normal fault which is segmented at the Earth's surface and which produced an earthquake of magnitude 6.4 along a single fault surface at depth. Numerical modeling of the segmented fault produces an asymmetric displacement profile and low shear stress at the relay zone consistent with the profile and fault interaction inferred from the held. We conclude that asymmetry of the displacement profile and large overlap/separation ratio of the unlinked relay zone can be indicators of the presence of a continuous fault surface underneath. These geometrical attributes are therefore important to consider in the probabilistic analysis of seismic hazard along segmented normal faults.
ISSN: 0191-8141

hal-00411161
https://hal.archives-ouvertes.fr/hal-00411161
DOI : 10.1016/j.jsg.2008.07.005

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
I use GPS interseismic velocities and classic 2D elastic half-space models with a screw dislocation to estimate the long-term fault slip rate, locking depth, and the offset between the surface fault trace and the location of the dislocation below the seismogenic zone for 13 segments along 8 major strike-slip faults. Using deduced strike-slip rates and the position of the dislocation to normalize the interseismic velocities to facilitate comparison of spatial patterns of deformation, I show that no substantial differences can be detected, ruling out a large asymmetry in interseismic velocities across the 8 faults used in this study. Only the Carrizo Plain segment of the San Andreas Fault shows a significant asymmetry that cannot be explained by shifting the position of the dislocation at depth relative to the fault trace. However, the resulting perturbation is less than 10% of total strike-slip rate. Fault traces are usually curved, defining a concave side. When the dislocation at depth is significantly offset from the fault trace, the shift is always toward the block on the concave side of the fault trace. This suggests that the fault zone in the lower crust may develop a simpler geometry more consistent with relative motion across the fault than its upper seismogenic part constrained by the structural complexity of the brittle crust. Since the faults used in this study are at different times in their interseismic period, comparing the interseismic velocity fields across them allows identification of possible variations of the interseismic velocities with time. When normalized by slip rate and dislocation location, all the faults show the same interseismic strain with no significant differences between deduced locking depths. These comparisons suggest that if temporal variations occur as suggested by some dynamic earthquake cycle models, they are small and below the accuracy of the available geodetic measurements.
ISSN: 0040-1951

hal-01171753
https://hal.archives-ouvertes.fr/hal-01171753
DOI : 10.1016/j.tecto.2015.01.013

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We use 3-D mechanical modeling representing faults as planar surfaces with frictional properties that obey Coulomb-failure process to explore the long-term slip rates of the Altyn Tagh fault and Kunlun faults in the north Tibetan plateau. Crustal theology is simplified as an elastoplastic upper crust and a viscoelastic lower crust. Far-field GPS velocities and late Quaternary fault slip rates are used to constrain the model results. Rheological tests show that effective fault friction lower than 0.1-0.08 leads to high slip rates that fit with geologically and geodetically determined slip rates of the Kunlun fault (10-11.7 +/- 1.5 mm/yr). Meanwhile, the modeled Altyn Tagh fault reaches slip rates similar to 13.7 mm/yr to similar to 7.8 mm/yr in its central portion, between ranges of the geological slip rates. Associated with high slip rates, our model predicts that central Tibet (similar to 84 degrees E-95 degrees E) from the Altyn Tagh fault to the north of the Himalayan arc accommodates north-south shortening and east-west extension rates of similar to 10-12 mm/yr and similar to 8-10 mm/yr, respectively. We also question the widely accepted idea that interseismic strain is driven at the base of the seismogenic zone by a screw dislocation. If this assumption fails, the presented model implies that interseismic strain around large strike-slip faults could be distributed in a much broader way if the lithosphere deforms as a thin elastic plate rather than an elastic half-space with an embedded dislocation. If this distributed deformation is ignored, and the instantaneous surface deformation field modeled as that resulting from slip on a dislocation below a specified depth embedded in an elastic half-space, the estimated slip rate will inevitably be lower than the true long-term slip rate. This appears to explain why geodetic slip rates proposed for the Altyn Tagh fault (510 mm/yr) are lower than some of the geological slip rates. (C) 2008 Elsevier B.V. All rights reserved.
ISSN: 0012-821X

hal-00411926
https://hal.archives-ouvertes.fr/hal-00411926
DOI : 10.1016/j.epsl.2008.06.049

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
We performed forward modelling of seismic anisotropy beneath a migrating strike-slip plate boundary to: (1) test if such a geodynamic context might explain teleseismic shear wave splitting data in the vicinity of the central part of the San Andreas Fault and (2) constrain the power of such data to unravel vertical and lateral variations in deformation patterns in the upper mantle. The modelling involves five steps: (1) thermomechanical modelling, using a finite-element code, of the deformation field, (2) viscoplastic self-consistent modelling of the resulting olivine and pyroxene crystal preferred orientations, (3) calculation of the elasticity tensors for different domains of the finite elements (FEs) model, (4) forward modelling of seismic wave propagation through the model using ray theory, finite-frequency theory and a full-wave approach and (5) performing splitting measurements on the synthetic seismograms. SKS splitting data in central California are best fitted by a model with a hotter geotherm within 60 km of the plate boundary accounting for the opening of an asthenospheric window due to the northward migration of the Mendocino Triple Junction. The westward motion of the plate boundary cannot however explain the rotation of fast polarization directions east of the San Andreas Fault in central California. Comparison between modelled and measured individual shear wave splitting also implies that the homogeneity of the two-layer models accounting for the observations in the vicinity of the San Andreas Fault indicates a sharp transition between lithospheric and asthenospheric deformations beneath this plate boundary. The ability of different synthetic approaches to localize horizontally and vertically the plate boundary-related deformation differs significantly. Splitting in ray theory synthetics closely follow variations in olivine crystallographic preferred orientations in the model. In contrast, splitting analysis on full-wave synthetics, which should be more representative of actual long-period SKS waves, results in smooth apparent lateral variations of anisotropy; the exact location and width of the plate boundary may only be retrieved by comparing fast-polarization profiles obtained using a multichannel analysis on waves with different periods.
ISSN: 0956-540X

hal-01236150
http://hal.univ-reunion.fr/hal-01236150
http://hal.univ-reunion.fr/hal-01236150/document
http://hal.univ-reunion.fr/hal-01236150/file/Bonnin_San_Andreas_modelling_GJI2012.pdf
DOI : 10.1111/j.1365-246X.2012.05650.x

É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 AGU and the Geochemical Society
Résumé : International audience
Analysis of short-term deformation along the southern part of Dead Sea Fault (DSF) provides a systematic view of kinematics this part of the continental transform. The southern DSF consists of two principal segments: the Wadi Araba and Jordan Valley faults. In addition to other regional continuous GPS data, this study uses new data from 25 survey sites and 4 continuous GPS stations in Jordan for improved near-field observations. Resulting velocities are reported with 1-sigma uncertainties ranging from 0.4-1.0 mm/yr. Application of elastic dislocation models yields estimates of slip rates for Wadi Araba and Jordan Valley faults are 4.9 +/- 0.4 mm/yr and 4.7 +/- 0.4 mm/yr, respectively. Modeling also suggests different depths of effective fault locking with 15 +/- 5 km and 8 +/- 5 km for the Wadi Araba and Jordan Valley faults, respectively. These slip rates are generally consistent with the upper end of the range of slip rates estimated from late Quaternary geology. Spatial variations in effective fault locking generally correspond with a heterogeneous mantle lithosphere. A similar observation can be observed along the southern San Andreas Fault, and this may reflect the influence of heterogeneity in the uppermost mantle on crustal faulting processes.
ISSN: 1525-2027

hal-00689031
https://hal.archives-ouvertes.fr/hal-00689031
DOI : 10.1029/2011GC003736

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : This paper presents an inversion method using the interseismic velocity field to determine effective rigidity of the lithosphere. The method is based on the minimization of a cost function defined as the quadratic measure of the difference between measured and modelled velocity fields on a discrete set of points. The continuous mapping of the rigidity is fulfilled with a limited set of parameters and the forward solution is achieved using a plane stress finite element code. The computation of the cost function gradient in the parameters' space allows one to iteratively find the best parameters set through a suitable optimization algorithm.;We first design a benchmark including an abrupt rigidity variation that cannot be described by a continuous function. For such a case, we show that increasing the number of parameters is a way to accurately describe sharp variations of the rigidity map. Then, we use a dense GPS velocity field over the southwestern United States to estimate the corresponding rigidity variations for different spatial resolutions of the parameters' grid. We analyse the conceptual and practical difficulties associated with our methodology. Finally, rigidity maps obtained by our inversion method in southwestern United States and particularly across the San Andreas Fault System are reviewed and compared to current plate rigidity estimates and geophysical data over this area.
ISSN: 0956-540X

hal-00669897
https://hal.archives-ouvertes.fr/hal-00669897
DOI : 10.1111/j.1365-246X.2011.05192.x