É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
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 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 Oxford University Press (OUP)
Résumé : Approximately 4 yr of campaign and continuous Global Positioning System (GPS) measurements across the Dead Sea fault system (DSFS) in Lebanon provide direct measurements of interseismic strain accumulation along a 200-km-long restraining bend in this continental transform fault. Late Cenozoic transpression within this restraining bend has maintained more than 3000 m of topography in the Mount Lebanon and Anti-Lebanon ranges. The GPS velocity field indicates 4-5 mm yr(-1) of relative plate motion is transferred through the restraining bend to the northern continuation of the DSFS in northwestern Syria. Near-field GPS velocities are generally parallel to the major, left-lateral strike-slip faults, suggesting that much of the expected convergence across the restraining bend is likely accommodated by different structures beyond the aperture of the GPS network (e.g. offshore Lebanon and, possibly, the Palmyride fold belt in SW Syria). Hence, these geodetic results suggest a partitioning of crustal deformation involving strike-slip displacements in the interior of the restraining bend, and crustal shortening in the outer part of the restraining bend. Within the uncertainties, the GPS-based rates of fault slip compare well with Holocene-averaged estimates of slip along the two principal strike-slip faults: the Yammouneh and Serghaya faults. Of these two faults, more slip occurs on the Yammouneh fault, which constitutes the primary plate boundary structure between the Arabia and Sinai plates. Hence, the Yammouneh fault is the structural linkage that transfers slip to the northern part of the transform in northwestern Syria. From the perspective of the regional earthquake hazard, the Yammouneh fault is presently locked and accumulating interseismic strain.
ISSN: 0956-540X

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

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
The Yakutat-St. Elias collision in SE Alaska and adjacent Canada represents a prime example of present-day tectonics associated with an indentor corner. Its eastern syntaxis is marked by high exhumation, a sharp structural bend, and strain concentration at the transition from shortening to oblique transpression. Here we present GPS velocity and strain rate fields that cover the syntaxis, including 11 new stations in the core of the St. Elias Mountains. These data are corrected for transient deformation (glacial isostatic adjustment and postseismic and interseismic loading) to produce residual velocities and strain rates representative of long-term tectonics. The main features of these velocity and strain rate fields are a peak in strain rates (strain knot) in the syntaxis at the junction between the main fault systems and a rapid rotation from convergence-parallel to convergence-normal orientations of the velocities and shortening axes around the syntaxis, leading to shortening across the southern Denali Fault. These features are consistent with the strain and tectonic patterns expected near an indentor corner at the transition from shortening to transpression, with a combination of diffuse and localized deformation. The GPS velocities and strain rates show diffuse deformation around the syntaxis, from pure convergence-parallel shortening in the orogenic wedge to oblique extension that accommodates the strain rotation at the front of the syntaxis. This indentor-corner model also results in a near-zero strike-slip rate on the southern Denali Fault and shows no clear evidence for a throughgoing fault hypothesized to link the Fairweather and Totschunda Faults.
ISSN: 0148-0227

hal-01243510
https://hal.archives-ouvertes.fr/hal-01243510
DOI : 10.1002/2014JB011842

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We develop a block model for the Aegean and surrounding areas, constrained by Global Positioning System (GPS), in order to investigate the degree of coupling on the Hellenic subduction interface (i.e., the fraction of the motion across the plate boundary accommodated by elastic strain accumulation). We use previously published models, and seismicity to define the geometry of the interface separating the down-going Nubian slab from the overriding Aegean. This model provides a good fit to the GPS observations; for the ∼200,000 km2 Aegean block the wrms of the residual velocities is 1.4 mm/yr for 80 GPS velocity estimates, approximately the 95% level of the GPS velocity uncertainties. We investigate the degree of coupling on the seismically active plate interface, the Hellenic trench splay fault (believed to be the source of the 365 AD Great Crete Earthquake and Tsunami), and the Kephalonia transform fault by comparing the modeled GPS residual velocity field for a range of coupling values. The GPS observations are almost insensitive to coupling on the Kephalonia transform fault, because of the vertical dip of the fault that creates interseismic deformation only close to the fault where few GPS sites exist. The absence of resolvable shortening of the leading edge of the Aegean Plate precludes coupling of more than 0.2 (20% of the full Nubia-Aegean convergence rate) on the modeled plate interface. Because of the shallow dip of the plate interface and trench splay fault, and high rate of convergence, if these boundaries were fully coupled, high elastic strain rates would be expected to extend well into the overriding Aegean plate. Based on our preferred value for the degree of coupling (0.1), and assuming characteristic earthquake behavior, we estimate a recurrence time for great earthquakes with slip similar to that for the 365 Crete event of 5700-8300 yr, consistent with the absence of subsequent great earthquakes on this segment of the subduction zone.
ISSN: 0012-821X

hal-01053486
https://hal.archives-ouvertes.fr/hal-01053486
DOI : 10.1016/j.epsl.2013.10.018

É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
P>The Strait of Hormuz area is a transition zone between the continental collision of the Zagros (west) and the subduction of an oceanic part of the Arabian Plate beneath the Makran wedge (east). Geology and recent GPS measurements indicate that about 15 mm yr-1 of relative motion in N10 degrees E direction is accommodated by two major fault systems: (1) the NNW-trending Minab-Zendan-Palami (MZP) fault system that connects the Main Zagros Thrust (MZT) to the inner Makran thrust system and the Frontal subduction thrust and (2) the N-trending Sabzevaran-Kahnuj-Jiroft (SKJ) fault system that bounds the Jazmurian depression to the west. We use dense GPS measurements along four transects across these fault systems in order to determine the strains spatial distribution. The northern GPS transect confirms the total fault slip rates for both fault systems estimated by the tectonic analyses (about 10 and 7.3 mm yr-1 in N10 degrees direction across the MZP and SKJ fault systems, respectively). For both fault systems, the elastic deformation spreads over shear zones that are several tens of kilometres wide. However, transects located close to latitude 27 degrees N reveal a much narrower shear zone (similar to 10 km) for the MZP fault system. Moreover, we confirm that most of the present-day strain is transferred towards the frontal subduction thrust rather than towards the inner Makran thrusts. In order to complement this new GPS velocity field with spatially dense measurements, we processed a set of ERS radar images by the radar interferometry (InSAR) technique. We used both a 'stacking' and a 'persistant-scatterers' approach to differentiate the ground deformation signal which spatial gradient is expected to be very low, from the atmospheric signal. Results from these interferograms appear to be relatively in agreement with the GPS-determined strain distribution. Nevertheless, they confirm the absence of any superficial creep behaviour since no sharp discontinuity on interferometric phase can be noted on any interferogram. Finally, we use a purely kinematic 'block model' inversion process to calculate slip rates and locking depths for each fault system from our GPS measurements. These models suggest that the relative quiescence over the last 200 yr has certainly produced a slip deficit as high as 2 m. So, we may wonder if the MZP fault system is not late in the interseismic phase of its earthquake cycle.
ISSN: 0956-540X

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