Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Laboratory models have been conducted to improve our understanding of the role that the resistance of the slab to bending and its coupling to the ambient mantle play in subduction dynamics over geological time scales. Our models are set up with a viscous plate of silicone (lithosphere) subducting under negative buoyancy in a viscous layer of glucose syrup (mantle). For our study, the lithosphere/upper mantle viscosity contrast has been systematically varied, from similar to 10 to similar to 10(5) in order to explore the parameter space between weak and strong slab dynamics. We found that subduction is characterized by a retreating mode for viscosity ratios > 10(4), by the coexistence of a retreating mode and an advancing mode for viscosity ratios between similar to 10(4) and similar to 10(2), and quasi-stationary, Rayleigh-Taylor like behaviour for ratios < 10(2). By combining our experimental results and kinematic data from current subduction zones in four reference frames which differ in the amount of net rotation, we infer that a lithosphere/upper mantle viscosity contrast of 150-500 is necessary to obtain realistic trench/subducting plate velocity ratios as well as the variability of subduction styles observed in nature.
ISSN: 0012-821X

hal-00412008
https://hal.archives-ouvertes.fr/hal-00412008
DOI : 10.1016/j.epsl.2008.04.006

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
1] Recent statistical analysis by Lallemand et al. (2008) of subduction zone parameters revealed that the back-arc deformation mode depends on the combination between the subducting (nu(sub)) and upper (nu(up)) plate velocities. No significant strain is recorded in the arc area if plate kinematics verifies nu(up) = 0.5 vsub - 2.3 (cm/a) in the HS3 reference frame. Arc spreading ( shortening) occurs if nu(up) is greater ( lower) than the preceding relationship. We test this statistical law with numerical models of subduction, by applying constant plate velocities far away from the subduction zone. The subducting lithosphere is free to deform at all depths. We quantify the force applied on the two converging plates to sustain constant surface velocities. The simulated rheology combined viscous (non-Newtonian) and brittle behaviors, and depends on water content. The influence of subduction rate vs is first studied for a fixed upper plate. After 950 km of convergence ( steady state slab pull), the transition from extensional to compressive stresses in the upper plate occurs for vs similar to 1.4 cm/a. The effect of upper plate velocity is then tested at constant subduction rate. Upper plate retreat ( advance) with respect to the trench increases extension ( compression) in the arc lithosphere and increases ( decreases) the subducting plate dip. Our modeling confirms the statistical kinematic relationship between vsub and nu(up) that describes the transition from extensional to compressive stresses in the arc lithosphere, even if the modeled law is shifted toward higher rates of upper plate retreat, using our set of physical parameters ( e. g., 100 km thick subducting oceanic plate) and short- term simulations. Our results make valid the choice of the HS3 reference frame for assessing plate velocity influence on arc tectonic regime. The subduction model suggests that friction along the interplate contact and the mantle Stokes reaction could be the two main forces competing against slab pull for upper mantle subductions. Besides, our simulations show that the arc deformation mode is strongly time dependent.
ISSN: 1525-2027

hal-00412328
https://hal.archives-ouvertes.fr/hal-00412328
https://hal.archives-ouvertes.fr/hal-00412328/document
https://hal.archives-ouvertes.fr/hal-00412328/file/arcayG32008.pdf
DOI : 10.1029/2007GC001875

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
This study presents microstructural evidence for deformation-controlled melt organization and for changes in olivine deformation associated with the presence of melt in an 80 m vertical section of the crust-mantle transition zone in the Oman ophiolite. This zone represents an 'end member' case for analyzing feedbacks between deformation and melt distribution in the upper mantle, since it experienced strong shear strains in presence of large melt fractions. It is characterized by a subhorizontal compositional layering at the mm to meter scale, from weakly impregnated dunites to olivine-rich gabbros, which parallels a pervasive foliation containing a strong stretching lineation. The parallelism between the compositional layering and the foliation, the diffuse limits of the layers, the alignment of elongated plagioclase-rich aggregates devoid of internal deformation structures with the elongation of olivine crystals in the dunitic layers, and the sharp compositional changes across some, but not all layer limits suggest deformation plays an essential role on the development of the layering. The variation on a mm-scale of the olivine crystal preferred orientation (CPO) symmetry as a function of the modal content: from axial-[100] symmetry in layers with <70% modal olivine to axial-[010] in more gabbroic levels (<40% olivine), which is repeated over the entire section, implies deformation in presence of variable melt fractions. Axial-[100] olivine CPO in olivine-rich layers is consistent with deformation by dislocation creep under high temperature, low pressure, dry conditions. Axial-[010] olivine CPO patterns imply additional sliding along preferentially wetted (010) grain boundaries, increase in the activity of [001] glide, or transpression localized in the melt-rich layers. Since the change in CPO symmetry is not accompanied by dispersion, instantaneous melt fractions must have remained <30-40% in all layers. The continuous variation in olivine CPO symmetry with decreasing olivine content implies therefore that the former depends on the cumulated strain in presence of melt rather than on the instantaneous melt fraction.
ISSN: 0012-821X

hal-00795508
https://hal.archives-ouvertes.fr/hal-00795508
DOI : 10.1016/j.epsl.2012.10.003

Éditeur(s) : HAL CCSD
Résumé : International audience
As subducting plates reach the base of the upper mantle, some appear to flatten and stagnate, while others seemingly go through unimpeded. This variable resistance to slab sinking has been proposed to affect long-term thermal and chemical mantle circulation. A review of observational constraints and dynamic models highlights that neither the increase in viscosity between upper and lower mantle (likely by a factor 20–50) nor the coincident endothermic phase transition in the main mantle silicates (with a likely Clapeyron slope of –1 to –2 MPa/K) suffice to stagnate slabs. However, together the two provide enough resistance to temporarily stagnate subducting plates, if they subduct accompanied by significant trench retreat. Older, stronger plates are more capable of inducing trench retreat, explaining why backarc spreading and flat slabs tend to be associated with old-plate subduction. Slab viscosities that are ∼2 orders of magnitude higher than background mantle (effective yield stresses of 100–300 MPa) lead to similar styles of deformation as those revealed by seismic tomography and slab earthquakes. None of the current transition-zone slabs seem to have stagnated there more than 60 m.y. Since modeled slab destabilization takes more than 100 m.y., lower-mantle entry is apparently usually triggered (e.g., by changes in plate buoyancy). Many of the complex morphologies of lower-mantle slabs can be the result of sinking and subsequent deformation of originally stagnated slabs, which can retain flat morphologies in the top of the lower mantle, fold as they sink deeper, and eventually form bulky shapes in the deep mantle.
ISSN: 1553-040X

hal-01571965
https://hal.archives-ouvertes.fr/hal-01571965
DOI : 10.1130/GES01476.1

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
Numerical models show that small-scale convection (SSC) occurring atop a mantle plume is a plausible mechanism to rejuvenate the lithosphere. The triggering of SSC depends on the density contrast and on the rheology of the unstable layer underlying the stagnant upper part of the thermal boundary layer (TBL). Partial melting may change both properties. We analyze, using 2-D numerical simulations, how partial melting influences the dynamics of time-dependent SSC instabilities and the resulting thermo-mechanical rejuvenation of an oceanic plate moving atop of a plume. Our simulations show a complex behavior, with acceleration, no change, or delay of the SSC onset, due to competing effects of the latent heat of partial melting, which cools the plume material, and of the buoyancy increase associated with both melt retention and depletion of residue following melt extraction. The melt-induced viscosity reduction is too localized to affect significantly SSC dynamics. Faster SSC triggering is promoted for low melting degrees (low plume temperature anomalies, thick lithosphere, or fast moving plates), which limit both the temperature reduction due to latent heat of melting and the accumulation of depleted buoyant residue in the upper part of the unstable layer. In contrast, high partial melting degrees lead to a strong temperate decrease due to latent heat of melting and development of a thick depleted layer within the sublithospheric convecting layer, which delay the development of gravitational instabilities. Despite differences in SSC dynamics, the thinning of the lithosphere is not significantly enhanced relatively to simulations that neglect partial melting.
ISSN: 1525-2027

hal-01277175
https://hal.archives-ouvertes.fr/hal-01277175
DOI : 10.1002/2015GC005967

É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 AGU and the Geochemical Society
Résumé : International audience
We image the lithospheric and upper asthenospheric structure of western continental Yemen with 24 broadband stations to evaluate the role of the Afar plume on the evolution of the continental margin and its extent eastwards along the Gulf of Aden. We use teleseismic tomography to compute relative P-wave velocity variations in southwestern Yemen down to 300 km depth. Published receiver function analysis suggest a dramatic and localized thinning of the crust in the vicinity of the Red Sea and the Gulf of Aden, consistent with the velocity structure that we retrieve in our model. The mantle part of the model is dominated by the presence of a low velocity anomaly in which we infer partial melting just below thick Oligocene flood basalts and recent off-axis volcanic events (from 15 Ma to present). This low velocity anomaly could correspond to an abnormally hot mantle and could be responsible for dynamic topography and recent magmatism in western Yemen. Our new P-wave velocity model beneath western Yemen suggests the young rift flank volcanoes beneath margins and on the flanks of the Red Sea rift are caused by focused small-scale diapiric upwelling from a broad region of hot mantle beneath the area. Our work shows that relatively hot mantle, along with partial melting of the mantle, can persist beneath rifted margins after breakup has occurred.
ISSN: 1525-2027

hal-01015496
https://hal.archives-ouvertes.fr/hal-01015496
https://hal.archives-ouvertes.fr/hal-01015496/document
https://hal.archives-ouvertes.fr/hal-01015496/file/Korostelev_etal_GGG_2014.pdf
DOI : 10.1002/2014GC005316

Éditeur(s) : HAL CCSD
Résumé : 38 pages
Rapport
The Chicxulub impact crater, México, is unique. It is the only known terrestrial impact structure that has been directly linked to a mass extinction event and the only terrestrial impact with a global ejecta layer. Of the three largest impact structures on Earth, Chicxulub is the best preserved. Chicxulub is also the only known terrestrial impact structure with an intact, unequivocal topographic peak ring. Chicxulub’s role in the Cretaceous/Paleogene (K-Pg) mass extinction and its exceptional state of preservation make it an important natural laboratory for the study of both large impact crater formation on Earth and other planets and the effects of large impacts on the Earth’s environment and ecology. Our understanding of the impact process is far from complete, and despite more than 30 years of intense debate, we are still striving to answer the question as to why this impact was so catastrophic.During International Ocean Discovery Program (IODP) Expedition 364, Paleogene sediments and lithologies that make up the Chicxulub peak ring were cored to investigate (1) the nature and formational mechanism of peak rings, (2) how rocks are weakened during large impacts, (3) the nature and extent of post-impact hydrothermal circulation, (4) the deep biosphere and habitability of the peak ring, and (5) the recovery of life in a sterile zone. Other key targets included sampling the transition through a rare midlatitude section that might include Eocene and Paleocene hyperthermals and/or the Paleocene/Eocene Thermal Maximum (PETM); the composition and character of the impact breccias, melt rocks, and peak-ring rocks; the sedimentology and stratigraphy of the Paleocene–Eocene Chicxulub impact basin infill; the chronology of the peak-ring rocks; and any observations from the core that may help us constrain the volume of dust and climatically active gases released into the stratosphere by this impact. Petrophysical property measurements on the core and wireline logs acquired during Expedition 364 will be used to calibrate geophysical models, including seismic reflection and potential field data, and the integration of all the data will calibrate impact crater models for crater formation and environmental effects. The proposed drilling directly contributes to IODP Science Plan goals:Climate and Ocean Change: How resilient is the ocean to chemical perturbations? The Chicxulub impact represents an external forcing event that caused a 75% level mass extinction. The impact basin may also record key hyperthermals within the Paleogene.Biosphere Frontiers: What are the origin, composition, and global significance of subseafloor communities? What are the limits of life in the subseafloor? How sensitive are ecosystems and biodiversity to environmental change? Impact craters can create habitats for subsurface life, and Chicxulub may provide information on potential habitats for life, including extremophiles, on the early Earth and other planetary bodies. Paleontological and geochemical studies at ground zero will document how large impacts affect ecosystems and effects on biodiversity.Earth Connections/Earth in Motion: What are the composition, structure and dynamics of Earth’s upper mantle? What mechanisms control the occurrence of destructive earthquakes, landslides, and tsunami? Mantle uplift in response to impacts provides insight into dynamics that differ between Earth and other rocky planets. Impacts generate earthquakes, landslides, and tsunami, and scales that generally exceed plate tectonic processes yield insight into effects, the geologic record, and potential hazards.IODP Expedition 364 was a Mission Specific Platform expedition to obtain subseabed samples and downhole logging measurements from the sedimentary cover sequence and peak ring of the Chicxulub impact crater. A single borehole was drilled into the Chicxulub impact crater on the Yucatán continental shelf, recovering core from 505.7 to 1334.73 m below seafloor with ~99% core recovery and acquiring downhole logs for the entire depth.
International Ocean Discovery Program

hal-01639520
https://hal.archives-ouvertes.fr/hal-01639520
DOI : 10.14379/iodp.pr.364.2017

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : We study the Africa-Iberia plate boundary in the vicinity of Gibraltar. Numerous models have been proposed for that region throughout the last decades, proposing mechanisms that range widely from continental delamination, convective removal, to subduction of oceanic lithosphere. To better constrain upper-mantle structure under the region, we study waveforms of P-waves that traverse the Alboran Sea region between Spain and Morocco. These show dispersive behavior, which, together with early arrival times, confirms the presence of an anomalous upper mantle structure under the Alboran Sea. The dispersion is consistent with that expected from subducted lithosphere. Waveforms of body waves therefore provide a way to better constrain the elusive mantle structure and dynamics of the Alboran Sea region
ISSN: 0094-8276

hal-00407670
https://hal.archives-ouvertes.fr/hal-00407670
DOI : 10.1029/2007GL030964

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The present study aims at better deciphering the different mechanisms involved in the functioning of the subduction interplate. A 2D thermo-mechanical model is used to simulate a subduction channel, made of oceanic crust, free to evolve. Convergence at constant rate is imposed under a 100 km thick upper plate. Pseudo-brittle and non-Newtonian behaviours are modelled. The influence of the subduction channel strength, parameterized by the difference in activation energy between crust and mantle () is investigated to examine in detail the variations in depth of the subduction plane down-dip extent, . First, simulations show that numerical resolution may be responsible for an artificial and significant shallowing of if the weak crustal layer is not correctly resolved. Second, if the age of the subducting plate is 100 Myr, subduction occurs for any . The stiffer the crust is, that is, the lower is, the shallower is (60 km depth if kJ/mol) and the hotter the fore-arc base is. Conversely, imposing a very weak subduction channel ( J/mol) leads there to an extreme mantle wedge cooling and inhibits mantle melting in wet conditions. Partial kinematic coupling at the fore-arc base occurs if kJ/mol. If the incoming plate is 20 Myr old, subduction can occur under the conditions that the crust is either stiff and denser than the mantle, or weak and buoyant. In the latter condition, cold crust plumes rise from the subduction channel and ascend through the upper lithosphere, triggering (1) partial kinematic coupling under the fore-arc, (2) fore-arc lithosphere cooling, and (3) partial or complete hindrance of wet mantle melting. then ranges from 50 to more than 250 km depth and is time-dependent if crust plumes form. Finally, subduction plane dynamics is intimately linked to the regime of subduction-induced corner flow. Two different intervals of are underlined: 80–120 kJ/mol to reproduce the range of slab surface temperature inferred from geothermometry, and 10–40 kJ/mol to reproduce the shallow hot mantle wedge core inferred from conditions of last equilibration of near-primary arc magmas and seismic tomographies. Therefore, an extra process controlling mantle wedge dynamics is needed to satisfy simultaneously the aforementioned observations. A mantle viscosity reduction, by a factor 4–20, caused by metasomatism in the mantle wedge is proposed. From these results, I conclude that the subduction channel down-dip extent, , should depend on the subduction setting, to be consistent with the observed variability of sub-arc depths of the subducting plate surface.
ISSN: 0031-9201

hal-01623155
https://hal.archives-ouvertes.fr/hal-01623155
DOI : 10.1016/j.pepi.2017.05.008

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
[1] Statistical analysis of modern oceanic subduction zone parameters, such as the age of a downgoing plate or the absolute plate motions, is performed in order to investigate which parameter controls the dip of a slab and, conversely, what the influence of slab geometry is on upper plate behavior. For that purpose, parameters have been determined from global databases along 159 transects from all subduction zones that are not perturbed by nearby collision or ridge/plateau/seamount subduction. On the basis of tomographic images, slabs that penetrate through, or lie on, the 670 km discontinuity are also identified. The results of the statistical analysis are as follows: (1) Back-arc stress correlates with slab dip, i.e., back-arc spreading is observed for deep dips (deeper than 125 km) larger than 50°, whereas back-arc shortening occurs only for deep dips less than 30°. (2) Slab dip correlates with absolute motion of the overriding plate. The correlation is even better when the slab lies on, or even more penetrates through, the 670 km discontinuity. (3) Slabs dip more steeply, by about 20° on average, beneath oceanic overriding plates than beneath continental ones. (4) Slabs dip more steeply on average by about 10° near edges. (5) Slab dip does not correlate with the magnitude of slab pull, the age of subducting lithosphere at the trench, the thermal regime of the subducting lithosphere, the convergence rate, or the subduction polarity (east versus west). The present study provides evidence that the upper plate absolute motion plays an important role on slab dip, as well as on upper plate strain. Retreating overriding plates are often oceanic ones and thus may partially explain the steeper slab dips beneath oceanic upper plates. One can infer that low slab dips correlate well with compression in continental advancing upper plates, whereas steep dips are often associated with extension in oceanic retreating upper plates. Excess weight of old slabs is often counterbalanced by other forces, probably asthenospheric in origin, such as lateral mantle flow near slab edges or anchor forces, to determine slab dip. Components: 12,676 words, 13 figures, 1 table.
ISSN: 1525-2027

hal-01261567
https://hal.archives-ouvertes.fr/hal-01261567
https://hal.archives-ouvertes.fr/hal-01261567/document
https://hal.archives-ouvertes.fr/hal-01261567/file/lallemandG32005.pdf
DOI : 10.1029/2005GC000917

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We discuss mantle structure and dynamics under the Gibraltar arc. For that region, a large variety of geodynamic models has been proposed, and there was, up to recently, no consensus on which of these models is the best one to describe the geodynamic situation of the region. The key to distinguishing the different models probably lies in the upper mantle, and we propose two specially adapted seismological techniques that help to distinguish each of the two models that are more frequently invoked. These are, on one hand, the subduction-rollback (SR) model, and on the other, continental delamination or convective removal (CDCR). The first technique that we use is based on observations of dispersed P-wavetrains that are expected for certain ray directions. Observations made on the Gibraltar arc suggest the presence of subducted oceanic lithosphere under the Alboran Sea. The second technique is based on shear-wave splitting of SKS phases, and it allows us to look for the radially-oriented flow that must be present in the mantle beneath the Alboran Sea, for a CDCR mechanism. There are however no indications of any radial flow up to this point. Both approaches therefore suggest that subduction-rollback is the best model to explain the complex geodynamic situation under the region.
ISSN: 0040-1951

hal-01388041
http://hal.univ-reunion.fr/hal-01388041
http://hal.univ-reunion.fr/hal-01388041/document
http://hal.univ-reunion.fr/hal-01388041/file/Bokelmann_Gibraltar_Tphy2011_hal.pdf
DOI : 10.1016/j.tecto.2010.08.004

Éditeur(s) : HAL CCSD
Résumé : International audience
In the frame of the Africa-Europe convergence, the Mediterranean tectonic system presents a complex interaction between subduction rollback and upper-plate deformation during the Tertiary. The western Mediterranean is characterized by the exhumation of the largest subcontinental mantle massif worldwide (the Ronda Peridotite) and a narrow arcuate geometryacross the Gibraltar arc within the Betic-Rif belt (the internal part being called the Alboran domain), where the relationship between slab dynamics and surface tectonics is not well understood. New structural and geochronological data are used to argue for 1/ hyperstrechting of the continental lithosphere allowing extensional mantle exhumation to shallow depths, followed by 2/ lower miocene thrusting. Two Lower Miocene E-W-trending strike-slip corridors played a major role in the deformation pattern of the Alboran Domain, in which E-W dextral strike-slip faults, N60°-trending thrusts and N140°-trending normal faults developed simultaneously during dextral strike-slip simple shear. The inferred continuous westward translation of the Alboran Domain is accommodated by a major E-W-trending lateral ramp (strike-slip) and a N60°-trending frontal thrust. At lithosphere-scale, we interpret the observed deformation pattern as the upper-plate expression of a lateral slab tear and of its westward propagation since Lower Miocene. The crustal emplacement of the Ronda Peridotites occurred at the onset of this westward motion.The Miocene tectonics of the western Alboran is therefore marked by the inversion of a continental rift, triggered by shortening of the upper continental plate and accommodated by E-W dextral strike-slip corridors. During thrusting and westward displacement of the Alboran domain with respect to Iberia, the hot upper plate, which involved the previously exhumed sub-continental mantle, underwent fast cooling.
AGU Fall Meeting 2015
San Francisco, United States

insu-01242043
https://hal-insu.archives-ouvertes.fr/insu-01242043

Éditeur(s) : HAL CCSD
Résumé : International audience
Understanding subduction rheology in both space and time has been a challenge since the advent of plate tectonics. We herein focus on "subduction infancy", that is the first ~1-5 My immediately following subduction nucleation, when a newly born slab penetrates into the upper plate mantle and heats up. The only remnants of this critical yet elusive geodynamic step are thin metamorphic soles, commonly found beneath pristine, 100-1000 km long portions of oceanic lithosphere emplaced on continents (i.e., ophiolites).Through the (i) worldwide compilation of pressure-temperature conditions of metamorphic sole formation augmented by pseudosection thermodynamic modeling, (ii) calculations of the viscosity of materials along the plate interface and (iii) generic numerical thermal models, we provide a conceptual model of dynamic plate interface processes during subduction infancy (and initiation s.l.).We show in particular how major rheological switches across the subduction interface control slab penetration, and the formation of metamorphic soles. Due to the downward progression of hydration and weakening of the mantle wedge with cooling, the lower plate (basalt, sediment) and the upper plate (mantle wedge) rheologies equalize and switch over a restricted temperature-time-depth interval (e.g., at ~800°C and ~1 GPa, during 0.1-2 My, for high-temperature metamorphic sole formation). These switches result in episodes of maximum interplate mechanical coupling, thereby slicing the top of the slab and welding pieces (basalt, sediment) to the base of the mantle wedge. Similar mechanical processes likely apply for the later, deeper accretion and exhumation of high-temperature oceanic eclogites in serpentinite mélanges, or for the accretion of larger tectonic slices.This model provides constraints on the effective rheologies of the crust and mantle and general understanding, at both rock and plate scale, for accretion processes and early slab dynamics.
AGU Fall Meeting 2015
San Francisco, United States

insu-01242103
https://hal-insu.archives-ouvertes.fr/insu-01242103

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
In the framework of the Africa-Europe convergence, the Mediterranean system presents a complex interaction between subduction rollback and upper-plate deformation during the Tertiary. The western end of the system shows a narrow arcuate geometry across the Gibraltar arc, the Betic-Rif belt, in which the relationship between slab dynamics and surface tectonics is not well understood. The present study focuses on the Western Betics, which is characterized by two major thrusts: 1/ the Internal/External Zone Boundary limits the metamorphic domain (Alboran Domain) from the fold-and-thrust belts in the External Zone; 2/ the Ronda Peridotites Thrust allows the juxtaposition of a strongly attenuated lithosphere section with large bodies of sub-continental mantle rocks on top of upper crustal rocks. New structural data show that two major E-W strike-slip corridors played a major role in the deformation pattern of the Alboran Domain, in which E-W dextral strike-slip faults, N60 thrusts and N140 normal faults developed simultaneously during dextral strike-slip simple shear. Olistostromic sediments of Lower Miocene age were deposited and deformed in this tectonic context and hence provide an age estimate for the inferred continuous westward translation of the Alboran Domain that is accommodated by an E-W lateral (strike-slip) ramp and a N60 frontal thrust. The crustal emplacement of large bodies of sub-continental mantle may occur at the onset of this westward thrusting in the Western Alboran domain. At lithosphere-scale, we interpret the observed deformation pattern as the subduction upper-plate expression of a lateral slab tear and its westward propagation since Lower Miocene.
ISSN: 0040-1951

insu-01163695
https://hal-insu.archives-ouvertes.fr/insu-01163695
https://hal-insu.archives-ouvertes.fr/insu-01163695/document
https://hal-insu.archives-ouvertes.fr/insu-01163695/file/Frasca-2015.pdf
DOI : 10.1016/j.tecto.2015.05.017

Éditeur(s) : HAL CCSD
Résumé : Joint inversions are now commonly used in Earth Sciences. They have been developed to better understand Earth structure since they provide more constraints on the inverted parameters. We propose a new process to simultaneously invert several data sets in order to better image 3-D crustal and upper mantle structures. Our inversion uses 3 kind of data which present a good complementarity : (1) P receiver functions to provide Moho depth variations, (2) teleseismic delay times of P waves to retrieve velocity anomalies in the crust and upper mantle, (3) gravity anomalies to image density variations at lithospheric scale. In our scheme, receiver functions are first inverted. The resulting Moho depths are then interpolated and incorporated as a priori information into the joint inversion of teleseismic delay times and gravity anomalies process. In our new approach, we perform a model space search for Moho variations, P-velocity and density structure to find acceptable fit to the three data sets. In order to preferentially sample the good data fit region, we chose the Neighborhood algorithm of Sambridge to optimistically survey the model space. We model the delay times with a 3-D raytracing using velocity nodes evenly spaced and linked with density nodes via a linear relationship. This relationship can evolve through depth, simulating its pressure dependance. We test our scheme on synthetic examples, and apply the method in a continental rifting region where deformation processes are complex and badly known. This joint inversion thus provides more constraints on the inverted parameters, and allows us to better image crust - mantle interactions in the case of the Hangai dome (Central Mongolia). Indeed, in this area, we expect an asthenospheric upwelling to dynamically support the topography, but we image a density - velocity anomaly in the lower crust. The effects on topography and the origin (composition vs. thermal) of this anomaly are still under debate.
Les inversions conjointes sont maintenant communément utilisées en Sciences de la Terre. Elles ont été développées afin d'améliorer notre connaissance et compréhension de la structure interne de la Terre puisqu'elles apportent de plus en plus de contraintes sur les paramètres inversés. Dans cette thèse, nous proposons un nouveau processus d'inversion conjointe qui prend en compte trois paramètres différents et qui mène à l'obtention d'un modèle lithosphérique tridimensionnel. Notre méthode utilise trois types de données gravimétriques et sismologiques qui présentent une bonne complémentarité : (1) Les fonctions récepteur P afin d'obtenir les variations de profondeur du Moho, (2) les délais de temps d'arrivées P des téléséismes pour retrouver les anomalies de vitesse dans la croûte et le manteau supérieur, (3) les anomalies gravimétriques qui donnent accès aux variations de densité à l'échelle lithosphérique. Dans le schéma d'inversion proposé, nous inversons d'abord les fonctions récepteur. Les variations de profondeur du Moho qui en résultent sont alors interpolées puis incorporées comme information a priori dans le processus d'inversion conjointe gravimétrie - tomographie télésismique. Que ce soit pour les fonctions récepteur ou pour la partie conjointe, notre méthode d'inversion se base sur un algorithme de recherche afin de trouver la structure qui minimise l'écart aux trois types de données. Nous avons opté pour l'algorithme de voisinage qui présente l'avantage de concentrer sa recherche dans les régions qui minimisent l'écart aux données sans toutefois abandonner les autres régions afin de limiter les risques de convergence vers un minimum local. Nous utilisons un modèle 3-D constitué de noeuds de vitesse afin de modéliser les délais de temps. Ces noeuds de vitesse sont liés aux noeuds de densité grâce à une relation linéaire entre vitesse et densité. Cette relation peut évoluer avec la profondeur de manière à simuler les effets dus aux variations de pression. Nous avons réalisé plusieurs tests synthétiques afin d'évaluer le comportement de notre nouvelle méthode. Ils ont permis de montrer que les variations du Moho ainsi que les anomalies de vitesse et densité sont bien retrouvées. Ils ont également montré l'importance de la prise en compte de ces trois types de données dans un même schéma d'inversion, notamment en ce qui concerne la faculté de discerner les variations d'interfaces des anomalies volumiques. Nous avons ensuite appliqué la méthode dans une région continentale en phase de rifting où les processus de déformation sont complexes et mal connus. Notre nouvelle méthode a apporté davantage de contraintes sur les paramètres inversés et a permis d'obtenir un modèle correctement résolu de la croûte et du manteau supérieur sous le dôme de Hangai (Mongolie Centrale). Une anomalie asthénosphérique sous le dôme est souvent invoquée pour expliquer sa topographie élevée mais c'est une anomalie (en vitesse et densité) crustale (croûte inférieure) qui s'exprime dans nos modèles. Ses effets sur la topographie et son origine (compositionnelle et/ou thermique) sont encore débattus et nécessiteraient des études plus poussées (paramétriques, VP =VS, ...).
https://tel.archives-ouvertes.fr/tel-00783524

tel-00783524
https://tel.archives-ouvertes.fr/tel-00783524
https://tel.archives-ouvertes.fr/tel-00783524/document
https://tel.archives-ouvertes.fr/tel-00783524/file/Basuyau2011.pdf

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
We performed forward modelling of seismic anisotropy beneath a migrating strike-slip plate boundary to: (1) test if such a geodynamic context might explain teleseismic shear wave splitting data in the vicinity of the central part of the San Andreas Fault and (2) constrain the power of such data to unravel vertical and lateral variations in deformation patterns in the upper mantle. The modelling involves five steps: (1) thermomechanical modelling, using a finite-element code, of the deformation field, (2) viscoplastic self-consistent modelling of the resulting olivine and pyroxene crystal preferred orientations, (3) calculation of the elasticity tensors for different domains of the finite elements (FEs) model, (4) forward modelling of seismic wave propagation through the model using ray theory, finite-frequency theory and a full-wave approach and (5) performing splitting measurements on the synthetic seismograms. SKS splitting data in central California are best fitted by a model with a hotter geotherm within 60 km of the plate boundary accounting for the opening of an asthenospheric window due to the northward migration of the Mendocino Triple Junction. The westward motion of the plate boundary cannot however explain the rotation of fast polarization directions east of the San Andreas Fault in central California. Comparison between modelled and measured individual shear wave splitting also implies that the homogeneity of the two-layer models accounting for the observations in the vicinity of the San Andreas Fault indicates a sharp transition between lithospheric and asthenospheric deformations beneath this plate boundary. The ability of different synthetic approaches to localize horizontally and vertically the plate boundary-related deformation differs significantly. Splitting in ray theory synthetics closely follow variations in olivine crystallographic preferred orientations in the model. In contrast, splitting analysis on full-wave synthetics, which should be more representative of actual long-period SKS waves, results in smooth apparent lateral variations of anisotropy; the exact location and width of the plate boundary may only be retrieved by comparing fast-polarization profiles obtained using a multichannel analysis on waves with different periods.
ISSN: 0956-540X

hal-01236150
http://hal.univ-reunion.fr/hal-01236150
http://hal.univ-reunion.fr/hal-01236150/document
http://hal.univ-reunion.fr/hal-01236150/file/Bonnin_San_Andreas_modelling_GJI2012.pdf
DOI : 10.1111/j.1365-246X.2012.05650.x

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
During ODP Leg 209, a magma-starved area of the Mid-Atlantic Ridge (MAR) was drilled in the vicinity of the Fifteen-Twenty Fracture Zone (FZ) that offsets one of the slowest portions of the spreading ridge. We present here the results of a bulk rock multi-elemental study of 27 peridotites drilled at Sites 1272 and 1274 (to the south and the north of the FZ, respectively). The peridotites comprise mainly of harzburgites with minor dunites. Clinopyroxene (Cpx), which is interstitial and interpreted as secondary, is observed in Site 1274 peridotites. Sites 1272 and 1274 peridotites have low Al2O3 contents (<= 1 anhydrous wt.%), high Mg# (>91.5), and bulk rock trace elements compositions mostly below 0.1 x primitive mantle (PM). These peridotites, and in particular Site 1272 peridotites, represent the most depleted peridotites yet sampled at a slow spreading ridge. Their compositions indicate high degrees of partial melting and melt extraction. A single open-system melting event (melting plus percolation of melts produced within upwelling mantle) can explain their highly depleted yet linear chondrite-normalized REE patterns, characterized by a steady depletion from HREE to LREE. Late melt-rock reactions and precipitation of Cpx explains the slightly less depleted compositions of Site 1274 peridotites. Hence, the differences in composition between Sites 1272 and 1274 peridotites do not provide evidence for regional variations in the degrees of partial melting from the south to the north of the FZ. The occurrence of highly refractory peridotites in the Fifteen-Twenty area suggests we sampled a more actively convecting mantle than generally supposed below slow spreading centers.
ISSN: 0012-821X

hal-00412001
https://hal.archives-ouvertes.fr/hal-00412001
DOI : 10.1016/j.epsl.2007.11.058

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
The nature and extent of the regional lithosphere-asthenosphere interaction beneath Ireland and Britain remains unclear. Although it has been established that ancient Caledonian signatures pervade the lithosphere, tertiary structure related to the Iceland plume has been inferred to dominate the asthenosphere. To address this apparent contradiction in the literature, we image the 3-D lithospheric and deeper upper-mantle structure beneath Ireland via non-linear, iterative joint teleseismic-gravity inversion using data from the ISLE (Irish Seismic Lithospheric Experiment), ISUME (Irish Seismic Upper Mantle Experiment) and GRACE (Gravity Recovery and Climate Experiment) experiments. The inversion combines teleseismic relative arrival time residuals with the GRACE long wavelength satellite derived gravity anomaly by assuming a depth-dependent quasilinear velocity-density relationship. We argue that anomalies imaged at lithospheric depths probably reflect compositional contrasts, either due to terrane accretion associated with Iapetus Ocean closure, frozen decompressional melt that was generated by plate stretching during the opening of the north Atlantic Ocean, frozen Iceland plume related magmatic intrusions, or a combination thereof. The continuation of the anomalous structure across the lithosphere-asthenosphere boundary is interpreted as possibly reflecting sub-lithospheric small-scale convection initiated by the lithospheric compositional contrasts. Our hypothesis thus reconciles the disparity which exists between lithospheric and asthenospheric structure beneath this region of the north Atlantic rifted margin.
ISSN: 0956-540X

hal-00617815
https://hal.archives-ouvertes.fr/hal-00617815
DOI : 10.1111/j.1365-246X.2011.04921.x

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Imaging deep active volcanic areas remains a challenge in our understanding of their activity and evolution, especially in subduction zones. Study of magnetic anomalies is appropriate to access such dynamics in depth. The magnetic anomaly pattern of the Lesser Antilles Arc (LAA) subduction is studied through Curie Point Depth (CPD), interpreted as the depth of the 580 °C isotherm, and developed to better assess the deep thermal structure of the arc. The depth of the estimated CPD exhibits a complex topography. Keeping in mind the overall uncertainty associated with this method, a main doming is evidenced below the Guadeloupe archipelago. Its apex is shifted towards the ancient arc, suggesting a very hot state of the fore-arc/arc domain. To better understand the LAA thermal state, we perform 2D thermo-mechanical simulations of the subduction zone. Recalling that magnetite is a serpentinization by-product, we simulate water transfer triggered by slab dehydration to test the assumption of fore-arc serpentinization suggested by the positive magnetic anomaly in the vicinity of the Guadeloupe archipelago. In this area, the subduction-induced arc lithosphere hydration and related weakening trigger a fast heating of the upper plate by basal convective removal. This process of fast arc lithosphere thinning may apply where simultaneously the volcanic arc is split in two and normal convergence is high enough. As serpentinization strongly decreases P-wave velocity, we propose a new interpretation of a published seismic profile below Guadeloupe. The seismic layer previously interpreted as the arc lower crust may rather be a layer of serpentinized mantle, as supported by spatial correlations between gravimetric and magnetic anomalies. Consequently, at the scale of Guadeloupe Island, the fore-arc Moho would be shallower than initially assumed, with a dome shape more consistent with both the extensive deformation active since the Oligocene in the inner fore-arc and the CPD doming.
ISSN: 0040-1951

Droits : http://creativecommons.org/licenses/by-nc-nd/
insu-01505300
https://hal-insu.archives-ouvertes.fr/insu-01505300
https://hal-insu.archives-ouvertes.fr/insu-01505300/document
https://hal-insu.archives-ouvertes.fr/insu-01505300/file/Gailler-2017.pdf
DOI : 10.1016/j.tecto.2017.03.014