@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} } @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{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} } @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{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} } @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} } @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} }