Patched removal of the mantle lithosphere under orogens: A systematic numerical study

Abstract

Delamination or convective thinning could cause large-scale and complete removal of the mantle lithosphere under orogens. However, geological and geophysical observations suggest that patched removal of the mantle lithosphere has occurred in some orogens, such as the northeastern Tibetan Plateau, the central Tianshan, and the central Andes. Dislocation-creep-induced strain localization cannot promote effective removal of the mantle lithosphere to the Moho on a small-scale. Recent rheological studies propose that dislocation-accommodated grain boundary sliding (DisGBS) may dominate upper mantle deformation. DisGBS could make the lower lithospheric mantle rheologically weaker than dry olivine. With 2-D high-resolution thermo-mechanical modeling, we systematically investigated the conditions for the initiation of small-scale lithospheric thinning under orogens and explored the minimum range of removal of the mantle lithosphere. The numerical results indicate that classic convective drip cannot effectively thin the mantle lithosphere to the Moho on a small-scale. In contrast, small-scale thinning can be induced by lithospheric heterogeneity with DisGBS and plasticity. The rheological heterogeneity can be verified by magmatism and metasomatism under the central Andes and orogens between terranes under the northeastern Tibetan Plateau or in Tianshan.

1 Introduction

Lithospheric mantle removal under orogens (e.g., the Himalayan–Tibetan Plateau, Central Anatolia, and Colorado) is often explained by delamination or convective thinning (Bird 1978; Bird 1979; Owens & Zandt 1997; Molnar et al., 1998; Gogus et al., 2017). Both delamination and convective thinning are driven by the negative buoyancy in the lithospheric mantle (and sometimes the eclogitized lower crust), even though the growth period of instability and the rheological control of them are different (Levander et al., 2011; Beall et al., 2017; Lei et al., 2019). The long wavelength of downwelling thickened lithosphere leads to mantle lithosphere thinning, and the range of the initial wavelength is commonly larger than the thickness of the mantle lithosphere (Conrad & Molnar 1997; Houseman & Molnar 1997). A laterally homogeneous rheological structure of the lithosphere with dislocation creep or diffusion creep is usually assumed in previous models of delamination or convective thinning, which predict large-scale peeling off of the mantle lithosphere (Bird, 1978; Bird 1979; Schott & Schmeling 1998; Morency & Doin 2004; Göğüş & Pysklywec 2008a; Göğüş & Pysklywec 2008b; Bajolet et al., 2012; Wang & Currie 2015; Beall et al., 2017; Lei et al., 2019). However, patched removal of the mantle lithosphere, with insignificant mafic magmatism, has been observed in the interior of collisional plates and the margin of subducting plates, such as the northeastern Tibetan Plateau, the Tianshan, and the central Andes (Beck & Zandt 2002; Schurr et al., 2006; Ducea 2011; Ducea et al., 2013; Li et al., 2022; Zhang et al., 2022). The width of detachment approaches the layer thickness of the mantle lithosphere under these orogens (∼100–200 km). Such patched removal of the mantle lithosphere under orogens cannot be readily explained by existing models of delamination or convective thinning because the convective dripping resulting from the short wavelength of perturbation tends to occur near the bottom of the mantle lithosphere and cannot remove the entire mantle lithosphere to the Moho (Conrad & Molnar 1997; Houseman & Molnar 1997; Currie et al., 2008; Lorinczi and Houseman, 2009; Ducea 2011; Beall et al., 2017; Liao et al., 2017).

Rheological weakening is the key factor for small-scale mantle lithosphere removal under orogens, that is, removal from the bottom of the mantle lithosphere to the Moho. However, dislocation creep and diffusion creep overestimate the upper mantle’s viscosity (van Hunen et al., 2005; Billen, 2008), which may not effectively promote the patched removal of the mantle lithosphere (Currie et al., 2008; Liao et al., 2017). Recent rheological studies indicate that dislocation-accommodated grain boundary sliding (DisGBS), which differs from dislocation creep or diffusion creep, may dominate the deformation of olivine in the Earth’s upper mantle (Tomohiro and Takaaki, 2015). The estimated viscosity of the upper mantle controlled by DisGBS is independent of depth and is lower than that of dislocation creep, indicating that the local mantle lithosphere may experience significant deformation via strain localization or hydration activity.

The mantle lithosphere in the northeastern Tibetan Plateau, the Tianshan, and the Central Andes manifests obvious rheological heterogeneities. Geological and geophysical observations indicate that the northeastern Tibetan Plateau consists of multiple terranes with contrasting lithologies, thermal states, and rheological structures. Patches of seismically fast and slow anomalies have been observed beneath the northeastern Tibetan Plateau (Deng et al., 2018). For example, the mantle lithosphere in the Songpan–Ganzi and southern Kunlun–Qaidam terranes has probably been removed, with traces of delaminated pieces. The entire mantle lithosphere in North Qilian has been removed, together with the base of the crust, leading to the sudden decrease of the Moho depth (Deng et al., 2018). This interpretation is consistent with the rapid uplift of the Qilian block at the Miocene-Quaternary from recent magnetostratigraphy and tectonosedimentology (Fang et al., 2013). Young volcanic rocks are exposed (<6 Myr) in the western North Qilian Mountains (Xia et al., 2011). Therefore, lithospheric heterogeneities may play a key role in the shortening of the northeastern Tibetan Plateau (Wang et al., 2008; Yin and Dang, 2008; Zhang and Wang, 2014; Deng et al., 2018). Compared with the Tarim Craton and Kazakh lithosphere (Li et al., 2022), the crust and mantle lithosphere of the sandwiched central Tianshan are assumed to be much rheologically weaker. The intracontinental deformation of the central Tianshan is controlled by rheological heterogeneity (Huangfu et al., 2021; Li et al., 2022). In contrast, the mantle lithosphere under the central Andes does not involve multiple terranes and sutures, thus differing from the northeastern TP and central Tianshan. The tectonic history of subduction has induced the migration of aqueous fluids from the subducting Nazca plate and may also lead to the weakening of the localized mantle lithosphere under the central Andes, which is related to variations in magmatism and metasomatism (Kay & Kay 1993; Beck & Zandt 2002; Jing et al., 2020; Contreras-Reyes and Diaz, 2021). Geochemical constraints show that the peridotite melting of mantle-derived magmatism may have occurred within a volcanic field over short time scales (1–5 Myr) in the Altiplano–Puna Plateau. The pattern of melting is consistent with the convective removal of the small-scale mantle lithosphere under the Altiplano–Puna Plateau (Ducea et al., 2013). Could such rheological heterogeneity in the mantle lithosphere account for the observed small-scale lithosphere thinning? Furthermore, with lithospheric heterogeneities, how does DisGBS regulate the removal of the whole lithospheric mantle? Can Rayleigh–Taylor instability developing in the perturbations of short wavelengths remove the mantle lithosphere with DisGBS to the Moho in the homogenous rheological model? In order to address these questions, we constructed a series of 2-D high-resolution thermo-mechanical models with homogeneous mantle lithosphere, including uniformly distributed perturbations, localized weak mantle lithosphere, and multiple terranes/blocks to systematically investigate the dynamics and rheological constraints of the small-scale removal of the mantle lithosphere under orogens.

2 Numerical model

2.1 Governing equations

We simulated lithospheric thinning in collisional orogens by numerically solving the governing equations of mass, momentum, and energy conservation in two-dimensional (2-D) finite difference models with a marker-in-cell technique (Gerya & Yuen 2003). The governing equations are as follows:where v is velocity, is the deviatoric stress tensor, P is pressure, ρ is density, g is gravity acceleration, c_p is heat capacity, T is temperature, k is thermal conductivity, H_r is the radioactive heating rate, H_a is the adiabatic heating rate (, α is thermal expansivity, and H_s is the shear heating rate (, strain rate tensor).

For the density of a specific rock type, depends on pressure (P) and temperature (T):where is the reference density under the conditions of = 0.1 MPa and = 298 K (i.e., pressure and temperature at Earth’s surface); and are the thermal expansion coefficient and the compressibility coefficient.

2.2 Viscoplastic rheology

1 Viscous rheology

where is the second invariant of the strain rate tensor; , the pre-exponential viscous factor; E, the activation energy; V, the activation volume; G, the grain size; , the water fugacity; and n (creep exponent) are experimentally determined flow law parameters (Table 1). The parameter R is the gas constant, p is the grain size exponent, and r is the water fugacity exponent.

Dislocation-accommodated grain boundary sliding is used in our numerical model. Recent studies have proposed that DisGBS may dominate upper mantle deformation (Hansen et al., 2011; Tomohiro and Takaaki, 2015). DisGBS is more effective at a larger strain rate and a smaller grain size, which is consistent with the upper mantle of orogens under convergence. We have applied the same rheological profile for DisGBS to the weak mantle lithosphere and asthenosphere, indicating that the mantle lithosphere of weak terrane is weakened from potential fluid/melt during previous oceanic subduction or terrane accretion. Consequently, the relatively low viscosity based on wet olivine rheology and low plastic strength are applied for the lithospheric mantle of the weak terrane (Tables 1, 2). We used the effective viscosity of “wet quartzite” for both the lower and upper crust of weak terrane (Tables 1, 2). In contrast, the flow law of “wet quartzite” is used for the upper continental crust and “Plagioclase ” for the lower continental crust in the strong terrane (Tables 1, 2). “Dry olivine” for the lithospheric mantle and the asthenosphere (Tables 1, 2) is generally used for the strong continental lithosphere of the strong terrane, where the fluid/melt weakening effects are neglected; thus, a high plastic effective friction coefficient () is used for the dry and strong lithospheric mantle.

2 Drucker–Prager plasticity

The extended Drucker–Prager yield criterion (e.g., Ranalli, 1995) is adopted in our model to simulate the viscoplastic behavior of the lithosphere:where is the viscosity of the Drucker–Prager plasticity, is the yield stress, P is the dynamic pressure, is the residual rock strength at p=0, and is the internal frictional angle of dry rocks. is the pore fluid/melt coefficient that controls the brittle strength of fluid/melt containing porous or fractured media: . Because the fluid/melt parameter is not directly calculated in the model, we use a range of “” values to represent the effective friction coefficient, based on previous systematic investigations (e.g., Gerya & Meilick 2011; Vogt et al., 2012; Li et al., 2016). The strain weakening effect is included in the plastic rheology, in which both the cohesion C₀ and effective friction coefficient decrease with increased strain, as shown in Table 2.

The minimum value of the viscous or plastic viscosity defines the effective viscosity in the model (Ranalli 1995):

2.3 Initial model configuration and boundary conditions

Large-scale models () were built to study the dynamics of the small-scale removal of the mantle lithosphere. Using a non-uniform rectangular numerical grid, the collision zone is represented by km high-resolution grids, while km grids are used for the rest of the model domain. More than 15 million active Lagrangian markers are used to trace and mark internal lithological boundaries, material properties, and temperature. Although the model length is 4,000 km, the deformation localizes in a relatively narrow (ca. 250 km and 1,000 km) region of interest. The effects of the Earth’s curvature are, therefore, neglected in this simplified Cartesian model.

The lithosphere includes a 20-km thick upper crust, a 15-km thick lower crust, and a 105-km thick lithospheric mantle, underlain by the asthenosphere. Some orogens, such as the central Andes, northeastern Tibetan Plateau, and Tianshan, could have experienced a series of subduction and lithospheric detachments (Ducea et al., 2013; Deng et al., 2018; Jing et al., 2020; Huangfu et al., 2021; Li et al., 2022) which could have reduced the viscosity of the mantle (Lei et al., 2019). Moreover, the lithospheric mantle under continental orogens that are adjacent to a subduction zone, such as South America and the Gibraltar Arc, may be modified by subducting plates, leading to the development of secondary downwelling in the continental interior (Levander and Berzada, 2014). Consequently, three groups of models were constructed to simulate the small-scale removing mantle under orogens, and a series of initial perturbations are integrated into the models of Type I∼Type III. In the model of Type I, the lithosphere is set to be weak and homogeneous, and a series of initial perturbations with a constant wavelength (u) and perturbation () are designed at the bottom of the mantle lithosphere (Figure 1A). In the model with heterogeneous rheology (Type II and Type III), the lithosphere is divided into strong and weak terranes, with the initial perturbation imposed on the bottom of the mantle lithosphere in the weak terrane (Figures 1B, C). The properties of various rock types are summarized in Tables 1 and 2. We assumed the same reference density for the mantle lithosphere and asthenosphere (Table 2). Thus, the gravitational instability of the mantle lithosphere arises from the thickening of the weak lithospheric mantle, which becomes denser than the asthenosphere because of the lower temperature. The initial thermal structure of the lithosphere (the white lines in Figure 1A) is laterally uniform with a linear gradient from at the surface to at the bottom of the lithosphere (Turcotte & Schubert 2002). The initial adiabatic thermal gradient in the asthenosphere is .

FIGURE 1

The velocity boundary conditions in the model of heterogeneity rheology include free slip for the left, right, and top boundaries and a permeable boundary for the lower boundary (Burg & Gerya 2005; Li et al., 2016). This infinity-like external free-slip condition along the lower boundary implies a free slip condition to be satisfied at about 100 km below the base of the model. The external free slip allows global conservation of mass in the computational domain and is implemented by using the following limitation for velocity components at the lower boundary: and , where is the vertical distance from the lower boundary to the external boundary where free slip (, ) is satisfied. The lithosphere is pushed from both sides with a constant convergence velocity () imposed on the left and the right parts of the models (Figure 1A).

The thermal boundary conditions contain a fixed temperature () at the upper boundary and zero horizontal heat flux across the vertical boundaries. For the lower thermal boundary, a constant temperature condition is imposed at a great depth (1,000 km) below the bottom of the model to allow both temperature and vertical heat flux to vary along the permeable lower boundary of the model domain and to be adjusted dynamically during the model’s evolution (Li et al., 2016).

3 Model results

3.1 Type I model with uniformly distributed perturbations

In the case of the homogeneous mantle lithosphere, models of classical delamination predict a large-scale peeling of the mantle lithosphere (Bird 1978, 1979; Lei et al., 2019), which may not adequately explain the small-scale lithospheric mantle removal under some orogens. Rayleigh–Taylor instability could occur in small-scale drips when viscosity is relatively low (Beall et al., 2017; Lei et al., 2019). However, will short wavelength drips remove the mantle lithosphere from its base to the Moho? In order to answer this question, we have performed a series of 2-D numerical experiments to systematically investigate the effects of DisGBS and plasticity on the removal of the mantle lithosphere with uniformly distributed perturbations.

3.1.1 Effect of perturbation wavelength on model evolution

In this case, we used dislocation-accommodated grain boundary sliding (DisGBS) and plasticity () to study the effects of perturbation wavelength on convective thinning. Figure 2 shows the model results of different perturbation wavelengths from 150 km to 800 km with constant amplitude () at 3.4 Myr, where convergent velocity . When the initial perturbation wavelength is small (), the drips cannot effectively remove the mantle lithosphere (Figures 2A, B). In contrast, with a larger wavelength initial perturbation (), the lower mantle lithosphere can be easily removed by the drips (Figures 2C, D).

FIGURE 2

3.1.2 Effect of perturbation amplitude on model evolution

In addition to various perturbation wavelengths, different perturbation amplitudes can also play an important role in convective thinning. Figure 3 shows the model results of different amplitudes from 1 km to 30 km with a constant of wavelength () at 3.4 Myr. With a small perturbation amplitude (Figure 3A, ), no convective dripping occurs. When the initial perturbation amplitude , the drips developed on the bottom of the mantle lithosphere; however, this process cannot effectively remove the lithospheric mantle to the Moho (Figure 3B). With a large amplitude (), obvious convective dripping developed on the bottom of the mantle lithosphere, leading to lower mantle lithosphere removal from ∼150 km to ∼100 km (Figure 3C). Furthermore, in the case of perturbation amplitude , the drips effectively induced a removal of the mantle lithosphere from ∼150 km to ∼50 km (Figure 3D).

FIGURE 3

3.2 Type II model with a weak block

Natural orogens often consist of numerous terranes with contrasting lithologies, thermal states, and rheological structures, or have experienced localized weakening in the mantle lithosphere through metasomatism induced by subduction or collision (Ducea et al., 2013; Deng et al., 2018; Jing et al., 2020). Could such rheological heterogeneity explain the observed small-scale lithospheric thinning? In order to solve these problems, we performed a series of experiments with localized weak lithospheric mantle to systematically investigate the effects of DisGBS and plasticity on the removal of the mantle lithosphere with a weak block.

3.2.1 Effect of weak crust and mantle

In the first set of models, we studied the effects of the weak crust and mantle lithosphere on the removal of the mantle lithosphere. Figure 4 shows the model results with DisGBS and significant plastic yielding (). At the beginning of the experiment, the localized weak mantle lithosphere is thickened by small scale drips (Figure 4B). With further convergence, the obvious fragmentary or small-scale thinning of the mantle lithosphere occurs in the section of weak lithosphere with weak crust (Figures 4C, D). Figures 4E–H show the corresponding effective viscosity and the second invariant of the strain rate. It shows that the strain rate localization goes through the entire mantle lithosphere by conjugate slip and has reduced viscosity in the same way.

FIGURE 4

The model in Figure 4 assumes a constant width of L = 300 km for the weak terrane. Furthermore, we tested the model with different widths of weak terrane (Figure 5). When the width of the weak block is small (L = 100 km), converging continents lead to pure shear thickening under weak terrane (Figure 5A). In contrast, the obvious removal of lithospheric mantle under weak terrane, with thickened crust, occurs at a width of L = 200 km (Figure 5B). The decreasing of plastic equivalent viscosity with high strain rates explains the quick development of instability in the lithospheric mantle on both sides of the weak terrane. This leads to the upwelling of the asthenosphere to the bottom of the crust to trigger delamination under weak terrane. With the increasing width of weak terrane, convective dripping can effectively remove the lithospheric mantle to the Moho without apparently thickened crust (Figures 5C–E).

FIGURE 5

3.2.2 Effect of localized weak mantle lithosphere

In this set of models, the localized weak mantle lithosphere, with a strong crust, is added into model setup to simulate the effects of the local weak mantle lithosphere on the small-scale thinning. Figure 6 shows the evolution results of the composition field with localized weak mantle, where the plastic yield When the width of the weak mantle lithosphere is small (L = 100–200 km), the drips could not cause effective removing under the weak mantle lithosphere (Figures 6A, B). In comparison with models of weak terrane (Figure 5), the removal of lithospheric mantle occurs in the critical width of the weak region (L = 300 km), which is larger than the weak terrane model (Figure 5). With the increasing width of weak regions, convective dripping can effectively remove the lithospheric mantle to the Moho without apparently thickened crust (Figures 6D, E).

FIGURE 6

3.3 Type III model with multiple rigid and weak blocks

We have further considered the effects of multiple weak terranes (Deng et al., 2018). Figure 7 shows the corresponding composition, effective viscosity, and the second invariant of the strain rate field. These results illustrate the feedback between strain rates and effective viscosity when DisGBS and plastic yielding are considered. Developed drips occur in the weak terrane and promote the mantle lithosphere peeling off to the crust (Figures 7A–D). As with the model results of a single section of weak terrane (Figure 4), strain weakening reduces both viscous and plastic viscosity, causing localized drips to develop in the bottom of the mantle lithosphere and the further removal of crust (Figures 7E–J). Like the model results of a single weak block (Figure 5), the effective thinning of the mantle lithosphere occurs in the weak terrane, when its width falls into the range of 200 km to 400 km (Figures 8B–D), except when L = 100 km (Figure 8A).

FIGURE 7

FIGURE 8

4 Discussion

4.1 Comparisons of the three types of models

Our study indicates that small-scale mantle lithosphere removal can develop under orogens, depending on the rheology of the mantle lithosphere. Here, we compare the effects of the rheological structure on the patched or small-scale removal of mantle lithosphere.

4.1.1 Homogeneous model (type I)

Given that dislocation-accommodated grain boundary sliding and plasticity is integrated into the mantle lithosphere of weak terranes, convective removing could occur with dripping; however, in that case it tends to occur in the model with a large perturbation wavelength and amplitude (Figures 2, 3). To quantify the effect of the initial perturbation wavelength and amplitude on the magnitude of the removal of the mantle lithosphere, we conducted an “available buoyancy” scaling analysis. The degree of convective instability depends on the thickness of a potentially unstable layer, which can be caused by the mechanical thickening of the layer. In this study, based on Conrad and Molnar (1999), the gravitational instability of the mantle lithosphere with thickness h can be represented by a Rayleigh number, :where n = 3 is the creep exponent of the accommodated grain boundary sliding (Table 1), is the temperature difference across the layer, g is the gravitational acceleration, is the thermal expansion coefficient, is the thermal diffusivity, is the density of mantle lithosphere, is the rheological strength parameter of the mantle lithosphere, h is the layer thickness of the mantle lithosphere, is the initial wavelength, and is the initial perturbation amplitude. This number is similar to the standard Rayleigh number with the additional factor , which is the “available buoyancy” of the layer. is calculated by the integral of the negative buoyancy divided by viscosity (Conrad & Molnar 1999). This integrates the variation of density, viscosity, and temperature with depth in an unstable layer; the dimensionless growth rate should depend only on the wavelength of the initial perturbation. Numerical experiments show that the layer will be unstable when is greater than 100 (Conrad & Molnar 1999; Conrad 2000).

In our models, the mantle lithosphere density profile depends on pressure and temperature, while the effective viscosity profile is given by both plasticity and the background strain rate associated with convergence. Figure 9 shows the calculated as a function of layer thickness at different initial perturbation wavelengths and amplitudes, respectively. At the small amplitude (), the analysis indicates that the whole mantle lithosphere remains stable, with for all thicknesses (Figure 9A). This agrees with the model result that indicates that the mantle lithosphere is stable and thickened during convergence when the small amplitude is applied (Figure 3A). When the initial perturbation amplitude falls in the range of 10 km to 30 km, the calculations show that the large wavelength of the initial perturbation becomes easier to super-exponentially grow and further remove part of the mantle lithosphere (Figures 9B–D). The thickness represents the minimum amount of mantle lithosphere that may be removed by gravitational instability, which indicates that the critical of 100 is obtained for layer thickness. In the models with the larger amplitude (), the calculated critical of 100 is obtained for the layer thickness of ∼45 km, ∼18 km, ∼8 km, and ∼5 km, corresponding to perturbation wavelengths of 300 km, 500 km, 800 km, and 1,000 km (Figure 9C). Because the initial amplitude () is larger than the critical thickness , corresponding to a wavelength range of 500 km–1000 km, the gravitational instability will develop at the base of mantle lithosphere. The results agree with the models, which indicate that the mantle lithosphere is unstable with dripping at a large wavelength of the initial perturbation (Figures 2C, D). For the same reason, in the models with constant wavelength (u=500 km), the calculated critical of 100 is obtained for a layer thickness of , ∼65 km, ∼18 km, and ∼10 km, corresponding to the respective perturbation amplitudes of 1 km, 10 km, 20 km, and 30 km (Figures 9A–D), where indicates that the mantle lithosphere is stable during convergence (Figure 3A). When the initial perturbation amplitude is 10 km, the minimum amount of removal of the mantle lithosphere (∼65 km) is larger than the initial amplitude (∼10 km) with a 500 km initial wavelength (Figure 9B); this indicates that the evolution of the model with a 10 km initial amplitude will have difficulty reaching super-exponential growth and the removal of the mantle lithosphere (Figure 3B). In contrast, with the initial amplitude , super-exponential growth develops at the base of the mantle lithosphere (Figure 3C) because the initial amplitude () is larger than the critical thickness (∼18 km), corresponding to a wavelength of 500 km (Figure 9C).

FIGURE 9

The thickness only estimates the minimum amount of the layer of removal by instability. Larger thicknesses are also potentially unstable ( >100 in Figures 9B, C). The rapid growth triggered by larger instability also depends on the convergence and weakening of strain rate localization. Thickening of the mantle lithosphere by horizontal shortening can lead to dripping to remove it in several obvious ways, which include the lithosphere weakening with increasing strain rate—as is expected for the mantle lithosphere with non-Newtonian viscosity (Molnar et al., 1998; Lei et al., 2019). Moreover, the effect of plasticity promotes dripping to remove the larger quantity of the mantle lithosphere than the model without plasticity. Figure 10 shows the amount of removed mantle in the models with an initial wavelength of 800 km and amplitude of 20 km, as well as a convergence velocity of 4.5 cm/yr. First, the models show that super-exponential growth occurs with a thickening of ∼30 km, which is calculated from the depth of the initial base of the thickened layer () to the depth at which the initial thickening develops super-exponential growth (Figures 10A, B). Second, the model with DisGBS (Figure 10A) shows that the thickness of the removed layer () is ∼30 km. In comparison, the larger amount of the removed layer (∼110 km) of occurs in the model with plasticity (Figure 10B), which indicates that plasticity promotes instability to remove the larger amount of potentially unstable mantle lithosphere ().

FIGURE 10

4.1.2 Heterogeneity model (type II and type III): Strain rate localization

Classical delamination and convective thinning are used to account for the removal of the mantle lithosphere with a homogeneous rheological structure, driven by dripping behaviors (Beall et al., 2017; Lei et al., 2019). Delamination would trigger large-scale peeling-off of the mantle lithosphere, which cannot explain the patched removal of the lithosphere under the Central Andes and the northeastern Tibetan Plateau (Ducea 2011; Lei et al., 2019). The feedback between the strain rate localization and reduced viscosity, based on dislocation creep, is concentrated at the base of the homogeneous lithosphere.

For orogens related to ocean-to-continent subduction, strong hydration processes can lead to local weakening of the mantle lithosphere and cause rheological heterogeneity. The feedback between the strain rate and effective viscosity can reduce the rheological strength of weak mantle lithosphere, causing patched removal of the mantle lithosphere under orogens (Figure 6). On the other hand, the model with multiple weak terranes and weak crust displays obviously patched removal of the mantle lithosphere, indicating a larger parameter range for lithospheric removal than the model of a local weakening mantle lithosphere (Figure 4 and Figure 5). However, the strain rate localization caused by dislocation creep (Ranalli 1995) cannot lead to conjugate slip, which instead promotes the small scale removal of mantle lithosphere (Lei et al., 2020). In contrast, dislocation-accommodated grain boundary sliding can promote strain rate localization in weak terranes, which causes the local weak mantle lithosphere to sink into the asthenosphere by dripping (Figure 11).

FIGURE 11

Lei et al., 2020). (B) Second invariants of the strain rate in the weak mantle lithosphere with dislocation-accommodated grain boundary sliding (Table 1).

Systematic numerical modeling suggests that the rheological range of the small-scale removal of the mantle lithosphere and the associated critical width depends on the heterogeneity in the rheological structure. The mantle lithosphere removal to the Moho in small-scale regions requires low plastic yielding stress () (Figure 5). In comparison, the model results with localized weak mantle lithosphere and strong crust indicate that obvious removal can occur with the condition of low plastic yielding stress () and larger critical width (∼300 km) than is the case with weak block/terranes (Figure 6).

4.2 Geological applications

The history of surface uplift and magmatism is key to discerning delamination or convective thinning under orogens (Göğüş & Pysklywec 2008a; Göğüş & Pysklywec 2008b; Li et al., 2016; Huangfu et al., 2018; Lei et al., 2019). However, it is difficult to identify the patched or fragmental removal of mantle under orogens by means of the two geological indicators because patched removal causes neither symmetric pattern of convective thinning nor migration of delamination.

Adiabatic upwelling of asthenospheric mantle triggered by dripping has been regarded as the most significant expected geological response (Kay & Kay 1993; Ducea & Saleeby 1998). However, it is difficult for small drip with short wavelength to remove the whole mantle lithosphere, which further restricts adiabatic upwelling of the asthenosphere in a few million years (Beall et al., 2017). Thus, it cannot account for the patched removal of lithospheric mantle (to the Moho) and rapid magmatism in the northeastern Tibetan Plateau, Tianshan, and central Andes (Ducea 2011; Ducea et al., 2013; Deng et al., 2018). The results of numerical modeling shed light on the tectonic evolution of orogens with localized weak mantle of metasomatism or the suture between terranes.

4.2.1 Northeastern Tibetan Plateau and Tianshan

The growth of the northeastern TP has been accompanied by continuous lithosphere-scale shortening since the Neogene, during which the inherited lithospheric heterogeneities have played a key role in this shortening (Wang et al., 2008; Yin and Dang, 2008; Zhang and Wang, 2014). The joint inversion of receiver functions reveals patches of seismically fast and slow anomalies in the mantle (Deng et al., 2018). This removal of the mantle lithosphere with patches under the northeastern Tibetan Plateau is reconciled in the models with weak lithospheric blocks (Figures 7A–D).

The central Tianshan together with the adjacent Kazakh Shield to the north and Tarim Craton to the south manifest obvious rheological heterogeneity (Bing et al., 2022; Li et al., 2022; Zhang et al., 2022). Previous studies attributed large-scale intracontinental subduction to the rapid uplift of the central Tianshan since ∼11 Myr, including the northward subduction of the Tarim Craton or the southward subduction of the Kazakh Shield (Gilligan et al., 2014; Zhang et al., 2022). However, the major shortening of the lithosphere is distributed throughout the Tianshan orogen (Thompson et al., 2002; Zubovich et al., 2010). Geophysical investigations reveal a weak block in the central Tianshan and a range of low-velocity anomalies in the lithosphere beneath the Tianshan at a depth of ∼180 km, indicating a local high temperature and weak rheology there (Lei and Zhao, 2007; Li et al., 2009). The weak mantle lithosphere and crust beneath the Tianshan accommodates strain induced by the distant India–Asia collision, which may have promoted the small-scale removal of the mantle lithosphere in the Tianshan since the Miocene. Our numerical model with localized weak crust and the mantle lithosphere shed light on the mantle lithosphere of the weak block, which is consistent with a high temperature mantle lithosphere presently under the Tianshan and can explain the onset of rapid uplift of the Tianshan since ∼11 Myr (Figure 12).

FIGURE 12

Li et al., 2022.

4.2.2 Central Andes

Another case with small-scale removal of the mantle lithosphere is the central Andes. In contrast to the northeastern TP, the mantle lithosphere under the central Andes does not include multiple terranes. The orogeny of the central Andes commenced in the Cenozoic (particularly in the past 30 Myrs). However, the Nazca Plate has been subducting along the western margin of the South American Plate for more than 200 Myrs (Isacks 1988; Allmendinger and Gubbels, 1996; Sobolev and Babeyko, 2005). This long history of subduction has caused a migration of aqueous fluids from the subducting Nazca plate, which may weaken the localized mantle lithosphere under the central Andes (Jing et al., 2020; Contreras-Reyes and Diaz, 2021; Wu et al., 2022). Geochemical constraints and geophysical data indicate a small-scale removal of the mantle lithosphere under the central Andes, with a pattern of melting in accordance with the process of small-scale foundering/dripping (<50 km diameter) of the thickened mantle lithosphere in the Altiplano–Puna Plateau, where mafic volcanic regions on the plateau manifest individual dripping (Drew et al., 2009; Ducea et al., 2013). This may be caused by the localized removal of the weak mantle lithosphere under the Altiplano–Puna Plateau. Geophysical observations also indicate a similar patched removal of lithospheric mantle under the southern Puna Plateau. This process is related to the activity of water released from the Nazca slab. Figure 13 shows the patched/small-scale removal of the mantle lithosphere under the southern Puna Plateau and that the low-velocity body beneath Cerro Galan (L2) reaches the greatest depth as it extends up to ∼50 km deep into the crust. This indicates that the localized weak mantle lithosphere may be removed, followed by adiabatic upwelling of the asthenosphere. Meanwhile, a high velocity body (H2) is located beneath the Moho, indicating patched or small-scale removal of the lithosphere (Jing et al., 2020). These processes may be well explained by the results of numerical modeling with a localized weak lithosphere (Figure 13C), in which the local weakness of the mantle lithosphere is dominated by DisGBS and strong plastic yielding ().

FIGURE 13

Jing et al., 2020), which causes the upwelling of fluids and melts into the crust, triggering the volcanic eruption of Cerro Galan (L2). (B) Interpretation cartoon under the southern Puna Plateau (modified from Jing et al., 2020). Water was released from the slab due to the dehydration reaction, which promotes the localized weak mantle lithosphere and further causes small-scale thinning under the southern Puna Plateau. (C) Model result with a weak mantle lithosphere, based on dislocation-accommodated grain boundary sliding and plastic yielding ().

5 Conclusion

References

• 1

  AllmendingerR. W.GubbelsT. (1996). Pure and simple shear plateau uplift, Altiplano-Puna, Argentina and Bolivia. Tectonophysics259, 1–13. 10.1016/0040-1951(96)00024-8

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  BajoletF.GaleanoJ.FunicielloF.MoroniM.NegredoA. M.FaccennaC. (2012). Continental delamination: Insights from laboratory models. Geochem. Geophys. Geosystems13, 1525–2027. 10.1029/2011gc003896

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  BeallA. P.MoresiL.SternT. (2017). Dripping or delamination? A range of mechanisms for removing the lower crust or lithosphere. Geophys. J. Int.210, 671–692. 10.1093/gji/ggx202

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  BeckS. L.ZandtG. (2002). The nature of orogenic crust in the central Andes. J. Geophys. Res. Solid Earth107 (10), ESE7-1–ESE7-16. 10.1029/2000jb000124

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  BillenM. I. (2008). Modeling the dynamics of subducting slabs. Annu. Rev. Earth Planet36, 325–356. 10.1146/annurev.earth.36.031207.124129

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  BingF. Z.XueW. B.YiX. X. (2022). Seismic anisotropy in the central Tien Shan unveils rheology-controlled deformation during intracontinental orogenesis. Geology50, 812–816. 10.1130/g49633.1

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  BirdP. (1979). Continental delamination and the Colorado plateau. J. Geophys. Res. Solid Earth84 (13), 7561–7571. 10.1029/jb084ib13p07561

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  BirdP. (1978). Initiation of intra continental subduction in the Himalaya. J. Geophys. Res. Solid Earth83 (10), 4975–4987. 10.1029/jb083ib10p04975

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  BittnerD.SchmelingH. (1995). Numerical modelling of melting processes and induced diapirism in the lower crust. Geophys. J. Int.123, 59–70. 10.1111/j.1365-246x.1995.tb06661.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  BurgJ. P.GeryaT. V. (2005). The role of viscous heating in barrovian metamorphism of collisional orogens: Thermomechanical models and application to the lepontine dome in the central alps. J. Metamorph. Geol.23, 75–95. 10.1111/j.1525-1314.2005.00563.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  ClauserC.HuengesE. (1995). “Thermal conductivity of rocks and minerals,” in Rock Physics and phase relations. Editor AhrensT. J. (Washington D. C.: AGU), 105–126. Reference Shelf 3.

  • Google Scholar
  • View reference in article

• 12

  ConradC. P. (2000). Convective instability of thickening mantle lithosphere. Geophys. J. Int.143, 52–70. 10.1046/j.1365-246x.2000.00214.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  ConradC. P.MolnarP. (1999). Convective instability of a boundary layer with temperature and strain-rate-dependent viscosity in terms of ‘available buoyancy. Geophys. J. Int.139, 51–68. 10.1046/j.1365-246x.1999.00896.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  ConradC. P.MolnarP. (1997). The growth of Rayleigh–Taylor type instabilities in the lithosphere for various rheological and density structures. Geophys. J. Int.129 (1), 95–112. 10.1111/j.1365-246x.1997.tb00939.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  Contreras-ReyesE.DiazD.Bello-GonzalezJ. P.SlezakK.PotinB.ComteD.et al (2021). Subduction zone fluids and arc magmas conducted by lithospheric deformed regions beneath the central Andes. Sci. Rep.11, 23078. 10.1038/s41598-021-02430-9

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  CurrieC. A.HuismansR. S.BeaumontC. (2008). Thinning of continental backarc lithosphere by flow-induced gravitational instability. Earth Planet. Sci. Lett.269 (3-4), 436–447. 10.1016/j.epsl.2008.02.037

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  DengY. F.LiJ. T.SongX. D.ZhuL. (2018). Joint inversion for lithospheric structures: Implications for the growth and deformation in northeastern Tibetan plateau. Geophys. Res. Lett.45 (9), 3951–3958. 10.1029/2018gl077486

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  DrewS. T.DuceaM. N.SchoenbohmL. M. (2009). Mafic volcanism on the Puna Plateau, NW Argentina: Implications for lithospheric composition and evolution with an emphasis on lithospheric foundering. Lithosphere1, 305–318. 10.1130/l54.1

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  DuceaM. N. (2011). Fingerprinting orogenic delamination. Geology39 (2), 191–192. 10.1130/focus022011.1

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  DuceaM. N.SeclamanA. C.MurrayK. E.JianuD.SchoenbohmL. (2013). Mantle-drip magmatism beneath the Altiplano-Puna plateau, central Andes. Geology41 (8), 915–918. 10.1130/g34509.1

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  DuceaM.SaleebyJ. B. (1998). A case for delamination of the deep batholithic crust beneath the Sierra Nevada, California. Int. Geol. Rev.40, 78–93. 10.1080/00206819809465199

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  FangX.LiuD.SongC.DaiS.MengQ. (2013). Oligocene slow and Miocene-quaternary rapid deformation and uplift of the Yumu Shan and north Qilian Shan: Evidence from high-resolution magnetostratigraphy and tectonosedimentology. Geol. Soc. Lond. Spec. Publ.373 (1), 149–171. 10.1144/SP373.5

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  GeryaT. V.MeilickF. I. (2011). Geodynamic regimes of subduction under an active margin: Effects of rheological weakening by fluids and melts. J. Metamorph. Geol.29 (1), 7–31. 10.1111/j.1525-1314.2010.00904.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  GeryaT. V.YuenD. A. (2003). Characteristics-based marker-in-cell method with conservative finite-differences schemes for modeling geological flows with strongly variable transport properties. Phys. Earth Planet. Interiors140 (4), 293–318. 10.1016/j.pepi.2003.09.006

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  GilliganA.RoeckerS. W.PriestleyK. F.NunnC. (2014). Shear velocity model for the Kyrgyz Tien Shan from joint inversion of receiver function and surface wave data. Geophys. J. Int.199, 480–498. 10.1093/gji/ggu225

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  GogusO. H.PysklywecR. N. (2008a). Mantle lithosphere delamination driving plateau uplift and synconvergent extension in eastern Anatolia. Geology36 (9), 723–726. 10.1130/g24982a.1

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  GogusO. H.PysklywecR. N. (2008b). Near-surface diagnostics of dripping or delaminating lithosphere. J. Geophys. Reserch113, B11404. 10.1029/2007jb005123

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  GogusO. H.PysklywecR. N.ŞengörA.GünE. (2017). Drip tectonics and the enigmatic uplift of the central anatolian plateau. Nat. Commun.8, 1538. 10.1038/s41467-017-01611-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  HansenL. N.ZimmermanM. E.KohlstedtD. L. (2011). Grain boundary sliding in San Carlos olivine: Flow law parameters and crystallographic-preferred orientation. J. Geophys. Reserch116, B08201. 10.1029/2011jb008220

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  HousemanG. A.MolnarP. (1997). Gravitational (Rayleigh-Taylor) instability of a layer with non-linear viscosity and convective thinning of continental lithosphere. Geophys. J. R. Astronomical Soc.128, 125–150. 10.1111/j.1365-246x.1997.tb04075.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  HuangfuP. P.LiZ. H.GeryaT.FanW.ZhangK. J.ZhangH.et al (2018). Multi-terrane structure controls the contrasting lithospheric evolution beneath the Western and central–eastern Tibetan plateau. Nat. Commun.9, 3780. 10.1038/s41467-018-06233-x

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  HuangfuP. P.LiZ. H.ZhangK. J.FanW.ZhaoJ.ShiY. (2021). India-Tarim lithospheric mantle collision beneath Western Tibet controls the Cenozoic building of Tian Shan. Geophys. Res. Lett.48, e2021GL094561. 10.1029/2021GL094561

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  IsacksB. L. (1988). Uplift of the Central Andean plateau and bending of the Bolivian orocline. J. Geophys. Res.93, 3211–3231. 10.1029/jb093ib04p03211

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  JingC.SofiaK. K.YuanX.HeitB.WuH.YangD.SchurrB.et al (2020). Lithospheric delamination beneath the southern Puna plateau Resolved by local earthquake tomography. J. Geophys. Res.125, e2019JB019040. 10.1029/2019JB019040

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  KayR. W.KayS. M. (1993). Delamination and delamination magmatism. Tectonophysics219, 177–189. 10.1016/0040-1951(93)90295-u

  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  KirbyS. H.KronenbergA. K. (1987). Rheology of the lithosphere: Selected topics. Rev. Geophys25, 1219–1244. 10.1029/rg025i006p01219

  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  KirbyS. H. (1983). Rheology of the lithosphere. Rev. Geophys. Space Phys.21, 1458–1487. 10.1029/rg021i006p01458

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  LeiJ.ZhaoD. (2007). Teleseismic P-wave tomography and the upper mantle structure of the central Tien Shan orogenic belt. Phys. Earth Planet. Interiors162, 165–185. 10.1016/j.pepi.2007.04.010

  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  LeiT.LiZ. H.LiuM. (2020). Dynamics of small-scale removing of mantle lithosphere under orogens. Chin. J. Geophys.63 (10), 3727–3739.

  • Google Scholar
  • View reference in article

• 40

  LeiT.LiZ. H.LiuM. (2019). Removing mantle lithosphere under orogens: Delamination versus convective thinning. Geophys. J. Int.219 (2), 877–896. 10.1093/gji/ggz329

  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  LevanderA.BerzadaM. J.NiuF.HumphreysE. D.PalomerasI.ThurnerS. M.et al (2014). Subduction-driven recycling of continental margin lithosphere. Nature515, 253–256. 10.1038/nature13878

  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  LevanderA.SchmandtB.MillerM.LiuK.KarlstromK.CrowR.et al (2011). Continuing Colorado plateau uplift by delamination-style convective lithospheric downwelling. Nature472 (7344), 461–465. 10.1038/nature10001

  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  LiW.ChenY.YuanX. H.XiaoW.WindleyB. F. (2022). Intracontinental deformation of the Tianshan Orogen in response to India-Asia collision[J]. Nat. Commun.13, 3738. 10.1038/s41467-022-30795-6

  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  LiZ. H.LiuM.GeryaT. V. (2016). Lithosphere delamination in continental collisional orogens: A systematic numerical study. J. Geophys. Res. Solid Earth121, 5186–5211. 10.1002/2016jb013106

  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  LiZ. W.RoeckerS.ZhihaiL.BinW.HaitaoW.SchelochkovG.et al (2009). Tomographic image of the crust and upper mantle beneath the Western Tien Shan from the MANAS broadband deployment: Possible evidence for lithospheric delamination. Tectonophysics477, 49–57. 10.1016/j.tecto.2009.05.007

  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  LiaoJ.WangQ.TarasG.BallmerM. D. (2017). Modeling Craton destruction by hydration-induced weakening of the upper mantle. J. Geophys. Res. Solid Earth122, 7449–7466. 10.1002/2017jb014157

  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  LorincziP.HousemanG. (2009). Lithospheric gravitational instability beneath the southeast Carpathians. Tectonophysics474 (1), 322–336. 10.1016/j.tecto.2008.05.024

  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  MolnarP.HousemanG. A.ConradC. P. (1998). Rayleigh—taylor instability and convective thinning of mechanically thickened lithosphere: Effects of non-linear viscosity decreasing exponentially with depth and of horizontal shortening of the layer. Geophys. J. Int.133 (3), 568–584. 10.1046/j.1365-246x.1998.00510.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  MorencyC.DoinM. P. (2004). Numerical simulations of the mantle lithosphere delamination. J. Geophys. Res. Solid Earth109 (3), B03410. 10.1029/2003jb002414

  • CrossRef
  • Google Scholar
  • View reference in article

• 50

  OwensT. J.ZandtG. (1997). Implications of crustal property variations for models of Tibetan plateau evolution. Nature387, 37–43. 10.1038/387037a0

  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  RanalliG.MurphyD. C. (1987). Rheological stratification of the lithosphere. Tectonophysics132 (4), 281–295. 10.1016/0040-1951(87)90348-9

  • CrossRef
  • Google Scholar
  • View reference in article

• 52

  RanalliG. (1995). Rheology of the Earth, deformation and flow process in Geophysics and geodynamics. 2nd ed.London: Chapman & Hall, 413.

  • Google Scholar
  • View reference in article

• 53

  SchottB.SchmelingH. (1998). Delamination and detachment of a lithospheric root. Tectonophysics296 (3-4), 225–247. 10.1016/s0040-1951(98)00154-1

  • CrossRef
  • Google Scholar
  • View reference in article

• 54

  SchurrB.RietbrockA.AschG.KindR.OnckenO. (2006). Evidence for lithospheric detachment in the central Andes from local earthquake tomography. Tectonophysics415 (1-4), 203–223. 10.1016/j.tecto.2005.12.007

  • CrossRef
  • Google Scholar
  • View reference in article

• 55

  SobolevS. V.BabeykoA. Y. (2005). What drives orogeny in the Andes?Geology33, 617–620. 10.1130/g21557ar.1

  • CrossRef
  • Google Scholar
  • View reference in article

• 56

  ThompsonS. C.WeldonR. J.RubinC. M.AbdrakhmatovK.MolnarP.BergerG. W. (2002). Late quaternary slip rates across the central Tien Shan, Kyrgyzstan, central Asia. J. Geophys. Res. Solid Earth107, ETG 7-1–ETG 7-32. 10.1029/2001jb000596

  • CrossRef
  • Google Scholar
  • View reference in article

• 57

  TomohiroO.TakaakiK.HigoY.FunakoshiK. I.SuzukiA.KikegawaT.et al (2015). Dislocation-accommodated grain boundary sliding as the major deformation mechanism of olivine in the Earth’s upper mantle. Sci. Adv.1 (9), e1500360. 10.1126/sciadv.1500360

  • CrossRef
  • Google Scholar
  • View reference in article

• 58

  TurcotteD. L.SchubertG. (2002). Geodynamics. 2nd ed.New York: Cambridge Univ, 456.

  • Google Scholar
  • View reference in article

• 59

  van HunenJ.ZhongS.ShapiroN. M.RitzwollerM. H. (2005). New evidence for dislocation creep from 3-D geodynamic modeling of the Pacific upper mantle structure. Earth Planet, Lett.238, 146–155. 10.1016/j.epsl.2005.07.006

  • CrossRef
  • Google Scholar
  • View reference in article

• 60

  VogtK.GeryaT.CastroA. (2012). Crustal growth at active continental margins: Numerical modeling. Phys. Earth Planet. Interiors192-193, 1–20. 10.1016/j.pepi.2011.12.003

  • CrossRef
  • Google Scholar
  • View reference in article

• 61

  WangH. L.CurrieC. A. (2015). Magmatic expressions of continental lithosphere removal. J. Geophys. Res. Solid Earth120, 7239–7260. 10.1002/2015jb012112

  • CrossRef
  • Google Scholar
  • View reference in article

• 62

  WangQ.WymanD. A.XuJ.DongY.VasconcelosP. M.PearsonN.et al (2008). Eocene melting of subducting continental crust and early uplifting of central Tibet: Evidence from central-Western Qiangtang high-K calc-alkaline andesites, dacites and rhyolites. Earth Planet. Sci. Lett.272 (1-2), 158–171. 10.1016/j.epsl.2008.04.034

  • CrossRef
  • Google Scholar
  • View reference in article

• 63

  WuY.LiaoJ.GuoF.WangX. C.ShenY. (2022). Styles of trench-parallel mid-ocean ridge subduction affect Cenozoic geological evolution in circum-Pacific continental margins. Geophys. Res. Lett.49. 10.1029/2022GL098428

  • CrossRef
  • Google Scholar
  • View reference in article

• 64

  XiaL. Q.LiX. M.MaZ. P.XuX.XiaZ. (2011). Cenozoic volcanism and tectonic evolution of the Tibetan Plateau. Gondwana Res.19 (4), 850–866. 10.1016/j.gr.2010.09.005

  • CrossRef
  • Google Scholar
  • View reference in article

• 65

  YinA.DangY. Q.ZhangM.ChenX. H.McRivetteM. W. (2008). Cenozoic tectonic evolution of the Qaidam basin and its surrounding regions (Part 3): Structural geology, sedimentation, and regional tectonic reconstruction. Geol. Soc. Am. Bull.120 (7-8), 847–876. 10.1130/B26232.1

  • CrossRef
  • Google Scholar
  • View reference in article

• 66

  ZhangB. F.BaoX. W.XuY. X. (2022). Seismic anisotropy in the central Tien Shan unveils rheology-controlled deformation during intracontinental orogenesis. Geology50, 812–816. 10.1130/g49633.1

  • CrossRef
  • Google Scholar
  • View reference in article

• 67

  ZhangZ.WangY.HousemanG. A.XuT.WuZ.YuanX.et al (2014). The Moho beneath Western Tibet: Shear zones and eclogitization in the lower crust. Earth Planet. Sci. Lett.408, 370–377. 10.1016/j.epsl.2014.10.022

  • CrossRef
  • Google Scholar
  • View reference in article

• 68

  ZubovichA. V.WangX. q.ScherbaY. G.SchelochkovG. G.ReilingerR.ReigberC.et al (2010). GPS velocity field for the Tien Shan and surrounding regions. Tectonics29, TC6014. 10.1029/2010tc002772

  • CrossRef
  • Google Scholar
  • View reference in article