@article{RodriguezPicedaScheckWenderothCacaceetal.2022,
  author    = {Rodriguez Piceda, Constanza and Scheck-Wenderoth, Magdalena and Cacace, Mauro and Bott, Judith and Strecker, Manfred},
  title     = {Long-Term Lithospheric Strength and Upper-Plate Seismicity in the Southern Central Andes, 29 degrees-39 degrees S},
  series = {Geochemistry, geophysics, geosystems},
  volume    = {23},
  journal   = {Geochemistry, geophysics, geosystems},
  number    = {3},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {1525-2027},
  doi       = {10.1029/2021GC010171},
  pages     = {22},
  year      = {2022},
  abstract  = {We examined the relationship between the mechanical strength of the lithosphere and the distribution of seismicity within the overriding continental plate of the southern Central Andes (SCA, 29 degrees-39 degrees S), where the oceanic Nazca Plate changes its subduction angle between 33 degrees S and 35 degrees S, from subhorizontal in the north (<5 degrees) to steep in the south (similar to 30 degrees). We computed the long-term lithospheric strength based on an existing 3D model describing variations in thickness, density, and temperature of the main geological units forming the lithosphere of the SCA and adjacent forearc and foreland regions. The comparison between our results and seismicity within the overriding plate (upper-plate seismicity) shows that most of the events occur within the modeled brittle domain of the lithosphere. The depth where the deformation mode switches from brittle frictional to thermally activated ductile creep provides a conservative lower bound to the seismogenic zone in the overriding plate of the study area. We also found that the majority of upper-plate earthquakes occurs within the realm of first-order contrasts in integrated strength (12.7-13.3 log Pam in the Andean orogen vs. 13.5-13.9 log Pam in the forearc and the foreland). Specific conditions characterize the mechanically strong northern foreland of the Andes, where seismicity is likely explained by the effects of slab steepening.},
  language  = {en}
}
@article{BaesSobolevGeryaetal.2020,
  author    = {Baes, Marzieh and Sobolev, Stephan and Gerya, Taras V. and Brune, Sascha},
  title     = {Plume-induced subduction initiation},
  series = {Geochemistry, geophysics, geosystems},
  volume    = {21},
  journal   = {Geochemistry, geophysics, geosystems},
  number    = {2},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {1525-2027},
  doi       = {10.1029/2019GC008663},
  pages     = {19},
  year      = {2020},
  abstract  = {Initiation of subduction following the impingement of a hot buoyant mantle plume is one of the few scenarios that allow breaking the lithosphere and recycling a stagnant lid without requiring any preexisting weak zones. Here, we investigate factors controlling the number and shape of retreating subducting slabs formed by plume-lithosphere interaction. Using 3-D thermomechanical models we show that the deformation regime, which defines formation of single-slab or multi-slab subduction, depends on several parameters such as age of oceanic lithosphere, thickness of the crust and large-scale lithospheric extension rate. Our model results indicate that on present-day Earth multi-slab plume-induced subduction is initiated only if the oceanic lithosphere is relatively young (<30-40 Myr, but >10 Myr), and the crust has a typical thickness of 8 km. In turn, development of single-slab subduction is facilitated by older lithosphere and pre-imposed extensional stresses. In early Earth, plume-lithosphere interaction could have led to formation of either episodic short-lived circular subduction when the oceanic lithosphere was young or to multi-slab subduction when the lithosphere was old.},
  language  = {en}
}
@misc{BaesSobolevGeryaetal.2020,
  author    = {Baes, Marzieh and Sobolev, Stephan Vladimir and Gerya, Taras V. and Brune, Sascha},
  title     = {Plume-induced subduction initiation},
  series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  number    = {2},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-52274},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-522742},
  pages     = {21},
  year      = {2020},
  abstract  = {Initiation of subduction following the impingement of a hot buoyant mantle plume is one of the few scenarios that allow breaking the lithosphere and recycling a stagnant lid without requiring any preexisting weak zones. Here, we investigate factors controlling the number and shape of retreating subducting slabs formed by plume-lithosphere interaction. Using 3-D thermomechanical models we show that the deformation regime, which defines formation of single-slab or multi-slab subduction, depends on several parameters such as age of oceanic lithosphere, thickness of the crust and large-scale lithospheric extension rate. Our model results indicate that on present-day Earth multi-slab plume-induced subduction is initiated only if the oceanic lithosphere is relatively young (<30-40 Myr, but >10 Myr), and the crust has a typical thickness of 8 km. In turn, development of single-slab subduction is facilitated by older lithosphere and pre-imposed extensional stresses. In early Earth, plume-lithosphere interaction could have led to formation of either episodic short-lived circular subduction when the oceanic lithosphere was young or to multi-slab subduction when the lithosphere was old.},
  language  = {en}
}
@article{BaesSobolevGeryaetal.2020,
  author    = {Baes, Marzieh and Sobolev, Stephan V. and Gerya, Taras V. and Brune, Sascha},
  title     = {Subduction initiation by Plume-Plateau interaction},
  series = {Geochemistry, geophysics, geosystems},
  volume    = {21},
  journal   = {Geochemistry, geophysics, geosystems},
  number    = {8},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {1525-2027},
  doi       = {10.1029/2020GC009119},
  pages     = {17},
  year      = {2020},
  abstract  = {It has recently been demonstrated that the interaction of a mantle plume with sufficiently old oceanic lithosphere can initiate subduction. However, the existence of large lithospheric heterogeneities, such as a buoyant plateau, in proximity to a rising plume head may potentially hinder the formation of a new subduction zone. Here, we investigate this scenario by means of 3-D numerical thermomechanical modeling. We explore how plume-lithosphere interaction is affected by lithospheric age, relative location of plume head and plateau border, and the strength of the oceanic crust. Our numerical experiments suggest four different geodynamic regimes: (a) oceanic trench formation, (b) circular oceanic-plateau trench formation, (c) plateau trench formation, and (d) no trench formation. We show that regardless of the age and crustal strength of the oceanic lithosphere, subduction can initiate when the plume head is either below the plateau border or at a distance less than the plume radius from the plateau edge. Crustal heterogeneity facilitates subduction initiation of old oceanic lithosphere. High crustal strength hampers the formation of a new subduction zone when the plume head is located below a young lithosphere containing a thick and strong plateau. We suggest that plume-plateau interaction in the western margin of the Caribbean could have resulted in subduction initiation when the plume head impinged onto the oceanic lithosphere close to the border between plateau and oceanic crust.},
  language  = {en}
}
@misc{WongMasonBruneetal.2019,
  author    = {Wong, Kevin and Mason, Emily and Brune, Sascha and East, Madison and Edmonds, Marie and Zahirovic, Sabin},
  title     = {Deep Carbon Cycling Over the Past 200 Million Years: A Review of Fluxes in Different Tectonic Settings},
  series = {Frontiers in Earth Science},
  volume    = {7},
  journal   = {Frontiers in Earth Science},
  publisher = {Frontiers Research Foundation},
  address   = {Lausanne},
  issn      = {2296-6463},
  doi       = {10.3389/feart.2019.00263},
  pages     = {22},
  year      = {2019},
  language  = {en}
}
@article{RajabiZieglerTingayetal.2016,
  author    = {Rajabi, Mojtaba and Ziegler, Moritz O. and Tingay, Mark and Heidbach, Oliver and Reynolds, Scott},
  title     = {Contemporary tectonic stress pattern of the Taranaki Basin, New Zealand},
  series = {Journal of geophysical research : Solid earth},
  volume    = {121},
  journal   = {Journal of geophysical research : Solid earth},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {2169-9313},
  doi       = {10.1002/2016JB013178},
  pages     = {6053 -- 6070},
  year      = {2016},
  abstract  = {The present-day stress state is a key parameter in numerous geoscientific research fields including geodynamics, seismic hazard assessment, and geomechanics of georeservoirs. The Taranaki Basin of New Zealand is located on the Australian Plate and forms the western boundary of tectonic deformation due to Pacific Plate subduction along the Hikurangi margin. This paper presents the first comprehensive wellbore-derived basin-scale in situ stress analysis in New Zealand. We analyze borehole image and oriented caliper data from 129 petroleum wells in the Taranaki Basin to interpret the shape of boreholes and determine the orientation of maximum horizontal stress (S-Hmax). We combine these data (151 S-Hmax data records) with 40 stress data records derived from individual earthquake focal mechanism solutions, 6 from stress inversions of focal mechanisms, and 1 data record using the average of several focal mechanism solutions. The resulting data set has 198 data records for the Taranaki Basin and suggests a regional S-Hmax orientation of N068 degrees E (22 degrees), which is in agreement with NW-SE extension suggested by geological data. Furthermore, this ENE-WSW average S-Hmax orientation is subparallel to the subduction trench and strike of the subducting slab (N50 degrees E) beneath the central western North Island. Hence, we suggest that the slab geometry and the associated forces due to slab rollback are the key control of crustal stress in the Taranaki Basin. In addition, we find stress perturbations with depth in the vicinity of faults in some of the studied wells, which highlight the impact of local stress sources on the present-day stress rotation.},
  language  = {en}
}
@phdthesis{JaraMunoz2016,
  author    = {Jara Mu{\~n}oz, Julius},
  title     = {Quantifying forearc deformation patterns using coastal geomorphic markers},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-102652},
  school      = {Universit{\"a}t Potsdam},
  pages     = {XXV, 213},
  year      = {2016},
  abstract  = {Rapidly uplifting coastlines are frequently associated with convergent tectonic boundaries, like subduction zones, which are repeatedly breached by giant megathrust earthquakes. The coastal relief along tectonically active realms is shaped by the effect of sea-level variations and heterogeneous patterns of permanent tectonic deformation, which are accumulated through several cycles of megathrust earthquakes. However, the correlation between earthquake deformation patterns and the sustained long-term segmentation of forearcs, particularly in Chile, remains poorly understood. Furthermore, the methods used to estimate permanent deformation from geomorphic markers, like marine terraces, have remained qualitative and are based on unrepeatable methods. This contrasts with the increasing resolution of digital elevation models, such as Light Detection and Ranging (LiDAR) and high-resolution bathymetric surveys. Throughout this thesis I study permanent deformation in a holistic manner: from the methods to assess deformation rates, to the processes involved in its accumulation. My research focuses particularly on two aspects: Developing methodologies to assess permanent deformation using marine terraces, and comparing permanent deformation with seismic cycle deformation patterns under different spatial scales along the M8.8 Maule earthquake (2010) rupture zone. Two methods are developed to determine deformation rates from wave-built and wave-cut terraces respectively. I selected an archetypal example of a wave-built terrace at Santa Maria Island studying its stratigraphy and recognizing sequences of reoccupation events tied with eleven radiocarbon sample ages (14C ages). I developed a method to link patterns of reoccupation with sea-level proxies by iterating relative sea level curves for a range of uplift rates. I find the best fit between relative sea-level and the stratigraphic patterns for an uplift rate of 1.5 +- 0.3 m/ka. A Graphical User Interface named TerraceM® was developed in Matlab®. This novel software tool determines shoreline angles in wave-cut terraces under different geomorphic scenarios. To validate the methods, I select test sites in areas of available high-resolution LiDAR topography along the Maule earthquake rupture zone and in California, USA. The software allows determining the 3D location of the shoreline angle, which is a proxy for the estimation of permanent deformation rates. The method is based on linear interpolations to define the paleo platform and cliff on swath profiles. The shoreline angle is then located by intersecting these interpolations. The accuracy and precision of TerraceM® was tested by comparing its results with previous assessments, and through an experiment with students in a computer lab setting at the University of Potsdam. I combined the methods developed to analyze wave-built and wave-cut terraces to assess regional patterns of permanent deformation along the (2010) Maule earthquake rupture. Wave-built terraces are tied using 12 Infra Red Stimulated luminescence ages (IRSL ages) and shoreline angles in wave-cut terraces are estimated from 170 aligned swath profiles. The comparison of coseismic slip, interseismic coupling, and permanent deformation, leads to three areas of high permanent uplift, terrace warping, and sharp fault offsets. These three areas correlate with regions of high slip and low coupling, as well as with the spatial limit of at least eight historical megathrust ruptures (M8-9.5). I propose that the zones of upwarping at Arauco and Topocalma reflect changes in frictional properties of the megathrust, which result in discrete boundaries for the propagation of mega earthquakes. To explore the application of geomorphic markers and quantitative morphology in offshore areas I performed a local study of patterns of permanent deformation inferred from hitherto unrecognized drowned shorelines at the Arauco Bay, at the southern part of the (2010) Maule earthquake rupture zone. A multidisciplinary approach, including morphometry, sedimentology, paleontology, 3D morphoscopy, and a landscape Evolution Model is used to recognize, map, and assess local rates and patterns of permanent deformation in submarine environments. Permanent deformation patterns are then reproduced using elastic models to assess deformation rates of an active submarine splay fault defined as Santa Maria Fault System. The best fit suggests a reverse structure with a slip rate of 3.7 m/ka for the last 30 ka. The register of land level changes during the earthquake cycle at Santa Maria Island suggest that most of the deformation may be accrued through splay fault reactivation during mega earthquakes, like the (2010) Maule event. Considering a recurrence time of 150 to 200 years, as determined from historical and geological observations, slip between 0.3 and 0.7 m per event would be required to account for the 3.7 m/ka millennial slip rate. However, if the SMFS slips only every ~1000 years, representing a few megathrust earthquakes, then a slip of ~3.5 m per event would be required to account for the long- term rate. Such event would be equivalent to a magnitude ~6.7 earthquake capable to generate a local tsunami. The results of this thesis provide novel and fundamental information regarding the amount of permanent deformation accrued in the crust, and the mechanisms responsible for this accumulation at millennial time-scales along the M8.8 Maule earthquake (2010) rupture zone. Furthermore, the results of this thesis highlight the application of quantitative geomorphology and the use of repeatable methods to determine permanent deformation, improve the accuracy of marine terrace assessments, and estimates of vertical deformation rates in tectonically active coastal areas. This is vital information for adequate coastal-hazard assessments and to anticipate realistic earthquake and tsunami scenarios.},
  language  = {en}
}
@article{LangeTilmannBarrientosetal.2012,
  author    = {Lange, Dietrich and Tilmann, Frederik and Barrientos, Sergio E. and Contreras-Reyes, Eduardo and Methe, Pascal and Moreno, Marcos and Heit, Ben and Agurto, Hans and Bernard, Pascal and Vilotte, Jean-Pierre and Beck, Susan},
  title     = {Aftershock seismicity of the 27 February 2010 Mw 8.8 Maule earthquake rupture zone},
  series = {Earth \& planetary science letters},
  volume    = {317},
  journal   = {Earth \& planetary science letters},
  number    = {2},
  publisher = {Elsevier},
  address   = {Amsterdam},
  issn      = {0012-821X},
  doi       = {10.1016/j.epsl.2011.11.034},
  pages     = {413 -- 425},
  year      = {2012},
  abstract  = {On 27 February 2010 the M-w 8.8 Maule earthquake in Central Chile ruptured a seismic gap where significant strain had accumulated since 1835. Shortly after the mainshock a dense network of temporary seismic stations was installed along the whole rupture zone in order to capture the aftershock activity. Here, we present the aftershock distribution and first motion polarity focal mechanisms based on automatic detection algorithms and picking engines. By processing the seismic data between 15 March and 30 September 2010 from stations from IRIS, IPGP, GFZ and University of Liverpool we determined 20,205 hypocentres with magnitudes M-w between 1 and 5.5. Seismic activity occurs in six groups: 1.) Normal faulting outer rise events 2.) A shallow group of plate interface seismicity apparent at 25-35 km depth and 50-120 km distance to the trench with some variations between profiles. Along strike, the aftershocks occur largely within the zone of coseismic slip but extend similar to 50 km further north, and with predominantly shallowly dipping thrust mechanisms. Along dip, the events are either within the zone of coseismic slip, or downdip from it, depending on the coseismic slip model used. 3.) A third band of seismicity is observed further downdip at 40-50 km depth and further inland at 150-160 km trench perpendicular distance, with mostly shallow dipping (similar to 28 degrees) thrust focal mechanisms indicating rupture of the plate interface significantly downdip of the coseismic rupture, and presumably above the intersection of the continental Moho with the plate interface. 4.) A deep group of intermediate depth events between 80 and 120 km depth is present north of 36 degrees S. Within the Maule segment, a large portion of events during the inter-seismic phase originated from this depth range. 5.) The magmatic arc exhibits a small amount of crustal seismicity but does not appear to show significantly enhanced activity after the M-w 8.8 Maule 2010 earthquake. 6.) Pronounced crustal aftershock activity with mainly normal faulting mechanisms is found in the region of Pichilemu (similar to 34.5 degrees S). These crustal events occur in a similar to 30 km wide region with sharp inclined boundaries and oriented oblique to the trench. The best-located events describe a plane dipping to the southwest, consistent with one of the focal planes of the large normal-faulting aftershock (M-w = 6.9) on 11 March 2010.},
  language  = {en}
}
@article{KonradSchmolkeZackO'Brienetal.2011,
  author    = {Konrad-Schmolke, Matthias and Zack, Thomas and O'Brien, Patrick J. and Barth, Matthias},
  title     = {Fluid migration above a subducted slab - Thermodynamic and trace element modelling of fluid-rock interaction in partially overprinted eclogite-facies rocks (Sesia Zone, Western Alps)},
  series = {Earth \& planetary science letters},
  volume    = {311},
  journal   = {Earth \& planetary science letters},
  number    = {3-4},
  publisher = {Elsevier},
  address   = {Amsterdam},
  issn      = {0012-821X},
  doi       = {10.1016/j.epsl.2011.09.025},
  pages     = {287 -- 298},
  year      = {2011},
  abstract  = {The amount and composition of subduction zone fluids and the effect of fluid-rock interaction at a slab-mantle interface have been constrained by thermodynamic and trace element modelling of partially overprinted blueschist-facies rocks from the Sesia Zone (Western Alps). Deformation-induced differences in fluid flux led to a partial preservation of pristine mineral cores in weakly deformed samples that were used to quantify Li, B, Stand Pb distribution during mineral growth, -breakdown and modification induced by fluid-rock interaction. Our results show that Li and 13 budgets are fluid-controlled, thus acting as tracers for fluid-rock interaction processes, whereas Stand Pb budgets are mainly controlled by the fluid-induced formation of epidote. Our calculations show that fluid-rock interaction caused significant Li and B depletion in the affected rocks due to leaching effects, which in turn can lead to a drastic enrichment of these elements in the percolating fluid. Depending on available fluid-mineral trace element distribution coefficients modelled fluid rock ratios were up to 0.06 in weakly deformed samples and at least 0.5 to 4 in shear zone mylonites. These amounts lead to time integrated fluid fluxes of up to 1.4-10(2) m(3) m(-2) in the weakly deformed rocks and 1-8-10(3) m(3) m(-2) in the mylonites. Combined thermodynamic and trace element models can be used to quantify metamorphic fluid fluxes and the associated element transfer in complex, reacting rock systems and help to better understand commonly observed fluid-induced trace element trends in rocks and minerals from different geodynamic environments.},
  language  = {en}
}
@phdthesis{Lange2008,
  author    = {Lange, Dietrich},
  title     = {The South Chilean subduction zone between 41° and 43.5°S : seismicity, structure and state of stress},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-18948},
  school      = {Universit{\"a}t Potsdam},
  year      = {2008},
  abstract  = {Die st{\"a}rksten Erdbeben treten an Subduktionszonen auf, so z.B. das st{\"a}rkste instrumentell jemals gemessene Erdbeben vom 22. Mai 1960 mit einer Magnitude von 9,5 Mw in S{\"u}d Chile. In dieser Arbeit werden lokal gewonnene seismologische Daten aus dem zentralen Bereich des 1960er-Bebens vorgestellt. Das seismologische Netzwerk umfasste den chilenischen Forearc zwischen Tiefseegraben und den vulkanischen Bogen zwischen 41,5°-43,5°S und {\"u}berdeckte sowohl die Insel Chilo{\´e} als auch die Nord-S{\"u}d-streichende Liqui{\~n}e-Ofqui St{\"o}rungszone (LOFZ). Zwischen November 2004 und Oktober 2005 konnten 364 lokale Ereignisse registriert werden. Die gewonnen Aufzeichnungen erlauben Aussagen sowohl {\"u}ber das aktuelle Spannungsfeld im Forearc als auch {\"u}ber das lokale Geschwindigkeitsmodell und die Geometrie der subduzierten Benioff-Zone. Mit einer Auswahl von P- und S-Laufzeiten von gut lokalisierbaren Erdbeben wurden ein Minimum 1-D Geschwindigkeitsmodell, Stationsresiduen und die Hypozentralparameter invertiert. Dieses Geschwindigkeitsmodell diente als Startmodell f{\"u}r die 2-D Tomographie. Das 2-D vp-Modell zeigt eine Zone erh{\"o}hter Geschwindigkeiten unterhalb des L{\"a}ngstals und des {\"o}stlichen Bereiches der Insel Chilo{\´e}, die als Mantelaufw{\"o}lbung interpretiert werden kann. Die Benioff-Zone wird als eine mit ca. 30° ostw{\"a}rts einfallende Struktur abgebildet. Die seismische Hauptaktivit{\"a}t findet parallel zur K{\"u}ste der Insel Chilo{\´e} in Tiefen zwischen 12 und 30 km statt; es handelt sich um Beben, die wahrscheinlich auf der Plattengrenzfl{\"a}che stattfinden. In Tiefen {\"u}ber 70 km l{\"a}sst die Seismizit{\"a}t bereits stark nach, die tiefsten Beben wurden in 120 km Tiefe registriert. Die Abwesenheit tieferer Seismizit{\"a}t wird auf das junge Alter (und eine damit verbundene hohe Temperatur) der ozeanischen Platte zur{\"u}ckgef{\"u}hrt. Neben der Seismizit{\"a}t in der Benioff-Zone treten flache, krustale Beben in verschiedenen H{\"a}ufungen entlang des magmatischen Bogens auf. Diese Bereiche erh{\"o}hter Seismizit{\"a}t sind r{\"a}umlich mit der LOFZ und den Vulkanen Chait{\´e}n, Michinmahuida und Corcovado verkn{\"u}pft. Beben bis zu einer Magnitude von 3,8 Mw zeigen die gegenw{\"a}rtige Aktivit{\"a}t der LOFZ. Herdfl{\"a}chen entlang der LOFZ wurden aus Momententensor-Inversion anhand von Amplitudenspektren von Raumwellen gewonnen. Ergebnisse einer Spannungsfeldinversion zeigen ein Blattverschiebungsregime f{\"u}r den magmatischen Bogen und ein {\"U}berschiebungsregime f{\"u}r Beben in der Benioff-Zone auf. Die hier gemachten seismologischen Beobachtungen, zusammen mit teleseismischen Erdbeben und geologischen Befunden, unterst{\"u}tzen die Modellvorstellung eines sich nordw{\"a}rts bewegenden kontinentalen Forearc-Blocks f{\"u}r S{\"u}d Chile.},
  language  = {en}
}
@phdthesis{Pilz2008,
  author    = {Pilz, Peter},
  title     = {Ein neues magmatisch-tektonisches Modell zur Asthenosph{\"a}rendynamik im Bereich der zentralandinen Subduktionszone S{\"u}damerikas},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-20206},
  school      = {Universit{\"a}t Potsdam},
  year      = {2008},
  abstract  = {Im Rahmen der Dissertation wurden an W{\"a}ssern und freien Gasen aus Thermalquellen sowie an weniger als 5 Millionen Jahre alten basischen Vulkaniten des zentralandinen Puna-Hochplateaus (NE-Argentinien) umfangreiche element- und isotopengeochemische Untersuchungen durchgef{\"u}hrt und die Edelgasgehalte und -isotopensignaturen in diesen Medien bestimmt. Damit soll ein Beitrag zum besseren Verst{\"a}ndnis der j{\"u}ngeren Subduktionsgeschichte im Bereich der s{\"u}dlichen Zentralanden geleistet, die Wechselwirkungen zwischen ozeanischer Unter- und kontinentaler Oberplatte sichtbar gemacht und die Edelgassystematik verbessert werden. Wie die Ergebnisse der Untersuchungen an Gasen aus den Thermalquellen der Puna-Region zeigen, ist der Anteil an Mantel-Helium in den Thermalquellen dieser Region mit bis zu 67 \% wesentlich h{\"o}her als in der westlich gelegenen vulkanisch aktiven Westkordillere und den anderen angrenzenden Gebieten. In einigen Quellen konnten sogar Anteile an Mantel-Neon nachgewiesen werden, was aufgrund von {\"U}berlagerungen mit Neon atmosph{\"a}rischen und krustalen Ursprungs weltweit bisher nur vereinzelt gelungen ist. F{\"u}r kontinentale Bereiche mit großer Krustendicke ist ein solch starker Mantelgasfluss {\"a}ußerst ungew{\"o}hnlich und bedeutet, dass Mantelschmelzen bis in die Kruste aufgedrungen sind und tief reichende Wegsamkeiten existieren, so dass die Mantelgase aufsteigen k{\"o}nnen, ohne stark krustal beeinflusst zu werden. Dass im Bereich der Puna rezent Mantelmaterial in die Kruste aufsteigt, zu diesem Ergebnis kommen auch aktuelle seismologische Untersuchungen. Zudem wurden junge, vorwiegend monogenetische Basalte bis basaltische Andesite geochemisch auf ihre Haupt-, Neben- und Spurenbestandteile sowie ihre Gehalte an Seltenenerdenelementen hin untersucht. Auch wurden die Isotopenverh{\"a}ltnisse von Sr, Nd und Pb in den Gesteinen bestimmt und petrographisch-mineralogische Analysen der darin enthaltenen Olivine und Pyroxene durchgef{\"u}hrt. Wie die Resultate belegen, haben die Magmen bei ihrem Aufstieg durch die Erdkruste insbesondere Material aus der Oberkruste assimiliert und sind zudem durch Fluide aus der abtauchenden Platte beeinflusst worden. Damit konnte gezeigt werden, dass einfache geochemische Methoden allein nicht ausreichen, um die Mantelquelle der Magmen ermitteln oder Aussagen {\"u}ber die Asthenosph{\"a}rendynamik in der Region machen zu k{\"o}nnen. Im Gegensatz dazu zeigen die Messungen der Edelgasisotopenverh{\"a}ltnisse in den Fluideinschl{\"u}ssen der Olivine und Pyroxene, dass deren Edelgaszusammensetzung nicht durch Krustenkontamination beeinflusst wurde, weil die Magmen erst nach der Olivin- bzw. Pyroxen-Kristallisation Schmelzen aus der Oberkruste assimiliert haben. Dar{\"u}ber hinaus konnten durch die Edelgasisotopenmessungen die bisher h{\"o}chsten magmatischen He- und Ne-Isotopenverh{\"a}ltnisse von ganz S{\"u}damerika nachgewiesen werden. Aus der unterschiedlichen H{\"o}he der Messwerte ist zu schließen, dass die im Osten der Puna vorkommenden {\"a}lteren Laven aus einem nichtkonvektiven (lithosph{\"a}rischen) Mantel stammen, w{\"a}hrend die am vulkanischen Bogen und Westrand der Puna gelegenen j{\"u}ngeren Laven, ihren Ursprung in einer konvektiven (asthenosph{\"a}rischen) Mantelquelle haben. Zudem konnte gezeigt werden, dass der Mantelgasfluss in der Region in den letzten 5 Millionen Jahren stark zunahm und sich die Eruption von mantelst{\"a}mmigen basischen Laven in dieser Zeit kontinuierlich in westliche Richtung zum aktiven Vulkanbogen hin verlagerte. Im daraus abgeleiteten Modell beruht dieser Prozess (1) auf einer an die Kontinentalverschiebung gekoppelten W-Drift des Kontinents und (2) auf einem mit der Versteilung der Unterplatte verbundenen Vordringen des subkontinentalen asthenosph{\"a}rischen Mantels nach W, nach dem Ende der Subduktion des unterseeischen aseismischen Juan Fern{\´a}ndez-R{\"u}ckens in der Region. Zudem gibt es starke Argumente daf{\"u}r, dass die asthenosph{\"a}rischen Magmen aus einer fluidreichen Zone in 500 - 600 km Tiefe parallel zur subduzierten Platte aufsteigen und nicht, wie bisher angenommen, durch Schmelzbildung in Bereichen unter 200 km Tiefe, allein durch Entw{\"a}sserung der abtauchenden Platte erzeugt werden. Zu diesem Resultat f{\"u}hrt vor allem die Kombination der He-Isotopenverh{\"a}ltnisse mit Ergebnissen seismologischer Untersuchungen.},
  language  = {de}
}
@phdthesis{Kellner2007,
  author    = {Kellner, Antje},
  title     = {Different styles of deformation of the fore-arc wedge along the Chilean convergent margin : insights from 3D numerical experiments},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-15898},
  school      = {Universit{\"a}t Potsdam},
  year      = {2007},
  abstract  = {The styles of deformation of the fore-arc wedges along the Chilean convergent margin are observed to vary significantly, despite similar plate kinematic conditions. Here, I focus on the analysis of fore-arc deformation on two regions along the Chilean convergent margin at 20°-24°S and 37°-42°S. Although both regions are subjected to the oblique subduction of the oceanic Nazca plate and backstopped by the Andes mountain chain; they display different patterns of deformation. The northern Chilean study area (20° - 24°S) is characterized by an exceptionally thick crust of about 60 km beneath the Altiplano - Puna plateau, lack of an accretionary wedge in the fore-arc due to hyperarid climate, and consequently a sediment starved trench. Two major margin parallel strike slip faults are observed in this area, the Atacama Fault Zone (AFZ) and the Precordilleran Fault System (PFS). Both strike-slip faults do not exhibit significant recent displacement. The southern study area (37° - 42°S), compared to the northern study area, is characterized by lower topography, high precipitation rates (~2000 mm/yr), and a younger subducted oceanic plate. An active strike-slip fault, the Liqui{\~n}e-Ofqui-Fault-Zone (LOFZ), shows ~1 cm/yr recent dextral movement and shapes the surface of this area. Thus, the southern Chilean study area exhibits localized strike-slip motion. Within this area the largest earthquake ever recorded, the 1960 Valdivia earthquake, occurred with a moment magnitude of MW=9.5. I have constructed 2D thermal models and 3D mechanical models for both Chilean study areas to study processes related to active subduction. The applied numerical method is the finite element technique by means of the commercial software package ABAQUS. The thermal models are focused on the thermal conditions along the plate interface. The thermal structure along the plate interface reveals the limits of coupling but also the type of transition from coupled to uncoupled and vice versa. The model results show that shear heating at the plate interface is an important mechanism that should be taken into account. The models also show that the thermal condition at the downdip limit of the coupling zone leads to a sharp decrease of friction along the interface. Due to the different geometries of the two Chilean study areas, such as the slab dip and the thickness of the continental crust, the downdip limit of the southern study area is slightly shallower than that of the northern study area. The results of the 2D thermal models are used to constrain the spatial extent of the coupling zone in the 3D mechanical models. 3D numerical simulations are used to investigate how geometry, rheology and mechanical parameters influence strain partitioning and styles of deformation in the Chilean fore-arc. The general outline of the models is based on the fore-arc geometry and boundary conditions as derived from geophysical and geological field data. I examined the influence of different rheological approaches and varying physical properties of the fore-arc to identify and constrain the parameters controlling the difference in surface deformation between the northern and southern study area. The results of numerical studies demonstrate that a small slab dip, a high coefficient of basal friction, a high obliquity of convergence, and a high Young's modulus favour localisation of deformation in the fore-arc wedge. This parameter study helped me to constrain preferred models for the two Chilean study areas that fit to first order observations. These preferred models explain the difference in styles of deformation as controlled by the angle of obliquity, the dip of subducting slab, and the strength of wedge material. The difference in styles can be even larger if I apply stronger coupling between plates within the southern area; however, several independent observations indicate opposite tendency showing southward decrease of intensity of coupling. The weaker wedge material of the preferred model for the northern study area is associated with advanced development of the adjacent orogen, the Central Andes. Analysis of world-wide examples of oblique subduction zones supports the conclusion that more mature subduction zones demonstrate less pronounced localization of strike-slip motion.},
  language  = {en}
}