@article{PourteauSchererSchornetal.2019, author = {Pourteau, Amaury and Scherer, Erik E. and Schorn, Simon and Bast, Rebecca and Schmidt, Alexander and Ebert, Lisa}, title = {Thermal evolution of an ancient subduction interface revealed by Lu-Hf garnet geochronology, Halilbagi Complex (Anatolia)}, series = {Geoscience Frontiers}, volume = {10}, journal = {Geoscience Frontiers}, number = {1}, publisher = {Elsevier}, address = {Amsterdam}, issn = {1674-9871}, doi = {10.1016/j.gsf.2018.03.004}, pages = {127 -- 148}, year = {2019}, abstract = {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.}, language = {en} } @article{DuesterhoeftQuinterosOberhaenslietal.2014, author = {Duesterhoeft, Erik and Quinteros, Javier and Oberh{\"a}nsli, Roland and Bousquet, Romain and de Capitani, Christian}, title = {Relative impact of mantle densification and eclogitization of slabs on subduction dynamics: A numerical thermodynamic/thermokinematic investigation of metamorphic density evolution}, series = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth}, volume = {637}, journal = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0040-1951}, doi = {10.1016/j.tecto.2014.09.009}, pages = {20 -- 29}, year = {2014}, abstract = {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.}, language = {en} } @misc{PourteauOberhaensliCandanetal.2016, author = {Pourteau, Amaury and Oberh{\"a}nsli, Roland and Candan, Osman and Barrier, Eric and Vrielynck, Bruno}, title = {Neotethyan closure history of western Anatolia: a geodynamic discussion}, series = {International journal of earth sciences}, volume = {105}, journal = {International journal of earth sciences}, publisher = {Springer}, address = {New York}, issn = {1437-3254}, doi = {10.1007/s00531-015-1226-7}, pages = {203 -- 224}, year = {2016}, abstract = {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.}, language = {en} } @article{RodriguezPicedaScheckWenderothGomezDacaletal.2020, author = {Rodriguez Piceda, Constanza and Scheck Wenderoth, Magdalena and Gomez Dacal, Maria Laura and Bott, Judith and Prezzi, Claudia Beatriz and Strecker, Manfred}, title = {Lithospheric density structure of the southern Central Andes constrained by 3D data-integrative gravity modelling}, series = {International journal of earth sciences}, volume = {110}, journal = {International journal of earth sciences}, number = {7}, publisher = {Springer}, address = {New York}, issn = {1437-3254}, doi = {10.1007/s00531-020-01962-1}, pages = {2333 -- 2359}, year = {2020}, abstract = {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.}, language = {en} } @misc{KonradSchmolkeHalama2014, author = {Konrad-Schmolke, Matthias and Halama, Ralf}, title = {Combined thermodynamic-geochemical modeling in metamorphic geology: Boron as tracer of fluid-rock interaction}, series = {Lithos : an international journal of mineralogy, petrology, and geochemistry}, volume = {208}, journal = {Lithos : an international journal of mineralogy, petrology, and geochemistry}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0024-4937}, doi = {10.1016/j.lithos.2014.09.021}, pages = {393 -- 414}, year = {2014}, abstract = {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.}, language = {en} }