Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
From outcrops located in Provence (South-East France), we describe the distribution, the microstructures, and the petrophysical properties of deformation band networks related to both contractional and extensional tectonic events. In contraction, pervasively distributed networks of reverse-sense compactional shear bands are observed in all folded sand units of the foreland, whereas localized networks of clustered reverse-sense shear bands are only observed close to a large-scale thrust. In extensional setting, networks of clustered normal-sense shear bands are generally observed adjacent to map-scale faults (100 m–10 km scale), although some randomly distributed bands are also observed between these faults. Normal-sense cataclastic faults, i.e. zone of deformation bands containing a localized slip-surface, are also observed to be restricted to sand units, suggesting that faults initiated in the sands during extension, but not during contraction. Shear bands and faults show cataclastic microstructures with high-permeability reduction whereas compactional shear bands show crush microbreccia or protocataclastic microstructures with moderate permeability reduction. This basin-scale analysis underlines the major role of tectonic settings (thrust-fault versus normal-fault andersonian-stress regime) and the influence of inherited large-scale faults on the formation/localization of low-permeability shear bands. We also provide a geometrical analysis of the band network properties (spacing, thickness, shear/compaction ratio, degree of cataclasis, petrophysical properties) with respect to the median grain size, porosity and grain sorting of host sand. This analysis suggests that grain size, although less important than stress-state conditions and the presence of large-scale faults, has a non-negligible effect on band network geometry. No correlations are observed between the grain sorting, porosity and band network geometry.
ISSN: 0264-8172

hal-01116421
https://hal.archives-ouvertes.fr/hal-01116421
DOI : 10.1016/j.marpetgeo.2013.12.020

Éditeur(s) : HAL CCSD
Résumé : Different applications of linear elasticity in structural geology are presented in this thesis through the development of three types of numerical computer codes. The first one uses forward modeling to study displacement and perturbed stress fields around complexly faulted regions. We show that incorporating inequality constraints, such as static Coulomb friction, enables one to explain the angle of initiation of jogs in extensional relays. Adding heterogeneous material properties and optimizations, such as parallelization on multicore architectures and complexity reduction, admits more complex models. The second type deals with inverse modeling, also called parameter estimation. Linear slip inversion on faults with complex geometry, as well as paleo-stress inversion using a geomechanical approach, are developed. The last type of numerical computer code is dedicated to restoration of complexly folded and faulted structures. It is shown that this technique enables one to check balanced cross-sections, and also to retrieve fault chronology. Finally, we show that this code allows one to smooth noisy 3D interpreted faulted and folded horizons using geomechanics.
Différentes applications de l'élasticité linéaire en géologie structurale sont présentées dans cette thèse à travers le développement de trois types de codes numériques. Le premier utilise la modélisation directe pour étudier les déplacements et champs de contraintes autour de zones faillées complexes. On montre que l'ajout de contraintes inégalitaires, telles que la friction de Coulomb, permet d'expliquer l'angle d'initiation des dominos dans les relais extensifs. L'ajout de matériaux hétérogènes et d'optimisations, telles la parallélisation sur processeurs multi-coeurs ainsi que la réduction de complexité des modèles, permettent l'étude de modèles beaucoup plus complexes. Le second type de code numérique utilise la modélisation inverse, aussi appelée estimation de paramètres. L'inversion linéaire de déplacements sur les failles ainsi que la détermination de paléo-contraintes utilisant une approche géomécanique sont développées. Le dernier type de code numérique concerne la restauration de structures complexes plissées et faillées. Il est notamment montré qu'une telle méthode permet de vérifier l'équilibre de coupes géologiques, ainsi que de retrouver la chronologie des failles. Finalement, nous montrons que ce même code permet de lisser des horizons 3D faillés, plissés et bruités en utilisant la géomécanique.
https://tel.archives-ouvertes.fr/tel-00537899

tel-00537899
https://tel.archives-ouvertes.fr/tel-00537899
https://tel.archives-ouvertes.fr/tel-00537899/document
https://tel.archives-ouvertes.fr/tel-00537899/file/thesemaerten.pdf

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Properly assessing the extent and magnitude of fault ruptures associated with large earthquakes is critical for understanding fault behavior and associated hazard. Submarine faults can trigger tsunamis, whose characteristics are defined by the geometry of seafloor displacement, studied primarily through indirect observations (e.g., seismic event parameters, seismic profiles, shipboard bathymetry, coring) rather than direct ones. Using deep-sea vehicles, we identify for the first time a marker of coseismic slip on a submarine fault plane along the Roseau Fault (Lesser Antilles), and measure its vertical displacement of ∼0.9 m in situ. We also map recent fissuring and faulting of sediments on the hangingwall, along ∼3 km of rupture in close proximity to the fault's base, and document the reactivation of erosion and sedimentation within and downslope of the scarp. These deformation structures were caused by the 2004 M w 6.3 Les Saintes earthquake, which triggered a subsequent tsunami. Their characterization informs estimates of earthquake recurrence on this fault and provides new constraints on the geometry of fault rupture, which is both shorter and displays locally larger coseismic displacements than available model predictions that lack field constraints. This methodology of detailed field observations coupled with near-bottom geophysical surveying can be readily applied to numerous submarine fault systems, and should prove useful in evaluating seismic and tsunamigenic hazard in all geodynamic contexts.
ISSN: 0012-821X

insu-01557932
https://hal-insu.archives-ouvertes.fr/insu-01557932
https://hal-insu.archives-ouvertes.fr/insu-01557932/document
https://hal-insu.archives-ouvertes.fr/insu-01557932/file/EPSL-2006-escartin.pdf
DOI : 10.1016/j.epsl.2016.06.024

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

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

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Wallace (1951) and Bott (1959) were the first to introduce the idea that the slip on each fault surface has the same direction and sense as the maximum shear stress resolved on that surface. However, this simplified hypothesis is questionable since fault mechanical interactions may induce slip reorientations. Earlier numerical geomechanical models confirmed that the slickenlines (slip vectors) are not necessarily parallel to the maximum resolved shear stress but are consistent with local stress perturbations. This leads us to ask as to what extent the Wallace and Bott simplifications are reliable as a basis hypothesis for stress inversion from fault slip data. Here, a geomechanical multi-parametric study using a 3D boundary element method, covering (i) fault geometries such as intersected faults or corrugated fault surfaces, (ii) the full range of Andersonian state of stress, (iii) fault friction, (iv) fault fluid pressure, (v) half space effect and (vi), rock properties, is performed in order to understand the effect of each parameter on the misfit angle between geomechanical slip vectors and the resolved shear stresses. It is shown that significant misfit angles can be found under specific configurations invalidating the Wallace and Bott assumptions, even though fault friction tends to minimize the misfit. We therefore conclude that in such cases, stress inversions based on fault slip data should be interpreted with care.
ISSN: 0040-1951

hal-01218418
https://hal.archives-ouvertes.fr/hal-01218418
DOI : 10.1016/j.tecto.2015.06.030

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The Parentis and Arzacq-Mauleon basins located in front of the V-shaped oceanic propagator in the Bay of Biscay present evidence for extreme crustal thinning. In this paper we investigate the 3D structure of these rift basins, based on field observations and the interpretation of seismic data. We compare these results with those obtained from two different and independent inversion methods: first a 3D gravity inversion and second the standard Euler deconvolution.;For the Mauleon Basin our results show that the positive gravimetric anomaly identified above its southern margin is the consequence of two shallower high density bodies that are separated by the Pamplona fault and a deeper high density body. The high density bodies can be explained by the presence of mid-crustal and mantle rocks that were exhumed or uplifted at shallower depth during Early Cretaceous rifting before they were reworked and integrated to the Pyrenean chain during compression phase. Also, during this reactivation phase, some slices of the exhumed mid-crustal and mantle rocks were sheared off and were integrated in the present-day thrust belt in the Mauleon basin. For the Parentis Basin we can demonstrate, based on seismic data and gravimetric inversion methods, a decrease in extension from west to east, which is compatible with the V-shape geometry of the overall basin. Along strike, a change in the fault geometry from downward concave top-basement detachment faults to upward concave high-angle faults can be observed eastwards, i.e. towards the termination of the basin. A key structure, controlling the evolution of the Parentis Basin, is the east-west trending Ibis fault. We interpret this fault to have initially formed as a strike slip fault before it was reactivated during later crustal thinning. At present, it forms the limit between an upper plate sag basin to the north and a lower plate sag basin, floored at least locally by a top-basement detachment faults to the south. The strong asymmetry of the basin is supported by the shape of the basin and the results of standard Euler deconvolution.
ISSN: 0040-1951

hal-00523444
https://hal.archives-ouvertes.fr/hal-00523444
DOI : 10.1016/j.tecto.2010.04.023

Éditeur(s) : HAL CCSD
Résumé : In order to identify and quantify the Rif (Morocco) present-day deformation, especially those related to active faults (location, geometry, kinematic), we perform a morphotectonic study. The studied area and faults are located in the Northern part (conjugate faults of Al-Hoceima Bay: Trougout, Boujibar and Rouadi), the Eastern part (Nekor fault) and the Southern part of the chain (Rif frontal thrust between Fes and Meknes cities).Morphotectonic tools (DEMs at different scales, satellite and aerial pictures), field studies and Quaternary dating methods (U/Th, OSL, 14C and in-situ cosmogenic isotopes 10Be, 3He, 36Cl) allow us to characterize transtensional deformation in the Northern Rif. Horizontal and vertical slip rates along the N-S-oriented Trougout major fault are ~1.6 mm/yr and ~2 mm/yr respectively. In the Southern part, the Rif frontal thrust shows a shortening rate of ~1.8 mm/yr into Fes and Meknes area. In the Eastern part, the deformation is focused on the Nekor left-lateral strike-slip major fault with ~1.5 mm/yr. Incision features and raised marine terraces involve an uplift of the chain about ~0.2 mm/yr. The overall deformation in the Rif seems to be continuous since the Messinian-Pliocene time.These new morphotectonic constraints are consistent with the GPS measurements showing southwestward overall motion of most of the Rif belt with respect to stable Africa, associated to clockwise rotation. Finally, we propose a geodynamic model in which the Rif, the Alboran sea and the Betics belong to a single microplate bounded by the strike-slip faults of the Betics on the North, the Rif frontal thrust on the South, the Nekor and TASZ (Trans Alboran Shear Zone) strike-slip faults on the East and North-East, and the western Rif folded structures on the West.
L'objectif de cette étude est d'identifier et de quantifier les déformations actuelles dans la chaine du Rif (Maroc), notamment celles associées aux failles actives (localisation, géométrie, cinématique). Pour ce faire, nous mettons en œuvre pour la première fois dans cette région, une approche morphotectonique. Les sites et failles sélectionnés et étudiés se situent dans la partie Nord (les failles conjuguées de la baie d'Al-Hoceima: Trougout, Boujibar et Rouadi), l'Est (la faille de Nékor) et le Sud (front chevauchant sud rifain entre Fes et Meknes).A partir des outils de la morphotectonique (MNTs à différentes résolutions, images satellitaires, photographies aériennes), d'analyses de terrain et grâce aux méthodes de datation du Quaternaire (U/Th, OSL, 14C et isotopes cosmogéniques produits in-situ 10Be, 3He, 36Cl), cette étude met en évidence et quantifie une déformation transtensive dans la partie Nord de la chaine. Les vitesses de glissement horizontales et verticales le long de la faille principale d'orientation méridienne de Trougout sont de ~1.6 mm/an et ~2 mm/an, respectivement. Dans la partie Sud, nous montrons que le front rifain entre Fez et Meknes est actif avec une vitesse de raccourcissement de ~1.8 mm/an. La déformation dans la région Est se localise principalement sur la faille du Nékor avec une vitesse décrochante senestre de ~1.5 mm/an. Par ailleurs, la présence de nombreuses figures d'incision et de surfaces marines soulevées au niveau du littoral Nord témoignent d'une surrection généralisée de la chaine de l'ordre de ~0.2 mm/an. L'ensemble de cette déformation semble être constante depuis la période du Messino-Pliocène. Ces résultats sont cohérents avec les données GPS qui montrent un déplacement de ~5mm/an vers le SW de l'ensemble du Rif associé à un mouvement rotatif horaire, dans le référentiel Afrique stable. D'un point de vue géodynamique, nous proposons que le Rif, la mer d'Alboran et les Bétiques appartiennent à une même microplaque limitée au Nord par les failles décrochantes des Bétiques, au Sud par le chevauchement rifain, à l'Est et au Nord-Est par les failles décrochantes de Nékor et de la TASZ (Trans Alboran Shear Zone) et à l'Ouest par les structures plissées bordières du Rif.
https://tel.archives-ouvertes.fr/tel-01170466

NNT : 2014MON20105
tel-01170466
https://tel.archives-ouvertes.fr/tel-01170466
https://tel.archives-ouvertes.fr/tel-01170466v2/document
https://tel.archives-ouvertes.fr/tel-01170466/file/41875_POUJOL_2014_archivage_cor.pdf

Éditeur(s) : HAL CCSD Elsevier Masson
Résumé : Normal faulting mechanisms observed in the northern foothills of the Central–Western Pyrenees are remarkable, since one expects thrust faults at a convergent plate boundary. To understand the mechanisms involved, we used numerical modeling and investigated the impact of the following processes: gravitational potential energy associated with topography and dense crustal blocks; isostatic compensation in response to denudation and/or sedimentation. To decipher the effects of each process, we designed three model geometries and added a pre-existing weak fault where most of the seismicity occurs. To evaluate our model results, we derived the fault slip rate from the focal mechanisms in the region where we have the fault in our model. We found a slip rate of ∼15 m/Ma, which is in agreement with our modelling results. We conclude that flexural rebound induced by surface processes is able to explain the seismicity in Central–Western Pyrenees.
ISSN: 1631-0713

hal-01346028
https://hal.archives-ouvertes.fr/hal-01346028
https://hal.archives-ouvertes.fr/hal-01346028/document
https://hal.archives-ouvertes.fr/hal-01346028/file/Genti.pdf
DOI : 10.1016/j.crte.2015.08.004

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The Cycladic Blueschist Unit (CBU) represents the northern passive margin of the Adria continental block and ophiolites that are the remnants of the Pindos Ocean, which were affected by high pressure-low temperature metamorphism in the blueschist and eclogite facies during the Eocene. Prior to the Lutetian, the ophiolitic and margin units were thrust towards the SW inside a NE-dipping subduction zone. Two subsequent exhumation stages are characterized by top-to-the-NE senses of shear: i) from mantle depths to lower crustal levels, before 37 Ma, which was accommodated by an extensional inversion of the Vardar suture zone, and ii) from deep crustal levels to the near-surface, between 30 and 20 Ma, which was accommodated by the North Cycladic Detachment. It is during this late stage of ductile exhumation that the central Cyclades Core Complex was exhumed. Segmentation of the CBU by normal faults started at 12 Ma, as recorded by the low temperature thermochronometers, and was coeval with block rotation and folding with NS-trending axes, which were predominantly controlled by the Myrthes-Ikaria strike-slip fault in the centre of the Cyclades. A map restoration of the CBU geometry, prior to segmentation, has been carried out using available structural and paleomagnetic data. The restoration shows that, following the CBU exhumation from mantle depths to crustal levels, a single core complex was exhumed along a single NE dipping detachment.
ISSN: 0040-1951

insu-00672320
https://hal-insu.archives-ouvertes.fr/insu-00672320
DOI : 10.1016/j.tecto.2011.12.025

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

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

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We review the evolution of concepts on and methods of estimating the state of stress from fault movements. Theories of failure in isotropic materials suggested a simple geometrical construction of optimal principal stress directions from a fault plane and its associated slip. These optimal directions align shear stress and slip directions and maximize the difference between shear stress and frictional resistance on the fault plane. Optimal stress directions for calcite twinning are obtained by a similar construction, with the difference that they maximize shear stress. Force representation of seismic sources independently introduced pressure, P, and tension, T, axes at positions that also maximize shear stress on both nodal planes. Frictional slip theory and the constraint that slip and shear stress directions be parallel allowed to address reactivation of pre-existing faults. This suggested that stress could also be inverted from reactivated fault and slip data or earthquake focal mechanisms. Early methods relied on geometrical constructions as a substitute for calculations, whereas later methods relied on software as these calculations became tractable with the help of computers. Similar methods were developed for the inversion of stress from crystal twin gliding with non-optimal geometry, with a different criterion that relies on a threshold of the component of shear stress along the gliding line. Even though these methods seek a common stress tensor compatible with fault and slip data, their main use is to separate polyphase data into homogeneous subsets and help deciphering complex tectonic histories. Fault and slip data can also be analyzed to constrain the strain rather than the stress tensor. In most cases this involves a summation and yields an average strain for the considered rock volume. Stress inversion thus appears better suited for differentiating heterogeneous data whereas strain analysis appears better suited for homogenizing them.
ISSN: 0040-1951

hal-00759043
https://hal.archives-ouvertes.fr/hal-00759043
DOI : 10.1016/j.tecto.2012.02.009

Éditeur(s) : HAL CCSD
Résumé : The present study addresses the question of increasing need in water supply. One line of research consists in investigating karst aquifers that represent significant water supplies for the areas around the Mediterranean. This work aims at developing a methodology of predictive evaluation of karst distribution and karst thickness within the carbonate massifs. It is based on the study of the Tethysian carbonate platform from the area of Montpellier, South of France. This area underwent a polyphased, complex, tectonic evolution, involving successive karst forming events. The approach is based on the relationships between eustacy and geodynamics, on one hand, and base-level variation recorded within palaeokarsts, on the other hand. Sedimentary fill of palaeokarsts allows to identify distinct karstification phases and sometimes to date them. It is a marker of base-level variations, which we attempt to quantify, that is responsible for karst formation and karst filling and sealing. For three palaeokarsts of regional extent, we analyze : karst sedimentary-fill, the distribution along upstream-downstream profile, the vertical extension of palaeokarsts, indicative of the base-level lowering. - The " Bauxite palaeokarst " results from the two-stage evolution of the " Isthme Durancien " : i) Bauxitic weathering of the Neocomian cover over the southern Massif Central , related to an inversion during Barremian ; ii) ≥1600 m uplift and exhumation during Aptian-Albian, inducing erosion of alterites, upstream, and deposition of detrital bauxites, downstream. - The " Early Paleocene palaeokarsts with marine sedimentary-fill " indicate base level variations of more than 350m in the upstream area. It is suggested here that such base-level variations, exceeding eustacy and recorded in tabular undeformed series, result from desiccation-flooding events in an endoreic silled basin. - The " Miocene palaeokarst with Cevennes-sourced sediment-fill " records a Serravalian-Tortonian base-level drop, depicted by canyon incision in the hinterland of the Gulf of Lion Margin. Uplift values of ≥ 400m upstream, about 250m in an intermediate area and 10's of m downstream, are measured. The present 3D geometry and distribution of carbonate massifs (exposed and covered) are modeled in the study area, from the Cevennes to the Mediterranean. The 3D geological model integrates seismic reflection, borehole and gravimetric data. Tectonic restoration of a N-S section, extracted from the 3D model, allows to reconstruct paleo-base-levels with respect with the carbonate massifs, at different stages of deformation and successive phases of karstification. The 3D paleo-flows through carbonates across the study area, are thus constrained, especially during the Messinian desiccation stage. Application of the results to hydrogeology leads to the following : - Uplift of upstream hinterland (Cretaceous, Miocene) produces large amplitude karsts, but they do not make significant reservoirs. - Karstification responsible for the development of important and deep karst reservoirs is related to desiccation of endoreic basins (early Paleocene and Messinian). - Theses two events have produced two main deep karst reservoirs, with a gravimetric signature : (1) Upstream of Montpellier Thrust, between Pic St Loup, Les Matelles Fault and Hérault Basin, a deep karsts developed over 900m thickness in the Malm limestone, during Paleocene. (2) Downstream (south) of Montpellier Thrust, the coastal karsts extend offshore over more than 10km, they are sealed by a thick Miocene to Quaternary cover, and they result from the Messinian event.
Cette étude s'inscrit dans un contexte de demande croissante en ressources en eau, conduisant à mieux contraindre et estimer la distribution des aquifères karstifiés qui représentent une ressource importante, notamment sur le pourtour méditerranéen. Le but de ce travail est d'élaborer une méthode prédictive de la distribution de la karstification et de l'épaisseur de massifs carbonatés karstifiés, notamment pour les karsts profonds. La méthodologie est développée à partir de l'étude des karsts de la plateforme carbonatée téthysienne du Languedoc montpelliérain; région qui a subi une histoire tectonique complexe et polyphasée, entrainant des phases de karstification successives. Elle est basée sur la compréhension des relations entre eustatisme, géodynamique et niveau de base, enregistrées par les paléokarsts de la région. Le remplissage sédimentaire des paléokarsts permet d'identifier les épisodes de karstification et de les dater ; c'est un marqueur des variations -que l'on cherche à quantifier- du niveau de base local responsable de la formation et du comblement des karsts. On analyse 3 paléokarsts d'extension régionale dont on peut étudier : le remplissage karstique, la distribution le long de profils amont-aval et l'extension verticale indicatrice de l'amplitude de l'abaissement du niveau de base dont il est issu. On en déduit les causes géodynamiques à l'origine de leur formation: - Le " paléokarst des bauxites " à remplissage caractéristique, résulte de l'évolution en 2 phases de " l'Isthme Durancien " : altération bauxitique de la couverture néocomienne sur le sud du Massif Central, liée à une inversion au Barrémien, suivi d'une surrection ≥ 1600m à l'Aptien-Albien, causant l'érosion des altérites en amont et le dépôt des bauxites détritiques en aval. - Le " paléokarst à remplissage marin paléocène inférieur " indique des variations du niveau de base local ≥350m dans la zone amont sud-cévenole. On propose que ces variations de niveau de base d'amplitude très supérieures aux variations eustatiques et dans une zone faiblement déformée, résultent d'assèchements et d'ennoiements d'un bassin endoréique à seuil. - Le " paléokarst à remplissage détritique d'origine cévenole ", enregistre la chute du niveau de base au Serravalien-Tortonien, qui se traduit par l'incision de canyons dans l'arrière-pays de la marge du Golfe du Lion. On détermine une surrection de plus de 400 mètres à l'amont, près de 250m dans les blocs faillés intermédiaires et quelques dizaines de mètres dans la zone avale. En parallèle, on établit la géométrie et la distribution actuelles des massifs carbonatés dans la zone d'étude (y compris sous couverture), en construisant un modèle géologique 3D, des Cévennes à la Méditerranée. On intègre les données de sismique réflexion, de forages et les données gravimétriques. La restauration tectonique séquentielle d'une coupe Nord-Sud extraite du modèle géologique permet d'appréhender la position des paléo-niveaux de base par rapport aux massifs carbonatés, pour chaque époque de karstification. On déduit l'organisation 3D des paléo-écoulements à travers la géométrie complexe des massifs carbonatés languedociens, notamment pendant l'évènement messinien. En terme d'évaluation de la ressource en eau, on déduit que: - Les surrections de l'arrière-pays (Crétacé " moyen ", Miocène) produisent des karsts de grande amplitude, mais ne constituent pas de réservoirs significatifs. - La karstification à l'origine de réservoirs exploitables importants et profonds est créée par l'assèchement de bassins endoréiques, qui ont affecté la région au Paléocène inférieur et au Messinien . - Ces deux évènements sont à l'origine de deux principaux réservoirs karstiques profonds qui ont une signature gravimétrique : (1) En amont du chevauchement de Montpellier, entre le Pic Saint Loup, la faille des Matelles et le bassin de l'Hérault, le karst profond se développe sur une épaisseur atteignant 900m, et est attribuée à l'événement du Paléocène. (2) Au Sud du chevauchement de Montpellier les karsts côtiers s'étendent jusqu'à plus de 10 km au large, scellés par une épaisse couverture Miocène et Plio-Quaternaire, et résultent de l'événement Messinien.
https://tel.archives-ouvertes.fr/tel-00936257

tel-00936257
https://tel.archives-ouvertes.fr/tel-00936257
https://tel.archives-ouvertes.fr/tel-00936257/document
https://tel.archives-ouvertes.fr/tel-00936257/file/thA_se_husson_couleurs.pdf

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
New magnetotelluric and receiver transfer function studies provide insights from the upper to the lower crust of the eastern Betic Cordillera, which is deformed by large folds, normal faults, and a major transcurrent left-lateral fault, the Carboneras Fault Zone (CFZ). Receiver function analysis determines a NNW dipping Moho reaching 20 degrees that increases in depth, from 20 km south of the CFZ up to 34 km in the Sierra de Los Filabres. In addition, seismic discontinuities determined in the upper crust are interpreted as major contacts between metamorphic complexes that are detached and folded. The MT inversion model reveals a conductive zone, also representing a crustal seismic discontinuity, associated with the Alpujarride/Nevado-Filabride contact and fitting the N vergent geometry of the Sierra Alhamilla antiform. A small flexure at Moho coincides with the CFZ, as revealed by the Bouguer anomaly trend, in agreement with the receiver function results. Moreover, the Bahr strike and tipper angle at the stations placed closest to the CFZ clearly reveal the continuity of the CFZ at least down to approximately 15 km in depth, crossing all the detected crustal discontinuities up to the Moho. The lack of a clear Moho offset associated with the Carboneras Fault supports the idea that some large strike-slip faults tend to accommodate the deformation by a broadening fault zone at lower crustal levels. Its nucleation could occur at the base of a thin crust, where melting processes critically reduced the lithospheric strength during the late Miocene, to then propagate upward, reaching the topographic surface. Northward, the lithosphere comprised moderately larger strength, and the crustal discontinuities favored the development of larger folds with kilometric amplitude instead of strike-slip faults since the late Miocene.
ISSN: 1525-2027

hal-00671024
https://hal.archives-ouvertes.fr/hal-00671024
DOI : 10.1029/2010GC003274

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Fault friction is a parameter that is difficult to assess along fault zones since its determination depends on the knowledge of any factor controlling the state of stress around faults. In brittle homogeneous rocks, a limited number of these factors, such as the shape of the fault surface, the vicinity of fault tips or the remote stress ratio, are crucial to constrain for this determination. In this paper, we propose to analyse a field example in which all these properties are met and where the nature of the slipped structure suggest differences in static friction. We compare the orientations of branching fractures at strike-slip relay zones between en echelon stylolites and en echelon joints both reactivated in shear. The field data are compared with both photoelastic and 3-D numerical models that consider the remote stress conditions and the role of the geometry of the strike-slip segments. Based on field observations, these analyses quantitatively demonstrate the significant role of fault friction on the local stress field orientation and subsequent fracture formation. This work points out that estimations of fault friction based on analyses of fracture patterns or in situ stresses must be accompanied with a thorough investigation of the 3-D fault shape, its segmentation and the remote stress state.
ISSN: 0191-8141

hal-00560999
https://hal.archives-ouvertes.fr/hal-00560999
DOI : 10.1016/j.jsg.2010.01.008

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
The outstanding exposure of deformation bands in a multilayer sequence of the Orange quarry, Provence, France, motivates a study investigating the development and controls of deformation band geometries by field mapping combined with finite element modeling. Field mapping yields new insights into angular relations between, and evolution of, the deformation bands as well as provides input parameters for the setup and boundary conditions of the numerical simulations using ADELI 2D. In particular, reverse-sense deformation bands are found to have developed within three of the four exposed layers of the multilayer stack with similar orientations to bedding, indicating layer-parallel contraction. Numerical simulations, carried out as a parametric study to investigate the influence on deformation band development and geometry, successfully reproduce the sequential development of deformation bands in the three layers as seen in the field as well as recreate the observed angular geometries of similar to 39 degrees to the maximum compression. This angular relationship is put into context with observations from other field studies and theoretical values to explore controls of deformation band orientations. Furthermore, band localization is impeded in the model by increased layer cohesive strength, demonstrating the importance of host-rock property variations on where localization of bands within the multilayer sequence. Results from both field investigation and numerical simulations are consistent with regional tectonics, where the strata underwent gentle folding due to thrust faulting at depth, and bring further insight into the development and geometry of deformation bands in both extensional and contractional tectonic regimes.
ISSN: 0148-0227

hal-00639452
https://hal.archives-ouvertes.fr/hal-00639452
DOI : 10.1029/2011JB008365

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Tectonically controlled landforms develop morphologic features that provide useful markers to investigate crustal deformation and relief growth dynamics. In this paper, we present results of morphotectonic experiments obtained with an innovative approach combining tectonic and surface processes (erosion, transport, and sedimentation), coupled with accurate model monitoring techniques. This approach allows for a qualitative and quantitative analysis of landscape evolution in response to active deformation in the three end-member geological settings: compression, extension, and strike–slip.Experimental results outline first that experimental morphologies evolve significantly at a short time scale. Numerous morphologic markers form continuously, but their lifetime is generally short because erosion and sedimentation processes tend to destroy or bury them. For the compressional setting, the formation of terraces above an active thrust appears mainly controlled by narrowing and incision of the main channel through the uplifting hanging-wall and by avulsion of deposits on fan-like bodies. Terrace formation is irregular even under steady tectonic rates and erosional conditions. Terrace deformation analysis allows retrieving the growth history of the structure and the fault slip rate evolution. For the extensional setting, the dynamics of hanging-wall sedimentary filling appears to control the position of the base level, which in turn controls footwall erosion. Two phases of relief evolution can be evidenced: the first is a phase of relief growth, and the second is a phase of upstream propagation of topographic equilibrium that is reached first in the sedimentary basin. During the phase of relief growth, the formation of triangular facets occurs by degradation of the fault scarp, and their geometry (height) becomes stationary during the phase of upstream propagation of the topographic equilibrium. For the strike–slip setting, the complex morphology of the wrench zone, composed of several interacting fault segments, enhances the interactions with the drainage network. Because of the widening of the main fault zone toward the surface, a significant amount of distributed deformation is observed along the wrench zone. Locally, where two terminations of fault segments interact, less than a quarter of the far field displacement can remain measurable using fault offsets, leading to a systematic underestimation of the real fault slip rate.These different experimental examples illustrate the great potential of the approach coupling deformation mechanisms and erosion–transport–sedimentation processes to investigate qualitatively and quantitatively the morphotectonic evolution of tectonically controlled landscapes.
ISSN: 0169-555X

hal-01243587
https://hal.archives-ouvertes.fr/hal-01243587
DOI : 10.1016/j.geomorph.2015.02.011

Éditeur(s) : HAL CCSD
Résumé : Today, one of the main challenges in the oil industry is the exploitation of new resources in naturally fractured reservoirs often located in structurally complex areas such as plate boundaries, mountain ranges or near salt structures. While numerical models of rock deformation based on continuum mechanics are becoming industry standard in providing efficient means for modelling natural fractures in reservoirs, the knowledge of past perturbed stress fields, which is referred as the type (normal, wrench or reverse), orientation and magnitude of the regional or local tectonic stresses that needs to be applied as boundary condition in numerical simulations, is often unknown or approximated.The essence of this thesis is to improve the recovery of the paleostress field needed to constrain numerical simulations. I first described and discuss through a comprehensive mechanically-based multi-parametric study how far we can rely on the commonly used paleostress inversion methods based on fault slip data and the Wallace and Bott assumptions. I then compared this method with a new generation technique based on geomechanics using field observations. Since slip markers on faults are hardly observed in core or image log, I used observed natural fracture data as main drivers for the inversion of the paleostress using geomechanically-based method. I demonstrated, through various outcrop and subsurface examples, how this can efficiently be done. Finally, I aimed at addressing the problem of polyphase fracture data sets with unknown mechanical type (joints, faults, stylolites ...) and expanded the mechanical stress inversion to the separation of tectonic phases.
Aujourd'hui, l'un des principaux défis dans l'industrie pétrolière est l'exploitation de nouvelles ressources dans les réservoirs naturellement fracturés souvent situés dans des zones structurellement complexes tels que les frontières de plaques, les chaînes de montagnes ou les structures de proches des dômes de sel. Alors que les modèles numériques de déformation des roches basée sur la mécanique des milieux continus fournissant des moyens efficaces pour la modélisation de fractures naturelles dans les réservoirs, sont devenus la norme dans l'industrie, les paléo-contraintes régionales ou locales appliquées comme condition aux limites dans les simulations numériques, définies par le régime (normal, décrochant ou inverse), l'orientation, et la magnitude des contraintes tectoniques, sont souvent inconnues ou approximées.L'essence de cette thèse est d'améliorer la prédiction du champ de paléo-contrainte nécessaire pour contraindre les simulations numériques. J’ai décrit dans un premier temps à travers une étude comparative multi-paramétrique complète basée sur la mécanique, dans quelle mesure nous pouvions compter sur les méthodes d'inversion des paléo-contraintes basées sur les hypothèses de Wallace et Bott et les données de glissement sur les failles. J’ai ensuite comparé ces méthodes avec une technique d’inversion de nouvelle génération basée sur la géomécanique en utilisant les observations de terrain. Puisque les marqueurs de glissement sur les failles sont rarement observés le long du forage, j’ai utilisé l’inversion géomécanique des contraintes pour les données de fractures naturelles observées. J’ai alors démontré, l’efficacité de cette technique à travers divers exemples d'affleurements et de subsurface. Enfin, j’ai proposé une solution afin de résoudre le problème des données de fractures polyphasées à type cinématique inconnu (joints, failles, stylolites ...) et étendu l’inversion de contrainte basée sur la géomécanique à la séparation des phases tectoniques.
https://hal.archives-ouvertes.fr/tel-01324110

tel-01324110
https://hal.archives-ouvertes.fr/tel-01324110
https://hal.archives-ouvertes.fr/tel-01324110/document
https://hal.archives-ouvertes.fr/tel-01324110/file/theselejri.pdf

Éditeur(s) : HAL CCSD
Résumé : The territory of Armenia was located in the collision zone between Arabia and Eurasia, and is the seat of active intercontinental deformations, which was attested by the many historical earthquakes that have occurred in this region. All these destructive earthquakes show a strong regional seismic activity, and therefore it is important to evaluate the seismic hazard associated with active structures that generates this seismicity.This study presents the results of the analysis of the active tectonics within the northwestern and southeastern extensions of the Pambak-Sevan-Syunik fault (PSSF), a major right-lateral strike-slip fault cutting through Armenia (NW - SE). Quantifying the deformations in terms of geometry, kinematics, slip rates and earthquake activity, using cosmogenic 3He, OSL/IRSL and radiocarbon dating techniques, reveal different behaviors between the two regions.Within the northwestern extension, in the region of Amassia, the PSSF bends to the west and splits into two main WNW-ESE trending reverse faults defining a compressional pop-up structure. We estimate an uplift rate and a shortening rate of 0.5 ± 0.1 mm/y and 1.4 ± 0.6 mm/y, respectively. This suggests that most of the ~2 mm/y right lateral movement of the PSSF seems to be absorbed within the Amasia pop-structure.Within the southeastern extension, in the region of Tsghuk volcano, the PSSF shows signs of dying out at the southernmost tip of the Syunik graben. A very slow NS trending normal faulting associated with a slight right-lateral movement characterizes the tectonic activity in the region of Tsghuk volcano. We estimate vertical slip rates, EW stretching rate, and right-lateral slip rate of ~0.2 mm/y, ~0.1 mm/y and ~0.05 mm/y, respectively. These results lead to the conclusion that the right lateral movement observed further north along the PSSF is mainly transferred within other active faults further west within the Karabagh (Hagari fault or other structures further northwestwards).
L’Arménie se situe dans la zone de collision entre Arabie et Eurasie et est le siège d’une déformation intracontinental active comme l’atteste les nombreux tremblements de terre historiques qui s’y sont produits. Tous ces séismes, destructeurs, attestent d’une forte activité sismique régionale, et il est par conséquent important d'évaluer l’aléa sismique associé aux structures actives qui génère cette sismicité.Le travail présenté ici, expose les résultats de l’analyse de la tectonique active au niveau des terminaisons nord-ouest et sud-est de la faille de Pambak-Sevan-Syunik (PSSF), une des failles décrochant majeure qui traverse l'Arménie du NW au SE. La quantification des déformations en termes de géométrie, cinématique, vitesse de glissement et paléosismicité, en utilisant les méthodes de datation 3He cosmogénique, OSL/IRSL et radiocarbone, révèlent des comportements différents entre les deux régions.Au niveau de la terminaison nord-ouest, dans la région d’Amassia, la faille PSSF s’incurve vers l'ouest et se subdivise en deux branches de direction WNW-ESE, et de cinématique inverse, définissant une structure en pop-up. Nous estimons une vitesse de soulèvement de 0.5 ± 0.1 mm/an et une vitesse de raccourcissement NNE-SSW de 1.4 ± 0.6 mm/an. Ces résultats suggère que l’essentiel des ~2 mm/an de mouvement dextre estimés le long de la faille de PSSF sont absorbés au niveau du pop-up d’Amassia.Au niveau de la terminaison sud-est, dans la région du volcan Tsghuk, la faille PSSF semble disparaitre. Le peu d’activité tectonique est caractérisée par des failles normales sub-méridiennes associées à une légère composante décrochant dextre. Nous déterminons une vitesse de glissement vertical de ~0,2 mm/an, une vitesse d’extension ~EW de ~0,1 mm/an associée à une composante dextre de ~0,05 mm/an. Ces résultats suggèrent que le mouvement dextre observé le long de la faille de PSSF plus au nord, a été transféré sur d'autres failles plus à l'ouest dans le Karabakh (faille Hagari ou autres structures situés encore plus au NW).
https://tel.archives-ouvertes.fr/tel-01621440

NNT : 2016MONTT152
tel-01621440
https://tel.archives-ouvertes.fr/tel-01621440
https://tel.archives-ouvertes.fr/tel-01621440v2/document
https://tel.archives-ouvertes.fr/tel-01621440/file/2016_MKRTCHYAN_archivage.pdf

Éditeur(s) : HAL CCSD
Résumé : In Taiwan , the subduction of the Chinese continental margin under the Luzon arc results in the progressive growth of an active orogenic wedge. It is one of the best places to study the complex relationships that occurs between the tectonic processes controlling deformation (plate rheology and kinematics, metamorphism and magmatism) and surface processes (erosion and sedimentation). In order to study the kinematics of the present day convergence of Philippine Sea plate to the Asian continent around Taiwan, we perform 2D and 3D sandbox modeling and we analyzed the models with PIV (Particle Image Velocimetry) software. In 2D experiments, a model Coulomb wedge is submitted to erosion and sedimentation under topography and flux steady state conditions. The 3D experiments are mainly used to study the escape tectonics occurring in the middle west and southwest Taiwan. A wedge shaped indenter and a basement relief were used to simulate the mechanical plate boundary between Philippine Sea plate and Eurasia plate and the basement topography to the west of Taiwan respectively. We examine the exhumation patterns, the mode of fault propagation and displacement patterns by strain partitioning of vertical vs. horizontal displacement and convergent parallel vs. orogen parallel components. The preliminary conclusions are: 1. The interactions between tectonics, erosion and sedimentation control geometry and kinematics of the doubly vergent transpressional orogen. 2. Exhumation processes are controlled by the balance between erosion, frontal accretion and underplating in different structural domains of the wedge. 3. Retrowedge erosion favorites backfolding and backthrusting, induces flower structures and fan shaped cleavage distribution. 4. Sedimentation promotes underplating and out-of sequences thrusting. 5. Continental subduction or subduction of topographic highs will enhance underthrusting/underplating, thick-skin tectonics and inversion of sedimentary basin.
Bolletino di Geofisica teorica ed applicata

hal-00411133
https://hal.archives-ouvertes.fr/hal-00411133

É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