Three-dimensional temperature variations in a fossil subduction zone resolved by RSCM thermometry (Tauern Window, Eastern Alps)

Author(s):  
Philip Groß ◽  
Jan Pleuger ◽  
Mark R. Handy ◽  
Timm John

<p>Knowledge on the thermal state of orogens and subduction zones is crucial in trying to understand the processes that take place in these zones, since temperature controls, e.g., rock strength, metamorphic reactions and fluid flow. These are all critical parameters for the dynamics of orogens and subduction zones and conversely, these parameters feed back on the thermal state in various ways. We investigated an example of a former subduction zone, exposed in the central Tauern Window (Eastern Alps), with the aim of reconstructing its three-dimensional temperature variations.</p><p>Structural and petrological observations in the central Tauern Window reveal a tens-of-kilometre-scale sheath fold that formed under high-pressure (HP) conditions (ca. 2 GPa). The fold is a composite structure that isoclinally folded the thrust of an oceanic nappe derived from Alpine Tethys onto a unit of the distal European continental margin, also affected by HP conditions. This structural assemblage is preserved between two younger domes at either end of the Tauern Window. The domes are associated with temperature-dominated Barrow-type metamorphism that overprints the HP-metamorphism partly preserved in the sheath fold.</p><p>Using Raman spectroscopy on carbonaceous material (RSCM) on 100 samples from this area, we were able to distinguish domains with the original, subduction-related peak temperature conditions from domains that were overprinted during later temperature-dominated (Barrovian) metamorphism. The distribution of RSCM-temperatures in the Barrovian domains indicates a decrease of peak temperature with increasing distance from the centres of the thermal domes, both in map view and cross section. This represents a geotherm where paleo-temperature increases downward, in line with previous studies using, e.g., oxygen isotope fractionation and calcite-dolomite equilibria. However, we observe the opposite temperature trend in the lower limb of the sheath fold, viz., tendentially an upward increase in paleo-temperature. We interpret this inverted temperature domain as the relic of a subduction-related temperature field. Towards the central part of the sheath fold’s upper limb, measured temperatures increase to a maximum of ca. 520°C. Further upsection in the hanging wall of the sheath fold, temperatures decrease to where they are indistinguishable from the peak temperatures of the overprinting Barrovian metamorphism. Isograds (i.e. contours of equal peak-temperature) are oriented roughly parallel to the nappe contacts and lithological layering, which results in an eye-shaped concentric isograd pattern in cross-section. This reflects a sheath-like three-dimensional geometry of the isograds. We propose the following hypothesis to explain the subduction-related peak-temperature pattern: The pattern reflects sheath folding of a subduction-related temperature field. Possibly, sheath folding occurred during exhumation, after the equilibration at peak pressure and temperature conditions. The preservation of the pattern implies fast exhumation and cooling of the rocks.</p>

2019 ◽  
Vol 5 (12) ◽  
pp. eaax6720 ◽  
Author(s):  
Jonathan R. Weiss ◽  
Qiang Qiu ◽  
Sylvain Barbot ◽  
Tim J. Wright ◽  
James H. Foster ◽  
...  

Deformation associated with plate convergence at subduction zones is accommodated by a complex system involving fault slip and viscoelastic flow. These processes have proven difficult to disentangle. The 2010 Mw 8.8 Maule earthquake occurred close to the Chilean coast within a dense network of continuously recording Global Positioning System stations, which provide a comprehensive history of surface strain. We use these data to assemble a detailed picture of a structurally controlled megathrust fault frictional patchwork and the three-dimensional rheological and time-dependent viscosity structure of the lower crust and upper mantle, all of which control the relative importance of afterslip and viscoelastic relaxation during postseismic deformation. These results enhance our understanding of subduction dynamics including the interplay of localized and distributed deformation during the subduction zone earthquake cycle.


2020 ◽  
Author(s):  
Nikolaus Froitzheim

<p>The Adula Nappe in the Central Alps and the Pohorje Nappe in the Eastern Alps are among the highest-pressure metamorphic complexes in the Alps. In both cases, Variscan continental crust containing post-Variscan intrusions was subducted, during the Cenomanian-Turonian in the case of Pohorje and during the Eocene in the case of Adula.</p><p>The Pohorje Nappe is exceptional in that ultrahigh pressures of 3.0 to 4.0 GPa are recorded by different rocks contrasting in rheology: competent lenses of kyanite eclogite and garnet peridotite as well as the surrounding incompetent matrix of diamond-bearing paragneiss. If pressure had been strongly non-lithostatic, rheologically different rock types would be expected to record different pressures. This is not the case, which rather suggests near-lithostatic pressure and, consequently, subduction to >100 km depth. Lu-Hf ages for UHP metamorphism in eclogite and garnet peridotite are similar (c. 96–92 Ma). Paragneiss yielded Permian to Triassic zircon cores and Cretaceous (c. 92 Ma) rims grown during UHP metamorphism. Hence, the rocks were subducted and exhumed together as a coherent, although strongly deformed unit.</p><p>The Adula Nappe originated from the southern passive continental margin of Europe. It was buried in and exhumed from a south-dipping subduction zone after Europe-Adria continent collision. Previous interpretations as a tectonic mélange were based on the mixture of gneiss with eclogite and garnet peridotite lenses. However, the eclogites also record an older (Variscan) metamorphism and thus do not represent Mesozoic oceanic crust but pre-Alpine continental basement, just like the gneisses. The Alpine subduction culminated around 37 Ma. Alpine metamorphic pressures show a strong gradient from c. 1.2 GPa at the front of the nappe in the North to >3 GPa in garnet peridotite and eclogite in the southernmost part (e.g. Alpe Arami), over a north-south distance of only c. 40 km. In contrast to Pohorje, indications of UHP metamorphism have not yet been found in the gneissic matrix surrounding eclogite and peridotite. During exhumation, the nappe was intensely sheared and folded but stayed coherent and did not mix with the surrounding units.  The exhumation of the Adula from deep in the subduction zone is recorded by mylonitic shearing in the gneissic matrix. Structures, strain, and textures indicate strongly three-dimensional, non-plane-strain flow. Differential loading, not buoyancy, is proposed to have caused the exhumation.</p><p>The main results from these two case studies are: (1) Subduction of continental crust to mantle depth is real and not a misinterpretation of non-lithostatic pressure; (2) not all subducted units are mélanges but some stay coherent during subduction and exhumation.</p>


2006 ◽  
Vol 73 (6) ◽  
pp. 984-994 ◽  
Author(s):  
Javier Quinteros ◽  
Pablo M. Jacovkis ◽  
Victor A. Ramos

The uplift and evolution of a noncollisional orogen developed along a subduction zone, such as the Andean system, is a direct consequence of the interrelation between plate tectonic stresses and erosion. Tectonic stresses are related to the convergence velocity and thermal state, among other causes. In this paper, a new model designed to investigate the evolution of the topography and the upper crustal deformation of noncollisional orogens in a subduction zone produced by the oceanic crust being subducted is presented. The mechanical behavior of the crust was modeled by means of finite elements methods to solve Stokes equations for a strain-rate-dependent viscoplastic rheology. The model takes into account erosion effects using interface-tracking methods to assign fictitious properties to nonmaterial elements.


Science ◽  
2018 ◽  
Vol 362 (6410) ◽  
pp. 58-61 ◽  
Author(s):  
Gavin P. Hayes ◽  
Ginevra L. Moore ◽  
Daniel E. Portner ◽  
Mike Hearne ◽  
Hanna Flamme ◽  
...  

Subduction zones are home to the most seismically active faults on the planet. The shallow megathrust interfaces of subduction zones host Earth’s largest earthquakes and are likely the only faults capable of magnitude 9+ ruptures. Despite these facts, our knowledge of subduction zone geometry—which likely plays a key role in determining the spatial extent and ultimately the size of subduction zone earthquakes—is incomplete. We calculated the three-dimensional geometries of all seismically active global subduction zones. The resulting model, called Slab2, provides a uniform geometrical analysis of all currently subducting slabs.


2012 ◽  
Vol 4 (2) ◽  
pp. 919-941 ◽  
Author(s):  
A. K. Bengtson ◽  
P. E. van Keken

Abstract. Quantifying the precise thermal structure of subduction zones is essential for understanding the nature of metamorphic dehydration reactions, arc volcanism, and intermediate depth seismicity. High resolution two-dimensional (2-D) models have shown that the rheology of the mantle wedge plays a critical role and establishes strong temperature gradients in the slab. The influence of three-dimensional (3-D) subduction zone geometry on thermal structure is however not yet well characterized. A common assumption for 2-D models is that the cross-section is taken normal to the strike of the trench with a corresponding velocity reduction in the case of oblique subduction, rather than taken parallel to velocity. A comparison between a full 3-D Cartesian model with oblique subduction and selected 2-D cross-sections demonstrates that the trench-normal cross-section provides a better reproduction of the slab thermal structure than the velocity-parallel cross-section. An exception is found in the case of strongly curved subduction, such as in the Marianas, where strong 3-D flow in the mantle wedge is generated. In this case it is shown that the full 3-D model should be evaluated for an accurate prediction of the slab thermal structure.


Tectonics ◽  
2020 ◽  
Vol 39 (2) ◽  
Author(s):  
Philip Groß ◽  
Mark R. Handy ◽  
Timm John ◽  
Gerhard Pestal ◽  
Jan Pleuger

2013 ◽  
Vol 50 (2) ◽  
pp. 33-41
Author(s):  
V. Barkans ◽  
P. Shipkovs ◽  
M. Vanags

In the work, the temperature field model is developed for the absorber of a round-pipe collector. As distinguished from previous models when the temperature of liquid was assumed to be constant over the entire pipe crosssection, the results obtained clearly show the temperature variations in the absorber’s cross-section. In the work, optimal values are found in the work for geometrical parameters of the collector (i.e. the plate thickness and the pipe diameter) that allow the highest possible temperature of liquid to be achieved.


2021 ◽  
Vol 7 (4) ◽  
pp. eabc9620 ◽  
Author(s):  
Zewei Wang ◽  
Dapeng Zhao

How mantle materials flow and how intraslab fabrics align in subduction zones are two essential issues for clarifying material recycling between Earth’s interior and surface. Investigating seismic anisotropy is one of a few viable technologies that can directly answer these questions. However, the detailed anisotropic structure of subduction zones is still unclear. Under a general hexagonal symmetry anisotropy assumption, we develop a tomographic method to determine a high-resolution three-dimensional (3D) P wave anisotropic model of the Japan subduction zone by inverting 1,184,018 travel time data of local and teleseismic events. As a result, the 3D anisotropic structure in and around the dipping Pacific slab is firstly revealed. Our results show that slab deformation plays an important role in both mantle flow and intraslab fabric, and the widely observed trench-parallel anisotropy in the forearc is related to the intraslab deformation during the outer-rise yielding of the subducting plate.


Solid Earth ◽  
2012 ◽  
Vol 3 (2) ◽  
pp. 365-373 ◽  
Author(s):  
A. K. Bengtson ◽  
P. E. van Keken

Abstract. Quantifying the precise thermal structure of subduction zones is essential for understanding the nature of metamorphic dehydration reactions, arc volcanism, and intermediate depth seismicity. High resolution two-dimensional (2-D) models have shown that the rheology of the mantle wedge plays a critical role and establishes strong temperature gradients in the slab. The influence of three-dimensional (3-D) subduction zone geometry on thermal structure is however not yet well characterized. A common assumption for 2-D models is that the cross-section is taken normal to the strike of the trench with a corresponding velocity reduction in the case of oblique subduction, rather than taken parallel to velocity. A comparison between a full 3-D Cartesian model with oblique subduction and selected 2-D cross-sections demonstrates that the trench-normal cross-section provides a better reproduction of the slab thermal structure than the velocity-parallel cross-section. An exception is found in the case of a strongly curved trench, such as in the Marianas, where strong 3-D flow in the mantle wedge is generated. In this case it is shown that the full 3-D model should be evaluated for an accurate prediction of the slab thermal structure. The models demonstrate that the use of a dynamic slab and wedge, separated by a kinematic boundary, yields good results for describing slab velocities in 3-D.


Sign in / Sign up

Export Citation Format

Share Document