An official website of the United States government
Abstract The thermal structure of subduction zones is fundamental to our understanding of the physical and chemical processes that occur at active convergent plate margins. These include magma generation and related arc volcanism, shallow and deep seismicity, and metamorphic reactions that can release fluids. Computational models can predict the thermal structure to great numerical precision when models are fully described but this does not guarantee accuracy or applicability. In a trio of companion papers, the construction of thermal subduction zone models, their use in subduction zone studies, and their link to geophysical and geochemical observations are explored. In this part II, the finite element techniques that can be used to predict thermal structure are discussed in an introductory fashion along with their verification and validation. Steady-state thermal structure for the updated subduction zone benchmark. a) Temperature predicted by TF for case 1; b) temperature difference between TF and Sepran using the penalty function (PF) method for case 1 at fm=1 where fmrepresents the smallest element sizes in the finite element grids near the coupling point; c) slab top temperature comparison for case 1; and d)–f) as a)–c) but now for case 2. The star indicates the position or temperature conditions at the coupling point.
more »
« less
Thermal modeling of subduction zones with prescribed and evolving 2D and 3D slab geometries
Abstract The determination of the temperature in and above the slab in subduction zones, using models where the top of the slab is precisely known, is important to test hypotheses regarding the causes of arc volcanism and intermediate-depth seismicity. While 2D and 3D models can predict the thermal structure with high precision for fixed slab geometries, a number of regions are characterized by relatively large geometrical changes over time. Examples include the flat slab segments in South America that evolved from more steeply dipping geometries to the present day flat slab geometry. We devise, implement, and test a numerical approach to model the thermal evolution of a subduction zone with prescribed changes in slab geometry over time. Our numerical model approximates the subduction zone geometry by employing time dependent deformation of a Bézier spline that is used as the slab interface in a finite element discretization of the Stokes and heat equations. We implement the numerical model using the FEniCS open source finite element suite and describe the means by which we compute approximations of the subduction zone velocity, temperature, and pressure fields. We compute and compare the 3D time evolving numerical model with its 2D analogy at cross-sections for slabs that evolve to the present-day structure of a flat segment of the subducting Nazca plate.
more »
« less
- Award ID(s):
- 2021027
- PAR ID:
- 10496584
- Publisher / Repository:
- Springer Science + Business Media
- Date Published:
- Journal Name:
- Progress in Earth and Planetary Science
- Volume:
- 11
- Issue:
- 1
- ISSN:
- 2197-4284
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The Tonga‐Kermadec subduction zone exhibits the fastest observed trench retreat and convergence near its northern end. However, a paradox exists: despite the rapid trench retreat, the Tonga slab maintains a relatively steep dip angle above 400 km depth. The slab turns flat around 400 km, then steepening again until encountering a stagnant segment near 670 km. Despite its significance for understanding slab dynamics, no existing numerical model has successfully demonstrated how such a distinct slab morphology can be generated under the fast convergence. Here we run subduction models that successfully reproduce the slab geometries while incorporating the observed subduction rate. We use a hybrid velocity boundary condition, imposing velocities on the arc and subducting plate while allowing the overriding plate to respond freely. This approach is crucial for achieving a good match between the modeled and observed Tonga slab. The results explain how the detailed slab structure is highly sensitive to physical parameters including the seafloor age and the mantle viscosity. Notably, a nonlinear rheology, where dislocation creep reduces upper mantle viscosity under strong mantle flow, is essential. The weakened upper mantle allows for a faster slab sinking rate, which explains the large dip angle. Our findings highlight the utilizing rheological parameters that lead to extreme viscosity variations within numerical models to achieve an accurate representation of complex subduction systems like the Tonga‐Kermadec zone. Our study opens new avenues for further study of ocean‐ocean subduction systems, advancing our understanding of their role in shaping regional and global tectonics.more » « less
-
Abstract Deep earthquakes require the cold temperatures found in sinking lithosphere to store elastic strain. It has also been proposed that sufficiently high rates of deformation are also required, regardless of the failure mechanism. However, this strain‐rate hypothesis is based on generic time‐dependent and visco‐plastic subduction models, positing a challenge for direct comparisons to present‐day earthquake observations. Here, we present a new numerical modeling approach incorporating location‐specific visco‐elasto‐plastic models to facilitate direct comparison with deep earthquake observations. We present a Proof‐of‐Concept Model using a 2D synthetic slab to demonstrate that this novel approach can reproduce stress and strain‐rate patterns and the stress orientations from a fully time‐dependent model. Applying this method to a 2D profile through the Tonga‐Kermadec subduction zone we find that variations in strain‐rate correlate with seismicity rate and regions of stress in the slab exceeding 500 MPa. Elasticity in the slab leads to formation of a clearly defined neutral plane extending into the transition zone and creating a narrow region of down‐dip compression along the top portion of the slab which broadens across the full width of the slab only within the deep transition zone. In addition, assuming that the strain‐rate hypothesis is correct, we show that peaks in strain‐rate, which are associated with bends in the slab, could be used to constrain the slab shape beyond the envelope of seismicity.more » « less
-
Abstract A global study of subduction zone dynamics indicates that the thermal structure of the overriding plate may control arc location. A fast convergence rate and a steep slab dip bring a hotter mantle further into the wedge corner, forming arc volcanoes closer to the trench. Separately, laboratory and numerical experiments showed that the development of a back‐arc spreading center (BASC) is driven by the migration of the subducting hinge, especially following changes in the slab geometry. As both arc location and the deformation regime of the overriding plate depend on slab kinematics and geometry, we investigate the possible correlations between BASC, the position of volcanic arcs, and slab dip at the scale of individual subduction zones. To do this, we compare the distance from trench to arc and trench to BASC at the Mariana, Scotia, Vanuatu, Tonga, and Kermadec subduction zones. In most cases, the arc and BASC are closer to the trench when the slab is dipping steeply. The correlation could result from an interplay between progressive changes in slab geometry and overriding plate deformation. This assumes, on the one hand, that the isotherm at the apex of which the arc forms is tied to a constant slab decoupling depth and, on the other hand, that back‐arc opening accommodates a change in slab dip. As slab dip decreases, both the BASC and the apex of the isotherm controlling the melt focusing move further from the trench. The observed trends are consistent with a slab anchored at 660 km depth.more » « less
-
Abstract Low‐angle subduction has been shown to have a profound impact on subduction processes. However, the mechanisms that initiate, drive, and sustain flat‐slab subduction are debated. Within all subduction zone systems, metamorphic dehydration reactions within the down‐going slab have been hypothesized to produce seismicity, and to produce water that fluxes melting of the asthenospheric wedge leading to arc magmatism. In this work, we examine the role hydration plays in influencing slab buoyancy and the geometry of the downgoing oceanic plate. When water is introduced to the oceanic lithosphere, it is incorporated into hydrous phases, which results in lowered rock densities. The net effect of this process is an increase in the buoyancy of the downgoing oceanic lithosphere. To better understand the role of water in low‐angle subduction settings, we model flat‐slab subduction in Alaska, where the thickened oceanic lithosphere of the Yakutat oceanic plateau is subducting beneath the continental lithosphere. In this work, we calculate the thermal conditions and stable mineral assemblages in the slab crust and mantle in order to assess the role that water plays in altering the density of the subducting slab. Our slab density results show that a moderate amount of hydration (1–1.5 wt% H2O) in the subducting crust and upper lithospheric mantle reduces slab density by 0.5%–0.8% relative to an anhydrous slab, and is sufficient to maintain slab buoyancy to 300–400 km from the trench. These models show that water is a viable factor in influencing the subduction geometry in Alaska, and is likely important globally.more » « less
