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