Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Superposed to ductile syn-metamorphic deformations, post-foliation deformations affect metamorphic units during their exhumation. Understanding the role of such deformations in the structuration of metamorphic units is key for understanding the tectonic evolution of convergence zones. We characterize post-foliations deformations using 3D modelling which is a first-order tool to describe complex geological structures, but a challenging task where based only on surface data. We propose a modelling procedure that combines fast draft models (interpolation of orientation data), with more complex ones where the structural context is better understood (implicit modelling), allowing us to build a 3D geometrical model of Syros Island blueschists (Cyclades), based on field data. With our approach, the 3D model is able to capture the complex present-day geometry of the island, mainly controlled by the su-perposition of three types of post-metamorphic deformations affecting the original metamorphic pile: i) a top-to-South ramp-flat extensional system that dominates the overall island structure, ii) large-scale folding of the metamorphic units associated with ramp-flat extensional system, and iii) steeply-dipping normal faults trending dominantly NNW-SSE and EW. The 3D surfaces produced by this method match outcrop data, are geologically consistent, and provide reasonable estimates of geological structures in poorly constrained areas.
ISSN: 0191-8141

insu-01188400
https://hal-insu.archives-ouvertes.fr/insu-01188400
DOI : 10.1016/j.jsg.2015.07.002

É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
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
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
Résumé : International audience
The Corinth Rift is the most seismically active zone in Europe. Since the year 2000, the Corinth Rift Laboratory (CRL, http://crlab.eu) consisting in a multidisciplinary natural observatory, aims at understanding the mechanics of faulting and earthquake nucleation in the Rift. Recent studies have improved our view about fault geometry and mechanics within CRL but there is still a critical need for a better knowledge of the structure at depth.In this project, we aim to analyze the complete seismological database (13 years of recordings) of CRL by using recently developed methodologies of structural imaging in order to determine, with high resolution, the local 3D structure and the earthquake locations. We perform an iterative joint determination of 3D velocity model and earthquake coordinates. P and S velocity models are estimated through a first arrival traveltime tomography by using the adjoint state method proposed by Taillandier et al. (2009). It consists in the minimization of the cost function between observed and theoretical arrival times,which is achieved by estimating the gradient of the cost function. Earthquake locations are also estimated with a linearized inversion approach by computing the gradient of the misfit function. Iterations are repeated until the cost function no longer decreases. We present preliminary results consisting in: (1) the adjustment of a 1D velocity model that is used as initial model of the 3D tomography and (2) a first attempt of the joint determination of 3D velocity model and earthquake locations.We also perform a preliminary receiver function analysis of teleseismic data recorded by the broadband stations of the CRL network. The multiscale approach (Gesret et. al., 2010) analysis of teleseismic P to S converted waves allows to analyze the signal toward the higher frequencies and thus to characterize the structure at depth with a high resolution. Our multiscale approach should give us reliable depth of interfaces beneath seismometers in the CRL area. In this first attempt, we adjust the 1D velocity model that produces a synthetic RF as similar as possible to the stack of observed RF for a subset of data. The identified interfaces could then be compared with heterogeneities imaged by the tomography.
AGU Fall Meeting
San Fransico, United States

hal-01417184
https://hal-mines-paristech.archives-ouvertes.fr/hal-01417184

Éditeur(s) : Université Nice Sophia Antipolis
Résumé : The Lesser Antilles is a case study of a very slow subduction (~2 cm/yr) of an old oceanic lithosphere (~84-100 Ma). The region presents a relatively low seismic activity, especially along the interplate contact, and the seismic hazard associated with a possible mega-thrust earthquake is still poorly known. This PhD thesis is a first step toward assessing the ability of the Lesser Antilles subduction zone to produce such a large subduction event. To do so, it aims at constraining the downdip width of the interplate's seismogenic zone. The lack of coverage of permanent seismological stations is a major limitation in the exploration of the Lesser Antilles subduction zone. It is due to the presence of only small aligned islands at far distances from the potentially seismogenic interplate area. Several oceanographic cruises were therefore planned that notably allowed the repeated deployment of ocean bottom seismometers; some of them being left for a few months of background seismicity recording. This thesis specifically focuses on two sets of wide-angle seismic data acquired offshore the Dominica and Martinique islands. From their analysis 3D and 2D tomographic models were produced respectively over the forearc region and across the whole subduction complex. These models constrain the plates' seismic structure as well as their geometry. They allow the discussion of how the imaged structures affect the subduction processes and give a first estimation of the downdip width of the seimogenic zone, defined as the segment of the interplate between the backstop and the upper plate's Moho. The joint interpretation of seismic models and earthquake localizations then refine this first assessment. Epicenter distribution from height months of seismic recording shows indeed that seismicity concentrates within the inner forearc region. The strong velocity gradient that characterize its basement suggests it is denser and more rigid than the more deformable outer forearc basement. The updip limit of the seismogenic zone could then lie arcward of the backstop at the contact of the interplate and the seaward end of the inner forearc crust. At depth, interplate earthquake mechanisms are observed between 35 and 45 km depth and interpreted to have occurred at the downdip limit of the seismogenic zone. The latter could reach a depth over 10,km deeper than the contact of the upper plate's Moho with the interplate, and therefore lie within the mantle wedge. All together, these results imply a large downdip width of the seismogenic zone (~70 km) offshore the Dominica and Martinique islands. Further work is, however, needed in order to fully comprehend the ability of the Lesser Antilles subduction zone to produce a possible mega-thrust earthquake. This would necessitate the evaluation of seismic coupling at the interplate contact and the possible segmentation of the seismogenic zone, for instance, due to the subduction of oceanic ridges.
Les Petites Antilles présentent un contexte géodynamique caractérisé par la subduction à très faible vitesse (~2 cm/an) d'une lithosphère océanique âgée (~84-100 Ma). L'activité sismique y est relativement faible, notamment à l'interplaque, où l'aléa sismique lié à un éventuel séisme de méga-chevauchement reste encore mal contraint. Cette thèse se veut être une première étape dans l'évaluation de la capacité de la zone de subduction des Petites Antilles à générer un tel évènement. Dans ce but, ces travaux tentent d'appréhender l'extension en profondeur du domaine sismogène de l'interplaque. Le manque de couverture des stations sismologiques permanentes est un inconvénient majeur dans l'exploration des Petites Antilles. Il s'explique en raison du peu de terres émergées et de leur éloignement de la zone potentiellement sismogène de l’interplaque. La région a donc fait l'objet de plusieurs campagnes océanographiques qui ont permis, notamment, le déploiement de sismomètres fond de mer (OBS); certains instruments étant restés immergés plusieurs mois afin de procéder à une écoute de la sismicité. Les travaux réalisés au cours de cette thèse se sont focalisés sur deux jeux de données de sismique grand-angle acquis au large des îles de la Dominique et de la Martinique. Leur analyse a permis la construction de modèles tomographiques 3D et 2D respectivement à l'échelle de l'avant-arc et de l'ensemble de la subduction. Ces modèles renseignent sur la structure sismique des plaques en convergence ainsi que sur leur géométrie. Ils permettent de discuter le rôle de la structure dans le fonctionnement de la subduction et d’obtenir une première estimation de l'extension en profondeur de la zone sismogène en considérant la portion de l'interplaque comprise entre le butoir et le Moho de la plaque supérieure. Dans un second temps, l'interprétation conjointe des modèles tomographiques et des localisations des séismes locaux a permis d'affiner cette estimation. Huit mois d'enregistrement de la sismicité montrent en effet une concentration des épicentres dans la région interne de l’avant-arc. Celle-ci présente un socle épais, à fort gradient de vitesse interprété comme plus dense et rigide que le socle de la région externe de l'avant-arc, plus déformable. La limite amont de la zone sismogène pourrait donc se situer en retrait de la position du butoir au contact de l'interplaque et de la limite entre ces deux zones de l'avant-arc. En profondeur, des mécanismes interplaques sont observés entre 35 et 45 km et interprétés comme des marqueurs de la limite aval de la zone sismogène. Cette dernière pourrait donc atteindre une profondeur jusqu'à 10 km supérieure à la limite précédemment évoquée et se trouver, par conséquent, au contact du manteau lithosphérique. L'ensemble de ces résultats suggèrent une extension en profondeur de la zone sismogène (i.e une largeur) de près de 70 km face aux îles de la Dominique et de la Martinique. Ces travaux doivent cependant être poursuivis afin d’évaluer pleinement la capacité de la zone de subduction des Petites Antilles à générer un éventuel séisme de méga-chevauchement. Le taux de couplage à l'interplaque doit être précisé ainsi que sa possible segmentation en lien, par exemple, avec à l'entrée en subduction des rides océaniques.

Droits : info:eu-repo/semantics/openAccess
http://archimer.ifremer.fr/doc/00172/28335/26628.pdf
http://archimer.ifremer.fr/doc/00172/28335/

Éditeur(s) : HAL CCSD Elsevier
Résumé : A synthesis of backscatter imagery coupled with a large 3D seismic dataset in the Lower Congo Basin (LCB) reveals a patchy distribution of features interpreted to be associated with fluid seepage from 300 m to 2500 m water depth. With the exception of one region of anomalous backscatter positive-relief mounds, all inferred seep sites occur in negative-relief pockmarks. The extensive 3D seismic dataset in the LCB offers a unique opportunity to study the plumbing system that is feeding surface cold seep systems, and in general, to reconstruct the relationship between tectonics and fluid flow in continental margins. The fluid seeps in the LCB are associated with morphologically, stratigraphically or tectonically controlled focused fluid flow. The integration of the geophysical datasets, backscatter imagery coupled to 3D seismic, clearly indicates that fluid seeps are not randomly distributed, but their seabed organization reflects 1) the location of the underlying structure (reservoir or trap) where the fluids are coming from, 2) the geometry and morphology of the reservoir/trap, and 3) the discontinuities in the sedimentary column along which fluids have migrated. In the LCB seafloor pockmarks are always associated with underlying tectonic structures (fault zones, salt diapirs, polygonal faults) or buried sedimentary bodies (turbiditic channels, erosional surfaces), whereas they never occur above sub-horizontal parallel-stratified fine-grained sediments. Even if triggering processes can not be clearly defined here, we propose a model of seafloor fluid seep organization, which represents a new tool for identifying the geometry of flow pathways and the underlying buried bodies where the fluids are originating from. This qualitative 3D model provides insight into the geohydrologic processes of continental margins.
ISSN: 0025-3227

hal-00406640
https://hal.archives-ouvertes.fr/hal-00406640
DOI : 10.1016/j.margeo.2007.06.003

Éditeur(s) : HAL CCSD Springer Verlag (Germany) Springer Verlag
Résumé : Production rates in horizontal shale gas wells display declines that are controlled by the low permeability and the intrinsic heterogeneity of the shale matrix. We present an original multi-continuum approach that yields a physical model able to reproduce the complexity of the decreasing gas rates. The model describes the dynamics of gas rate as function of the physical reservoir parameters and geometry, while the shale matrix heterogeneity is accounted for by a stochastic description of transmissivity field. From the 3D (Dimensional) problem setting, including the heterogeneous shale matrix, the fractures generated by the hydrofracking operations, as well as the production well characteristics, we establish an effective upscaled 1D model for the gas pressures in fracture and matrix as well as the volumetric flux. We analyse the decline curves behaviour, and we identify the time scales that characterize the dynamics of the gas rate decline using explicit analytical Laplace space solutions of the upscaled process model. Asymptotically, the flux curves decrease exponentially, while in an intermediate regime we find a power-law behaviour, in which the flux scales with a power law in time as (t^{-beta }), where (beta ) reflects the medium heterogeneity. We use this solution to fit a set of real data displaying distinctly different decline trends and study the sensitivity of the model to the reservoir parameters in order to identify their respective controls at the different stages of the decline curve dynamics. Results indicate that the initial value of the gas rate is determined by the transmissivity of the fractures and the initial pressure of the gas in the shale matrix. The latter causes mainly a shift in the entire decline curve. The early time of decline curve shape is primarily controlled by the fracture properties (compressibility and transmissivity). During the main part of the economically valuable production times, i.e. before the production rate drops exponentially, the decline curve is strongly controlled by the properties of the shale rocks including their heterogeneity, which is modelled by two parameters describing the non-Fickian pressure diffusion effects in a stochastic framework.
ISSN: 0169-3913

hal-01217247
https://hal.archives-ouvertes.fr/hal-01217247
DOI : 10.1007/s11242-015-0504-y

Éditeur(s) : HAL CCSD
Résumé : International audience
The Corinth Rift is the most seismically active zone in Europe. The area is characterized by very localized NSextension at a rate of1.5cm/year, the occurrence of frequent and intensive microseismic crises and occasionalmoderate to large earthquakes like in 1995 (Mw=6.1).Since the year 2000, the Corinth Rift Laboratory (CRL, http://crlab.eu) consisting in a multidisciplinary naturalobservatory, aims at understanding the mechanics of faulting and earthquake nucleation in the Rift. Recent studieshave improved our view about fault geometry and mechanics within CRL, but there is still a critical need for a betterknowledge of the structure at depth both for the accuracy of earthquake locations and for mechanical interpretationof the seismicity.In this project, we aim to analyze the complete seismological database (13 years of recordings) of CRL by usingrecently developed methodologies of structural imaging, in order to determine at the same time and with highresolution, the local 3D structure and the earthquake locations.We perform an iterative joint determination of 3D velocity model and earthquake coordinates. In a first step, Pand S velocity models are determined using first arrival time tomography method proposed by Taillandier et al.(2009). It consists in the minimization of the cost function between observed and theoretical arrival times which isachieved by the steepest descent method (e.g. Tarantola 1987). This latter requires computing the gradient of thecost function by using the adjoint state method (Chavent 1974). In a second step, earthquakes are located in thenew velocity model with a non-linear inversion method based on a Bayesian formulation (Gesret et al. 2015). Step1 and 2 are repeated until the cost function no longer decreases.We present preliminary results consisting in: (1) the adjustement of a 1D velocity model that is used as initialmodel of the 3D tomography and (2) a first attempt of the joint determination of 3D velocity model and earthquakelocation.In addition to the tomographic imaging, we perform a preliminary receiver function analysis of teleseismic datarecorded by the broadband stations of the CRL network. The RF analysis should provide the interface depthsbeneath seismometers and increase the imaging resolution of the upper crustal structures provided by the 3Dtomography. In this first attempt, we adjust the 1D velocity model that produces a synthetic RF as similar aspossible to the observed RF for a subset of data. We compare the identified interfaces with structures imaged bythe tomography.
EGU general assembly 2016
Vienne, Austria

hal-01306224
https://hal-mines-paristech.archives-ouvertes.fr/hal-01306224

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
An exploration 3D seismic data set from the Gjallar Ridge off mid-Norway images a giant fluid seep structure 3 × 5 km wide, which connects to late Palaeocene magmatic sills at depth. Two of the pipes that have developed as hydrothermal vents reach all the way to the modern seafloor implying that they either were active much longer than the original hydrothermal activity or have been reactivated. We combine detailed seismic analysis of the northern pipe and sandbox modeling to constrain pipe initiation and propagation. Although both the seismic data and the sandbox models suggest that fluids at depth are focused through a vertical conduit, sandbox models show that fluids ascend and reach a critical depth migration where focused migration abruptly transforms into distributed fluid flow through unconsolidated sediments. This indicates that at this level the sediments are intensely deformed during pipe propagation, creating a V-shaped structure, i.e. an inverted cone at depth and a positive relief anomaly, 5 to 10 m high, at the seafloor, which is clearly identified on 3D seismic data. Comparison of the geometries observed in sandbox modeling with the seismically observed geometries of the Giant Gjallar Vent suggests that the Giant Gjallar Vent may be a proto-fluid seep at an early stage of its development, preceding the future collapse of the structure forming a seafloor depression. Our results imply that the Gjallar Giant Vent can be used as a window into the geological processes active in the deep parts of the Vøring Basin.
ISSN: 0025-3227

hal-00772465
https://hal.archives-ouvertes.fr/hal-00772465
DOI : 10.1016/j.margeo.2012.08.010

Éditeur(s) : HAL CCSD Elsevier
Résumé : A combination of statistical studies on present-day subduction zones and three-dimensional (3D) laboratory models is performed with the aim to clarify the way that plate kinematics control the geometry of the slab and the overriding plate deformation in subduction zones. In 3D laboratory models, the analogue of a two layer linearly viscous lithosphere–upper mantle system is achieved by means of silicon putty glucose syrup tank experiment. The subducting and overriding plate velocities are systematically changed by exploring the variability field of natural plate kinematics. Both statistical and modelling approaches recognize the importance of overriding plate motion on subduction process behavior: (1) trenches migrate at a rate close to the overriding plate motion, but always move slower than the overriding plates. The mechanism at work is a direct consequence of “slab anchoring” opposed by both lithosphere and mantle viscous resistance and is responsible for overriding plate deformation and slab geometry variability. (2) An overriding plate shortens when the overriding plate moves toward the trench and conditions that are favourable for overriding plate extension are created when the overriding plate moves away from the trench. (3) Shallow and steep dips are found if the overriding plate moves toward and away from the trench, respectively.
ISSN: 0012-821X

hal-00406598
https://hal.archives-ouvertes.fr/hal-00406598
DOI : 10.1016/j.epsl.2007.02.004

Éditeur(s) : HAL CCSD Société géologique de France
Résumé : In the western Mediterranean Sea, the Liguro-Provençal Basin (LPB) is a key area for studying passive margins because of its recent formation and abundance of onshore and offshore data. The Nerthe area located in the northern margin of LPB provides the unique continuous Oligo-Miocene deposits contemporaneous of the transition rifting to drifting. However, the age of the deposits remains debated and the link between outcrops and offshore seismic data is poorly constrained. The purpose of this paper is double. First, we intend to propose a new chronostratigraphic frame based on bio- (planktonic foraminifera, calcareous nannofossils) and magneto-stratigraphy. Second, we aim to make, through the integration of new highly time-resolved seismic data and field works, a coherent onshore-offshore link concretized by a 3D geological model. The new temporal and spatial data presented in this paper allow correlating the Oligo-Miocene sequences, defining their geometry and specifying precisely the timing of syn- and post-rift stages. The first marine transgression is now precisely dated latest Chattian within the syn-rift deposits and appears to be synchronous with the first marine deposits in the offshore wells and other marginal basins. The transition from syn-rift to post-rift appears to last 3.3 Ma at maximum, between 21.8 and 18.5 Ma (late Aquitanian to early Burdigalian). It is underlined by two major erosional unconformities bearing a hiatus of around 1 Ma. The post-rift started with a major marine transgression that is now dated from middle Burdigalian, at around 18.5 Ma, as elsewhere in the LPB. Contrarily to recent proposals, the post-rift deposits are widely represented on the northeastern margin of the "Golfe du Lion". There, the subsidence of the margin was low during the syn-rift and the transitional periods and high during the post-rift. The onset of this high post-rift subsidence appears to be synchronous with the slowdown of the Corsica-Sardinia block (CSb) motion.
EISSN: 1777-5817

hal-00558497
https://hal.archives-ouvertes.fr/hal-00558497
DOI : 10.2113/gssgfbull.181.6.591

Éditeur(s) : HAL CCSD National Research Council Canada
Résumé : International audience
Back-arc extension in the Aegean, which was driven by slab rollback since 45 Ma, is described here for the first time in two stages. From Middle Eocene to Middle Miocene, deformation was localized leading to i) the exhumation of high-pressure metamorphic rocks to crustal depths, ii) the exhumation of high-temperature metamorphic rocks in core complexes and iii) the deposition of sedimentary basins. Since Middle Miocene, extension distributed over the whole Aegean domain controlled the deposition of onshore and offshore Neogene sedimentary basins. We reconstructed this two-stage evolution in 3D and four steps at Aegean scale by using available ages of metamorphic and sedimentary processes, geometry and kinematics of ductile deformation, paleomagnetic data and available tomographic models. The restoration model shows that the rate of trench retreat was around 0.6 cm/y during the first 30 My and then accelerated up to 3.2 cm/y during the last 15 My. The sharp transition observed in the mode of extension, localized versus distributed, in Middle Miocene correlates with the acceleration of trench retreat and is likely a consequence of the Hellenic slab tearing documented by mantle tomography. The development of large dextral NE-SW strike-slip faults, since Middle Miocene, is illustrated by the 450 Km-long fault zone, offshore from Myrthes to Ikaria and onshore from Izmir to Balikeshir, in western Anatolia. Therefore, the interaction between the Hellenic trench retreat and the westward displacement of Anatolia started in Middle Miocene, almost 10 Ma before the propagation of the North Anatolian Fault in the North Aegean.
ISSN: 0008-4077

insu-01271296
https://hal-insu.archives-ouvertes.fr/insu-01271296
https://hal-insu.archives-ouvertes.fr/insu-01271296/document
https://hal-insu.archives-ouvertes.fr/insu-01271296/file/Brun%20etal%20CJES%202016.pdf
DOI : 10.1139/cjes-2015-0203

Éditeur(s) : HAL CCSD
Résumé : International audience
Aegean extension is a process driven by slab rollback that, since 45 Ma, shows a twostage evolution. From Middle Eocene to Middle Miocene it is accommodated by localized deformation leading to i) the exhumation of high-pressure metamorphic rocks from mantle to crustal depths, ii) the exhumation of high-temperature rocks in core complexes and iii) the deposition of Paleogene sedimentary basins. Since Middle Miocene, extension is distributed over the whole Aegean domain giving a widespread development of onshore and offshore Neogene sedimentary basins. We reconstructed this two-stage evolution in 3D at Aegean scale by using available ages of metamorphic and sedimentary processes, geometry and kinematics of ductile deformation, paleomagnetic data and available tomographic models. The restorationmodel shows that the rate of trench retreat was around 0.6 cm/y during the first 30 My and then accelerated up to 3.2 cm/y during the last 15 My. The sharp transition observed in the mode of extension, localized versus distributed, which occurred in Middle Miocene correlates with the acceleration of trench retreat and is more likely a consequence of the Hellenic slab tearing documented by mantle tomography. The development of large dextral NE-SW strike-slip faults during the second stage of Aegean extension, since Middle Miocene, is illustrated by the 450 Km-long fault, recently put in evidence, offshore from Myrthes to Ikaria and onshore from Izmir to Balikeshir, in western Anatolia. Therefore, the interaction between the Hellenic trench retreat and the westward displacement of Anatolia started in Middle Miocene,almost 10 Ma before the propagation of the North Anatolian Fault in the North Aegean. This raises a fundamental issue concerning the dynamic relationship between slab tearing and Anatolia displacement
ISSN: 0438-9557

insu-01574659
https://hal-insu.archives-ouvertes.fr/insu-01574659
https://hal-insu.archives-ouvertes.fr/insu-01574659/document
https://hal-insu.archives-ouvertes.fr/insu-01574659/file/EFFECTS_OF_SLAB_ROLLBACK_ACCELERATION_ON_AEGEAN_EX.pdf
DOI : 10.12681/bgsg.11697

É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
Résumé : Gravitropism is a major process in vertical growth as it allows the positioning of plant in the gravity field and thus, for trees, it takes part in the light foraging strategies. Without any gravitropism, trees would adopt a weeping habit or fall on the ground due to the increasing with growth supported weight, especially at the advance regeneration stage (sapling stage) in the understorey, where tree are very slender and poorly stiff. Most of the gravitropic studies are made by physiologists in order to understand the process and thus are associated to complex laboratory methods (genetic engineering, anatomical studies, accurate 2D ou 3D measurements of form changes with time). In order to study the ecological relevance and the specific diversity of gravitropic performances in natural forests, we developed a framework that uses different kinds of datas. Synchronical analysis has been done i) at the population level in large sample of trees and ii) at the tree level from destructive measurement of the dissymmetry maturation strains known to be the main motor of the posture control reaction of trees. In both cases, biomechanical models allows to estimate some aspects of the dynamic gravitropic process. Diachronical analysis has been done in natural conditions on permanent plots, or in a more traditional way, by observing in greenhouse experiments the righting movements of plants artificially tilted. In the first case, the problem is the high variability of natural disturbances and the slowness of growth and movements in the dark understorey. The second choice allows a measurement of a capacity (i.e. the ability of the species to react in quite extreme conditions of verticality disturbance but good for growth light and nutrient conditions) that is more easy to analyse, but is maybe not so relevant for estimating the ecological strategy (some species could have developed a high efficience of gravitropic reaction in natural conditions by avoiding the mechanical disturbance or adapting the reaction to slow growth and changing environments in natural conditions). All approaches use the same basic biomechanical modelling, i) to analyse simultaneously the reaction and the disturbance (due to the weight increase) during growth, or ii) to combine different parameters (geometry, growth, weight, maturation strains) involved in shape changes and posture control. They give a classification of the species according to their gravitropic performances. Results on 15 species are discussed, with relation to shade-tolerance or avoidance.
6th Plant Biomechanics Conference
Cayenne, France

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