Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
The seismic behavior of the Himalayan arc between central Nepal and Arunachal Pradesh remains poorly understood due to the lack of observations concerning the timing and size of past major and great earthquakes in Bhutan. We present here the first paleoseismic study along the Himalayan topographic front conducted at two sites in southern central Bhutan. Paleoseismological excavations and related OxCal modeling reveal that Bhutan experienced at least two great earthquakes in the last millennium: one between the seventeenth and eighteenth century and one during medieval times, producing a total cumulative vertical offset greater than 10 m. Along with previous studies that reported similar medieval events in Central Nepal, Sikkim, and Assam, our investigations support the occurrence of either (i) a series of great earthquakes between A.D. 1025 and A.D. 1520 or (ii) a single giant earthquake between A.D. 1090 and A.D. 1145. In the latter case, the surface rupture may have reached a total length of ~800 km and could be associated with an earthquake of magnitude Mw = 8.7–9.1.
ISSN: 0148-0227

hal-01622276
https://hal.archives-ouvertes.fr/hal-01622276
DOI : 10.1002/2015JB012733

Éditeur(s) : HAL CCSD
Résumé : Ocean Drilling Program Hole 1256D reached for the first time the transition zone between the sheeted dike complex and the uppermost gabbros. The recovered crustal section offers a unique opportunity to study the deepest part of the hydrothermal system in present-day oceanic crust. We present a structural analysis of electrical borehole wall images. We identified, and measured the orientations of four categories of structures: major faults, minor fractures, possibly hydrothermal veins, and dikes. All structures tend to strike parallel to the paleo-ridge axis. Three major fault zones (meter thick) and dikes are steeply dipping (~ 75° on average) outward the ridge. Centimeter-thick moderately conductive planar features are interpreted as hydrothermal veins, are organized in arrays of consistent spacing, thickness, and orientation, and are dipping about 15-20° toward the ridge. This structural pattern is interpreted as an on-axis paleohydrothermal circulation system, with vertical, dike-parallel fractures, and sub-horizontal high-temperature hydrothermal veins at the base of the sheeted dike, which was subsequently rotated ~ 15° westward around a ridge-parallel, sub-horizontal axis. This rotation can be caused by upper-crustal block rotation along a listric normal fault, and/or subsidence at the ridge axis.
Ofioliti

hal-00757592
https://hal.archives-ouvertes.fr/hal-00757592
DOI : 10.4454/ofioliti.v37i1.402

Éditeur(s) : HAL CCSD Springer Verlag (Germany) Springer Verlag
Résumé : Direct numerical simulations of flow through two millimeter-scale rock samples of limestone and sandstone are performed using three diverse fluid dynamic simulators. The resulting steady-state velocity fields are compared in terms of the associated empirical probability density functions (PDFs) and key statistics of the velocity fields. The pore space geometry of each sample is imaged at 5.06−μm voxel size resolution using X-ray microtomography. The samples offer contrasting characteristics in terms of total connected porosity (about 0.31 for the limestone and 0.07 for the sandstone) and are typical of several applications in hydrogeology and petroleum engineering. The three-dimensional fluid velocity fields within the explicit pore spaces are simulated using ANSYS® FLUENT® ANSYS Inc. (2009), EULAG Prusa et al. (Comput. Fluids 37, 1193–1207 2008), and SSTOKES Sarkar et al. (2002). These computational approaches are highly disperse in terms of algorithmic complexity, differ in terms of their governing equations, the adopted numerical methodologies, the enforcement of internal no-slip boundary conditions at the fluid-solid interface, and the computational mesh structure. As metrics of comparison to probe in a statistical sense the internal similarities/differences across sample populations of velocities obtained through the computational systems, we consider (i) integral quantities, such as the Darcy flux and (ii) main statistical moments of local velocity distributions including local correlations between velocity fields. Comparison of simulation results indicates that mutually consistent estimates of the state of flow are obtained in the analyzed samples of natural pore spaces despite the considerable differences associated with the three computational approaches. We note that in the higher porosity limestone sample, the structures of the velocity fields obtained using ANSYS FLUENT and EULAG are more alike than either compared against the results obtained using SSTOKES. In the low-porosity sample, the structures of the velocity fields obtained by EULAG and SSTOKES are more similar than either is to the fields obtained using ANSYS FLUENT. With respect to macroscopic quantities, ANSYS FLUENT and SSTOKES provide similar results in terms of the average vertical velocity for both of the complex microscale geometries considered, while EULAG tends to render the largest velocity values. The influence of the pore space structure on fluid velocity field characteristics is also discussed.
ISSN: 1420-0597

hal-01174152
https://hal.archives-ouvertes.fr/hal-01174152
DOI : 10.1007/s10596-015-9486-7

É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 Oxford University Press (OUP)
Résumé : International audience
Advances in groundwater storage monitoring are crucial for water resource management and hydrological processes understanding. The evaluation of water storage changes (WSC) often involve point measurements (observation wells, moisture probes, etc.), which may be inappropriate in heterogeneous media. Over the past few years, there has been an increasing interest in the use of gravimetry for hydrological studies. In the framework of the GHYRAF (Gravity and Hydrology inAfrica) project, 3 yr of repeated absolute gravity measurements using a FG5-type gravimeter have been undertaken at Nalohou, a Sudanian site in northern Benin. Hydrological data are collected within the long-term observing system AMMA-Catch. Once corrected for solid earth tides, ocean loading, air pressure effects, polar motion contribution and non-local hydrology, seasonal gravity variations reach up to 11 μGal, equivalent to a WSC of 260-mm thick infinite layer of water. Absolute temporal gravity data are compared to WSC deduced from neutron probe and water-table variations through a direct modelling approach. First, we use neutronic measurements available for the whole vertical profile where WSC occur (the vadose zone and a shallow unconfined aquifer). The RMSD between observed and modelled gravity variations is 1.61 μGal, which falls within the error bars of the absolute gravity data. Second, to acknowledge for the spatial variability of aquifer properties, we use a 2-D model for specific yield (Sy) derived from resistivity mapping and Magnetic Resonance Soundings (MRS). The latter provides a water content (θMRS) known to be higher than the specific yield. Hence, we scaled the 2-D model of θMRS with a single factor (α). WSC are calculated from water-table monitoring in the aquifer layer and neutronic measurements in the vadose layer. The value of α is obtained with aMonte-Carlo sampling approach, minimizing the RMSD between modelled and observed gravity variations. This leads to α = Sy/θMRS = 0.63 ± 0.15, close to what is found in the literature on the basis of pumping tests experiments, with a RMSD value of 0.94 μGal. This hydrogeophysical experiment is a first step towards the use of time-lapse gravity data as an integrative tool to monitor interannual WSC even in complicated subsurface distribution.
ISSN: 0956-540X

hal-00821181
https://hal.archives-ouvertes.fr/hal-00821181
DOI : 10.1093/gji/ggt146

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : International audience
Over the years, a number of different models and techniques have been proposed to both quantify and explain the glacial isostatic adjustment (GIA) process. There are serious challenges, however, to obtaining accurate results from measurements, due to noise in the data and the long periods of time necessary to identify the relatively small-magnitude signal in certain regions. The primary difficulty, in general, is that most of the geophysical signals that occur in addition to GIA are nonstationary in nature. These signals are also corrupted by random as well as correlated noise added during data acquisition. The nonstationary characteristic of the data makes it difficult for traditional frequency-domain denoising approaches to be effective. Time-frequency filters present a more robust and reliable alternative to deal with this problem. This paper proposes an extended S transform filtering approach to separate the various signals and isolate that associated with GIA. Continuous global positioning system (GPS) data from eastern Canada for the period from June 2001 to June 2006 are analyzed here, and the vertical velocities computed after filtering are consistent with the GIA models put forward by other researchers.
ISSN: 0033-4553

hal-00795451
https://hal.archives-ouvertes.fr/hal-00795451
DOI : 10.1007/s00024-011-0404-1

Éditeur(s) : HAL CCSD
Résumé : In SW Yukon – SE Alaska, the boundary between the Pacific and North America plates is characterized by a syntaxis at the transition between the Aleutian subduction to the W and the Fairweather – Queen Charlotte strike-slip faults to the SE. The relative motion is oblique to the main fault structures, and the area is marked by the Yakutat block collision. From the Chugach – Saint Elias mountains in the plate boundary zone (up to 6 000 m high) to the intraplate strike-slip faults, markers of the present-day deformation give information on its partition in the system.During my PhD, I first measure surface deformation using a dense GPS network, deployed up to 500 km inland the North America plate. After precise processing and corrections of transient effects in the area (postseismic and glacial isostatic rebound), a new residual velocity field is produced for the syntaxis area, from which I derive strain rates. Those data allow me to quantify the fault slip rates for Fairweather and southern Denali strike-slip faults, and to characterize a bi-modal deformation pattern: Along the plate boundary, the deformation is localized on large-scale structures (accretionary prism to the W, Fairweather to the E); In the syntaxis area, strain rates are the highest and the GPS data shows a diffuse intraplate deformation, similar to an indentor pattern. The Yakutat block seems to strongly drive the North America plate deformation.This indentor pattern induces strong lateral variations on the large intraplate faults: the Denali – Totschunda – Duke River system. In a second part, I realize a regional geomorphological study to characterize the role and slip rate of those faults. From very high-resolution Digital Elevation Models (~ 1 m), a detailed cartography is done. On the basis of fieldwork observations, I measure offsets of fluvial and glacial markers, which are sampled for dating. A dextral cumulative deformation is highlighted on the northern Denali Fault, where as all southern Denali is marked by vertical deformation. This study allows me to quantify new slip rates for the system Denali – Totschunda – Duke River, and to show the leading role of the Totschunda (~ 14 mm/a) and Duke River (~ 6 mm/a) faults, contrary to the Denali Fault (~ 1 mm/a) North of the syntaxis.The new tectonic model for the Yakutat collision provides an important case study for the understanding of indentor systems. The concomitance of a rigid-block deformation (to the West) and diffuse deformation (to the East), as well as the near-zero slip rate on the lithospheric-scale southern Denali Fault highlight the major control of boundary conditions and the structural heritage on the orogen deformation.
Au SO Yukon – SE Alaska, la frontière entre les plaques Pacifique et Nord-Amérique est marquée par une syntaxe à la transition entre la subduction des Aléoutiennes à l’O et le décrochement de Fairweather – Queen Charlotte au SE. Le mouvement relatif est oblique et la région est marquée par la collision du bloc Yakutat. De la chaine des Chugach – Saint Elias en frontière de plaque (jusqu’à 6 000 m d’altitude) aux grands décrochements intraplaque, les marqueurs de la déformation actuelle apportent des informations sur sa partition dans ce système.Au cours de cette thèse, je m’intéresse dans un premier temps à mesurer la déformation de surface à l’aide d’un réseau GPS dense, déployé jusqu’à 500 km à l’intérieur de la plaque Nord-Amérique. Après un travail minutieux de traitement des données et de correction des effets transitoires (rebond post-glaciaire et post-sismique), un champ de vitesses résiduelles robuste et inédit est produit au niveau de la syntaxe, dont je dérive des taux de déformation. Ces données me permettent de quantifier les vitesses de glissement les failles décrochantes de Fairweather et Denali au Sud, et également de caractériser une déformation bimodale : En frontière de plaques, la déformation est localisée sur les grandes structures (prisme d’accrétion à l’O, Fairweather à l’E) ; Au niveau de la syntaxe, les taux de déformation sont les plus importants et les données GPS mettent en évidence une déformation intraplaque diffuse, similaire à un champ attendu à l’aplomb d’un indenteur. Le bloc Yakutat semble donc contrôler fortement la déformation de la plaque Nord-Amérique.Ce champ de déformation induit des variations latérales fortes sur les grandes failles intraplaque : le système Denali-Totschunda-Duke River. Dans une seconde partie, je réalise une étude géomorphologique régionale pour caractériser le rôle et la vitesse de ces failles. À partir de Modèles Numériques de Terrain très haute résolution (~ 1 m), une cartographie de détail est réalisée. Puis une mission de terrain me permet de mesurer des décalages sur des marqueurs fluviatiles et glaciaires et collecter des échantillons pour les dater. Je mets en évidence une déformation cumulée dextre sur le segment Nord de la faille de Denali, alors que toute la partie Sud déforme verticalement la surface. Cette étude me permet de quantifier de nouvelles vitesses de failles sur le système Denali-Totschunda-Duke River, et de montrer le rôle prépondérant des failles de Totschunda (~ 14 mm/a) et Duke River (~ 6 mm/a) contrairement à la faille de Denali (~ 1 mm/a) au Nord de la syntaxe.Le nouveau modèle cinématique proposé pour la région de collision Yakutat permet d’apporter un nouvel exemple à la compréhension des systèmes d’indenteur. La concomitance d’une déformation de blocs (à l’Ouest) et diffuse (à l’Est), ainsi que l’absence de déformation sur la faille lithosphérique de Denali Sud met en évidence le contrôle majeur des conditions aux limites et de l’héritage structural dans la déformation des orogènes.
https://hal.archives-ouvertes.fr/tel-01241688

tel-01241688
https://hal.archives-ouvertes.fr/tel-01241688
https://hal.archives-ouvertes.fr/tel-01241688/document
https://hal.archives-ouvertes.fr/tel-01241688/file/manuscrit_anais_v2_optimise.pdf

É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 Nature Publishing Group
Résumé : International audience
Collisional mountain belts grow as a consequence of continental plate convergence and eventuallydisappear under the combined effects of gravitational collapse and erosion. Using a decade ofGPS data, we show that the western Alps are currently characterized by zero horizontal velocityboundary conditions, offering the opportunity to investigate orogen evolution at the time ofcessation of plate convergence. We find no significant horizontal motion within the belt, but GPS andlevelling measurements independently show a regional pattern of uplift reaching ~2.5 mm/yr in thenorthwestern Alps. Unless a low viscosity crustal root under the northwestern Alps locally enhancesthe vertical response to surface unloading, the summed effects of isostatic responses to erosion andglaciation explain at most 60% of the observed uplift rates. Rock-uplift rates corrected from transientglacial isostatic adjustment contributions likely exceed erosion rates in the northwestern Alps. In theabsence of active convergence, the observed surface uplift must result from deep-seated processes.
ISSN: 2045-2322

hal-01357450
https://hal.archives-ouvertes.fr/hal-01357450
https://hal.archives-ouvertes.fr/hal-01357450/document
https://hal.archives-ouvertes.fr/hal-01357450/file/noquet2016.pdf
DOI : 10.1038/srep28404

É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 Springer Verlag
Résumé : International audience
Continuous Global Positioning System (CGPS) position time series from eastern North America constrain the pattern and magnitude of regional crustal deformation. Initial analysis delineates consistent uplift patterns, as expected from glacial isostatic adjustment (GIA) predictions, but the associated horizontal deformation is not definitive, in part due to the short time periods for a significant number of the available stations. We employ an eigenpattern decomposition in order to define a unique, finite set of deformation patterns for this continuous GPS data. Similar in nature to the empirical orthogonal functions historically employed in the analysis of atmospheric and oceanographic phenomena, the method derives the eigenvalues and eigenstates from the diagonalization of the correlation matrix using a Karhunen-Loeve expansion (KLE). The KLE technique is used to identify the important modes in both time and space for the CGPS data, modes that potentially include signals such as horizontal and vertical GIA, tectonic strain, and seasonal effects. Here we filter both the vertical and horizontal velocity patterns on different spatiotemporal scales in order to study the potential geophysical sources, after the removal of correlated and random noise. The method is successful in disaggregating the linear vertical signal from more variable and less spatially correlated signals. The vertical and horizontal results are compared to the predictions of the ICE-3G GIA loading model with a number of plausible mantle viscosity profiles. The horizontal velocity analysis allows for qualitative differentiation between several potential GIA models and suggests that, with longer time series, this technique can be employed to remove correlated noise and improve estimates of crustal strain patterns and their sources.
ISSN: 0033-4553

hal-00795454
https://hal.archives-ouvertes.fr/hal-00795454
DOI : 10.1007/s00024-011-0420-1

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Detailed structural and petrological mapping in the Beni Bousera peridotite (Rif Belt, northern Morocco) shows that this orogenic peridotite massif is composed of four tectono-metamorphic domains with consistent kinematics, marked by a pervasive, shallowly dipping foliation with a NW-SE stretching lineation that progressively rotates towards a NNE-SSW orientation in the lowermost part of the massif. From top to bottom, these domains are garnet-spinel mylonites, Ariègite subfacies fine-grained porphyroclastic spinel peridotites, Ariègite-Seiland subfacies porphyroclastic, and Seiland subfacies coarse-porphyroclastic to coarse-granular spinel peridotites. Microstructures and crystal preferred orientations point to deformation dominantly by dislocation creep in all domains, but the continuous increase in average olivine grain size indicates decreasing plastic work rates from top to bottom. This evolution in deformation conditions is consistent with the change in synkinematic pressure and temperature conditions, from 900°C at 2*0 GPa in garnet-spinel mylonites to 1150°C at 1*8 GPa in the Seiland domain. A pervasive diffuse dunitic-websteritic layering subparallel to the foliation suggests deformation in the presence of melt in the Seiland domain. Gravitational instabilities owing to local melt accumulation may account for <200 m wide areas exhibiting a vertical lineation in this domain. To account for the consistent kinematics and the tectono-metamorphic evolution, which implies a temperature gradient of c. 125°C km−1 preserved across the Beni Bousera massif, we propose that the entire massif records the functioning of a low-angle shear zone, a few kilometres wide, which accommodated exhumation of the base of the lithosphere from ∼90 to ∼60 km depth. Partial melting in the Seiland domain may be explained by fast decompression of the footwall, without the need for exotic heat sources. Moreover, if the present-day orientation of the shear zone is similar to that when it was active in the mantle, the stretching lineations at high angle to the metamorphic gradient imply that shearing parallel to the trend of the belt accompanied thinning; that is, a transtensional deformation of the margin.
ISSN: 0022-3530

hal-00950056
https://hal.archives-ouvertes.fr/hal-00950056
DOI : 10.1093/petrology/egt067

É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
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 American Geophysical Union
Résumé : International audience
More than 10 years of geodetic measurements demonstrate an uplift rate of 1–3 mm/yr of the high topography region of the Western Alps. By contrast, no significant horizontal motion has been detected. Two uplift mechanisms have been proposed: (1) the isostatic response to denudation responsible for only a fraction of the observed uplift and (2) the rebound induced by the Wurmian ice cap melting which predicts a broader uplifting region than the one evidenced by geodetic observations. Using a numerical model to fit the geodetic data, we show that a crustal viscosity contrast between the foreland and the central part of the Alps, the latter being weaker with a viscosity of 1021 Pa s, is needed. The vertical rates are enhanced if the strong uppermost mantle beneath the Moho is interrupted across the Alps, therefore allowing a weak vertical rheological anomaly over the entire lithosphere.
ISSN: 0094-8276

hal-01332757
https://hal.archives-ouvertes.fr/hal-01332757
DOI : 10.1002/2016GL067821

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The island of Taiwan is located on the convergent boundary between the Philippine Sea plate and the Chinese continental margin. It offers very active mountain building and collapsing processes well illustrated by the rugged topography, rapid uplift and denudation, young tectonic landforms, active faulting and numerous earthquakes. In this paper, using simple models, we have estimated vertical movements and associated absolute gravity variations which can be expected along a profile crossing the southern part of the island and probably suffering the highest rates of rising. The two different tectonic styles proposed for the island, thin-skinned and thick-skinned, were taken into account. Horizontal and vertical movements were modeled by an elastic deformation code. Gravity variations due to these deformations are then modeled at a second step. They are dominated by plate and free-air effects, i.e. elevation of the topography, with several mu Gal yr(-1). By comparison, gravity changes generated by mass transfers are weak: maximum 0.1 mu Gal yr(-1) with the thin-skinned tectonic and 0.3 mu Gal yr(-1) with the thick-skinned tectonic. Though elastic rheology has limitations, this modeling offers interesting results on what gravity signal can be expected from the AGTO project (Absolute Gravity in the Taiwanese Orogen), which proposes to study the dynamic of these mountain ranges using absolute gravimetry (AG) and also including relative gravimetry (RG) and GPS measurements.
ISSN: 0264-3707

hal-00456035
https://hal.archives-ouvertes.fr/hal-00456035
DOI : 10.1016/j.jog.2009.09.004

Éditeur(s) : HAL CCSD Berlin: A. Asher & Co.
Résumé : International audience
In this work, we compare IASI-retrieved vertical water vapour profiles and related precipitable water over a complex region, namely the Rhine Valley area, during the pre-operational period of IASI exploitation (June?August 2007). Both IASI water vapour mixing ratio profiles and integrated water vapour content are retrieved from L1C radiances spectra using two techniques and compared with water vapour related observations acquired during the Convective and Orographically-induced Precipitation Study (COPS) field campaign that took place in this area at that time (i.e. lidars, radiosoundings and a global positioning system - GPS - station network). This work addresses the issue of IASI vertical spatial resolution and its capability to detect two-layer water vapour structures such as those observed in a mountainous area and which play an important role in convective initiation or inhibition. We found that this capability mostly relies on the type of a-priori background vector (climatology or space-time colocated ECMWF analysis), which is used within the retrieval scheme. Systematic comparison of water vapour products derived from 71 IASI spectra confirms that IASI can retrieve water vapour amounts in 2 km width layers, in the lower troposphere, with an accuracy of approximately 10%.
ISSN: 0369-1845

hal-00923710
https://hal.archives-ouvertes.fr/hal-00923710
https://hal.archives-ouvertes.fr/hal-00923710/document
https://hal.archives-ouvertes.fr/hal-00923710/file/masiello.pdf
DOI : 10.1127/0941-2948/2013/0430

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Coastal aquifers often display seawater intrusion resulting in the formation of a salty water wedge progressing inland. This study investigates the mass transfers in the mixing zone at the freshwater-seawater interface where the water is out of equilibrium with the rock-forming carbonates. Investigations were conducted in two boreholes, separated by 5 m, at the Ses Sitjoles test site (Mallorca Island, Spain) where repeated electrical conductivity logs of the formation and the saturating fluid, as well as regular pore-water sampling and permanent downhole multi-parameter monitoring of the water were performed over a period of 9 a. In the mixing zone, the significant acidification, the calcite saturation index profile and the Ca concentration profile cannot be explained by conservative mixing nor by dissolution-precipitation reactions only. Conversely, the analysis of organic C content and of the distinctly different time-resolved pH profiles measured in the two boreholes suggests the development of perennial biomass that triggers calcite dissolution. Moreover, the presence of biomass seems to be correlated with the permeability and vertical connectivity at the meter-scale. It is speculated that the mechanism could be self-activated because the microbiological activity induces calcite dissolution and tends to increase porosity and permeability that favors biomass development.
ISSN: 0883-2927

hal-00937665
https://hal.archives-ouvertes.fr/hal-00937665
DOI : 10.1016/j.apgeochem.2013.09.011

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
The reactivation of faults with near-optimal orientations is commonly considered to control the state of stress in the crust. Near the surface, where a principal stress direction is vertical, the attitude of such faults is explained by Anderson's theory. This raises the questions of how prevalent this type of faulting actually is in current seismicity, down to what depth it frequently occurs, and what range of friction angles explains it best. The Global Centroid Moment Tensor catalog is analyzed to address these questions. Dip-slip and strike-slip mechanisms are dominant, and oblique slips are relatively rare for well-constrained events with depths shallower than 30 km. Preferred Andersonian faulting is the simplest, but not the unique, explanation for this dominance. Isolating reverse, strike- slip, and normal events reveals an asymmetry in the distribution of nodal plane dips and plunges of the P, B, and T axes between reverse and normal faulting. Assuming that the most frequent events correspond to reactivation near optimal orientations yields 40-60 degrees and 0-20 degrees friction angles for reverse and normal faults, respectively. This indicates that reverse and normal faulting mechanics are not symmetrical with respect to stress configuration as predicted by Anderson's theory.
ISSN: 8755-1209

hal-00412159
https://hal.archives-ouvertes.fr/hal-00412159
DOI : 10.1029/2007RG000240