@misc{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 = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {7}, issn = {1866-8372}, doi = {10.25932/publishup-56262}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-562628}, pages = {29}, 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} } @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} } @phdthesis{Rehak2008, author = {Rehak, Katrin}, title = {Pliocene-Pleistocene landscape evolution in south-central Chile : interactions between tectonic, geomorphic, and climatic processes}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-19793}, school = {Universit{\"a}t Potsdam}, year = {2008}, abstract = {Landscapes evolve in a complex interplay between climate and tectonics. Thus, the geomorphic characteristics of a landscape can only be understood if both, climatic and tectonic signals of past and ongoing processes can be identified. In order to evaluate the impact of both forcing factors it is crucial to quantify the evolution of geomorphic markers in natural environments. The Cenozoic Andes are an ideal setting to evaluate tectonic and climatic aspects of landscape evolution at different time and length scales in different natural compartments. The Andean Cordillera constitutes the type subduction orogen and is associated with the subduction of the oceanic Nazca Plate beneath the South American continent since at least 200 million years. In Chile and the adjacent regions this convergent margin is characterized by active tectonics, volcanism, and mountain building. Importantly, along the coast of Chile megathrust earthquakes occur frequently and influence landscape evolution. In fact, the largest earthquake ever recorded occurred in south-central Chile in 1960 and comprised a rupture zone of ~ 1000 km length. However, on longer time scales beyond historic documentation of seismicity it is not well known, how such seismotectonic segments have behaved and how they influence the geomorphic evolution of the coastal realms. With several semi-independent morphotectonic segments, recurrent megathrust earthquakes, and a plethora of geomorphic features indicating sustained tectonism, the margin of Chile is thus a key area to study relationships between surface processes and tectonics. In this study, I combined geomorphology, geochronology, sedimentology, and morphometry to quantify the Pliocene-Pleistocene landscape evolution of the tectonically active south-central Chile forearc. Thereby, I provide (1) new results about the influence of seismotectonic forearc segmentation on the geomorphic evolution and (2) new insights in the interaction between climate and tectonics with respect to the morphology of the Chilean forearc region. In particular, I show that the forearc is characterized by three long-term segments that are not correlated with short-lived earthquake-rupture zones that may. These segments are the Nahuelbuta, Tolt{\´e}n, and Bueno segments, each recording a distinct geomorphic and tectonic evolution. The Nahuelbuta and Bueno segments are undergoing active tectonic uplift. The long-term behavior of these two segments is manifested in form of two doubly plunging, growing antiforms that constitute an integral part of the Coastal Cordillera and record the uplift of marine and river terraces. In addition, these uplifting areas have caused major changes in flow directions or rivers. In contrast, the Tolt{\´e}n segment, situated between the two other segments, appears to be quasi-stable. In order to further quantify uplift and incision in the actively deforming Nahuelbuta segment, I dated an erosion surface and fluvial terraces in the Coastal Cordillera with cosmogenic 10Be and 26Al and optically stimulated luminescence, respectively. According to my results, late Pleistocene uplift rates corresponding to 0.88 mm a-1 are faster than surface-uplift rates averaging over the last 5 Ma, which are in the range of 0.21 mm a-1. This discrepancy suggests that surface uplift is highly variable in time and space and might preferably concentrate along reverse faults as indicated by a late Pleistocene flow reversal. In addition, the results of exposure dating with cosmogenic 10Be and 26Al indicate that the morphotectonic segmentation of this region of the forearc has been established in Pliocene time, coeval with the initiation of uplift of the Coastal Cordillera about 5 Ma ago, inferred to be related to a shift in subduction mode from erosion to accretion. Finally, I dated volcanic clasts obtained from alluvial surfaces in the Central Depression, a low-relief sector separating the Coastal from the Main Cordillera, with stable cosmogenic 3He and 21Ne, in order to reveal the controls of sediment accumulation in the forearc. My results document that these gently sloping surfaces have been deposited 150 to 300 ka ago. This deposition may be related to changes in the erosional regime during glacial episodes. Taken together, the data indicates that the overall geomorphic expression of the forearc is of post-Miocene age and may be intimately related to a climatic overprint of the tectonic system. This climatic forcing is also reflected in the topography and local relief of the Central and Southern Andes that vary considerably along the margin, determined by the dominant surface process that in turn is eventually controlled by climate. However, relief also partly reflects surface processes that have taken place under past climatic conditions. This emphasizes that due care has to be exercised when interpreting landscapes as mirrors of modern climates.}, language = {en} } @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} } @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} } @phdthesis{Pons2023, author = {Pons, Micha{\"e}l}, title = {The Nature of the tectonic shortening in Central Andes}, doi = {10.25932/publishup-60089}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-600892}, school = {Universit{\"a}t Potsdam}, pages = {160}, year = {2023}, abstract = {The Andean Cordillera is a mountain range located at the western South American margin and is part of the Eastern- Circum-Pacific orogenic Belt. The ~7000 km long mountain range is one of the longest on Earth and hosts the second largest orogenic plateau in the world, the Altiplano-Puna plateau. The Andes are known as a non-collisional subduction-type orogen which developed as a result of the interaction between the subducted oceanic Nazca plate and the South American continental plate. The different Andean segments exhibit along-strike variations of morphotectonic provinces characterized by different elevations, volcanic activity, deformation styles, crustal thickness, shortening magnitude and oceanic plate geometry. Most of the present-day elevation can be explained by crustal shortening in the last ~50 Ma, with the shortening magnitude decreasing from ~300 km in the central (15°S-30°S) segment to less than half that in the southern part (30°S-40°S). Several factors were proposed that might control the magnitude and acceleration of shortening of the Central Andes in the last 15 Ma. One important factor is likely the slab geometry. At 27-33°S, the slab dips horizontally at ~100 km depth due to the subduction of the buoyant Juan Fernandez Ridge, forming the Pampean flat-slab. This horizontal subduction is thought to influence the thermo-mechanical state of the Sierras Pampeanas foreland, for instance, by strengthening the lithosphere and promoting the thick-skinned propagation of deformation to the east, resulting in the uplift of the Sierras Pampeanas basement blocks. The flat-slab has migrated southwards from the Altiplano latitude at ~30 Ma to its present-day position and the processes and consequences associated to its passage on the contemporaneous acceleration of the shortening rate in Central Andes remain unclear. Although the passage of the flat-slab could offer an explanation to the acceleration of the shortening, the timing does not explain the two pulses of shortening at about 15 Ma and 4 Ma that are suggested from geological observations. I hypothesize that deformation in the Central Andes is controlled by a complex interaction between the subduction dynamics of the Nazca plate and the dynamic strengthening and weakening of the South American plate due to several upper plate processes. To test this hypothesis, a detailed investigation into the role of the flat-slab, the structural inheritance of the continental plate, and the subduction dynamics in the Andes is needed. Therefore, I have built two classes of numerical thermo-mechanical models: (i) The first class of models are a series of generic E-W-oriented high-resolution 2D subduction models thatinclude flat subduction in order to investigate the role of the subduction dynamics on the temporal variability of the shortening rate in the Central Andes at Altiplano latitudes (~21°S). The shortening rate from the models was then validated with the observed tectonic shortening rate in the Central Andes. (ii) The second class of models are a series of 3D data-driven models of the present-day Pampean flat-slab configuration and the Sierras Pampeanas (26-42°S). The models aim to investigate the relative contribution of the present-day flat subduction and inherited structures in the continental lithosphere on the strain localization. Both model classes were built using the advanced finite element geodynamic code ASPECT. The first main finding of this work is to suggest that the temporal variability of shortening in the Central Andes is primarily controlled by the subduction dynamics of the Nazca plate while it penetrates into the mantle transition zone. These dynamics depends on the westward velocity of the South American plate that provides the main crustal shortening force to the Andes and forces the trench to retreat. When the subducting plate reaches the lower mantle, it buckles on it-self until the forced trench retreat causes the slab to steepen in the upper mantle in contrast with the classical slab-anchoring model. The steepening of the slab hinders the trench causing it to resist the advancing South American plate, resulting in the pulsatile shortening. This buckling and steepening subduction regime could have been initiated because of the overall decrease in the westwards velocity of the South American plate. In addition, the passage of the flat-slab is required to promote the shortening of the continental plate because flat subduction scrapes the mantle lithosphere, thus weakening the continental plate. This process contributes to the efficient shortening when the trench is hindered, followed by mantle lithosphere delamination at ~20 Ma. Finally, the underthrusting of the Brazilian cratonic shield beneath the orogen occurs at ~11 Ma due to the mechanical weakening of the thick sediments covered the shield margin, and due to the decreasing resistance of the weakened lithosphere of the orogen. The second main finding of this work is to suggest that the cold flat-slab strengthens the overriding continental lithosphere and prevents strain localization. Therefore, the deformation is transmitted to the eastern front of the flat-slab segment by the shear stress operating at the subduction interface, thus the flat-slab acts like an indenter that "bulldozes" the mantle-keel of the continental lithosphere. The offset in the propagation of deformation to the east between the flat and steeper slab segments in the south causes the formation of a transpressive dextral shear zone. Here, inherited faults of past tectonic events are reactivated and further localize the deformation in an en-echelon strike-slip shear zone, through a mechanism that I refer to as "flat-slab conveyor". Specifically, the shallowing of the flat-slab causes the lateral deformation, which explains the timing of multiple geological events preceding the arrival of the flat-slab at 33°S. These include the onset of the compression and of the transition between thin to thick-skinned deformation styles resulting from the crustal contraction of the crust in the Sierras Pampeanas some 10 and 6 Myr before the Juan Fernandez Ridge collision at that latitude, respectively.}, language = {en} } @misc{LipkeKruegerRoessler2008, author = {Lipke, Katrin and Kr{\"u}ger, Frank and R{\"o}ßler, Dirk}, title = {Subduction zone structure along Sumatra from receiver functions}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-18260}, year = {2008}, abstract = {Receiver functions are a good tool to investigate the seismotectonic structure beneath the a seismic station. In this study we apply the method to stations situated on or near Sumatra to find constraints on a more detailed velocity model which should improve earthquake localisation. We estimate shallow Moho-depths (~ 21 km) close to the trench and depths of ~30 km at greater distances. First evidences for the dip direction of the slab of ~60° are provided. Receiver functions were calculated for 20 stations for altogether 110 earthquakes in the distance range between 30° and 95° from the receiver. However the number of receiver functions per station is strongly variable as it depends on the installation date, the signal-to-noise-ratio of the station and the reliability of the acquisition.}, 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} } @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} }