TY - JOUR A1 - Socquet, Anne A1 - Valdes, Jesus Pina A1 - Jara, Jorge A1 - Cotton, Fabrice Pierre A1 - Walpersdorf, Andrea A1 - Cotte, Nathalie A1 - von Specht, Sebastian A1 - Ortega-Culaciati, Francisco A1 - Carrizo, Daniel A1 - Norabuena, Edmundo T1 - An 8month slow slip event triggers progressive nucleation of the 2014 Chile megathrust JF - Geophysical research letters N2 - The mechanisms leading to large earthquakes are poorly understood and documented. Here we characterize the long-term precursory phase of the 1 April 2014 M(w)8.1 North Chile megathrust. We show that a group of coastal GPS stations accelerated westward 8months before the main shock, corresponding to a M(w)6.5 slow slip event on the subduction interface, 80% of which was aseismic. Concurrent interface foreshocks underwent a diminution of their radiation at high frequency, as shown by the temporal evolution of Fourier spectra and residuals with respect to ground motions predicted by recent subduction models. Such ground motions change suggests that in response to the slow sliding of the subduction interface, seismic ruptures are progressively becoming smoother and/or slower. The gradual propagation of seismic ruptures beyond seismic asperities into surrounding metastable areas could explain these observations and might be the precursory mechanism eventually leading to the main shock. KW - seismology KW - GPS KW - subduction KW - precursor Y1 - 2017 U6 - https://doi.org/10.1002/2017GL073023 SN - 0094-8276 SN - 1944-8007 VL - 44 SP - 4046 EP - 4053 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Jozi Najafabadi, Azam A1 - Haberland, Christian A1 - Le Breton, Eline A1 - Handy, Mark R. A1 - Verwater, Vincent F. A1 - Heit, Benjamin A1 - Weber, Michael T1 - Constraints on crustal structure in the vicinity of the adriatic indenter (European Alps) from Vp and Vp/Vs local earthquake tomography JF - Journal of geophysical research : Solid earth N2 - In this study, 3-D models of P-wave velocity (Vp) and P-wave and S-wave ratio (Vp/Vs) of the crust and upper mantle in the Eastern and eastern Southern Alps (northern Italy and southern Austria) were calculated using local earthquake tomography (LET). The data set includes high-quality arrival times from well-constrained hypocenters observed by the dense, temporary seismic networks of the AlpArray AASN and SWATH-D. The resolution of the LET was checked by synthetic tests and analysis of the model resolution matrix. The small inter-station spacing (average of similar to 15 km within the SWATH-D network) allowed us to image crustal structure at unprecedented resolution across a key part of the Alps. The derived P velocity model revealed a highly heterogeneous crustal structure in the target area. One of the main findings is that the lower crust is thickened, forming a bulge at 30-50 km depth just south of and beneath the Periadriatic Fault and the Tauern Window. This indicates that the lower crust decoupled both from its mantle substratum as well as from its upper crust. The Moho, taken to be the iso-velocity contour of Vp = 7.25 km/s, agrees with the Moho depth from previous studies in the European and Adriatic forelands. It is shallower on the Adriatic side than on the European side. This is interpreted to indicate that the European Plate is subducted beneath the Adriatic Plate in the Eastern and eastern Southern Alps. KW - European Alps KW - crustal structure KW - subduction KW - seismic tomography KW - body waves Y1 - 2022 U6 - https://doi.org/10.1029/2021JB023160 SN - 2169-9313 SN - 2169-9356 VL - 127 IS - 2 PB - American Geophysical Union CY - Washington ER - TY - THES A1 - Freisleben, Roland T1 - Deciphering the mechanisms of permanent forearc deformation based on marine terraces T1 - Untersuchung der Mechanismen dauerhafter Forearc-Deformation anhand von marinen Terrassen N2 - The Andes reflect Cenozoic deformation and uplift along the South American margin in the context of regional shortening associated with the interaction between the subducting Nazca plate and the overriding continental South American plate. Simultaneously, multiple levels of uplifted marine terraces constitute laterally continuous geomorphic features related to the accumulation of permanent forearc deformation in the coastal realm. However, the mechanisms responsible for permanent coastal uplift and the persistency of current/decadal deformation patterns over millennial timescales are still not fully understood. This dissertation presents a continental-scale database of last interglacial terrace elevations and uplift rates along the South American coast that provides the basis for an analysis of a variety of mechanisms that are possibly responsible for the accumulation of permanent coastal uplift. Regional-scale mapping and analysis of multiple, late Pleistocene terrace levels in central Chile furthermore provide valuable insights regarding the persistency of current seismic asperities, the role of upper-plate faulting, and the impact of bathymetric ridges on permanent forearc deformation. The database of last interglacial terrace elevations reveals an almost continuous signal of background-uplift rates along the South American coast at ~0.22 mm/yr that is modified by various short- to long-wavelength changes. Spatial correlations with crustal faults and subducted bathymetric ridges suggest long-term deformation to be affected by these features, while the latitudinal variability of climate forcing factors has a profound impact on the generation and preservation of marine terraces. Systematic wavelength analyses and comparisons of the terrace-uplift rate signal with different tectonic parameters reveal short-wavelength deformation to result from crustal faulting, while intermediate- to long-wavelength deformation might indicate various extents of long-term seismotectonic segments on the megathrust, which are at least partially controlled by the subduction of bathymetric anomalies. The observed signal of background-uplift rate is likely accumulated by moderate earthquakes near the Moho, suggesting multiple, spatiotemporally distinct phases of uplift that manifest as a continuous uplift signal over millennial timescales. Various levels of late Pleistocene marine terraces in the 2015 M8.3 Illapel-earthquake area reveal a range of uplift rates between 0.1 and 0.6 mm/yr and indicate decreasing uplift rates since ~400 ka. These glacial-cycle uplift rates do not correlate with current or decadal estimates of coastal deformation suggesting seismic asperities not to be persistent features on the megathrust that control the accumulation of permanent forearc deformation over long timescales of 105 years. Trench-parallel, crustal normal faults modulate the characteristics of permanent forearc-deformation; upper-plate extension likely represents a second-order phenomenon resulting from subduction erosion and subsequent underplating that lead to regional tectonic uplift and local gravitational collapse of the forearc. In addition, variable activity with respect to the subduction of the Juan Fernández Ridge can be detected in the upper plate over the course of multiple interglacial periods, emphasizing the role of bathymetric anomalies in causing local increases in terrace-uplift rate. This thesis therefore provides new insights into the current understanding of subduction-zone processes and the dynamics of coastal forearc deformation, whose different interacting forcing factors impact the topographic and geomorphic evolution of the western South American coast. N2 - Die känozoischen Anden resultieren aus der regionalen Verkürzung und Hebung der kontinentalen Kruste entlang des südamerikanischen Kontinentalrandes infolge der Interaktion zwischen der subduzierenden Nazca-Platte und der südamerikanischen Platte. Zahlreiche, durch tektonische Prozesse angehobene marine Abrasionsplattformen entlang der Küsten am westlichen Kontinentalrand Südamerikas bilden lateral kontinuierliche Stufen, welche die Akkumulation dauerhafter Deformation im küstennahen Forearc-Bereich widerspiegeln. Die Mechanismen, welche für die dauerhafte Hebung der Küste und die Beständigkeit derzeitiger/dekadischer Deformationsmuster auch über Jahrtausende hinweg verantwortlich sind, sind jedoch noch nicht vollständig geklärt. Vor diesem Hintergrund wird in dieser Dissertation zunächst eine kontinentale Datenbank der gehobenen marinen Terrassen des letzten Interglazials (125,000 Jahre), ihrer Höhenverteilung sowie der daraus abgeleiteten Hebungsraten für die Westküste von Südamerika vorgestellt. Diese Datenbank bildet die Grundlage für die Analyse einer Vielzahl von Mechanismen, welche möglicherweise für die Akkumulation dauerhafter Küstenhebung und Deformation verantwortlich sind. In einem weiteren Schritt wurden zusätzlich spätpleistozäne Terrassenniveaus in Zentralchile regionalmaßstäblich kartiert und hinsichtlich ihrer räumlich-zeitlichen Entwicklung analysiert. Diese Informationen liefern wertvolle Erkenntnisse über die Beständigkeit von seismischen Asperitäten, den Einfluss bathymetrischer Rücken auf die dauerhafte Deformation der Oberplatte sowie die Rolle von krustalen Störungen hinsichtlich der lokalen topographischen Differenzierung seismotektonischer Segmente. Die Datenbank der im letzten Interglazial entstandenen marinen Terrassen und der abgeleiteten Hebungsraten zeigt ein nahezu kontinuierliches Signal von Hintergrund-Hebungsraten entlang der südamerikanischen Küste von ~0,22 mm/Jahr, das durch verschiedene kurz- bis langwellige Modifikationen gekennzeichnet ist. Räumliche Übereinstimmungen mit krustalen Störungen und subduzierten bathymetrischen Rücken deuten darauf hin, dass die langfristige Deformation von diesen Parametern beeinflusst wird, während die breitengradabhängige Variabilität von Klimafaktoren tiefgreifende Auswirkungen auf die Entstehung und Erhaltung von marinen Terrassen hat. Systematische Wellenlängenanalysen und Vergleiche des Signals der Terrassenhebungsrate mit verschiedenen tektonischen Parametern zeigen, dass Deformation im kurzwelligen Bereich auf krustale Störungen zurückzuführen ist. Die Deformation im mittel- bis langwelligen Bereich könnte hingegen auf verschiedene Ausdehnungen langfristig aktiver seismotektonischer Segmente der Subduktionszone hindeuten, die zumindest teilweise durch die Subduktion bathymetrischer Anomalien kontrolliert werden. Das beobachtete Signal der Hintergrund-Hebungsrate wird wahrscheinlich durch Erdbeben mittlerer Magnituden in der Nähe der Moho generiert, was auf mehrere, räumlich und zeitlich getrennte Hebungsphasen hindeutet, die sich über Jahrtausende hinweg als kontinuierliches Hebungssignal manifestieren. Gehobene spätpleistozäne marine Terrassen im Illapel-Erdbebengebiet von 2015 weisen eine Hebungsrate zwischen 0,1 und 0,6 mm/Jahr auf und deuten insgesamt auf abnehmende tektonische Hebungsraten seit ~400,000 Jahren hin. Diese langfristigen Hebungsraten auf glazialen Zeitskalen korrelieren nicht mit den Hebungsraten auf kurzen Zeitskalen, was darauf schließen lässt, dass seismische Asperitäten im Bereich der Subduktionszone keine dauerhaften Merkmale sind, welche die Deformation und Hebung der Forearc-Regionen über lange Zeiträume von mehreren 105 Jahren steuern. Wie auch an anderen aktiven konvergenten Plattenrändern existieren in Zentralchile krustale Abschiebungen parallel zum Plattenrand, die das Erscheinungsbild dauerhafter Forearc-Deformation stark beeinflussen. Die Extensionsprozesse der Oberplatte sind dabei wahrscheinlich ein Phänomen zweiter Ordnung, welches durch Subduktionserosion und anschließende Akkretion hervorgerufen wird und sich in regionaler tektonischer Hebung sowie einem gravitativen Kollaps der Forearc-Region äußert. Darüber hinaus lassen sich über verschiedene Glazial- und Interglazialzyklen hinweg Änderungen in der Bedeutung des Juan-Fernández-Rückens hinsichtlich der Deformations-prozesse in der Oberplatte nachweisen; dies unterstreicht die transiente Rolle bathymetrischer Anomalien bei der Beeinflussung lokaler Änderungen der Terrassenhebungsrate. Die vorgelegte Arbeit liefert daher neue Einblicke für ein besseres Verständnis der Deformationsprozesse in Subduktionszonen und den küstennahen Forearc-Regionen sowie ihrer unterschiedlich interagierenden Antriebsfaktoren, welche die topographische und geomorphologische Entwicklung der westlichen südamerikanischen Küste beeinflussen. KW - South America KW - marine terrace KW - tectonic geomorphology KW - subduction KW - Südamerika KW - marine Terrassen KW - Subduktion KW - tektonische Geomorphologie Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-610359 ER - TY - JOUR A1 - Konrad-Schmolke, Matthias A1 - O'Brien, Patrick J. A1 - Zack, Thomas T1 - Fluid Migration above a Subducted Slab-Constraints on Amount, Pathways and Major Element Mobility from Partially Overprinted Eclogite-facies Rocks (Sesia Zone, Western Alps) JF - Journal of petrology N2 - The Western Alpine Sesia-Lanzo Zone (SLZ) is a sliver of eclogite-facies continental crust exhumed from mantle depths in the hanging wall of a subducted oceanic slab. Eclogite-facies felsic and basic rocks sampled across the internal SLZ show different degrees of retrograde metamorphic overprint associated with fluid influx. The weakly deformed samples preserve relict eclogite-facies mineral assemblages that show partial fluid-induced compositional re-equilibration along grain boundaries, brittle fractures and other fluid pathways. Multiple fluid influx stages are indicated by replacement of primary omphacite by phengite, albitic plagioclase and epidote as well as partial re-equilibration and/or overgrowths in phengite and sodic amphibole, producing characteristic step-like compositional zoning patterns. The observed textures, together with the map-scale distribution of the samples, suggest open-system, pervasive and reactive fluid flux across large rock volumes above the subducted slab. Thermodynamic modelling indicates a minimum amount of fluid of 0 center dot 1-0 center dot 5 wt % interacting with the wall-rocks. Phase relations and reaction textures indicate mobility of K, Ca, Fe and Mg, whereas Al is relatively immobile in these medium-temperature-high-pressure fluids. Furthermore, the thermodynamic models show that recycling of previously fractionated material, such as in the cores of garnet porphyroblasts, largely controls the compositional re-equilibration of the exhumed rock body. KW - fluid migration KW - subduction KW - fluid-rock interaction KW - Sesia Zone Y1 - 2011 U6 - https://doi.org/10.1093/petrology/egq087 SN - 0022-3530 VL - 52 IS - 3 SP - 457 EP - 486 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Saki, Adel A1 - Miri, Mirmohammad A1 - Oberhänsli, Roland T1 - High temperature - low pressure metamorphism during subduction of Neo-Tethys beneath the Iranian plate BT - evidence for mafic migmatite formation in the Alvand complex (western Iran) JF - Mineralogy and petrology N2 - Subduction of Neo-Tethys oceanic lithosphere beneath the Iranian plate during the Mesozoic formed several igneous bodies of ultramafic to intermediate and felsic composition. Intrusion of these magmas into a regional metamorphic sequence (the Sanandaj-Sirjan Zone) caused partial melting and formation of migmatites with meta-pelitic protoliths. The Alvand complex (west Iran) is a unique area comprising migmatites of both mafic and pelitic protoliths. In this area, the gabbroic rocks contain veins of leucosome at their contact with pyroxenite and olivine gabbro. These leucosomes are geochemically and mineralogically different from leucosomes of the meta-pelitic migmatites and clearly show properties of I-type granites. Microscopic observations and whole rock compositions of the mafic migmatite leucosomes show that migmatization occurred through partial melting of biotite, hornblende and plagioclase. Thermobarometric calculations indicate 800 degrees C and 3.7 kbar for partial melting, although phase diagram modeling demonstrates that the presence of water could decrease the solidus temperature by about 40 degrees C. Our results suggest an asthenospheric magma upwelling as the source of heat for partial melting of the gabbroic rock during subduction of Neo-Tethys oceanic crust under the western edge of the Iranian plate. The present study also reveals relationships between migmatization and formation of S- and I -type granites in the area. KW - partial melting KW - mafic migmatite KW - pelitic migmatite KW - Neo-Tethys KW - subduction KW - Gabbroic rocks Y1 - 2020 U6 - https://doi.org/10.1007/s00710-020-00721-z SN - 0930-0708 SN - 1438-1168 VL - 114 IS - 6 SP - 539 EP - 557 PB - Springer CY - Wien [u.a.] ER - TY - THES A1 - Muldashev, Iskander T1 - Modeling of the great earthquake seismic cycles T1 - Modellierung der seismischen Zyklen von Mega-Erdbeben N2 - The timing and location of the two largest earthquakes of the 21st century (Sumatra, 2004 and Tohoku 2011, events) greatly surprised the scientific community, indicating that the deformation processes that precede and follow great megathrust earthquakes remain enigmatic. During these phases before and after the earthquake a combination of multi-scale complex processes are acting simultaneously: Stresses built up by long-term tectonic motions are modified by sudden jerky deformations during earthquakes, before being restored by multiple ensuing relaxation processes. This thesis details a cross-scale thermomechanical model developed with the aim of simulating the entire subduction process from earthquake (1 minute) to million years’ time scale, excluding only rupture propagation. The model employs elasticity, non-linear transient viscous rheology, and rate-and-state friction. It generates spontaneous earthquake sequences, and, by using an adaptive time-step algorithm, recreates the deformation process as observed naturally over single and multiple seismic cycles. The model is thoroughly tested by comparing results to those from known high- resolution solutions of generic modeling setups widely used in modeling of rupture propagation. It is demonstrated, that while not modeling rupture propagation explicitly, the modeling procedure correctly recognizes the appearance of instability (earthquake) and correctly simulates the cumulative slip at a fault during great earthquake by means of a quasi-dynamic approximation. A set of 2D models is used to study the effects of non-linear transient rheology on the postseismic processes following great earthquakes. Our models predict that the viscosity in the mantle wedge drops by 3 to 4 orders of magnitude during a great earthquake with magnitude above 9. This drop in viscosity results in spatial scales and timings of the relaxation processes following the earthquakes that are significantly different to previous estimates. These models replicate centuries long seismic cycles exhibited by the greatest earthquakes (like the Great Chile 1960 Earthquake) and are consistent with the major features of postseismic surface displacements recorded after the Great Tohoku Earthquake. The 2D models are also applied to study key factors controlling maximum magnitudes of earthquakes in subduction zones. Even though methods of instrumentally observing earthquakes at subduction zones have rapidly improved in recent decades, the characteristic recurrence interval of giant earthquakes (Mw>8.5) is much larger than the currently available observational record and therefore the necessary conditions for giant earthquakes are not clear. Statistical studies have recognized the importance of the slab shape and its surface roughness, state of the strain of the upper plate and thickness of sediments filling the trenches. In this thesis we attempt to explain these observations and to identify key controlling parameters. We test a set of 2D models representing great earthquake seismic cycles at known subduction zones with various known geometries, megathrust friction coefficients, and convergence rates implemented. We found that low-angle subduction (large effect) and thick sediments in the subduction channel (smaller effect) are the fundamental necessary conditions for generating giant earthquakes, while the change of subduction velocity from 10 to 3.5 cm/yr has a lower effect. Modeling results also suggest that having thick sediments in the subduction channel causes low static friction, resulting in neutral or slightly compressive deformation in the overriding plate for low-angle subduction zones. These modeling results agree well with observations for the largest earthquakes. The model predicts the largest possible earthquakes for subduction zones of given dipping angles. The predicted maximum magnitudes exactly threshold magnitudes of all known giant earthquakes of 20th and 21st centuries. The clear limitation of most of the models developed in the thesis is their 2D nature. Development of 3D models with comparable resolution and complexity will require significant advances in numerical techniques. Nevertheless, we conducted a series of low-resolution 3D models to study the interaction between two large asperities at a subduction interface separated by an aseismic gap of varying width. The novelty of the model is that it considers behavior of the asperities during multiple seismic cycles. As expected, models show that an aseismic gap with a narrow width could not prevent rupture propagation from one asperity to another, and that rupture always crosses the entire model. When the gap becomes too wide, asperities do not interact anymore and rupture independently. However, an interesting mode of interaction was observed in the model with an intermediate width of the aseismic gap: In this model the asperities began to stably rupture in anti-phase following multiple seismic cycles. These 3D modeling results, while insightful, must be considered preliminary because of the limitations in resolution. The technique developed in this thesis for cross-scale modeling of seismic cycles can be used to study the effects of multiple seismic cycles on the long-term deformation of the upper plate. The technique can be also extended to the case of continental transform faults and for the advanced 3D modeling of specific subduction zones. This will require further development of numerical techniques and adaptation of the existing advanced highly scalable parallel codes like LAMEM and ASPECT. N2 - Zeitpunkt und Ort der zwei größten Erdbeben des 21. Jahrhunderts (Sumatra 2004 und Tohoku 2011) überraschten die wissenschaftliche Gemeinschaft, da sie darauf hindeuten, dass die einem Megathrust-Erdbeben vorangehenden und folgenden Deformationsprozesse weiterhin rästelhaft bleiben. Ein Problem ist die komplexe Art der Subduktionsprozesse auf unterschiedlichen Skalen. Spannungen werden durch langzeitliche, tektonische Bewegungen aufgebaut, von plötzlicher, ruckartiger Deformation während Erdbeben modifiziert und anschließend durch verschiedene Entspannungsprozesse wiederhergestellt. In dieser Arbeit wird ein skalen-übergreifendes thermomechanisches Modell entwickelt mit dem Ziel den vollständigen Subduktionsprozess von kleiner Skala (Minuten, z.B. Erdbeben) bis zu Millionen Jahren zu simulieren. Dabei bleibt nur das dynamische Prozess der unberücksichtigt. Das Modell nutzt Elastizität, nicht-lineare transient viskose Rheologie und „rate-and-state“ Reibungsgesetz. Es erzeugt spontane Erdbeben-Sequenzen und stellt durch das Nutzen eines adaptiven Zeitschritt Algorithmus den Deformationsprozess nach, wie er in der Natur während einzelner und mehrerer seismischer Zyklen beobachtet wird. Anhand der 2D-Modell Reihe werden die Effekte der nicht-linearen transient (viskosen) Rheologie auf postseismische Prozesse nach großen Erdbeben untersucht. Die Modelle sagen eine Verringerung der Mantelkeil-Viskosität um 3 bis 4 Größenordnungen während der großen Erdbeben (Mw > 9) vorher. Im Gegensatz zur momentanen Meinung, ergeben sich andere räumliche und zeitliche Verteilungen des Entspannungsprozesses nach großen Erdbeben. Jahrhunderte lange seismische Zyklen, wie zum Beispiel das große Erdbeben in Chile 1960, werden reproduziert und stimmen mit den Hauptmerkmalen der postseismischen Oberflächenverschiebungen überein, die nach dem großen Tohoku Erdbeben aufgezeichnet wurden. Mithilfe der 2D-Modelle werden außerdem Schlüsselfaktoren untersucht, die die maximale Magnitude von Erdbeben bestimmen. Obwohl die instrumentellen Methoden zur Beobachtung von Erdbeben in Subduktionszonen in den letzten Jahrzehnten stetig verbessert wurden, kann das Auftreten der größten Erdbeben (Mw > 8.5) und die dafür vorrauszusetzenden Bedingungen nicht genau bestimmt werden, da die charackeristische Zeit für ein Wiederauftreten deutlich größer als der Beobachtungszeitraum ist. Wir versuchen in dieser Arbeit diese Beobachtungen zu erklären und die kontrollierenden Schlüsselfaktoren zu bestimmen. Eine 2D-Modell Reihe mit großen, seismischen Erdbebenzyklen untersucht unterschiedliche Geometrien von Subduktionszonen sowie Reibungskoeffizieten an der Subduktions Platten-Grenzoberfläche und Konvergenzraten. Subduktion mit geringem Winkel (großer Effekt) und mächtige Sedimente im Subduktionskanal (kleiner Effekt) sind grundlegend notwendige Bedingungen für riesige Erdbeben, während eine Änderung des Reibungsparameters nur geringe Auswirkungen hat. Die Modellierungsergebnisse sind in viiÜbereinstimmung mit den Beobachtungen von den größten Erdbeben. Die maximale vohergesagte Magnituden liegen an der oberen Grenze alle historische Beben der 20. und 21. Jahrhunderden. Die größte Einschränkung der entwickelten Modelle ist ihre 2-dimensionale Natur. Um die Interaktion zwischen zwei großen Unebenheiten („asperities“) auf der Subduktions Platten-Grenzoberfläche, die von einer aseismischen Lücke („gap“) voneinander getrennt sind, zu untersuchen, wurde eine 3D-Modell Reihe mit geringer Auflösung durchgeführt. Neu an diesem Modell ist, dass das Verhalten der Unebenheiten während mehrerer seismischer Zyklen berücksichtigt wird. Wie erwartet zeigt das Modell, dass eine schmale, aseismische Lücke ein Propagieren von Brüchen von einer Unebenheit zur nächsten nicht verhindern kann und diese Brüche das gesamte Modell durchkreuzen. Wird die Lücke zu groß interagieren die Unebenheiten nicht mehr miteinander und brechen unabhängig voneinander. Allerdings wurde in dem Modell mit mittlerer Breite der aseismischen Lücke eine interessante Art von Interaktion beobachtet: Nach einigen seismischen Zyklen beginnen die Unebenheiten kontinuerlich und zwar in Gegen-Phasen zu brechen. Diese Ergebnisse sind andeutend, aber können aufgrund der geringen Auflösung des Models nur als vorläufig angesehen werden und erfordern weitere Bestätigung mit hoch-auflösenden Modellen. KW - earthquake modeling KW - seismic cycle modeling KW - subduction KW - numerical modeling KW - Erdbeben Modellierung KW - Modellierung des seismischen Zyklus KW - Subduktion KW - numerische Modellierung Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-398926 ER - TY - THES A1 - Martin, Sebastian T1 - Subduction zone wave guides : deciphering slab structure using intraslab seismicity at the Chile-Peru subduction zone T1 - Wellenleiter in Subduktionszonen : Bestimmung der Struktur der Chile-Peru Subduktionszone mit Hilfe von intra-platten Seismizität N2 - Subduction zones are regions of intense earthquake activity up to great depth. Sources are located inside the subducting lithosphere and, as a consequence, seismic radiation from subduction zone earthquakes is strongly affected by the interior slab structure. The wave field of these intraslab events observed in the forearc region is profoundly influenced by a seismically slow layer atop the slab surface. This several kilometer thick low-velocity channel (wave guide) causes the entrapment of seismic energy producing strong guided wave phases that appear in P onsets in certain regions of the forearc. Observations at the Chile-Peru subduction zone presented here, as well as observations at several other circum-pacific subduction zones show such signals. Guided wave analysis contributes details of immense value regarding the processes near the slab surface, such as layering of subducted lithosphere, source locations of intraslab seismicity and most of all, range and manner of mineralogical phase transitions. Seismological data stem from intermediate depth events (depth range 70 km - 300 km) recorded in northern Chile near 21 Grad S during the collaborative research initiative " Deformation Processes in the Andes" (SFB 267). A subset of stations - all located within a slab-parallel transect close to 69 Grad W - show low-frequency first arrivals (2 Hz), sometimes followed by a second high-frequency phase. We employ 2-dimensional finite-difference simulations of complete P-SV wave propagation to explore the parameter space of subduction zone wave guides and explain the observations. Key processes underlying the guided wave propagation are studied: Two distinct mechanisms of decoupling of trapped energy from the wave guide are analyzed - a prerequisite to observe the phases at stations located at large distances from the wave guide (up to 100 km). Variations of guided wave effects perpendicular to the strike of the subduction zone are investigated, such as the influence of phases traveling in the fast slab. Further, the merits and limits of guided wave analysis are assessed. Frequency spectra of the guided wave onsets prove to be a robust quantity that captures guided wave characteristics at subduction zones including higher mode excitation. They facilitate the inference of wave guide structure and source positioning: The peak frequency of the guided wave fundamental mode is associated with a certain combination of layer width and velocity contrast. The excitation strength of the guided wave fundamental mode and higher modes is associated with source position and orientation relative to the low-velocity layer. The guided wave signals at the Chile-Peru subduction zone are caused by energy that leaks from the subduction zone wave guide. On the one hand, the bend shape of the slab allows for leakage at a depth of 100 km. On the other, equalization of velocities between the wave guide and the host rocks causes further energy leakage at the contact zone between continental and oceanic crust (70 km depth). Guided waves bearing information on deep slab structure can therefore be recorded at specific regions in the forearc. These regions are determined based on slab geometry, and their locations coincide with the observations. A number of strong constraints on the structure of the Chile-Peru slab are inferred: The deep wave guide for intraslab events is formed by a layer of 2 km average width that remains seismically slow (7 percent velocity reduction compared to surrounding mantle). This low-velocity layer at the top of the Chile-Peru slab is imaged from a depth of 100 km down to at least 160 km. Intermediate depth events causing the observed phases are located inside the layer or directly beneath it in the slab mantle. The layer is interpreted as partially eclogized lower oceanic crust persisting to depth beyond the volcanic arc. N2 - Subduktionszonen sind bis in große Tiefen von intensiver Erdbebentätigkeit geprägt. Die Erdbebenquellen befinden sich in der subduzierten Lithosphäre (Slab), ihr Wellenfeld wird deshalb stark von der internen Slab-Struktur beeinflusst. Eine Schicht mit reduzierter seismischer Geschwindigkeit im oberen Bereich der Platte kann als Wellenleiter für diese Signale fungieren. In der nur wenige Kilometer dicken Schicht entstehen sogenannte geführte Wellen, die in Teilen des Forearc beobachtet werden. Diese Phasen bergen wertvolle Informationen über die Struktur nahe der Slab-Oberfläche, wie zum Beispiel Dicke der Schichtung, Herdlokationen und vor allem Tiefe und Art mineralogischer Umsetzungen. Die Beobachtungen stammen von mitteltiefen Beben (70 km - 300 km) im Untersuchungsgebiet in Nord-Chile und wurden im Rahmen des Sonderforschungsbereich 267 " Deformationsprozesse in den Anden" aufgezeichnet. Stationen in einem Streifen um 69 Grad W, der sich parallel zum Streichen der Subduktionszone erstreckt, zeigen niederfrequente Ersteinsätze, denen teilweise höherfrequente Phasen folgen. Mit Hilfe eines 2-dimensionalen Finite-Differenzen-Algorithmus werden die P-SV Wellenausbreitung simuliert, und die Beobachtungen erklärt. Zentrale Fragestellungen zu Wellenleitern in Subduktionszonen werden untersucht: Es werden zwei Mechanismen, die das Auskoppeln seismischer Energie aus dem Wellenleiter ermöglichen beschrieben - eine Grundvoraussetzung für das Auftreten von geführten Wellen in großen Entfernungen vom Wellenleiter (bis zu 100 km). Des weiteren werden Stärken und Grenzen der Analyse von geführten Wellen erörtert. Die Spektren der geführten Wellenzüge erweisen sich als robuste Messgröße, um die Charakteristika des Wellenleiters zu bestimmt. Struktur des Wellenleiters und Quellpositionen können so für festgelegte Quell-Empfänger-Geometrien abgeleitet werden. Die Peak-Frequenz der Grundmode wird durch eine Kombination aus Dicke der Schicht und Geschwindigkeitskontrast bestimmt. Die Stärke der Anregung der Grundmode und höherer Moden lässt auf die Lage und Orientierung der Erdbebenquelle relativ zur Schicht schließen. Geschwindigkeitskontrast, Schichtdicke und Quellposition sind von herausragender Bedeutung, um mineralogische Interpretationen des Wellenleiters zu überprüfen. Aufbauend auf die Simulationen werden die Beobachtungen interpretiert und Auskunft über die Struktur der Chile-Peru Subduktionszone erhalten: Eine dünne Schicht an der Slab-Oberfläche (durchschnittlich 2 km dick) trägt geringere seismische Geschwindigkeiten als der umgebende Mantel und fungiert als Wellenleiter für intra-platten Ereignisse in Tiefen von 100 bis mindestens 160 km. Ereignisse, die geführte Wellen hervorrufen, liegen in dieser Schicht oder direkt darunter im subduzierten Mantel. Um zu den Stationen in der Forearc-Region zu gelangen, entkoppelt ein Teil der geführten Wellen in einer Tiefe von circa 100 km aus der Niedergeschwindigkeitsschicht. Die Krümmung des Slab erlaubt das Austreten der Wellen und nimmt auch Einfluss auf die Pulsformen. Der Wellenleiter in der Chile-Peru Subduktionszone ergibt sich als unregelmäßige Schicht mit reduzierter seismischer Geschwindigkeit, in der geführte Wellen entstehen, in unterschiedlichen Tiefen wieder austreten, und an die freie Oberfläche gelangen. Die Beobachtungsgebiete befinden sich im Forearc und werden durch die Geometrie und Struktur der subduzierten Platte festgelegt. Die nur wenige Kilometer dicke, seismisch langsame Schicht an der Oberfläche des Chile-Peru Slab legt nahe, dass die Unterkruste der subduzierten Platte bis in große Tiefen besteht und nicht vollständig eklogitisiert ist. Abgeleitete Schichtdicke, Geschwindigkeitskontrast KW - Anden KW - Subduktion KW - Wellenleiter KW - Erdbeben KW - ozeanische Kruste KW - Gabbro-Eklogit KW - Transformation KW - geführte Wellen KW - subduction KW - guided waves KW - oceanic crust KW - Andes KW - finite difference simulation Y1 - 2005 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-5820 ER - TY - JOUR A1 - Ohuchi, Tomohiro A1 - Nishihara, Yu A1 - Kawazoe, Takaaki A1 - Spengler, Dirk A1 - Shiraishi, Rei A1 - Suzuki, Akio A1 - Kikegawa, Takumi A1 - Ohtani, Eiji T1 - Superplasticity in hydrous melt-bearing dunite Implications for shear localization in Earth's upper mantle JF - Earth & planetary science letters N2 - Deformation experiments on hydrous melt-bearing dunite (olivine+4 vol% orthopyroxene+4 vol% clinopyroxene with less than 2.5 vol% of the melt phase) were conducted at pressures of 1.3-5.7 GPa and temperatures of 1270-1490 K in order to explore the effect of intergranular fluids on the plastic flow of olivine in Earth's upper mantle. The strain rate was proportional to steady-state creep strength to the 2.1 power, and the creep strength markedly increased with increase in grain size. Developments of the crystallographic preferred orientation of olivine and flattening of olivine grains were hardly observed even after 33-55% shortening of the samples. These observations show that grain boundary sliding (GBS) dominated the deformation of olivine (i.e., superplasticity). The creep strength of hydrous melt-bearing dunite was 2-5 times lower than that of melt-free dunite. The dependence of creep rate on melt fraction is known to be expressed empirically as (epsilon) over dot(phi) = (epsilon) over dot(0) exp(alpha phi), where alpha is a constant and phi is the melt fraction. The experimentally obtained value of alpha was in the range of 150-230, corresponding to 5-7 times the reported values for the olivine-basalt system at 0.3 GPa (i.e., creep strength of dunite was efficiently reduced by the hydrous melt). Superplasticity is the dominant creep mechanism of olivine in fluid-bearing fine-grained peridotites under low-temperature and high-stress conditions (i.e., peridotite shear zones in the upper mantle). Superplasticity induced by geological fluids would play an important role in the shear localization (and thus initiation of subduction) in the upper mantle. KW - olivine KW - hydrous melt KW - grain boundary sliding KW - superplasticity KW - shear localization KW - subduction Y1 - 2012 U6 - https://doi.org/10.1016/j.epsl.2012.04.032 SN - 0012-821X VL - 335 IS - 12 SP - 59 EP - 71 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Konrad-Schmolke, Matthias A1 - Babist, Jochen A1 - Handy, Mark R. A1 - O'brien, Patrick J. T1 - The physico-chemical properties of a subducted slab from garnet zonation patterns (Sesia Zone, Western Alps) JF - Journal of petrology N2 - Garnets in continentally derived high-pressure (HP) rocks of the Sesia Zone (Western Alps) exhibit three different chemical zonation patterns, depending on sample locality. Comparison of observed garnet zonation patterns with thermodynamically modelled patterns shows that the different patterns are caused by differences in the water content of the subducted protoliths during prograde metamorphism. Zonation patterns of garnets in water-saturated host rocks show typical prograde chemical zonations with steadily increasing pyrope content and increasing XMg, together with bell-shaped spessartine patterns. In contrast, garnets in water-undersaturated rocks have more complex zonation patterns with a characteristic decrease in pyrope and XMg between core and inner rim. In some cases, garnets show an abrupt compositional change in core-to-rim profiles, possibly due to water-undersaturation prior to HP metamorphism. Garnets from both water-saturated and water-undersaturated rocks show signs of intervening growth interruptions and core resorption. This growth interruption results from bulk-rock depletion caused by fractional garnet crystallization. The water content during burial influences significantly the physical properties of the subducted rocks. Due to enhanced garnet crystallization, water-undersaturated rocks, i.e. those lacking a free fluid phase, become denser than their water-saturated equivalents, facilitating the subduction of continental material. Although water-bearing phases such as phengite and epidote are stable up to eclogite-facies conditions in these rocks, dehydration reactions during subduction are lacking in water-undersaturated rocks up to the transition to the eclogite facies, due to the thermodynamic stability of such hydrous phases at high P-T conditions. Our calculations show that garnet zonation patterns strongly depend on the mineral parageneses stable during garnet growth and that certain co-genetic mineral assemblages cause distinct garnet zonation patterns. This observation enables interpretation of complex garnet growth zonation patterns in terms of garnet-forming reactions and water content during HP metamorphism, as well determination of detailed P-T paths. KW - dehydration KW - high-pressure metamorphism KW - Sesia Zone KW - subduction KW - thermodynamic modelling Y1 - 2006 U6 - https://doi.org/10.1093/petrology/egl039 SN - 0022-3530 VL - 47 IS - 11 SP - 2123 EP - 2148 PB - Oxford Univ. Press CY - Oxford ER - TY - THES A1 - Barrionuevo, Matías T1 - The role of the upper plate in the Andean tectonic evolution (33-36°S): insights from structural geology and numerical modeling T1 - El rol de la placa superior en la evolución tectónica andina (33-36°S): aportes desde la geología estructural y el modelado numérico T1 - Die Rolle der oberen Platte in der tektonischen Entwicklung der Anden (33-36°S): Erkenntnisse aus der Strukturgeologie und der numerischen Modellierung N2 - Los Andes Centrales del Sur (33-36°S) son un gran laboratorio para el estudio de los procesos de deformación orogénica, donde las condiciones de borde, como la geometría de la placa subductada, imponen un importante control sobre la deformación andina. Por otro lado, la Placa Sudamericana presenta una serie de heterogeneidades que también imparten un control sobre el modo de deformación. El objetivo de esta tesis es probar el control de este último factor sobre la construcción del sistema orogénico andino. A partir de la integración de la información superficial y de subsuelo en el área sur (34°-36°S), se estudió la evolución de la deformación andina sobre el segmento de subducción normal. Se desarrolló un modelo estructural que evalúa el estado de esfuerzos desde el Mioceno hasta la actualidad, el rol de estructuras previas y su influencia en la migración de fluidos. Con estos datos y publicaciones previas de la zona norte del área de estudio (33°-34ºS), se realizó un modelado numérico geodinámico para probar la hipótesis del papel de las heterogeneidades de la placa superior en la evolución andina. Se utilizaron dos códigos (LAPEX-2D y ASPECT) basados en elementos finitos/diferencias finitas, que simulan el comportamiento de materiales con reologías elastoviscoplásticas bajo deformación. Los resultados del modelado sugieren que la deformación contraccional de la placa superior está significativamente controlada por la resistencia de la litósfera, que está definida por la composición de la corteza superior e inferior y por la proporción del manto litosférico, que a su vez está definida por eventos tectónicos previos. Estos eventos previos también definieron la composición de la corteza y su geometría, que es otro factor que controla la localización de la deformación. Con una composición de corteza inferior más félsica, la deformación sigue un modo de cizalla pura mientras que las composiciones más máficas provocan un modo de deformación tipo cizalla simple. Por otro lado, observamos que el espesor inicial de la litósfera controla la localización de la deformación, donde zonas con litósfera más fina es propensa a concentrar la deformación. Un límite litósfera-astenósfera asimétrico, como resultado del flujo de la cuña mantélica tiende a generar despegues vergentes al E. N2 - The Southern Central Andes (33°-36°S) are an excellent natural laboratory to study orogenic deformation processes, where boundary conditions, such as the geometry of the subducted plate, impose an important control on the evolution of the orogen. On the other hand, the South American plate presents a series of heterogeneities that additionally impart control on the mode of deformation. This thesis aims to test the control of this last factor over the construction of the Cenozoic Andean orogenic system. From the integration of surface and subsurface information in the southern area (34-36°S), the evolution of Andean deformation over the steeply dipping subduction segment was studied. A structural model was developed evaluating the stress state from the Miocene to the present-day and its influence in the migration of magmatic fluids and hydrocarbons. Based on these data, together with the data generated by other researchers in the northern zone of the study area (33-34°S), geodynamic numerical modeling was performed to test the hypothesis of the decisive role of upper-plate heterogeneities in the Andean evolution. Geodynamic codes (LAPEX-2D and ASPECT) which simulate the behavior of materials with elasto-visco-plastic rheologies under deformation, were used. The model results suggest that upper-plate contractional deformation is significantly controlled by the strength of the lithosphere, which is defined by the composition of the upper and lower crust, and by the proportion of lithospheric mantle, which in turn is determined by previous tectonic events. In addition, the previous regional tectono-magmatic events also defined the composition of the crust and its geometry, which is another factor that controls the localization of deformation. Accordingly, with more felsic lower crustal composition, the deformation follows a pure-shear mode, while more mafic compositions induce a simple-shear deformation mode. On the other hand, it was observed that initial lithospheric thickness may fundamentally control the location of deformation, with zones characterized by thin lithosphere are prone to concentrate it. Finally, it was found that an asymmetric lithosphere-astenosphere boundary resulting from corner flow in the mantle wedge of the eastward-directed subduction zone tends to generate east-vergent detachments. N2 - Die südlichen Zentralanden (33°-36°S) sind eine ausgezeichnete, natürliche Forschungsumgebung zur Untersuchung gebirgsbildender Deformationsprozesse, in der Randbedingungen, wie die Geometrie der subduzierten Platte, einen starken Einfluss auf die Evolution des Gebirges besitzen. Anderseits sind die Deformationsmechanismen geprägt von der Heterogenität der Südamerikanischen Platte. In dieser Arbeit wird die Bedeutung dieses Mechanismus für die Herausbildung der Anden während des Känozoikums untersucht. Im südlichen Teil (34-36°S), in dem die subduzierte Platte in einem steileren Winkel in den Erdmantel absinkt, wird die Entwicklung der Andendeformation mithilfe von oberflächlich aufgezeichneten und in tiefere Erdschichten reichenden Daten untersucht. Das darauf aufbauende Strukturmodell ermöglicht die Abschätzung der tektonischen Spannungen vom Miozän bis in die Neuzeit und den Einfluss der Bewegungen von magmatischen Fluiden, sowie Kohlenwasserstoffen. Auf Grundlage dieser Daten und solcher, die von Wissenschaftlern im nördlichen Bereich des Untersuchungsgebietes (33-34°S) erfasst wurden, wurde eine geodynamische, numerische Modellierung durchgeführt, um die Hypothese des Einflusses der Heterogenität der oberen Platten auf die Gebirgsbildung der Anden zu überprüfen. Die genutzte geodynamische Softwares (LAPEX-2D und ASPECT) simulieren das Verhalten von elasto-viskoplastischen Materialien, wenn diese unter Spannung stehen. Die Modellierungsergebnisse zeigen, dass die Kontraktionsprozesse hauptsächlich durch die Stärke der Lithosphäre beeinflusst werden. Diese Kenngröße wird aus der Zusammensetzung von Ober- und Unterkruste und dem Anteil des lithosphärischen Mantels, der durch vorhergehende tektonische Vorgänge überprägt ist, bestimmt. Diese räumlich begrenzten tektono-magmatischen Events definieren ebenfalls die Zusammensetzung und die Geometrie der Erdkruste, welche einen großen Einfluss auf das räumliche Auftreten von Deformationsprozessen hat. Eine eher felsische Unterkruste führt vorrangig zu pure-shear, während eine eher mafisch zusammengesetzte Unterkruste primär zu einem Deformationsmechanismus führt, der simple-shear genannt wird. Weiterhing wurde beobachtet, dass die Dicke der Lithosphäre vor der Deformation einen fundamentalen Einfluss auf die räumliche Eingrenzung von Deformation hat, wobei Regionen mit einer dünnen Lithosphärenschicht verstärkt Deformation aufweisen. Eine asymmetrische Grenzschicht zwischen Lithosphäre und Asthenosphäre ist das Resultat von Fließprozessen im Erdmantel, im Keil zwischen der obenliegenden Platte und der sich ostwärts absinkenden Subduktionszone, und verstärkt die Herausbildung von nach Osten gerichteten Abscherungen in der Erdkruste. KW - structural geology KW - tectonics KW - subduction KW - geodynamic modeling KW - geodynamische Modellierung KW - Strukturgeologie KW - Subduktion KW - Tektonik Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-515909 ER - TY - GEN 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, Halilbağı Complex (Anatolia) T2 - Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 632 KW - subduction KW - Lu/Hf dating of garnet KW - metamorphic sole KW - eclogite KW - blueschist KW - lawsonite Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-424651 SN - 1866-8372 IS - 632 ER - TY - THES A1 - Rodriguez Piceda, Constanza T1 - Thermomechanical state of the southern Central Andes T1 - Thermomechanischer Zustand der südlichen Zentral Anden BT - implications for active deformation patterns in the transition from flat to steep subduction BT - Implikationen für aktive Deformationsmuster beim Übergang von flacher zu steiler Subduktion N2 - The Andes are a ~7000 km long N-S trending mountain range developed along the South American western continental margin. Driven by the subduction of the oceanic Nazca plate beneath the continental South American plate, the formation of the northern and central parts of the orogen is a type case for a non-collisional orogeny. In the southern Central Andes (SCA, 29°S-39°S), the oceanic plate changes the subduction angle between 33°S and 35°S from almost horizontal (< 5° dip) in the north to a steeper angle (~30° dip) in the south. This sector of the Andes also displays remarkable along- and across- strike variations of the tectonic deformation patterns. These include a systematic decrease of topographic elevation, of crustal shortening and foreland and orogenic width, as well as an alternation of the foreland deformation style between thick-skinned and thin-skinned recorded along- and across the strike of the subduction zone. Moreover, the SCA are a very seismically active region. The continental plate is characterized by a relatively shallow seismicity (< 30 km depth) which is mainly focussed at the transition from the orogen to the lowland areas of the foreland and the forearc; in contrast, deeper seismicity occurs below the interiors of the northern foreland. Additionally, frequent seismicity is also recorded in the shallow parts of the oceanic plate and in a sector of the flat slab segment between 31°S and 33°S. The observed spatial heterogeneity in tectonic and seismic deformation in the SCA has been attributed to multiple causes, including variations in sediment thickness, the presence of inherited structures and changes in the subduction angle of the oceanic slab. However, there is no study that inquired the relationship between the long-term rheological configuration of the SCA and the spatial deformation patterns. Moreover, the effects of the density and thickness configuration of the continental plate and of variations in the slab dip angle in the rheological state of the lithosphere have been not thoroughly investigated yet. Since rheology depends on composition, pressure and temperature, a detailed characterization of the compositional, structural and thermal fields of the lithosphere is needed. Therefore, by using multiple geophysical approaches and data sources, I constructed the following 3D models of the SCA lithosphere: (i) a seismically-constrained structural and density model that was tested against the gravity field; (ii) a thermal model integrating the conversion of mantle shear-wave velocities to temperature with steady-state conductive calculations in the uppermost lithosphere (< 50 km depth), validated by temperature and heat-flow measurements; and (iii) a rheological model of the long-term lithospheric strength using as input the previously-generated models. The results of this dissertation indicate that the present-day thermal and rheological fields of the SCA are controlled by different mechanisms at different depths. At shallow depths (< 50 km), the thermomechanical field is modulated by the heterogeneous composition of the continental lithosphere. The overprint of the oceanic slab is detectable where the oceanic plate is shallow (< 85 km depth) and the radiogenic crust is thin, resulting in overall lower temperatures and higher strength compared to regions where the slab is steep and the radiogenic crust is thick. At depths > 50 km, largest temperatures variations occur where the descending slab is detected, which implies that the deep thermal field is mainly affected by the slab dip geometry. The outcomes of this thesis suggests that long-term thermomechanical state of the lithosphere influences the spatial distribution of seismic deformation. Most of the seismicity within the continental plate occurs above the modelled transition from brittle to ductile conditions. Additionally, there is a spatial correlation between the location of these events and the transition from the mechanically strong domains of the forearc and foreland to the weak domain of the orogen. In contrast, seismicity within the oceanic plate is also detected where long-term ductile conditions are expected. I therefore analysed the possible influence of additional mechanisms triggering these earthquakes, including the compaction of sediments in the subduction interface and dehydration reactions in the slab. To that aim, I carried out a qualitative analysis of the state of hydration in the mantle using the ratio between compressional- and shear-wave velocity (vp/vs ratio) from a previous seismic tomography. The results from this analysis indicate that the majority of the seismicity spatially correlates with hydrated areas of the slab and overlying continental mantle, with the exception of the cluster within the flat slab segment. In this region, earthquakes are likely triggered by flexural processes where the slab changes from a flat to a steep subduction angle. First-order variations in the observed tectonic patterns also seem to be influenced by the thermomechanical configuration of the lithosphere. The mechanically strong domains of the forearc and foreland, due to their resistance to deformation, display smaller amounts of shortening than the relatively weak orogenic domain. In addition, the structural and thermomechanical characteristics modelled in this dissertation confirm previous analyses from geodynamic models pointing to the control of the observed heterogeneities in the orogen and foreland deformation style. These characteristics include the lithospheric and crustal thickness, the presence of weak sediments and the variations in gravitational potential energy. Specific conditions occur in the cold and strong northern foreland, which is characterized by active seismicity and thick-skinned structures, although the modelled crustal strength exceeds the typical values of externally-applied tectonic stresses. The additional mechanisms that could explain the strain localization in a region that should resist deformation are: (i) increased tectonic forces coming from the steepening of the slab and (ii) enhanced weakening along inherited structures from pre-Andean deformation events. Finally, the thermomechanical conditions of this sector of the foreland could be a key factor influencing the preservation of the flat subduction angle at these latitudes of the SCA. N2 - Die Anden sind eine ~7000 km lange N-S-verlaufende Hochgebirgskette, die entlang des westlichen südamerikanischen Kontinentalrandes entstanden ist. Aufgrund der Subduktion der ozeanischen Nazca-Platte unter die kontinentale südamerikanische Platte ist die Bildung des nördlichen und zentralen Teils des Gebirges typisch für eine nicht-kollisionale Orogenese. In den südlichen Zentralanden (SZA, 29-39° S) verändert sich der Subduktionswinkel der ozeanischen Platte zwischen 33 ° S und 35 ° S von fast horizontal (< 5° Einfallen) im Norden zu einem steileren Winkel (~ 30 ° Einfallen) im Süden. Begleitet wird dieser Trend von systematischen, Süd-gerichteten Abnahmen der topographischen Erhebung, der Krusteneinengung und der Vorland- und Orogenbreite, sowie von Variationen im Deformationsstil des Vorlandes, wo die Einengung des Deckgebirges in unterschiedlichem Maße von einer entsprechenden Deformation des Grundgebirges begleitet wird. . Darüber hinaus sind die SZA eine seismisch sehr aktive Region. Die Kontinentalplatte zeichnet sich durch eine relativ flache Seismizität (< 30 km Tiefe) aus, die sich hauptsächlich auf die Übergänge vom Orogen zu den Vorlandbereichen konzentriert; im Gegensatz dazu tritt tiefere Seismizität in den zentralen Bereichen des nördlichen Vorlandes auf. Darüber hinaus ist häufig auftretende Seismizität auch in den flachen Teilen der ozeanischen Platte und im Plattensegment mit flach einfallender Subduktion zwischen 31 ° S und 33 ° S festzustellen. Die beobachtete räumliche Heterogenität der tektonischen und seismischen Deformation in den SZA wurde auf mehrere Ursachen zurückgeführt, darunter Schwankungen der Sedimentmächtigkeit, das Vorhandensein vererbter Strukturen und Veränderungen des Subduktionswinkels der ozeanischen Platte. Es gibt jedoch bislang keine Studie, die den Zusammenhang zwischen der langfristigen rheologischen Konfiguration der SZA und den räumlichen Deformationsmustern untersucht hat. Darüber hinaus wurden die Auswirkungen der Dichte- und Mächtigkeitsvariationen in der kontinentalen Oberplatte und der verschiedenen Subduktionswinkel auf den rheologischen Zustand der Lithosphäre noch nicht grundlegend untersucht. Da die Rheologie von der Gesteinsart, dem Druck und der Temperatur abhängt, ist eine detaillierte Charakterisierung der Zusammensetzung, Struktur und des thermischen Feldes der Lithosphäre erforderlich. Daher habe ich unter Verwendung kombinierter Modellierungsansätze und geophysikalischer Daten die folgenden 3D Modelle für die Lithosphäre der SZA konstruiert: (i) ein auf seismischen Daten basierendes Struktur- und Dichtemodell, das anhand des beobachteten Schwerefeldes validiert wurde; (ii) ein thermisches Modell, das die Umwandlung von Mantelscherwellengeschwindigkeiten in Temperaturen mit Berechnungen des konduktiven Wärmetransports für stationäre Bedingungen in der obersten Lithosphäre (<50 km Tiefe) integriert und durch Temperatur- und Wärmeflussmessungen validiert wurde; und (iii) ein rheologisches Modell der langfristig bedingten Lithosphärenfestigkeit, das auf den zuvor erzeugten Modellen gründet. Die Ergebnisse dieser Dissertation zeigen, dass die thermischen und rheologischen Bedingungen in den heutigen SZA durch verschiedene Mechanismen in unterschiedlichen Tiefen gesteuert werden. In flachen Tiefen (< 50 km) wird das thermomechanische Feld durch die heterogene Zusammensetzung der kontinentalen Lithosphäre differenziert. Eine Überprägung durch die ozeanische Platte ist dort nachweisbar, wo die ozeanische Platte flach (< 85 km tief) und die radiogene Kruste dünn ist, was insgesamt zu niedrigeren Temperaturen und einer höheren Festigkeit im Vergleich zu Bereichen führt, in denen die Platte steil einfällt und die radiogene Kruste dick ist. In Tiefen > 50 km treten die größten Temperaturschwankungen dort auf, wo die subduzierten Platte nachgewiesen wurde, was bedeutet, dass das tiefe thermische Feld den Subduktionswinkel gesteuert wird. Die Ergebnisse dieser Doktorarbeit legen nahe, dass der langfristige thermomechanische Zustand der Lithosphäre die räumliche Verteilung rezenter Seismizität beeinflusst. Der größte Anteil innerhalb der Kontinentalplatte registrierter Erdbebentätigkeit tritt oberhalb des modellierten Übergangs von spröden zu duktilen Bedingungen auf. Außerdem besteht eine räumliche Korrelation zwischen Erdbebenclustern und den Übergängen von den mechanisch rigideren Vorlandbereichen (Forearc und Foreland) zum mechanisch schwächeren Orogen. Demgegenüber wird vermehrte Seismizität innerhalb der ozeanischen Platte auch dort nachgewiesen, wo entsprechend der Modellierung langfristig duktile Bedingungen erwartet werden. Ich habe daher den möglichen Einfluss zusätzlicher Mechanismen untersucht, die ein Auslösen dieser Erdbeben begünstigen könnten, darunter die Kompaktion von Sedimenten an der Subduktionsgrenzfläche und Dehydrationsreaktionen innerhalb der Platte. Dazu habe ich eine qualitative Analyse des Hydratationszustandes des Mantels unter Verwendung des Verhältnisses zwischen Kompressions- und Scherwellengeschwindigkeit (Vp/Vs-Verhältnis aus einemseismischen Tomographiemodell) durchgeführt. Die Ergebnisse dieser Analyse zeigen, dass der Großteil der Seismizität räumlich mit hydratisierten Bereichen in der subduzierten Platte und im darüber liegenden kontinentalen Mantel korreliert, mit Ausnahme eines Erdbebenclusters, das innerhalb des flachen Plattensegments auftritt. In diesem Bereich wechselt die subduzierte Platte von einem flachen in einen steilen Subduktionswinkel und Erdbeben werden wahrscheinlich durch Biegevorgänge in der Platte ausgelöst. Auch die wichtigsten Variationen in den beobachteten tektonischen Mustern scheinen durch die thermomechanische Konfiguration der Lithosphäre beeinflusst zu sein. Die mechanisch starken Bereiche von Forearc und Foreland zeigen aufgrund ihrer Verformungsbeständigkeit geringere Verkürzungsraten als der relativ schwache Bereich des Orogens. Darüber hinaus bestätigen die in dieser Dissertation modellierten strukturellen und thermomechanischen Eigenschaften der Lithosphäre auch frühere Analysen geodynamischer Simulationen, denen zufolge der Deformationsstil im Orogen- und Vorlandbereich jeweils von Variationen in der Lithosphären- und Krustendicke, im Vorhandensein schwacher Sedimente und in der gravitativen potentiellen Energie kontrolliert wird. Eine Sonderstellung nimmt der nordöstliche Vorlandbereich der SZA ein, wo eine verstärkte Seismizität und eine das Deck-und Grundgebirge erfassende Deformation zu beobachten sind, obwohl die modellierte Krustenfestigkeit dort Werte übersteigt, die für die in diesem Gebiet anzunehmenden tektonischen Spannungen typisch wären. . Mechanismen zur Lokalisierung verstärkter Deformation in einem Gebiet beitragen können, das nach den vorliegenden Modellen einer tektonischen Verformung widerstehen sollte, sind: (i) erhöhte tektonische Kräfte durch ein steileres Abtauchen der Platte und (ii) Schwächezonen in der Kruste, die auf prä-andine Deformationsereignisse zurückgehen. Schließlich könnten die thermomechanischen Bedingungen in diesem Teil des Vorlands einchlüsselfaktor für die Erhaltung des flachen Subduktionswinkels in diesen Breiten der SZA sein. KW - Andes KW - Anden KW - subduction KW - Subduktion KW - lithosphere KW - Lithosphäre KW - earthquakes KW - Erdbeben KW - modelling KW - Modellierung Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-549275 ER - TY - JOUR A1 - Muldashev, Iskander A. A1 - Sobolev, Stephan T1 - What controls maximum magnitudes of giant subduction earthquakes? JF - Geochemistry, geophysics, geosystems N2 - Giant earthquakes with magnitudes above 8.5 occur only in subduction zones. Despite the developments made in observing large subduction zone earthquakes with geophysical instruments, the factors controlling the maximum size of these earthquakes are still poorly understood. Previous studies have suggested the importance of slab shape, roughness of the plate interface contact, state of the strain in the upper plate, thickness of sediments filling the trenches, and subduction rate. Here, we present 2-D cross-scale numerical models of seismic cycles for subduction zones with various geometries, subduction channel friction configurations, and subduction rates. We found that low-angle subduction and thick sediments in the subduction channel are the necessary conditions for generating giant earthquakes, while the subduction rate has a negligible effect. We suggest that these key parameters determine the maximum magnitude of a subduction earthquake by controlling the seismogenic zone width and smoothness of the subduction interface. This interpretation supports previous studies that are based upon observations and scaling laws. Our modeling results also suggest that low static friction in the sediment-filled subduction channel results in neutral or moderate compressive deformation in the overriding plate for low-angle subduction zones hosting giant earthquakes. These modeling results agree well with observations for the largest earthquakes. Based on our models we predict maximum magnitudes of subduction earthquakes worldwide, demonstrating the fit to magnitudes of all giant earthquakes of the 20th and 21st centuries and good agreement with the predictions based on statistical analyses of observations. KW - giant earthquakes KW - earthquake modeling KW - subduction Y1 - 2020 U6 - https://doi.org/10.1029/2020GC009145 SN - 1525-2027 VL - 21 IS - 9 PB - American Geophysical Union CY - Washington ER -