TY - JOUR A1 - Pourteau, Amaury A1 - Scherer, Erik E. A1 - Schorn, Simon A1 - Bast, Rebecca A1 - Schmidt, Alexander A1 - Ebert, Lisa T1 - Thermal evolution of an ancient subduction interface revealed by Lu-Hf garnet geochronology, Halilbagi Complex (Anatolia) JF - Geoscience Frontiers N2 - The thermal structure of subduction zones exerts a major influence on deep-seated mechanical and chemical processes controlling arc magmatism, seismicity, and global element cycles. Accretionary complexes exposed inland may comprise tectonic blocks with contrasting pressure-temperature (P-T) histories, making it possible to investigate the dynamics and thermal evolution of former subduction interfaces. With this aim, we present new Lu-Hf geochronological results for mafic rocks of the Halilbagi Complex (Anatolia) that evolved along different thermal gradients. Samples include a lawsonite-epidote blueschist, a lawsonite-epidote eclogite, and an epidote eclogite (all with counter-clockwise P-T paths), a prograde lawsonite blueschist with a "hairpin"-type P-T path, and a garnet amphibolite from the overlying sub-ophiolitic metamorphic sole. Equilibrium phase diagrams suggest that the garnet amphibolite formed at similar to 0.6-0.7 GPa and 800-850 degrees C, whereas the prograde lawsonite blueschist records burial from 2.1 GPa and 420 degrees C to 2.6 GPa and 520 degrees C. Well-defined Lu-Hf isochrons were obtained for the epidote eclogite (92.38 +/- 0.22 Ma) and the lawsonite-epidote blueschist (90.19 +/- 0.54 Ma), suggesting rapid garnet growth. The lawsonite-epidote eclogite (87.30 +/- 0.39 Ma) and the prograde lawsonite blueschist (ca. 86 Ma) are younger, whereas the garnet amphibolite (104.5 +/- 3.5 Ma) is older. Our data reveal a consistent trend of progressively decreasing geothermal gradient from granulite-facies conditions at similar to 104 Ma to the epidote-eclogite facies around 92 Ma, and the lawsonite blueschist-facies between 90 Ma and 86 Ma. Three Lu-Hf garnet dates (between 92 Ma and 87 Ma) weighted toward the growth of post-peak rims (as indicated by Lu distribution in garnet) suggest that the HP/LT rocks were exhumed continuously and not episodically. We infer that HP/LT metamorphic rocks within the Halilbagi Complex were subjected to continuous return flow, with "warm" rocks being exhumed during the tectonic burial of "cold" ones. Our results, combined with regional geological constraints, allow us to speculate that subduction started at a transform fault near a mid-oceanic spreading centre. Following its formation, this ancient subduction interface evolved thermally over more than 15 Myr, most likely as a result of heat dissipation rather than crustal underplating. (C) 2018, China University of Geosciences (Beijing) and Peking University. Production and hosting by Elsevier B.V. KW - Subduction KW - Lu/Hf dating of garnet KW - Metamorphic sole KW - Eclogite KW - Blueschist KW - Lawsonite Y1 - 2019 U6 - https://doi.org/10.1016/j.gsf.2018.03.004 SN - 1674-9871 VL - 10 IS - 1 SP - 127 EP - 148 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Duesterhoeft, Erik A1 - Quinteros, Javier A1 - Oberhänsli, Roland A1 - Bousquet, Romain A1 - de Capitani, Christian T1 - Relative impact of mantle densification and eclogitization of slabs on subduction dynamics: A numerical thermodynamic/thermokinematic investigation of metamorphic density evolution JF - Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth N2 - Understanding the relationships between density and spatio-thermal variations at convergent plate boundaries is important for deciphering the present-day dynamics and evolution of subduction zones. In particular, the interaction between densification due to mineralogical phase transitions and slab pull forces is subject to ongoing investigations. We have developed a two-dimensional subduction zone model that is based on thermodynamic equilibrium assemblage calculations and includes the effects of melting processes on the density distribution in the lithosphere. Our model calculates the "metamorphic density" of rocks as a function of pressure, temperature and chemical composition in a subduction zone down to 250 km. We have used this model to show how the hydration, dehydration, partial melting and fractionation processes of rocks all influence the metamorphic density and greatly depend on the temperature field within the subduction system. These processes are largely neglected by other approaches that reproduce the density distribution within this complex tectonic setting. Our model demonstrates that the initiation of edogitization (i.e., when crustal rocks reach higher densities than the ambient mantle) of the slab is not the only significant process that makes the descending slab denser and generates the slab pull force. Instead, the densification of the lithospheric mantle of the sinking slab starts earlier than eclogitization and contributes significantly to slab pull in the early stages of subduction. Accordingly, the complex metamorphic structure of the slab and the mantle wedge has an important impact on the development of subduction zones. (C) 2014 Elsevier B.V. All rights reserved. KW - Density KW - Melt KW - Metamorphism KW - Subduction KW - Thermodynamic modeling KW - Thermo-mechanical modeling Y1 - 2014 U6 - https://doi.org/10.1016/j.tecto.2014.09.009 SN - 0040-1951 SN - 1879-3266 VL - 637 SP - 20 EP - 29 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Pourteau, Amaury A1 - Oberhänsli, Roland A1 - Candan, Osman A1 - Barrier, Eric A1 - Vrielynck, Bruno T1 - Neotethyan closure history of western Anatolia: a geodynamic discussion JF - International journal of earth sciences N2 - This paper addresses the lithosphere-scale subduction-collision history of the eastern termination of the Aegean retreating subduction system, i.e. western Anatolia. Although there is some general consensus on the protracted subduction evolution of the Aegean since the early Cenozoic at least, correlation with western Anatolia has been widely debated for more than several decades. In western Anatolia, three main tectonic configurations have been envisaged in the past years to reconstruct slab dynamics during the closure of the Neotethyan oceanic realm since the Late Cretaceous. Some authors have suggested an Aegean-type scenario, with the continuous subduction of a single lithospheric slab, punctuated by episodic slab roll-back and trench retreat, whereas others assumed a discontinuous subduction history marked by intermittent slab break-off during either the Campanian (ca. 75 Ma) or the Early Eocene (ca. 55-50 Ma). The third view implies three partly contemporaneous subduction zones. Our review of these models points to key debated aspects that can be re-evaluated in the light of multidisciplinary constraints from the literature. Our discussion leads us to address the timing of subduction initiation, the existence of hypothetical ocean basins, the number of intervening subduction zones between the Taurides and the Pontides, the palaeogeographic origin of tectonic units and the possibility for slab break-off during either the Campanian or the Early Eocene. Thence, we put forward a favoured tectonic scenario featuring two successive phases of subduction of a single lithospheric slab and episodic accretion of two continental domains separated by a continental trough, representing the eastern end of the Cycladic Ocean of the Aegean. The lack of univocal evidence for slab break-off in western Anatolia and southward-younging HP/LT metamorphism in continental tectonic units (from similar to 85, 70 to 50 Ma) in the Late Cretaceous-Palaeogene period suggests continuous subduction since similar to 110 Ma, marked by roll-back episodes in the Palaeocene and the Oligo-Miocene, and slab tearing below western Anatolia during the Miocene. KW - Subduction KW - Anatolia KW - Aegean KW - Neotethys KW - Slab break-off KW - Slab roll-back Y1 - 2016 U6 - https://doi.org/10.1007/s00531-015-1226-7 SN - 1437-3254 SN - 1437-3262 VL - 105 SP - 203 EP - 224 PB - Springer CY - New York ER - TY - JOUR A1 - Rodriguez Piceda, Constanza A1 - Scheck Wenderoth, Magdalena A1 - Gomez Dacal, Maria Laura A1 - Bott, Judith A1 - Prezzi, Claudia Beatriz A1 - Strecker, Manfred T1 - Lithospheric density structure of the southern Central Andes constrained by 3D data-integrative gravity modelling JF - International journal of earth sciences N2 - The southern Central Andes (SCA) (between 27 degrees S and 40 degrees S) is bordered to the west by the convergent margin between the continental South American Plate and the oceanic Nazca Plate. The subduction angle along this margin is variable, as is the deformation of the upper plate. Between 33 degrees S and 35 degrees S, the subduction angle of the Nazca plate increases from sub-horizontal (< 5 degrees) in the north to relatively steep (similar to 30 degrees) in the south. The SCA contain inherited lithological and structural heterogeneities within the crust that have been reactivated and overprinted since the onset of subduction and associated Cenozoic deformation within the Andean orogen. The distribution of the deformation within the SCA has often been attributed to the variations in the subduction angle and the reactivation of these inherited heterogeneities. However, the possible influence that the thickness and composition of the continental crust have had on both short-term and long-term deformation of the SCA is yet to be thoroughly investigated. For our investigations, we have derived density distributions and thicknesses for various layers that make up the lithosphere and evaluated their relationships with tectonic events that occurred over the history of the Andean orogeny and, in particular, investigated the short- and long-term nature of the present-day deformation processes. We established a 3D model of lithosphere beneath the orogen and its foreland (29 degrees S-39 degrees S) that is consistent with currently available geological and geophysical data, including the gravity data. The modelled crustal configuration and density distribution reveal spatial relationships with different tectonic domains: the crystalline crust in the orogen (the magmatic arc and the main orogenic wedge) is thicker (similar to 55 km) and less dense (similar to 2900 kg/m(3)) than in the forearc (similar to 35 km, similar to 2975 kg/m(3)) and foreland (similar to 30 km, similar to 3000 kg/m(3)). Crustal thickening in the orogen probably occurred as a result of stacking of low-density domains, while density and thickness variations beneath the forearc and foreland most likely reflect differences in the tectonic evolution of each area following crustal accretion. No clear spatial relationship exists between the density distribution within the lithosphere and previously proposed boundaries of crustal terranes accreted during the early Paleozoic. Areas with ongoing deformation show a spatial correlation with those areas that have the highest topographic gradients and where there are abrupt changes in the average crustal-density contrast. This suggests that the short-term deformation within the interior of the Andean orogen and its foreland is fundamentally influenced by the crustal composition and the relative thickness of different crustal layers. A thicker, denser, and potentially stronger lithosphere beneath the northern part of the SCA foreland is interpreted to have favoured a strong coupling between the Nazca and South American plates, facilitating the development of a sub-horizontal slab. KW - Central andes KW - Lithospheric structure KW - Crustal density KW - Gravity KW - modelling KW - Subduction Y1 - 2020 U6 - https://doi.org/10.1007/s00531-020-01962-1 SN - 1437-3254 SN - 1437-3262 VL - 110 IS - 7 SP - 2333 EP - 2359 PB - Springer CY - New York ER - TY - JOUR A1 - Konrad-Schmolke, Matthias A1 - Halama, Ralf T1 - Combined thermodynamic-geochemical modeling in metamorphic geology: Boron as tracer of fluid-rock interaction JF - Lithos : an international journal of mineralogy, petrology, and geochemistry N2 - Quantitative geochemical modeling is today applied in a variety of geological environments from the petrogenesis of igneous rocks to radioactive waste disposal. In addition, the development of thermodynamic databases and computer programs to calculate equilibrium phase diagrams has greatly advanced our ability to model geodynamic processes. Combined with experimental data on elemental partitioning and isotopic fractionation, thermodynamic forward modeling unfolds enormous capacities that are far from exhausted. In metamorphic petrology the combination of thermodynamic and trace element forward modeling can be used to study and to quantify processes at spatial scales from mu m to km. The thermodynamic forward models utilize Gibbs energy minimization to quantify mineralogical changes along a reaction path of a chemically open fluid/rock system. These results are combined with mass balanced trace element calculations to determine the trace element distribution between rock and melt/fluid during the metamorphic evolution. Thus, effects of mineral reactions, fluid-rock interaction and element transport in metamorphic rocks on the trace element and isotopic composition of minerals, rocks and percolating fluids or melts can be predicted. Here we illustrate the capacities of combined thermodynamic-geochemical modeling based on two examples relevant to mass transfer during metamorphism. The first example focuses on fluid-rock interaction in and around a blueschist-facies shear zone in felsic gneisses, where fluid-induced mineral reactions and their effects on boron (B) concentrations and isotopic compositions in white mica are modeled. In the second example, fluid release from a subducted slab, the associated transport of B as well as variations in B concentrations and isotopic compositions in liberated fluids and residual rocks are modeled. We compare the modeled results of both examples to geochemical data of natural minerals and rocks and demonstrate that the combination of thermodynamic and geochemical models enables quantification of metamorphic processes and insights into element cycling that would have been unattainable if only one model approach was chosen. (C) 2014 Elsevier B.V. All rights reserved. KW - Thermodynamic-geochemical modeling KW - Fluid-rock interaction KW - Subduction KW - Dehydration KW - Boron isotopes Y1 - 2014 U6 - https://doi.org/10.1016/j.lithos.2014.09.021 SN - 0024-4937 SN - 1872-6143 VL - 208 SP - 393 EP - 414 PB - Elsevier CY - Amsterdam ER -