É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 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 Oxford University Press (OUP)
Résumé : International audience
Two-dimensional simulations using a thermomechanical model based on a finite-difference method on a staggered grid and a marker in cell method are performed to study the plume-lithosphere interaction beneath moving plates. The plate and the convective mantle are modelled as a homogeneous peridotite with a Newtonian temperature- and pressure-dependent viscosity. A constant velocity, ranging from 5 to 12.5 cm yr−1, is imposed at the top of the plate. Plumes are generated by imposing a thermal anomaly of 150 to 350 K on a 50 km wide domain at the base of the model (700 km depth); the plate atop this thermal anomaly is 40 Myr old. We analyse (1) the kinematics of the plume as it impacts the moving plate, (2) the dynamics of time-dependent small-scale convection (SSC) instabilities developing in the low-viscosity layer formed by spreading of hot plume material at the base of the lithosphere and (3) the resulting thermal rejuvenation of the lithosphere. The spreading of the plume material at the base of the lithosphere, characterized by the ratio between the maximum down- and upstream horizontal (dimensionless) velocities in the plume-fed sublithospheric layer, Peup/Pedown depends on the ratio between the maximum plume upwelling velocity and the plate velocity, Peplume/Peplate. For fast plate velocities and sluggish plumes (low Peplume/Peplate), plate motion drags most plume material and downstream flow is dominant. As Peplume/Peplate increases, an increasing part of the plume material flows upstream. SSC systematically develops in the plume-fed sublithospheric layer, downstream from the plume. Onset time of SSC decreases with the Rayleigh number. For vigorous plumes, it does not depend on plate velocity. For more sluggish plumes, however, variations in the plume spreading behaviour at the base of the lithosphere result in a decrease in the onset time of SSCs with increasing plate velocity. In any case, SSC results in uplift of the isotherm 1573 K by up to 20 km relative to its initial equilibrium depth at the impact point.
ISSN: 0956-540X

hal-00856763
https://hal.archives-ouvertes.fr/hal-00856763
DOI : 10.1093/gji/ggt128

Éditeur(s) : HAL CCSD European Geosciences Union
Résumé : The properties of the subduction interplate domain are likely to affect not only the seismogenic potential of the subduction area but also the overall subduction process, as it influences its viability. Numerical simulations are performed to model the long-term equilibrium state of the subduction interplate when the diving lithosphere interacts with both the overriding plate and the surrounding convective mantle. The thermomechanical model combines a non-Newtonian viscous rheology and a pseudo-brittle rheology. Rock strength here depends on depth, temperature and stress, for both oceanic crust and mantle rocks. I study the evolution through time of, on one hand, the brittle-ductile transition (BDT) depth, zBDT, and, on the other hand, of the kinematic decoupling depth, zdec, simulated along the subduction interplate. The results show that both a high friction and a low ductile strength at the asthenospheric wedge tip shallow zBDT. The influence of the weak material activation energy is of second order but not negligible. zBDT becomes dependent on the ductile strength increase with depth (activation volume) if the BDT occurs at the interplate decoupling depth. Regarding the interplate decoupling depth, it is shallowed (1) significantly if mantle viscosity at asthenospheric wedge tip is low, (2) if the difference in mantle and interplate activation energy is weak, and (3) if the activation volume is increased. Very low friction coefficients and/or low asthenospheric viscosities promote zBDT = zdec. I then present how the subducting lithosphere age affects the brittle-ductile transition depth and the kinematic decoupling depth in this model. Simulations show that a rheological model in which the respective activation energies of mantle and interplate material are too close hinders the mechanical decoupling at the down-dip extent of the interplate, and eventually jams the subduction process during incipient subduction of a young (20-Myr-old) and soft lithosphere under a thick upper plate. Finally, both the BDT depth and the decoupling depth are a function of the subducting plate age, but are not influenced in the same fashion: cool and old subducting plates deepen the BDT but shallow the interplate decoupling depth. Even if BDT and kinematic decoupling are intrinsically related to different mechanisms of deformation, this work shows that they are able to interact closely. Comparison between modelling results and observations suggests a minimum friction coefficient of 0.045 for the interplate plane, even 0.069 in some cases, to model realistic BDT depths. The modelled zdec is a bit deeper than suggested by geophysical observations. Eventually, the better way to improve the adjustment to observations may rely on a moderate to strong asthenosphere viscosity reduction in the metasomatised mantle wedge.
ISSN: 1869-9510

hal-00795559
https://hal.archives-ouvertes.fr/hal-00795559
DOI : 10.5194/se-3-467-2012, 2012

Éditeur(s) : HAL CCSD CNRS éditions
Résumé : International audience
Die antike Siedlung von Mandeure und Mathay (Epomanduodurum), welche durch ihre Größe und ihre monumentale Ausstattung als die zweitwichtigste Stadt der Sequaner hinter der Civitashauptstadt Besançon betrachtet wird, ist Gegenstand eines interdisziplinarischen Forschungsprogramms. Dieser Artikel stellt eine Bilanz der Ergebnisse der zwischen 2001 und 2004 durchgeführten Forschungen dar. Im Mittelpunkt der Forschungen der ersten Phase stand der zentralen Teil der Siedlung mit dem Bereich Heiligtum-Theater-spätrömische Befestigung, der für die Problematik der Entstehung, der Entwicklung und dem Ende der römischen Besiedlung besonders interessant ist. Verschiedene, sich ergänzende Projekte wurden eingeleitet: Erforschung des alten Dokumentarbestandes, geophysikalische Prospektionen, geoarchäologische Sondierungen, Plangrabungen, architektonische Untersuchunge. Durch diese Projekte war es möglich, die Kenntnisse über die drei wichtigen Großbauten (Heiligtum-Theater-spätrömische Befestigung) grundsätzlich zu vertiefen und insbesondere die geophysikalischen Prospektionsergebnisse lassen in bemerkenswerter Weise einen dichten Baubestand dieses Stadtteils erkennen. Schließlich konnte die Funktion und die Entwicklung der Flußebene des Doubs in diesem Bereich besser verstanden werden als auch die Art und Weise in welcher die antike Monumentalanlage in diesen einzigartigen natürlichen Rahmen einbezogen wurde.
Since 2001, a collective pluridisciplinary research programme is being undertaken on the antique town of Epomanduodurum in the communes of Mandeure and Mathay. This site is considered like the second one of the Sequani, for its dimensions and the monuments, after the civitas capital, Besançon. This article is an assessment of the investigations conducted between 2001 and 2004. During this period, the central part of the site was studied (sanctuary, theatre, Late fortification area), considered as the most significant for the question of the birth, development and decline of the town. Various complementary approaches have been carried out: study of archives, geophysical and geoarchaeological survey, extent excavations, architectural analyses. They lead to a better understanding of the three main buildings (the Celtic and Gallo-Roman sanctuary, the theatre of the Principate, the Late Antique fortification) and they add considerable informations about the pattern of settlement of this area; so it is possible to detect more precisely the spatial organization and finally to understand how the alluvial plain of the Doubs river was functioning and how the monumental complex has been integrated in such a peculiar landscape.
L'agglomération antique d'Epomanduodurum, située sur les communes de Mandeure et Mathay, et considérée comme la seconde du pays séquane par la taille et la parure monumentale, derrière la capitale de cité, Besançon, fait l'objet d'un programme collectif de recherche, pluridisciplinaire, depuis 2001. Cet article dresse un bilan des résultats des fouilles réalisées entre 2001 et 2004. Dans cette première phase, les investigations ont porté principalement sur la partie centrale du site (secteur sanctuaire-théâtrefortification tardive), perçue comme particulièrement intéressante pour la problématique de la naissance, du développement et du déclin de l'agglomération. Différentes approches complémentaires ont été lancées : étude du fond documentaire ancien, prospections géophysiques, sondages géoarchéologiques, fouilles extensives, études architecturales. Elles ont permis d'approfondir la connaissance des trois principaux édifices (sanctuaire laténien et gallo-romain, théâtre du Haut-Empire, fortification du Bas-Empire), de compléter de façon très notable la trame de l'occupation de ce secteur, et donc d'en mieux percevoir l'organisation spatiale, enfin de mieux comprendre le fonctionnement de la plaine alluviale du Doubs et l'insertion dans ce contexte naturel particulier d'un complexe monumental antique.
ISSN: 0016-4119

halshs-00718356
https://halshs.archives-ouvertes.fr/halshs-00718356

É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

É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