Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We use a strike-slip fault analog model to study experimentally the role played by along-fault non-uniform and asymmetric applied normal stress on both coseismic slip and long-term fault behavior. Our model is based on a visco-elasto-plastic multi-layered rheology that allows to produce several hundreds of scaled analog microquakes and associated seismic cycles. Uniform or heterogeneous applied normal stress along the fault plane is imposed and maintained constant during the whole experiment durations. Our results suggest that coseismic slip patterns are strongly controlled by spatial normal stress variations and subsequent accumulated shear stress along fault strike. Major microquakes occur preferentially in zones of major shear stress asperities. Coseismic slip distributions exhibit a pattern similar to the along-fault applied normal stress distribution. The occurrence of isolated low to moderate microquakes where residual stresses persist around secondary stress asperities, indicates that stress conditions along the fault also control the whole variability of fault slip events. Moreover, when fault slip stability conditions are modulated by normal stress distribution, our experiments suggest that the along-fault stress heterogeneity influences the seismic cycle regularity and, consequently, long-term fault slip behavior. Uniform applied normal stress favors irregular seismic cycles and the occurrence of earthquakes clustering, whereas non-uniform normal stress with a single high amplitude stress asperity generates strong characteristic microquake events with stable return periods. Together our results strengthen the assumption that coseismic slip distribution and earthquake variability along an active fault may provide relevant information on long term tectonic stress and could thus improve seismic hazard assessment.
ISSN: 0012-821X

hal-01666379
https://hal.archives-ouvertes.fr/hal-01666379
DOI : 10.1016/j.epsl.2017.10.009

É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 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

É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 Oxford University Press (OUP)
Résumé : International audience
Iran is an active continental domain accommodating the convergence between the Arabia and Eurasia plates. In northwestern Iran, deformation between the Central Iranian block and the Caucasus domain is mainly accommodated by right lateral strike-slip on the Tabriz fault. Cities and villages, including the city of Tabriz, have been destroyed by several strong historical earthquakes (M similar to 7). In this study, we compare the slip-rates estimated from geodetic measurements (radar interferometry and UPS) with those determined by dating a geomorphological offset of an alluvial fan along the Tabriz fault. The GPS measurements along two profiles normal to the Tabriz fault suggest a slip-rate of 7.3 +/- 1.3 mm yr(-1). The persistent scatterer radar interferometry analysis of Envisat satellite archives from 2003 to 2010 shows a velocity gradient (6 +/- 3 mm yr(-1)) across the Tabriz fault in agreement with GPS results. Moreover, it reveals that most of the area located south of the Tabriz fault is affected by subsidence, and that some sections of the fault probably act as barriers to fluid migration which may have an impact on its mechanical behaviour. West of Tabriz morphotectonic investigations on an alluvial fan surface show a right-lateral cumulative offset of 320 +/- 40 m Luminescence analyses of the coarse matrix alluvial fan deposits provide an age of 46 +/- 3 ka. This yields a slip-rate comprised between 6.5 and 7.3 mm yr(-1) along this segment. These results suggest that the Late Quaternary slip-rate is in agreement with the present-day slip-rate estimated by geodetic measurements, showing no slip-rate changes during the past 45 000 yr. Short-term variations within the 45 000 yr related to temporal earthquake clustering over few seismic cycles cannot be ruled out, but if they exist, they do not affect the geodetic and the geomorphological estimates. This study is in agreement with previous ones suggesting that long-term slip-rates (i.e. averaged over several tens of seismic cycles) are consistent with geodetic estimated slip-rates (i.e. extrapolated from few years of interseismic observations), and suggests that perturbations of fault slip-rates are related to variations over few seismic cycles.
ISSN: 0956-540X

hal-00844695
https://hal.archives-ouvertes.fr/hal-00844695
https://hal.archives-ouvertes.fr/hal-00844695v2/document
https://hal.archives-ouvertes.fr/hal-00844695/file/rizza2013.pdf
DOI : 10.1093/gji/ggt041

Éditeur(s) : HAL CCSD
Résumé : International audience
Optimizing the resource extraction and limiting the environmental constraint motivate to accurately characterize reservoirs. The extraction of resources induces volume change of the reservoir causing displacements in the surrounding medium up to the surface. The resulting surface deformation is usually significant enough to be measurable by geodetic instruments (GPS, satellites, tiltmeters, strainmeters). However, the data needs to be associated to a reservoir model using an inverse scheme to allow the estimation of physical properties of the reservoir. Complete geodetic data including InSAR, GPS, tilt, strain and gravimetric data are uncommon for conventional, unconventional or mining exploitations. In this study, we work on the one site that gathers a wide range of geodetic data: the salt exploitation of Vauvert, France, operated by KemOne. Our approach aims at estimating the volume changes of the cavities created by the extraction combining all the available geodetic data. To do so, we minimize a global functional including weighting factors for each data type. The weighting factors are optimally estimated for our study and can be adapted depending on the issue. Hence, an optimal distribution of data can be suggested for a specific reservoir characterization.
79th EAGE Conference
Paris, France

Droits : http://hal.archives-ouvertes.fr/licences/copyright/
hal-01575426
https://hal.archives-ouvertes.fr/hal-01575426
https://hal.archives-ouvertes.fr/hal-01575426/document
https://hal.archives-ouvertes.fr/hal-01575426/file/reservoir-characterization-inversion.pdf
DOI : 10.3997/2214-4609.201700987

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We present the 1998-2009 GPS-derived velocity field for the Armenia region based on a survey-mode observation network of 31 GPS sites.We combine our results with previous GPS studies of the region to better assess the deformation of the Lesser Caucasus and Kura basin region. The results show that the Kura basin and the Lesser Caucasus regions are two different blocks, and that the main fault (Pambak-Sevan-Sunik) between these two regions has a right-lateral slip rate of 2±1 mm/yr. This is consistentwithmorphotectonic estimates and suggests a fairly constant slip rate over the last 120-300 ka. The right-lateral slip rate on one of the southern branch of the Pambak-Sevan-Sunik fault is lower than 1 mm/yr and the good agreementwith a geologically estimated slip rate suggests a constant slip rate over the last 1.4 Myr. The Sardarapat and Akhurian faults experience some shortening. This shortening is consistentwith someindependent geological estimates and shows the Arabian push.However, NNW-SSE-orientated faults have an extensional fault normal component instead of the expected shortening due to the Arabia-Eurasia convergence. This substantial extensive strain, and the sharp azimuth change of the velocity vectors between the Arabia promontory and the Lesser Caucasus suggest that processes other than "extrusion", possibly related to old subduction or delamination, contribute to the geodynamics of the region.
ISSN: 0040-1951

hal-00812623
https://hal.archives-ouvertes.fr/hal-00812623
DOI : 10.1016/j.tecto.2013.02.002

Éditeur(s) : HAL CCSD Wiley-Blackwell
Résumé : We use InSAR to measure deformation and kinematics of the Mw = 4.9 Borujerd (2005/05/03) and Mw = 6.1 Chalan-Chulan (2006/03/31) earthquakes that occurred in the Zagros fold-and-thrust belt. The focal mechanism of the 2006 event is consistent with right lateral strike-slip motion and the event ruptured the Dorud-Borujerd segment of the Main Recent Fault. An Envisat interferogram spanning the 2006 event shows peak ground deformation of 9 cm in the satellite line-of-sight along a 10 km long fault portion. The interferogram spanning the 2005 earthquake is rather related to atmospheric artefact than to ground deformation. Dislocation models of the 2006 Chalan-Chulan event indicate dextral slip amounting to a maximum of 90 cm at a depth of 4 km. The predicted vertical displacements are in good agreement with differential levelling data. The 2006 event filled only a small part of the seismic gap located between large M = 7 events that occurred in 1909 and 1957.
ISSN: 0954-4879

hal-00411313
https://hal.archives-ouvertes.fr/hal-00411313
DOI : 10.1111/j.1365-3121.2008.00797.x

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : International audience
In this study, we quantify and analyze the spatial and temporal evolution of the surface displacement of Kahrod landslide located in the center of Alborz range (Iran) within the Haraz valley. This landslide represents a threat for this main drainage axis and its numerous infrastructures. We present three sets of displacement vectors based on GPS technique. An 8-benchmark network has been surveyed four times on a 1-year period basis. It provides accurate information on the rate of displacement within the landslide, and addresses the problem of the mechanical resistance of a small hillock, down slope, under the stress imposed by the landslide. Then. this network is densified (57 marks) and measured twice in 6 months using a rapidstatic approach. This yields to a dense description of surface deformation over the whole landslide. Finally, a 1-year time series of permanent GPS recordings is presented and compared to rainfall. Furthermore, we analyze Envisat radar differential interferograms (DInSAR) spanning the same period as permanent GPS. These geodetic data allow to precisely determine the limits of the current sliding zone and to describe the spatial and temporal evolution of surface displacement. The combination of geodesy and field observations leads to a precise description of the past and present kinematics behavior of Kahrod landslide. The chaotic nature of the sliding mass suggests a first catastrophic landslide in a first episode. During the period of observation, the landslide appears to deform quite steadily, and the evidence of short-term correlation between rainfall and deformation amplitude needs to be confirmed by future measurements. Carrying on the acquisition of GPS and InSAR data within the sliding mass but also within adjacent bedrock should give fundamental information with regards to major activation processes (river sapping, water seeping, earthquakes, or failure within the frontal hill of bedrock) and their potential consequences.
ISSN: 1435-9529

hal-00411318
https://hal.archives-ouvertes.fr/hal-00411318
DOI : 10.1016/j.enggeo.2008.02.013

Éditeur(s) : HAL CCSD Oxford University Press (OUP): Policy P - Oxford Open Option A Oxford University Press (OUP)
Résumé : International audience
While the kinematics of Anatolia plate and the North Anatolian Fault System (NAFS) has been studied extensively, the slip rate and locking depth along the NAFS are usually assumed constant in the analyses due to the lack of sufficient data. This is also partly due to the reasonably good fit of Euler small circle and partly due to the lack of spatial resolution of observations to determine slip rates independently from locking depths. On the other hand, recent geodetic studies show a contrast for locking depth between Marmara and other parts of the NAFS, implying a non-uniform locking depth across the NAFS. In this study, we analyse new GPS data and homogenously combine available data sets covering the eastern part of the NAFS to form the most complete data set. In particular, we incorporate the first results of Turkish Real-Time Kinematic GPS Network (CORS-TR) into our data set. A detailed analysis of three profiles within the NAFS reveals an increase of locking depth in the middle profile to 19.1 ± 3.4 km from 11.9 ± 3.5 km in the easternmost profile while the slip rate is nearly constant (20–22 mm yr−1), which implies a variation of strain rate of ∼100 nanostrain yr−1. Assuming a constant locking depth throughout whole NAFS gives an average locking depth of 14.3 ± 1.7 km. Our best estimates of slip rates in block modelling which takes the variation of locking depths into account are in the range between 22.5 and 22.8 mm yr−1 over eastern part of the NAFS.
ISSN: 0956-540X

hal-01218392
https://hal.archives-ouvertes.fr/hal-01218392
DOI : 10.1093/gji/ggv274

É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 American Geophysical Union
Résumé : International audience
In Taiwan, about one third of the lithospheric plate convergence between Eurasia and the Philippine Sea plate is accommodated on the eastern coast across the narrow Longitudinal Valley (LV). The Longitudinal Valley Fault (LVF) is the main seismically active fault zone in this region. However, the spatial distribution of ground deformation due to interseismic loading on locked fault segments remains widely unknown. To address this problem, we use a "permanent scatterer" (PS) interferometric synthetic aperture radar (InSAR) approach for processing the ERS satellite archives. The PS mean velocity maps allow accurately mapping and quantifying the ground deformation all along the LV. Assuming that, close to the surface, the fault plane along the central LVF creeping segment dips 55 degrees to the east, we derive a slip vector at the surface of similar to 25 mm/yr with a mean rake of 70 degrees. In agreement with our PS observations, shallow seismicity indicates a clear change in the fault mechanical behavior of the creeping segment around Fuli town. Finally, we propose a combination of the PS velocity fields for estimating a 3-D map of ground deformation. In the southern part of the LV, the vertical component reveals an uplift of the Coastal Range with respect to the LV of similar to 10 mm/yr. Moreover, between Taitung and Luyeh, the vertical deformation is essentially accommodated along the Luyeh Strand. Finally, north of Rueisuei, the Coastal Range uniformly subsides at similar to 15 mm/yr relative to the western side of the LV. However, a local uplift pattern is detected around latitude 23 degrees 42'.
ISSN: 0148-0227

hal-00595532
https://hal.archives-ouvertes.fr/hal-00595532
DOI : 10.1029/2010JB007898