Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
The interpretation of coseismic surface deformation measurements through inversion techniques is of major importance to understand the mechanical behaviour of a seismic fault. Dense geodetic data sets in the vicinity of the ruptured fault provide unique constraints on detailed fault geometry and slip distribution at depth, making them complementary to seismological data. Bam earthquake (Mw 6.6, 2003 December 26) induced surface deformation has been precisely mapped by Envisat ASAR interferometry and by subpixel correlation techniques applied to Spot-5 and ASAR amplitude images. These oblique and horizontal estimations of deformation have been completed with one levelling profile along the main road crossing the rupture from west to east. We process these data (separately and jointly) in a two-step inversion technique, within the elastic half-space theory framework. Our objective is to determine the dislocation model at depth that satisfies simultaneously all the geodetic constraints. Also, we estimate the relative contribution of each geodetic data set to this inversion process.We first use a stochastic direct approach called neighbourhood algorithm in order to estimate the average characteristics of the rupture, and their relative uncertainty. Constraining in this way the geometry of the ruptured fault, we then linearize the inverse problem and compute the slip distribution on the fault using a standard weighted least-square technique, assuming the solution is smooth to some degree. At each step, we discuss the optimal models, their stability as well as the relative influence of each data set on the derived models parameters. Our preferred model reveals a shallow dislocation on a quasi-vertical fault, slightly dipping towards east. The slip vector has a strike-slip component as high as 2 m, while the dip-slip component seems negligible.However, the estimation of the resolution matrices emphasizes the fact that the details of deep fault slip distribution remain out of the scope of this ill-conditioned inverse problem. Yet, our preferred model suggests a main dislocation limited at depth between 1 and 6 km. By contrast, the aftershocks observed in the months following the earthquake are located just beneath the estimated main shock.
ISSN: 0956-540X

hal-00314173
https://hal.archives-ouvertes.fr/hal-00314173
https://hal.archives-ouvertes.fr/hal-00314173/document
https://hal.archives-ouvertes.fr/hal-00314173/file/GJI-2007-169-849-briole.pdf
DOI : 10.1111/j.1365-246X.2007.03358.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 Oxford University Press (OUP)
Résumé : International audience
Ocean loading effects cause 3-D displacements large enough to affect space geodesy measurements either at the subdiurnal periods or at longer time scales by the means of spurious signals. GPS measurements, in turn, could provide local improvements of the models in coastal areas if their ability to precisely monitor such rapid displacements is assessed. In this paper, we use 105 days of continuous GPS measurements collected in 2004 in the French Brittany and Cotentin region to investigate: (1) the precision achieved by the GPS analysis on measuring 3-D subdiurnal displacements and amplitude and phase of the tidal constituents, (2) the quality of the most recent ocean tide model FES2004 in such a complex coastal context. Indeed, in this area, tide amplitudes are among the largest in the world (up to 16 cm of loading displacements on the vertical component) and are believed to show strong shallow-water tides. From a state of the art GPS analysis using the scientific GAMIT software over 2-h sessions, we test two independent strategies for the realization of the reference frame. The position time series are then compared with the displacements predicted by the FES2004 model applied on an elastic Earth model. The two sets of results are consistent with each other at the same level of agreement than with the predicted displacements, namely 3–5 mm on the horizontal components, 10 mm on the vertical. This assesses the capability of this technique for measuring 3-D ocean tide loading deformation. We validate the FES2004 model in the Brittany area, even though it slightly (2–7 mm) underestimates the three components amplitudes. With a harmonic analysis of the observed position time series, we obtain nevertheless an agreement at a submillimetre level for the M2, N2, O1, Q1 tidal constituents and at a millimetre level for the K1 and S2 tidal constituents. Moreover, we can extract a significant M4 load signal at the 95 per cent confidence level from the GPS time series at the stations located in the Mont St Michel area. The detection of other shallow-water constituents would be helpful to understand the amplitude deficit between the FES2004 predicted and GPS observed displacements.
ISSN: 0956-540X

insu-00335909
https://hal-insu.archives-ouvertes.fr/insu-00335909
https://hal-insu.archives-ouvertes.fr/insu-00335909/document
https://hal-insu.archives-ouvertes.fr/insu-00335909/file/Vergnolle_GJI_08_revised.pdf
DOI : 10.1111/j.1365-246X.2008.03734.x

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The Betic-Rif orocline surrounding the Alboran Sea, the westernmost tip of the Mediterranean Sea, accommodates the NW–SE convergence between the Nubia and Eurasia plates. Recent GPS observations indicate a ∼4 mm/yr SW motion of the Rif Mountains, relative to stable Nubia, incompatible with a simple two-plate model. New gravity data acquired in this study define a pronounced negative Bouguer anomaly south of the Rif, interpreted as a ∼40 km-thick crust in a state of non-isostatic equilibrium. We study the correlation between these present-day kinematic and geodynamic processes using a finite-element code to model in 2-D the first-order behavior of a lithosphere affected by a downward normal traction (representing the pull of a high-density body in the upper mantle). We show that intermediate viscosities for the lower crust and uppermost mantle (View the MathML source1021–1022Pas) allow an efficient coupling between the mantle and the base of the brittle crust, thus enabling (1) the conversion of vertical movement, resulting from the downward traction, to horizontal movement and (2) shortening in the brittle upper crust. Our results show that incipient delamination of the Nubian continental lithosphere, linked to slab pull, can explain the present-day abnormal tectonics, contribute to the gravity anomaly observed in northern Morocco, and give insight into recent tectonics in the Western Mediterranean region.
ISSN: 0012-821X

hal-01467130
https://hal.archives-ouvertes.fr/hal-01467130
DOI : 10.1016/j.epsl.2016.08.041

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Tilt fluctuations can potentially reflect the response of hydrosystems to important rainfall. In this context, long baseline tiltmeters have been installed in an underground tunnel penetrating the Fontaine de Vaucluse karst to study the medium deformation related to solicitations exerted by water infiltrating the hydrosystem. The instruments monitor the tilt as well as its spatial variation. Northward tilts reaching a 1 μrad amplitude are observed consecutively to rainfalls. The tilt amplitude is highly correlated with the Fontaine de Vaucluse outlet flow fluctuations. The measured tilt signal is also relatively homogeneous over a 150 m length. Different types of structure likely to produce such observations are tested in order to identify their location with respect to the tiltmeters, their dimension as well as the amount of water level variation in the structure. Following rainfalls, the infiltration of water modifies the pore pressure, inducing a medium deformation. The hypothesis of an homogeneous surface loading on the Vaucluse plateau is first refuted since the related tilt is much lower than the one measured. The water supplied by rainfalls has to accumulate in discontinuities in order to generate a higher tilt. So, the deformation related to a pressure exerted on a fracture filled by water is assessed. A first study reveals the interest of the tilt homogeneity information that constrains strongly the fracture properties. Thus, the fracture must be located at a distance more than a few hundreds metres from the tiltmeters in order to produce a tilt homogeneous in space. If the fracture is initially dry, it must also be filled on a height higher than 150 m consecutive to a rainfall in order to generate a tilt amplitude in the same magnitude as the one measured. Then, we explore the influence of water level variations on the tilt produced by a fracture located at the interface between the saturated and unsaturated zones, which are thereby permanently flooded. Since several parameters of that model satisfactorily explain the field observations, we discuss how simultaneous geodetical observations could provide complementary information that would further constrain the geometry of the structure at the origin of the medium deformation.
ISSN: 0956-540X

insu-01408482
https://hal-insu.archives-ouvertes.fr/insu-01408482
https://hal-insu.archives-ouvertes.fr/insu-01408482/document
https://hal-insu.archives-ouvertes.fr/insu-01408482/file/Lesparre-2017.pdf
DOI : 10.1093/gji/ggw446

Éditeur(s) : HAL CCSD
Résumé : When the upper crust below the mountain ranges shows crustal extension with a direction collinear to the convergence direction, it is generally accepted that it is the result of gravitational collapse. However, recent studies challenge this paradigm by showing that erosion induces uplift and extension in the central part of the low convergent mountain ranges. In this study, we investigate the impact of denudation on the seismotectonic regime of mountain ranges.In order to identify a relationship between seismotectonic regime and erosion, the first part of this work presents a compilation of data in the mountain ranges. Based on these observations, a simple kinematic model is proposed to predict the seismotectonic regime of the study ranges.For low convergence rate and moderate mean elevation ranges, this model predicts an extension regime when the denudation rate is 15% higher than the convergence rate.The second part is devoted to the development of thermomechanical 2D finite element models to study the impact of surface processes on the deformation of the Pyrenees. The results show that the isostatic response to erosion reactivates pre-existing structures. The kinematics of an inherited fault plane can be predicted due to the gradient of the horizontal surface velocity profile. Thus, a plane located in the eroded zone is reactivated in normal fault when near the border of the erosion area, this same plane is reactivated in reverse fault. These results suggest that the current deformation in the North Pyrenean Zone could be the result of surface processes.Given the small number of studies quantifying erosion rates in the Pyrenees, the models developed in the second part suffer from high uncertainty. To remedy this, we study the central Pyrenees through a study that combines two types of data : watershed denudation rates from cosmogenic isotopes concentration, and incision rates from sediments buried in the karst. These results are presented in chapter 3. Denudation profiles obtained are consistent with reactivation of a normal fault plane located in the North Zone of the Central Pyrenees.From our study in the Pyrenees, a good correlation appears between the value of the rate of erosion and the glacier occupation. Given that the glaciers where more developed in the Alps, it raises the question of the impact of the late-würmian deglaciation on the present deformation of this range. A numerical model detailing this relationship is presented in the fourth chapter. Theresults show that deglaciation of the western Alps is controlled by the rheological heterogeneity of the crust. Some of our models predict uplift rates consistent with those highlighted by geodesy.
Lorsque l’extension de la croûte sous les parties hautes des chaînes de montagnes est colinéaire à la direction de convergence, il est traditionnellement admis que le moteur est l’effondrement gravitaire. Pourtant, des études récentes remettent en cause ce paradigme en montrant que l’érosion induit un soulèvement et de l’extension dans la partie centrale des chaînes de montagne à faible taux de convergence. L’objectif de notre étude est d’étudier l’impact de la dénudation de la topographie sur le régime sismo-tectonique des chaînes de montagnes.La première partie de ce travail présente une compilation de données dans les chaînes de montagnes afin de dégager des relations entre régime sismo-tectonique et érosion. Sur la base de ces observations, un modèle cinématique simple permettant de prévoir le régime de la chaîne est proposé. Ainsi, pour les chaînes à faible taux de convergence et d’élévation moyenne, ce modèleprédit de l’extension lorsque le taux de dénudation est 15% plus élevé que le taux de convergence.La deuxième partie est consacrée au développement d’un modèle thermo-mécanique 2D en éléments finis pour étudier l’impact des processus de surface sur la déformation des Pyrénées.Les résultats montrent que la réponse isostatique à l’érosion permet de réactiver des structurespré-existantes. La cinématique d’un plan de faille hérité peut être prédite grâce au gradient du profil des vitesses de surface horizontales. Ainsi, un plan situé dans la zone d’érosion est réactivé en faille normale alors qu’en bordure de cette zone une faille est réactivée en régime inverse. Ces résultats suggèrent que la déformation actuelle des Pyrénées pourrait être la conséquence d’unprocessus d’érosion.Compte tenu du faible nombre d’études quantifiant les taux d’érosion dans les Pyrénées, les modèles développés dans la deuxième partie souffrent d’une forte incertitude. Pour y remédier, nous avons cherché à les quantifier dans les Pyrénées Centrales grâce à une étude qui combine deux types de données : taux de dénudation des bassins versants à partir des isotopes cosmogéniques, et vitesses d’incision à partir des sédiments piégés dans les karsts. Ces résultats sontprésentés dans la troisième partie. Les profils de dénudation obtenus sont compatibles avec unrejeu en faille normale d’un plan situé dans la Zone Nord Pyrénéenne des Pyrénées Centrales.Notre étude dans les Pyrénées, montre une bonne corrélation entre la valeur du taux d’érosion et la cartographie de l’occupation glaciaire. Les glaciers ayant été plus développés dans les Alpes, cela pose la question de l’impact de la déglaciation tardi-Wurmienne sur la déformation actuelle de cette chaîne. Un modèle numérique détaillant cette relation est présenté dans le quatrième chapitre. Les résultats montrent que la déglaciation des Alpes occidentales est contrôléepar l’hétérogénéité rhéologique de la croute. Certains de nos modèles prédisent des vitesses de surrection compatibles avec celles mises en évidence par la géodésie.
https://hal.archives-ouvertes.fr/tel-01262333

tel-01262333
https://hal.archives-ouvertes.fr/tel-01262333
https://hal.archives-ouvertes.fr/tel-01262333/document
https://hal.archives-ouvertes.fr/tel-01262333/file/these-Genti-2015-postjury.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 Elsevier
Résumé : This paper is the first presentation of a project called GHYRAF (Gravity and Hydrology in Africa) devoted to the detailed comparison between models and multidisciplinary observations (ground and satellite gravity, geodesy, hydrology, meteorology) of the variations of water storage in Africa from the Sahara and part to the monsoon equatorial part. We describe the various actions planned in this project. We first detail the actions planned in gravimetry which consist in two main surface gravity experiments: on the one hand the periodic repetition of absolute gravity measurements along a north-south monsoonal gradient of rainfall in West Africa, going from Tamanrasset (20 mm/year) in southern Algeria to Djougou (1200 mm/year) in central Benin; on the other hand the continuous measurements at Djougou (Benin) with a superconducting gravimeter to monitor with a higher sampling rate the gravity changes related to an extreme hydrological cycle. Another section describes the actions planned in GPS which will maintain and develop the present-day existing network in West Africa. The third type of actions deals with hydrology and we review the three sites that will be investigated in this joint hydrogeophysics project namely Wankama (near Niamey) and Bagara (near Diffa) in the Niger Sahelian zone and Nalohou (near Djougou) in the Benin monsoon zone. We also address the question of the ground truth of satellite-derived missions: in this context the GHYRAF project will lead to test the hydrology models by comparison both with in situ and satellite data such as GRACE, as well as to an important increase of our knowledge of the seasonal water cycle in Africa. We finally present preliminary results in GPS based on the analysis of the vertical motion of the Djougou site. The resulting absolute gravity changes related to the 2008 monsoon are finally given.
ISSN: 0264-3707

hal-00475660
https://hal.archives-ouvertes.fr/hal-00475660
DOI : 10.1016/j.jog.2009.09.014

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
No major earthquake occurred in North Chile since the 1877 Mw 8.6 subduction earthquake that produced a huge tsunami. However, geodetic measurements conducted over the last decade in this area show that the upper plate is actually deforming, which reveals some degree of locking on the subduction interface. This accumulation of elastic deformation is likely to be released in a future earthquake. Because of the long elapsed time since 1877 and the rapid accumulation of deformation (thought to be 6-7 cm yr−1), many consider this area is a mature seismic gap where a major earthquake is due and seismic hazard is high. We present a new Global Positioning System (GPS) velocity field, acquired between 2008 and 2012, that describes in some detail the interseismic deformation between 18°S and 24°S. We invert for coupling distribution on the Nazca-South America subduction interface using elastic modelling. Our measurements require that, at these latitudes, 10 to 12 mm yr−1 (i.e. 15 per cent of the whole convergence rate) are accommodated by the clockwise rotation of an Andean block bounded to the East by the subandean fold-and-thrust belt. This reduces the accumulation rate on the subduction interface to 56 mm yr−1 in this area. Coupling variations on the subduction interface both along-strike and along-dip are described. We find that the North Chile seismic gap is segmented in at least two highly locked segments bounded by narrow areas of weak coupling. This coupling segmentation is consistent with our knowledge of the historical ruptures and of the instrumental seismicity of the region. Intersegment zones (Iquique, Mejillones) correlate with high background seismic rate and local tectonic complexities on the upper or downgoing plates. The rupture of either the Paranal or the Loa segment alone could easily produce a Mw 8.0-8.3 rupture, and we propose that the Loa segment (from 22.5°S to 20.8°S) may be the one that ruptured in 1877.
ISSN: 0956-540X

hal-00941279
https://hal.archives-ouvertes.fr/hal-00941279
DOI : 10.1093/gji/ggt183

É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

Éditeur(s) : HAL CCSD European Geosciences Union
Résumé : International audience
Earth deformation is a multi-scale process ranging from seconds (seismic deformation) to millions of years (tectonic deformation). Bridging short- and long-term deformation and developing seismotectonic models has been a challenge in experimental tectonics for more than a century. Since the formulation of Reid's elastic rebound theory 100 years ago, laboratory mechanical models combining frictional and elastic elements have been used to study the dynamics of earthquakes. In the last decade, with the advent of high-resolution monitoring techniques and new rock analogue materials, laboratory earthquake experiments have evolved from simple spring-slider models to scaled analogue models. This evolution was accomplished by advances in seismology and geodesy along with relatively frequent occurrences of large earthquakes in the past decade. This coincidence has significantly increased the quality and quantity of relevant observations in nature and triggered a new understanding of earthquake dynamics. We review here the developments in analogue earthquake modelling with a focus on those seismotectonic scale models that are directly comparable to observational data on short to long timescales. We lay out the basics of analogue modelling, namely scaling, materials and monitoring, as applied in seismotectonic modelling. An overview of applications highlights the contributions of analogue earthquake models in bridging timescales of observations including earthquake statistics, rupture dynamics, ground motion, and seismic-cycle deformation up to seismotectonic evolution.
ISSN: 1869-9510

hal-01553120
https://hal.archives-ouvertes.fr/hal-01553120
https://hal.archives-ouvertes.fr/hal-01553120/document
https://hal.archives-ouvertes.fr/hal-01553120/file/rosenau2017.pdf
DOI : 10.5194/se-8-597-2017

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
A network of 54 survey GPS sites, 28 continuous GPS stations and three absolute gravity (AG) observation sites have been set up in the Alborz mountain range to quantify the present-day kinematics of the range. Our results allow us to accurately estimate the motion of the South Caspian block (SCB) for the first time, and indicate rotation of the SCB relative to Eurasia, accounting for the left lateral motion in the Alborz range. In light of these new results, it clearly appears that deformation rates vary along the range, the eastern part accommodating mainly left lateral strike slip (2 mm yr(-1) south of the range and 5 mm yr(-1) north of the range) with a very low range normal shortening rate on the Khazar thrust fault (similar to 2 mm yr(-1)), and the western part accommodating range normal shortening (similar to 6 mm yr-1) on the Khazar thrust fault with a left lateral component of similar to 2 mm yr(-1) north of the range and 1 mm yr(-1) south of the range. These present-day kinematics agree with geomorphologic estimated slip rates, but not the long-term deformation, corroborating the idea that the kinematics of the range have changed recently due to the change of SCB motion.;Modelling of the interseismic deformation suggests a deep locking depth on the central-western segment of the Khazar fault (similar to 30 km) in agreement with the Baladeh earthquake rupture and aftershock ranging between 10 and 30 km. Given this unusual deep locking depth and the 34 degrees dip of the thrust, a large part of the Alborz range is located above the seismically coupled part of the fault. Based on our AG measurements this part of the range seems to uplift at a rate of 1-5 mm yr(-1), in agreement with terrace uplift.
ISSN: 0956-540X

hal-00617146
https://hal.archives-ouvertes.fr/hal-00617146
https://hal.archives-ouvertes.fr/hal-00617146v2/document
https://hal.archives-ouvertes.fr/hal-00617146/file/djamour2010.pdf
DOI : 10.1111/j.1365-246X.2010.04811.x

Éditeur(s) : HAL CCSD American Geophysical Union (AGU)
Résumé : International audience
Nowadays, technological advances in satellite imagery measurements as well as the development of dense geodetic and seismologic networks allow for a detailed analysis of surface deformation associated with active fault seismic cycle. However, the study of earthquake dynamics faces several limiting factors related to the difficulty to access the deep source of earthquake and to integrate the characteristic time scales of deformation processes that extend from seconds to thousands of years. To overcome part of these limitations and better constrain the role and couplings between kinematic and mechanical parameters, we have developed a new experimental approach allowing for the simulation of strike-slip fault earthquakes and analyze in detail hundreds of successive seismic cycle. Model rheology is made of multilayered visco-elasto-plastic analog materials to account for the mechanical behavior of the upper and lower crust and to allow simulating brittle/ductile coupling, postseismic deformation phase and far-field stress transfers. The kinematic evolution of the model surface is monitored using an optical system, based on subpixel spectral correlation of high-resolution digital images. First, results show that the model succeed in reproducing the deformation mechanisms and surface kinematics associated to the main phases of the seismic cycle indicating that model scaling is satisfactory. These results are comforted by using numerical algorithms to study the strain and stress distribution at the surface and at depth, along the fault plane. Our analog modeling approach appears, then, as an efficient complementary approach to investigate earthquake dynamics.
ISSN: 0278-7407

hal-01174176
https://hal.archives-ouvertes.fr/hal-01174176
DOI : 10.1002/2014TC003701

É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 American Geophysical Union
Résumé : International audience
GPS measurements adjacent to the southern Red Sea and Afar Triple Junction, indicate that the Red Sea Rift bifurcates south of 17 degrees N latitude with one branch following a continuation of the main Red Sea Rift (similar to 150 degrees Az.) and the other oriented more N-S, traversing the Danakil Depression. These two rift branches account for the full Arabia-Nubia relative motion. The partitioning of extension between rift branches varies approximately linearly along strike; north of similar to 16 degrees N latitude, extension (similar to 15 mm/yr) is all on the main Red Sea Rift while at similar to 13 degrees N, extension (similar to 20 mm/yr) has transferred completely to the Danakil Depression. The Danakil Block separates the two rifts and rotates in a counterclockwise sense with respect to Nubia at a present-day rate of 1.9 +/- 0.1 degrees/Myr around a pole located at 17.0 +/- 0.2 degrees N, 39.7 +/- 0.2 degrees E, accommodating extension along the rifts and developing the roughly triangular geometry of the Danakil Depression. Rotating the Danakil Block back in time to close the Danakil Depression, and assuming that the rotation rate with respect to Nubia has been roughly constant, the present width of the Danakil Depression is consistent with initiation of block rotation at 9.3 +/- 4 Ma, approximately coincident with the initiation of ocean spreading in the Gulf of Aden, and a concomitant similar to 70% increase in the rate of Nubia-Arabia relative motion.
ISSN: 0094-8276

hal-00475536
https://hal.archives-ouvertes.fr/hal-00475536
DOI : 10.1029/2009GL041127

Éditeur(s) : HAL CCSD American Association for the Advancement of Science
Résumé : International audience
Large earthquakes produce crustal deformation that can be quantified by geodetic measurements, allowing for the determination of the slip distribution on the fault. We used data from Global Positioning System (GPS) networks in Central Chile to infer the static deformation and the kinematics of the 2010 moment magnitude (M-w) 8.8 Maule megathrust earthquake. From elastic modeling, we found a total rupture length of similar to 500 kilometers where slip (up to 15 meters) concentrated on two main asperities situated on both sides of the epicenter. We found that rupture reached shallow depths, probably extending up to the trench. Resolvable afterslip occurred in regions of low coseismic slip. The low-frequency hypocenter is relocated 40 kilometers southwest of initial estimates. Rupture propagated bilaterally at about 3.1 kilometers per second, with possible but not fully resolved velocity variations.
ISSN: 0036-8075

hal-00783213
https://hal.archives-ouvertes.fr/hal-00783213
DOI : 10.1126/science.1204132

Éditeur(s) : HAL CCSD Elsevier
Résumé : TAG LITHO
In this study we attempt to understand the water storage variations in a karst aquifer on the Larzac Plateau (South of France) using absolute gravimetry. On this karst system, water input is exclusively rainfall and draining occurs at the Durzon perennial spring in a karstic valley. Our basis assumption is that water storage in the vadose zone (epikarst and infiltration zone) has a significant effect on surface gravity through Newtonian attraction. The karst aquifer is hence being monitored since January 2006 with monthly absolute gravity measurements at three sites. The gravity measurements are corrected for regional scale gravity using water storage models in order to obtain gravity variations related exclusively to the local water storage variations. The gravity variations exhibit an important seasonal component (10-15 mu gal corresponding to a variation of equivalent water stab thickness of similar to 24-36 cm). Water storage variation seen by gravimetry is not spatially uniform on the studied karst system, corroborating current geomorphologic observations and interpretations. Finally, we use a global mass balance consideration linking rainfall, evapotranspiration and spring discharge to provide a conceptual framework for the understanding of observed gravity variations.
ISSN: 0022-1694

hal-00411884
https://hal.archives-ouvertes.fr/hal-00411884
DOI : 10.1016/j.jhydrol.2008.06.020