É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é : 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 Seismological Society of America
Résumé : International audience
The objective of this paper is to provide a better understanding of the relationship between the microseismicity, active tectonics, and crustal structures in the southwest Yukon Territory, Canada, in order to improve seismic‐hazard assessments in this region. We utilize data from a new dense seismic array that was deployed in the southwest Yukon in the summer of 2010. Data from the new seismic array significantly improve the magnitude completeness level in the region from ML 3.0 to 1.0. We analyze 980 events ranging in magnitudes from ML 0.2 to 4.7, at depths from 0 to 35 km. Relocation analysis using the progressive multiple event location shows that seismicity is concentrated in four main areas: (1) Yakutat block northern boundary-Fairweather fault, (2) Duke River fault, (3) southern Denali fault, and (4) a previously unrecognized northeast trending cluster that may highlight a previously unknown active fault. This cluster may contribute to stress and strain transfer inland from the Yakutat block region.
ISSN: 0037-1106

hal-00940935
https://hal.archives-ouvertes.fr/hal-00940935
DOI : 10.1785/0120130068

Éditeur(s) : HAL CCSD
Résumé : This work provides new constraints on the interactions between deformation and percolation of melt or fluids, as well as on the implications of these processes for the seismic properties of the mantle wedge. It is based on the analysis of spinel peridotites from the massif of Ronda (Spain) and two xenolith suites from active subduction zones (Kamchatka, Papua New Guinea). A coupled structural, petrological and geochemical study of these samples shows that they underwent a reactive percolation of melts or Si-rich fluids, which was synchronous to a deformation event occurring under high temperature and low stress conditions, consistent with the PT conditions at the base of the lithosphere or in the asthenosphere. This reactive percolation is responsible for enrichment in pyroxenes, mainly orthopyroxene, which is often localized in bands parallel to the high-temperature foliation. This enrichment is associated with grain size reduction and dispersion of the crystallographic orientation of olivine. The dominant slip system in olivine is {0kl}[100], which results in fast S-wave polarization parallel to the flow direction in the mantle. The decrease in the intensity of the olivine crystal preferred orientations associated with the enrichment in pyroxene results in significant decrease of the anisotropy that may induce an underestimation of the thickness of the anisotropic layer by up to 33%. The observed orthopyroxene enrichment also lowers the Vp/Vs ratio, but cannot explain Vp/Vs < 1.7 mapped locally in the fore-arc mantle in Japan and the Andes. Such low Vp/Vs ratios may however be explained by considering the intrinsic anisotropy of the peridotites, which is generally ignored in large-scale Vp/Vs ratio mapping of the mantle wedge. Infrared analyses show that olivine from the both xenolith collections contains less water than the theoretical saturation calculated for their estimated equilibrium temperature in the spinel stability field. These low water content are similar to those observed in spinel peridotites from other subduction zones and probably record both the low solubility of water in olivine at relatively low pressure and dehydration during exhumation of the xenoliths. These measured water contents as well as theoretical saturation estimations are not sufficient to change the dominant slip direction of in olivine from [100] to [001].
Ce travail apporte de nouvelles contraintes sur les interactions entre déformation et processus d'hydratation et de percolation de magma ou de fluides dans le manteau à l'aplomb d'une zone de subduction et leurs implications sur les propriétés sismiques dans le coin mantellique. Il se base sur l'analyse de péridotites à spinelle provenant du massif de Ronda (Espagne) et deux séries de xénolites issues de zones de subduction actives (Kamchatka, Papouasie-Nouvelle-Guinée). L'étude structurale, pétrologique et géochimique de ces échantillons montrent qu'ils ont subit une percolation réactive de magma ou de fluide synchrone d'une déformation de haute température basse contrainte cohérente avec les condition PT de la base de la lithosphère ou de l'asthénosphère. Cette percolation réactive est responsable d'un enrichissement en pyroxènes, qui est souvent localisé dans des bandes parallèles à la foliation. Cet enrichissement est associé à une dispersion de l'orientation cristallographique de l'olivine. Les systèmes de glissement dominant dans l'olivine sont dans tous les cas {0kl}[100], ce qui implique que la direction de polarisation rapide des ondes S dans la partie supérieure du coin mantellique est parallèle à la direction d'écoulement dans le manteau. La décroissance de l'intensité des OPR de l'olivine associée à l'enrichissement en pyroxènes a pour conséquence une réduction non négligeable de l'anisotropie sismique qui peut induire jusqu'à 33% d'erreur sur l'interprétation de l'épaisseur de la couche anisotrope. Un enrichissement en orthopyroxène peut aussi entraîner une baisse du rapport Vp/Vs, mais ne peut expliquer les rapports Vp/Vs <1,7 cartés dans certains avant-arcs. Cependant de telles valeurs peuvent être expliquées si l'anisotropie des péridotites du coin mantellique est prise en compte. Les analyses par spectroscopie infra-rouge montrent que les olivines des deux séries de xénolites étudiées contiennent moins d'eau que la saturation théorique calculée pour leurs températures d'équilibre dans le domaine du spinelle. Ces faibles concentrations en eau sont similaires à celles observées dans les olivines des péridotites à spinelle des autres zones de subduction. Elles enregistrent probablement la faible solubilité de l'eau dans l'olivine à des pressions relativement basses et la déshydratation au cours de l'exhumation des xénolites. Ces teneurs en eau mesurées ainsi que les estimations de la saturation théorique ne sont pas suffisantes pour changer la direction de glissement dominante dans l'olivine de [100] à [001].
https://tel.archives-ouvertes.fr/tel-00689832

tel-00689832
https://tel.archives-ouvertes.fr/tel-00689832
https://tel.archives-ouvertes.fr/tel-00689832/document
https://tel.archives-ouvertes.fr/tel-00689832/file/These-Soustelle2010.pdf

Éditeur(s) : HAL CCSD Geological Society of London
Résumé : International audience
Low-velocity zones ('bright spots') imaged by the INDEPTH seismic experiment in southern Tibet are extensively interpreted as widespread partial melt within the crust, which has given a strong support for the channel flow model. These suggest that a continuous seismic low-velocity zone underlies Tibet on the large scale. Here we take advantage of the Hi-CLIMB seismic experiment which includes a dense south-north profile and a lateral 2D seismic network to assess the vertical and the horizontal extension of low-velocity zones in southern Tibet. Several approaches including migration, amplitude analysis and waveform inversion of receiver functions are performed to detect crustal low-velocity zones using this new seismological dataset. Our results reveal localized and discontinuous low-velocity zones in Tibet. They indicate that the vertical extension of the low-velocity zones is about 10 km, and their maximum horizontal length appears to be c. 50 km. Our study suggests a partial correlation between the location of these low-velocity zones and the spatial distribution of Tibetan grabens. These results, especially the non-continuity of low-velocity zones, together with the observed regular value of mean crustal VP/VS ratio, question the existence of widespread partial melt of the southern Tibetan crust and, therefore, the viability of the channel flow model.
Growth and Collapse of the Tibetan Plateau

hal-00853971
https://hal.archives-ouvertes.fr/hal-00853971
DOI : 10.1144/SP353.6

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Mongolia represents the northernmost area affected by the India–Asia collision, and it is actively deformed along transpressive belts closely associated with large-scale strike-slip faults. The active and past mantle flow beneath this region is, however, poorly known. In order to investigate deep mantle deformation beneath central Mongolia and its relation with the surrounding major structures such as the Siberian craton, the Gobi–Altay belt and the Baikal rift, a NS-trending profile of broadband seismic stations has been deployed in the summer 2003 from the southern Siberian craton to the Gobi–Altay range, crossing the entire Hangay dome. Mantle flow is deduced from the splitting of teleseismic shear waves such as SKS phases. In eastern Mongolia, the permanent station ULN in Ulaanbaatar reveals the presence of two anisotropic layers, the upper one being oriented NE–SW, close to the trend of lithospheric structures and the lower one NW–SE, close to the trend of Eurasia absolute plate motion. Along the NS profile in central Mongolia, seismic anisotropy deduced from SKS splitting reveals a homogeneous NW–SE trending structure, fully consistent with the observations made in the Altay–Sayan in western Mongolia. The observed delay times of 1.5 to more than 2.0 s favor consistent mantle flow over large mantle thicknesses. Since the lithosphere is less than 100 km thick beneath central Mongolia and since the observed fast directions are parallel to the trend of the lithospheric structures but also close to the trend of the absolute plate motion, we propose that both the lithosphere and the asthenosphere may join their anisotropic effects beneath central Mongolia to explain the large delay times. Although GPS vectors represent the instantaneous displacement of the Earth's surface and SKS splitting the time and vertical integration of finite strain at depth, we use the opportunity of the dense geodetic measurements available in this region to discuss the anisotropy pattern in term of present-day deformation. In the Eurasia-fixed reference frame, GPS and SKS both depict a similar trend beneath central Mongolia, suggesting a lithospheric block “escaping” toward the east that could orient olivine a-axes in the upper mantle, within a transpressive tectonic regime. A different behaviour is observed in western Mongolia: the GPS vectors trend NS, close to the regional compression direction, whereas the fast SKS directions trend EW, suggesting a tectonic regime close to a mode of axial shortening, generating the development of an EW-trending foliation at depth. We therefore propose that Mongolia is a place where active and frozen lithospheric deformation may add their effects, together with the sublithospheric flow. In the three sources of anisotropy inferred, a fundamental role is played by the Siberian craton that acted as an undeformable core of the continent through time: the frozen Paleozoic and Mesozoic structures wrap around the craton, building up the fast anisotropic direction pattern; the present-day sublithospheric flow induced by the plate motion is likely deflecting around its deep roots; finally, the present-day tectonic regime appears to be controlled by the presence of the craton to the north.
ISSN: 0012-821X

hal-00407853
https://hal.archives-ouvertes.fr/hal-00407853
https://hal.archives-ouvertes.fr/hal-00407853/document
https://hal.archives-ouvertes.fr/hal-00407853/file/Barruol_SKS_Mongolia_EPSL2008_hal.pdf
DOI : 10.1016/j.epsl.2008.07.027

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The Mosha and North Tehran faults correspond to the nearest seismic sources for the northern part of the Tehran megacity. The present-day structural relationships and the kinematics of these two faults, especially at their junction in Lavasanat region, is still a matter of debate. In this paper, we present the results of a morphotectonic analysis (aerial photos and field investigations) within the central part of the Mosha and eastern part of the North Tehran faults between the Mosha valley and Tehran City. Our investigations show that, generally, the traces of activity do not follow the older traces corresponding to previous long-term dip-slip thrusting movements. The recent faulting mainly occurs on new traces trending E-W to ENE-WSW affecting Quaternary features (streams, ridges, risers, and young glacial markers) and cutting straight through the topography. Often defining en-echelon patterns (right- and left-stepping), these new traces correspond to steep faults with either north- or south-dipping directions, along which clear evidences for left-lateral strike-slip motion are found. At their junction zone, the two sinistral faults display a left-stepping en-echelon pattern defining a positive flower structure system clearly visible near Ira village. Further west, the left-lateral strike-slip motion is transferred along the ENE-WSW trending Niavaran fault and other faults. The cumulative offsets associated with this left-lateral deformation is small compared with the topography associated with the previous Late Tertiary thrusting motion, showing that it corresponds to a recent change of kinematics.
ISSN: 0040-1951

hal-00617778
https://hal.archives-ouvertes.fr/hal-00617778
DOI : 10.1016/j.tecto.2010.09.013

É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 American Geophysical Union
Résumé : International audience
Giant earthquake (moment magnitude M-w >= 8.5) forecasts for subduction zones have been empirically related to both tectonic stresses and geometrical irregularities along the subduction interface. Both of these controls have been suggested as able to tune the ability of rupture to propagate laterally and, in turn, exert an important control on giant earthquake generation. Here we test these hypotheses, and their combined influence, by compiling a dataset of trench fill thickness (a proxy for smoothing of subducting plate relief by sediment input into the subduction channel) and upper plate strain (a proxy for the tectonic stresses applied to the subduction interface) for 44 segments of the global subduction network. We statistically compare relationships between upper plate strain, trench sediment thickness and maximal earthquake magnitude. We find that the combination of both large trench fill (>= 1 km) and neutral upper plate strain explains spatial patterns of giant earthquake occurrence to a statistically significant degree. In fact, the concert of these two factors is more highly correlated with giant earthquake occurrence than either factor on its own. Less frequent giant earthquakes of lower magnitude are also possible at subduction zones with thinner trench fill and compressive upper plate strain. Extensional upper plate strain and trench fill < 0.5 km appear to be unfavorable conditions, as giant earthquakes have not been observed in these geodynamical environments during the last 111 years. Citation: Heuret, A., C. P. Conrad, F. Funiciello, S. Lallemand, and L. Sandri (2012), Relation between subduction megathrust earthquakes, trench sediment thickness and upper plate strain, Geophys. Res. Lett., 39, L05304, doi: 10.1029/2011GL050712.
ISSN: 0094-8276

hal-00745636
https://hal.archives-ouvertes.fr/hal-00745636
https://hal.archives-ouvertes.fr/hal-00745636/document
https://hal.archives-ouvertes.fr/hal-00745636/file/heuretGRL2012.pdf
DOI : 10.1029/2011GL050712

Éditeur(s) : HAL CCSD Springer-Verlag
Résumé : The seismology is the most effective method to explore the Structure of subduction zones to great depth. The distinguishing feature of the mantle in the subduction regions is the presence of hydrated phases, which transport water into the Earth's interior and release it with dramatic local consequences, triggering earthquakes and melting. The seismological detection of these hydrous phases and geo-dynamic interpretation of flow in the hydrous mantle depend on knowledge of the anisotropic elastic properties and the characteristics of the wave propagation in anisotropic media. In this paper we briefly recall the distinguishing features of anisotropic wave propagation and the observable parameters. We suggest the Vp/Vs ratio is a physically sound parameter than can be observed by seismology in anisotropic regions of the Earth, whereas the Poisson's ratio, which is often quoted, is not directly observable and does not correspond to the characteristics of wave propagation in an anisotropic or isotropic medium. We report for the first time the ratios of Vp/Vs1 and Vp/Vs2, where Vs1 is the fastest and Vs2 slowest S-wave velocities of an anisotropic media. We present the current knowledge of the anisotropic seismic properties of hydrous minerals in the upper mantle, transition zone, and lower mantle that are stable along low temperature geotherms associated with subduction, and identify which minerals are likely to influence seismological observations because they have very high volume fractions, or very high anisotropies, or both of these.;In the upper mantle antigorite and talc are exceptionally anisotropic (Vp 71% Vs 68% and Vp 65%, Vs 68%, respectively) and chlorite is also very anisotropic for S-waves (Vp 35%, Vs 76%). Comparatively less anisotropic are hornblende (Vp 27%, Vs 31%), used as proxy for tremolite in subduction zones, and clinohumite (21.8%, Vs 15.9%). Brucite at 4 GPa (Vp 26.5%, Vs 30.9%) and the dense hydrogen magnesium silicate (DHMS) phase A at 9GPa (Vp 9.3%, Vs 17.6%) are the only hydrous minerals stable in the upper mantle that have had their elastic properties measured at in situ mantle pressures. Except for the phase A, all these minerals are more anisotropic than olivine for at least one parameter (Vp, Vs, Vp/Vs1, or Vp/Vs2). In the transition zone the major phases hydrous wadsleyite and ringwoodite C have moderate to weak anisotropies (Vp 16.3%, Vs 16.5%, and Vp 1.9%, Vs 4.4%, respectively). The DHMS Superhydrous B has a moderate anisotropy (Vp 6.9%, Vs 11.6%). At greater depth the DHMS phase D Lit 24 GPa (Vp 10.8%, Vs 18.0%) is the only hydrous phase that can transport hydrogen from the transition zone into the lower mantle. The phase D is more anisotropic for S than P waves like many of the hydrous phases. From these data it is clear that hydrous phases are in general very anisotropic. However, pressure can play a strong role in reducing anisotropy. It is the case for brucite and talc, in which increasing pressure from ambient to 4 GPa reduces the anisotropy by about 50% for both P and S waves. In contrast, the anisotropy of the DHMS phase A does not change significantly with pressure.;Our picture of the seismic anisotropy of hydrated minerals remains incomplete; the elastic properties of many have not been measured even at ambient conditions (e.g., 10 angstrom phase and phase E) or not measured in their true elastic symmetry (e.g., clinochlore). The majority of hydrous minerals have not been measured at high pressure. and none related to hydrated mantle have been measured at elevated temperature. We have shown in the few cases where hydrous minerals have been measured as a function of pressure, that this variable has an important effect oil the velocity distribution and in most cases reduces the degree of anisotropy, hence we Would expect seismic anisotropy to play a key role in the determination of the shallow structure of subduction zones in the upper mantle.
Subduction zone geodynamics

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

É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 Elsevier
Résumé : International audience
In this paper, we present the first estimate of the Holocene deformation along the southern front of Gibraltar arc (Morocco) and the first field constraints on the local 1755 CE Fes-Meknes surface rupturing earthquake which could be associated to the “Great Lisbon Earthquake” (M > 8.5) in November 1st, 1755. Using satellite imagery, aerial photographs and field investigations, we carried out a morphotectonic study along the ~ 150 km-long Southern Rif Front (SRF) to identify the most recent evidences of tectonic activity. Analyzed offset alluvial deposits confirm that (i) the last ~ 5 ka cumulative deformation leading to a slip rate of ~ 3.5 ± 1 mm/yr for this segment of the SRF is consistent with the GPS derived horizontal shortening rate of 2–4 mm/yr and (ii) a recent major earthquake ruptured a ~ 30 km-long segment along the SRF. Based on deposits dating and historical seismicity we propose that this seismic event occurred in 1755 as a local earthquake. Even though this 1755 local event cannot be considered as a strong aftershock of the main Lisbon seismic event (M > 8.5), their temporal closeness, their occurrence under the same convergent stress regime (~ NNW-SSE-oriented compression) and the fact that Fes-Meknes area was strongly shaken during the Lisbon earthquake, raises the question of the possible triggering of the Fes earthquake. Anyway, our new results suggest that most of the Nubia-Rif belt convergence is accommodated by the SRF, making it potentially the most destructive structure of the Rif.
ISSN: 0040-1951

hal-01622048
https://hal.archives-ouvertes.fr/hal-01622048
DOI : 10.1016/j.tecto.2017.05.036

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We investigated the response of the crystal structure, lattice parameters, and unit-cell volume of hydrous layered silicate phlogopite at conditions relevant to subduction zone settings. We have used first principles simulation based on density functional theory to calculate the equation of state and full elastic constant tensor. Based on the generalized gradient approximation, the full single crystal elastic constant tensor with monoclinic symmetry shows significant anisotropy with the compressional elastic constants: c11 = 181 GPa, c22 = 185 GPa, c33 = 62 GPa, the shear elastic constants c44 = 14 GPa, c55 = 20 GPa, c66 = 68 Ga, and c46 = −6 GPa; the off diagonal elastic constants c12 = 48 GPa, c13 = 12 GPa, c23 = 12 GPa, c15 = −16 GPa, c25 = −5 GPa and c35 = −1 GPa at zero pressure. The elastic anisotropy of phlogopite is larger than most of the layered hydrous phases relevant in the subduction zone conditions. The shear anisotropy, AVS for phlogopite is ∼77% at zero pressure condition and although it decreases upon compression it remains relatively high compared to other hydrous phases relevant in the subduction zone settings. We also note that the shear elastic constants for phlogopite are relatively low. Phlogopite also has a high isotropic bulk VP/VS ratio ∼2.0. However, the VP/VS ratio also exhibits significant anisotropy with values as low as 1.49. Thus, phlogopite bearing metasomatized mantle could readily explain unusual VP/VS ratio as observed from seismological studies from the mantle wedge regions of the subduction zone.
ISSN: 0031-9201

hal-01115670
https://hal.archives-ouvertes.fr/hal-01115670
DOI : 10.1016/j.pepi.2014.05.004

Éditeur(s) : HAL CCSD Springer Verlag (Germany)
Résumé : International audience
Surface displacements due to temporal changes in environmental mass redistributions are observable in the coordinate time series of many Global Navigation Satellite System (GNSS) sites. In this study, we investigated the effect of loading on estimates of tectonic velocity computed from campaign-style GNSS observations. The study region is in the Pyrenees mountain range between France and Spain (ResPyr campaigns). In this area, seismic activity is continuous and moderate and the expected amplitude of the horizontal tectonic velocity is less than 0.5 mm/yr. In order to determine the velocity, 4 sparse GNSS campaigns were carried out from 1995 to 2010. Considering this small rate of deformation, loading phenomena can contribute a non-negligible artifact to the velocity computation that could affect our geodynamical interpretation. In this investigation, we specifically considered the atmospheric, hydrological, and non-tidal ocean loading phenomena. The computed loading deformations for this region show the horizontal displacements are dominated by the non-tidal ocean loading (maximum 4 mm for the North and 3.1 mm for the East components); the main vertical contributions are due to the atmospheric and continental water storage loading (maximum 14.3 for the atmosphere and 8.1 mm for the hydrology, respectively). We have found that the dominant loading effect on the horizontal velocity is the non-tidal ocean loading (mean of 0.11 mm/yr), whereas the vertical component is dominated by the hydrological loading (mean of 0.21 mm/yr). Since the study area is in a mountainous region, we also analyzed the difference between the atmospheric and the topography dependent atmospheric loading models at our GNSS campaign sites. We did not find any significant difference between the two atmospheric loading models in terms of horizontal velocity. Finally, we performed simulations to identify the optimum timing and frequency of future GNSS campaigns in this area that would minimize the loading effects on tectonic velocity estimates.
ISSN: 0039-3169

hal-00963629
https://hal.archives-ouvertes.fr/hal-00963629
DOI : 10.1007/s11200-012-0458-2

Éditeur(s) : HAL CCSD Société géologique de France
Résumé : International audience
Guadeloupe island (West French Indies) is one of the twenty islands that compose the Lesser Antilles arc, which results from the subduction of the Atlantic ocean plate beneath the Caribbean one. The island lies in a complex volcano-tectonic system and the need to understand its geological context has led to numerous on- and offshore geophysical investigations. This work presents a compilation and the processing of available, on-land, airborne and marine, gravity and magnetic data acquired during the last 40 years on Guadeloupe islands and at the scale of the Lesser Antilles arc. The overall dataset provides new Bouguer and reduced to the pole magnetic anomaly maps at the highest achievable resolution. Regionally, the main central negative gravity trend of the arc allows defining two subsident areas. The first one is parallel to the arc direction (~N160°E) to the north, whereas the second unexpected southern one is oriented parallel to oceanic ridges (N130°E). Along the Outer arc, the long wavelength positive anomaly is interpreted, at least along the Karukera spur, as an up-rise of the volcanic basement in agreement with the seismic studies. To the NE of Guadeloupe, the detailed analysis of the geophysical anomalies outlines a series of structural discontinuities consistent with the main bathymetric morphologies, and in continuity of the main fault systems already reported in this area. Based on geophysical evidences, this large scale deformation and faulting of the Outer arc presumably primarily affects the Atlantic subducting plate and secondarily deforms the upper Caribbean plate and the accretion prism, as evidenced in bathymetry as well as on the islands. At the scale of Guadeloupe island, combined gravity and magnetic modeling has been initiated based on existing interpretation of old seismic refraction profiles, with a general structure in three main layers. According to our geophysical anomalies, additional local structures are also modeled in agreement with geological observations: i) the gravity and magnetic signals confirm an up-rise of the volcanic basement below the limestone platforms outcropping on Grande-Terre island ; ii) the ancient volcanic complexes of Basse-Terre island are modeled with high density and reverse magnetized formations; iii) the recent volcanic centre is associated with formations consistent with the low measured density and the underlying hydrothermal system. The E-W models coherently image a NNW-SSE depression structure in half-graben beneath Basse-Terre island, its western scarp following the arc direction in agreement with bathymetric and seismic studies to the north of the island. The so-defined depressed area, and particularly its opening in half-graben toward the SW, is interpreted as the present-day front of deformation of the upper plate, associated with the recent volcanic activity on and around Guadeloupe. Based on this regional deformation model, perspectives are given for further integrated investigation of key targets to address the internal structure and evolution of the Lesser Antilles arc and Guadeloupe volcanic system.
EISSN: 1777-5817

hal-00853352
https://hal.archives-ouvertes.fr/hal-00853352
DOI : 10.2113/gssgfbull.184.1-2.77

Éditeur(s) : HAL CCSD
Résumé : International audience
Seismic anisotropy is a ubiquitous feature of the subcontinental mantle. This can be inferred both from direct seismic observations of shear wave splitting from teleseismic shear waves, as well as the petrofabric analyses of mantle nodules from kimberlite pipes. The anisotropy is principally due to the strain-induced lattice preferred orientation (LPO) of olivine. The combined use of these mantle samples, deformation experiments on olivine, and numerical modeling of LPO, provides a critical framework for making inferences about mantle deformation from observed seismic anisotropy. In most cases there is a close correspondence between mantle deformation derived from seismic observations of anisotropy, and crustal deformation, from the Archean to the present. This implies that the mantle plays a major, if not dominant role in continental deformation. No clear evidence is found for a continental asthenospheric decoupling zone, suggesting that continents are probably coupled to general mantle circulation.
Mantle Petrology: Field Observations and high Pressure experimentations: a Tribute to Francis R. (Joe) Boyd

hal-01391514
http://hal.univ-reunion.fr/hal-01391514
http://hal.univ-reunion.fr/hal-01391514/document
http://hal.univ-reunion.fr/hal-01391514/file/Silver_et_al_BoydVolume 1999_hal.pdf

Éditeur(s) : HAL CCSD Société Géologique de France
Résumé : Based on a High-Resolution 3D seismic block acquired in the Gulf of Lion in 2004–2005 we investigated fluid pipes and pockmarks on the top of the interfluve between the Hérault canyon and the Bourcart canyon both created by turbidity currents and gravity flows from the shelf to the deep basin in the north-western Mediterranean Sea. Combining the geometry of the potential fluid pipes with the induced deformation of surrounding sediments leads then to the ability to differentiate between potential fluid sources (root vs source) and to better estimate the triggering mechanisms (allochtonous vs. autochtonous cause). We linked together a set of derived attributes, such as Chaos and RMS amplitude, to a 3D description of pipes along which fluids may migrate. As previously shown in other basins, the induced deformation, creating cone in cone or V-shaped structures, may develop in response to the fluid pipe propagation in unconsolidated sediments in the near surface. The level at the top of a cone structure is diachronous. It means that stratigraphic levels over this surface are deformed at the end of the migration. They collapse forming a depression called a pockmark. These pipes are the result of repeated cycles of fluid expulsion that might be correlated with rapid sea-level rise instead of sediment loading. The most recent event (MIS 2.2 stage) has led to the formation of a pockmark on the modern seafloor. It has been used as a reference for calculating the effect of a rapid sea-level rise on fluid expulsion. As all physical and geometrical parameters are constrained, we were able to define that a + 34 m of sea level rise may account for triggering fluid expulsion from a very shallow silty-sandy layer at 9 m below seafloor since the last glacial stage. This value is consistent with a sea level rise of about 102 m during this period. This study shows that the episodic nature of fluid release resulted from hydromechanical processes during sea-level rise due to the interactivity between high pressure regimes and principal in situ stresses.
ISSN: 0037-9409

hal-01667241
https://hal.archives-ouvertes.fr/hal-01667241
DOI : 10.1051/bsgf/2017190

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The seismic anisotropy of the deep Earth is reviewed as a profile from the upper mantle to the solid inner core at the centre of the Earth. The upper mantle is by far the most anisotropic region of the Earth, followed the D” layer above the core mantle boundary. In contrast it is shown that many minerals that are present in the upper mantle, transition zone, lower mantle, D” layer and solid inner core are elastically anisotropic to different degrees. The anisotropy of hydrous phases present in subduction zones is briefly introduced. The basic concepts of single crystal and polycrystalline elasticity and the extrapolation of elastic properties to high temperature and pressure are presented. The specific features of elastic wave propagation in an anisotropic medium of any arbitrary symmetry are illustrated using mantle and inner core phases. The roles of crystal preferred orientation, water and melt in producing seismic anisotropy in the upper mantle are discussed.
Treatese of Geophysics, vol.2

hal-00408321
https://hal.archives-ouvertes.fr/hal-00408321
DOI : 10.1016/B978-044452748-6.00045-6

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : GPS velocities and seismicity across the Betic-Rif Arc structural domains (Morocco and Iberia) provide a basis to evaluate present-day seismotectonic processes between different deformation belts. The results show asymmetric movements in the complex Alboran system accommodating the convergence between the African (Nubian) and Eurasian plates. While the Betic Mountains are attached to Iberia, moving toward the southeast with respect to Africa, the Rif is divided into three blocks with distinct displacements relative to Nubia: (1) the Tangier block moving southeastward, (2) the Central Rif block moving SSW, and (3) the Oriental Rif block undergoing clockwise rotation. GPS-derived motions decrease in rate from the Rif nappes complex to the foreland.
ISSN: 0010-7999

hal-00411144
https://hal.archives-ouvertes.fr/hal-00411144
DOI : 10.1016/j.crte.2008.02.003

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

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