Overriding plate topography provides constraints on subduction zone geodynamics. We investigate its evolution using fully dynamic laboratory models of subduction with techniques of stereoscopic photogrammetry and particle image velocimetry. Model results show that the topography is characterized by an area of forearc dynamic subsidence, with a magnitude scaling to 1.44–3.97 km in nature, and a local topographic high between the forearc subsided region and the trench. These topographic features rapidly develop during the slab free-sinking phase and gradually decrease during the steady state slab rollback phase. We propose that they result from the variation of the vertical component of the trench suction force along the subduction zone interface, which gradually increases with depth and results from the gradual slab steepening during the initial transient slab sinking phase. The downward mantle flow in the nose of the mantle wedge plays a minor role in driving forearc subsidence.
- Three-dimensional dynamic analogue subduction models are used to study time-evolving topography of overriding plate
- Forearc topographic subsidence rapidly develops during slab free-sinking phase and decreases during steady state slab rollback phase
- Trench suction at subduction interface is predominantly responsible for the transient forearc subsidence
The origin, evolution, and spatial variability of topography at subduction zones have been studied with geodynamic models for more than two decades (Bonnardot, Hassani, & Tric, 2008; Buiter et al., 2001; Guillaume et al., 2010; Gvirtzman & Nur, 2001; Hampel & Pfiffner, 2006; Hassani et al., 1997; Husson et al., 2012; Martinod et al., 2013, 2015; Schellart & Spakman, 2015; Yang et al., 2016; Zhong & Gurnis, 1992, 1994). Topography shapes the free top surface of the lithosphere and contains important information about the dynamics of the tectonic plates and the sublithospheric mantle. For example, dynamic topography is caused by the vertical movement of the lithosphere in response to its viscous coupling with the underlying sublithospheric mantle. It is different from isostatically supported topography generated by isostatic equilibrium (e.g., crustal thickness variations) and has instead been ascribed to the vertical component of mantle flow below the plates (Boschi et al., 2010; Braun, 2010; Coblentz & Karlstrom, 2011; Flament et al., 2013; Gurnis et al., 1998; Hager et al., 1985; Lithgow-Bertelloni & Silver, 1998; Pysklywec & Mitrovica, 1998; Saleeby & Foster, 2004). Another important form of topography results from the tectonic forces in the plates and at the plate boundaries, such as subduction zones. Investigating topography around subduction zones can provide quantitative and conceptual insights into the interaction between the plates, the slabs, mantle flow, and the associated stresses. As the present-day topography is only transient, geodynamic modeling can be an effective tool to study the evolution of topography back in time and to gain insights into the driving mechanisms of topography formation and evolution.
A number of geodynamic modeling studies have previously been conducted to investigate the evolution of subduction-associated topography during progressive subduction (Cattin et al., 1997; Hampel & Pfiffner, 2006; Hassani et al., 1997; Husson et al., 2012; Zhong & Gurnis, 1994), during slab detachment (Bonnardot, Hassani, & Tric, 2008; Buiter et al., 2001; Gvirtzman & Nur, 2001; Pysklywec & Mitrovica, 1998), or during progressive slab sinking after slab detachment (Braun, 2010; Flament et al., 2013; Gurnis et al., 1998; Hager et al., 1985; Schellart & Spakman, 2015; Yang et al., 2016). In earlier analogue subduction models investigating topography, an overriding plate was either excluded (Husson et al., 2012) or included, but in that case there was either a large physical gap between the two plates (Guillaume et al., 2013, 2010) or convergence in the models was kinematically imposed with a piston to sustain subduction (Luth et al., 2010; Martinod et al., 2013; Sokoutis & Willingshofer, 2011). Such models provide new insights that can help understand the geodynamics of topography associated with subduction. However, the lack of an overriding plate or the large physical gap between the plates prevented the existence of a subduction zone interface with proper dimensions and rheological scaling. Therefore, to study the effect of forces acting at the subduction zone interface on topography, we used a modeling approach that is based on recently developed buoyancy-driven analogue models of subduction incorporating an overriding plate and a realistic interplate coupling (Chen et al., 2015, 2016; Duarte et al., 2013, 2015; Meyer & Schellart, 2013). This modeling approach allows one to investigate the evolution of topography in the forearc region of subduction zones in a fully dynamic framework.
In this work, we have conducted fully dynamic analogue models of time-evolving subduction in three-dimensional space (no externally imposed velocity or force boundary condition). In order to capture and record the progressive subduction process, we have used a stereoscopic particle image velocimetry (sPIV) system that can simultaneously track surface topography, surface deformation, and flow in the sublithospheric mantle. In our models, topographic subsidence was observed in the forearc lithosphere and the mantle flow velocity field in the mantle wedge was mapped to determine if any causal link could be established. Based on an analysis of forces acting on the forearc lithosphere and a comparison of our model results with previous geodynamic modeling studies, we provide new insights into the mechanism that drives subsidence in the forearc region of the overriding lithosphere at subduction zones.
The experiments were carried out in a transparent square tank (Figure 1). The tank was filled with low-viscosity glucose syrup to a depth of 13.3 cm to simulate the viscous sublithospheric upper mantle, such that the rigid tank bottom represents the 670 km discontinuity. Two mixtures of high-viscosity silicone putty and iron powder in different proportions were used to represent the subducting oceanic lithosphere and an overriding lithosphere and were placed on the surface of the glucose syrup. The overriding plate was neutrally buoyant with respect to the upper mantle, whereas a density difference of 100 kg/m3 between the subducting plate and the sublithospheric upper mantle was used to simulate the negative buoyancy of the slab (Cloos, 1993). The viscosity ratio between the subducting plate and the upper mantle was 188–203, which falls in the suggested range of ~100–500 as proposed in earlier subduction modeling work (Funiciello et al., 2008; Ribe, 2010; Schellart, 2008). To achieve a weak coupling at the subduction zone interface, the top surface of the subducting plate was lubricated with a mixture of petrolatum and paraffin oil (Duarte et al., 2013), which has a relatively low effective viscosity of ~0.8–1.5 Pa s and low yield stress. It results in a reasonable viscosity ratio (191–383) between the sublithospheric upper mantle (286–306 Pa s) and the weak lubricant material. The weak lubricant has proven realistic in simulating the subduction zone interface from a comparison between dynamical subduction experiments and natural observations (Duarte et al., 2015).
Both the trailing edges and the lateral edges of the subducting plate were free, representing mid-oceanic ridges and strike-slip faults that offer negligible resistance to plate motion. The lateral edges of the overriding plate were free, while the trailing edge was either free or fixed, representing a small and relatively mobile plate or a large, relatively immobile plate, respectively. Subduction initiation was generated by manually pushing down an ~3–5 cm long slab segment at the tip of the subducting plate. Thereafter, the subduction process was self-controlled due to the negative buoyancy of the subducted slab, which is the only driving force in our experiments. Our experiments were scaled to nature assuming that the sinking velocity of the slab can be approximated by the Stokes velocity (please refer to Duarte et al., 2013 and Jacoby, 1973 for the scaling procedure and factors used in this type of experiments). In the models, 1 cm corresponds to 50 km in nature, 1 s corresponds to 8300 years in nature, and 0.01 mm/s is equivalent to 0.6 cm/yr in nature.
To monitor simultaneously the overriding plate topography and the upper mantle flow occurring in the mantle wedge during progressive subduction, a stereoscopic particle image velocimetry (sPIV) technique was used (Figure 1). The sPIV technique can monitor the subduction-induced mantle flow in the mantle wedge during progressive subduction by using a laser sheet that illuminates fluorescent particles that are homogeneously distributed throughout the mantle fluid. Please refer to Strak and Schellart (2014) for detailed technical information about this type of recording. Furthermore, with a stereoscopic photogrammetry technique we can compare top view photographs to calculate the overriding plate topography. Two high-resolution PIV cameras were installed above the model with an angle of ~30° between them and very light white powder was seeded on the top surface of the overriding plate to allow mapping the overriding plate topography. Another high-resolution PIV camera was located on the side to map the cross-sectional velocity field of the subduction-induced mantle flow occurring in the mantle wedge of the subduction zone along the central plane. To allow computing the cross-sectional mantle flow velocity field we used an interval time, defined as one loop (L). Two normal cameras, one on top and another on the side, were also used to track the kinematics of the subducting and overriding plates using white passive tracers.
3.1 Topography of Overriding Plate
During the early free sinking phase of subduction (defined as the early subduction phase during which the negatively buoyant slab is sinking through the sublithospheric upper mantle before the slab tip starts interacting with the upper-lower mantle discontinuity), the topography of the overriding plate was characterized by an elongated area of subsidence striking parallel to the trench and separated from the trench by a distance of 1–5 cm (Figures 2a and 2b and 2e and 2f). The magnitude of the maximum depression increased with progressive subduction, reaching 1.35–1.79 mm (scaling to 3.00–3.97 km) in the center of the subduction zone (see Figures 3a and 3b, 3e and 3f, and 4b and 4c). Closer to the trench, the topography was still negative but characterized by a trench-parallel ridge with a local relative high. After a phase of interaction between the slab tip and the rigid bottom discontinuity, subduction reached a steady state rollback phase during which the slab geometry, the trench kinematics, and plate kinematics remained relatively constant (Figures 3c and 3d, 3g and 3h, and 4b and 4c). During this phase the magnitude of the depression progressively decreased to 0.65–0.80 mm (scaling to 1.44–1.78 km) (Figures 2c and 2d, 2g and 2h, and 4b and 4c). During the entire subduction process, the rest of the overriding plate, including the backarc region, had a relative flat topography. The horizontal position of the maximum depression in the forearc region was always located within 5 mm (scaling to 25 km) of the deepest point of contact between the subducting and overriding plates at the subduction zone interface (Figure 3). Furthermore, the local ridge in between the depression and the trench became progressively less pronounced during the steady state rollback phase and almost disappeared in the final stages of subduction (Figures 3c and 3d and 3g and 3h).
3.2 Subduction-Induced Mantle Flow in the Mantle Wedge
In the center of the subduction zone, the slab sinking and rolling back in the sublithospheric upper mantle forced the mantle wedge material to flow in a poloidal fashion (Figure 3). This poloidal cell had a downwelling component in the mantle wedge just above the top surface of the slab, as a result of viscous drag by the sinking slab. The downwelling velocity reached ~7 mm/min (scaling to 7 cm/yr) during the free sinking phase and ~3 mm/min (scaling to 3 cm/yr) during the steady state phase (Figure 3), which correlates well with the velocities attained by the subducting plate during these phases (Figures 4b and 4c). We also observed upward mantle flow in the far backarc mantle wedge. This upwelling was stronger during the free sinking phase with a maximum velocity of up to ~2 mm/min (scaling to 2 cm/yr), and it was located at ~75–175 mm (scaling to 375–875 km) from the trench (Figure 3a). The velocity magnitude thereafter decreased during the steady state phase, and the upwelling migrated farther away from the trench as the poloidal cell became elongated horizontally as a response to the slab draping over the 670 km discontinuity (Figures 3c and 3d and 3g and 3h). We note that the vertical velocities of the mantle wedge poloidal flow were higher during the free sinking phase and lower during the steady state phase, generally, in good correlation with the subducting plate velocity (Figures 4b and 4c).
4.1 Origin of the Topographic Subsidence in the Forearc Region
An important observation concerns the horizontal position of the maximum depression, which corresponds with that of the deepest contact point between the subducting and overriding plates at the subduction zone interface (Figure 3). This observation will help us to differentiate between the different possible driving forces responsible for the forearc subsidence. We propose that there are three potential candidates to explain the subsidence formation. Our subduction models are buoyancy driven as in nature and have a realistic subduction interface. Therefore, the three forces acting on the forearc region and likely leading to the forearc topographic subsidence in the models should be comparable in natural subduction settings. These three forces are the shear force at the subduction zone interface, the trench suction force (normal to the subduction zone interface), and the viscous drag force induced by the vertical movement of the mantle flow in the mantle wedge region (Figure 4a).
In several previous works using two-dimensional numerical subduction models, the influence of the subduction zone interface shear force on the forearc topography has been investigated, through varying the coefficient of friction at the subduction zone interface (Buiter et al., 2001; Cattin et al., 1997; Hampel & Pfiffner, 2006; Hassani et al., 1997; Zhong & Gurnis, 1994). Some modeling studies show that increasing the coefficient of friction leads to an increase in forearc depression (Buiter et al., 2001; Hampel & Pfiffner, 2006; Zhong & Gurnis, 1994), whereas others show an opposite relationship (Cattin et al., 1997; Hassani et al., 1997). In our models, an internal subsidence area, located at 1–5 cm away from the trench, has been observed particularly during the slab free-sinking phase. However, this transient subsidence cannot be explained by the shear force at the subduction zone interface.
We have quantified the stress that drives dynamic subsidence in our models. To get the maximum vertical stress for subsidence, the forearc subsidence map has been used, for which the maximum depth of the depression has been quantified, which is ~1.5 mm (Figure 3). The density contrast between the sublithospheric upper mantle (1,428 kg/m3) and air (1.225 kg/m3) is 1,426.775 kg/m3, while the gravitational acceleration is 9.8 m/s2. From the product of the maximum depression depth, density contrast, and gravitational acceleration, we calculate the maximum vertical stress for forearc dynamic subsidence in our models, which is of the order 20.97 Pa. We can also estimate the shear stress at the subduction zone interface. The effective flow stress of the lubrication material filled in the subduction channel of our models is 1.0–1.5 Pa at a shear strain rate of 0.1 s−1 (Duarte et al., 2014). We assume that the thickness of the subduction channel in our models is ~1 mm (Duarte et al., 2013). The average subduction velocity in our models is ~0.06–0.13 mm/s (Figure 4), giving an average shear rate in the subduction channel of 0.06–0.13 s−1, and the subduction dip angle is approximately 30–45°. From this we can calculate that the magnitude of the vertical component of the shear stress at the subduction zone interface is 0.54–1.25 Pa, which is much less than the estimated vertical stress for forearc dynamic subsidence (20.97 Pa). This comparison indicates that the shear force at the subduction zone interface plays only a minor role in driving the internal transient forearc subsidence. Nevertheless, the shear force at the subduction zone interface applies a bending moment to the overriding plate, which would induce a more gradual increase in depression in the overriding plate toward the trench, as observed in the experiments (Figures 3d and 3h).
The second force possibly causing the forearc topographic subsidence is trench suction, which is normal to the subduction zone interface. This force has been proposed by Shemenda (1993) using analogue models of subduction, in which it was described as hydrostatic suction that prevents the plates from being separated, and it was suggested to provide a driver for backarc extension. In our models, trench suction also plays a role in maintaining the two plates in contact. Furthermore, the magnitude of this force would correlate with the level of resistance to translation of the overriding plate. Specifically, the magnitude of this force would be related to the velocity of the slab in the direction perpendicular to the subduction zone interface during progressive slab rollback (Figures 4d and 4e). Moving a longer distance implies larger trench suction. During the slab free-sinking phase, we observe a progressive increase of the slab dip angle with progressive subduction (Figures 3a and 3b and 3e and 3f). An increase in the slab dip angle results in a gradual increase of the vertical component of the trench suction toward the lower part of the interface (Figure 4d), thereby causing a relatively deep topographic depression in the forearc region, with the maximum depth corresponding to the lowest contact point of the interface (Figure 3). In contrast, during the steady state slab rollback phase the vertical component of the trench suction is nearly the same along the subduction zone interface because of the steady state slab rollback velocity along the interface and the relatively constant slab dip angle (Figure 4e). This results in minor variations in the vertical component of the trench suction along the interface, thereby generating a topographic subsidence in the forearc region that is spatially more constant with respect to the slab free-sinking phase. Therefore, during the steady state phase the trench suction force, like the shear force at the subduction zone interface, will mostly bend the overriding plate downward at the trench, thereby promoting a forearc topographic subsidence that gradually increases toward the trench.
We have confirmed that the suction force at the subduction zone interface is predominantly responsible for the transient forearc subsidence. This indicates that the vertical force for the dynamic subsidence and vertical component of trench suction at the subduction zone interface have approximately equal magnitudes. We can roughly estimate the trench suction force from quantifying the force that drives dynamic subsidence. In our models, the width (perpendicular to trench) of the depression is ~0.05 m (Figure 2), and the length of the depression is ~0.15 m. The average depth of the maximum depression is ~0.001 m. So the total volume of the depression is ~3.75 × 10−6 m3. By multiplying the volume of the depression by the density contrast of 1426.775 kg/m3 and the gravitational acceleration, we get an average vertical force that drives the depression in our models of 0.052 N, which is also the approximate magnitude of the suction force. Furthermore, in our models the width and thickness of the overriding plate are 0.15 m and 0.015 m, respectively, and we assume that the slab dip angle is approximately 30–45°. This gives a subduction zone interface surface area of 3.18–4.50 × 10−3 m2. The total shear force results from multiplying the shear stress (1.0–1.5 Pa, given above) by the surface area of the subduction zone interface. The vertical component of the total shear force at the interface is then 0.0016–0.0048 N, which is about an order of magnitude less than the trench suction force.
We can also estimate the negative buoyancy force of the subducting slab and then compare it with the suction force. The density contrast between the slab and lithospheric upper mantle in our models is 100 kg/m3. The width and thickness of the slab are 0.15 m and 0.02 m, respectively. We will assume that there is a 10 cm long subducting slab. This gives a total negative buoyancy force of the subducting slab of 0.29 N, compared with the estimated vertical suction force (0.052 N). Considering that the driving force in our models only comes from the negatively buoyant subducting slab, our model results indicate that some 17.8% of the negative buoyancy of the slab is used to drive the forearc dynamic subsidence during the transient slab free-sinking phase.
The last force to be discussed as being a potential cause for the formation of overriding plate topography is the force generated by vertical movement of mantle flow in the mantle wedge. Previous natural observations and modeling studies have demonstrated that the vertical mantle flow can cause a long-wavelength, low-amplitude topography, which is referred to as dynamic topography (Gurnis et al., 1998; Hager et al., 1985; Lithgow-Bertelloni & Silver, 1998; Pysklywec & Mitrovica, 1998). In our models, an area of maximum downward mantle velocity has always been observed in the mantle wedge very close to the downgoing slab within 10 cm, scaling to 500 km, from the trench (Figure 3). Previous analogue and numerical subduction models have shown that the corner flow induced by slab rollback in the mantle wedge can generate a pressure gradient with higher suction in the nose of the mantle wedge and progressively lower away from it (Hall et al., 2000; Kneller & van Keken, 2007, 2008; MacDougall et al., 2014). Such suction produced by the vertical component of the subduction-induced mantle flow should also exist in our models and would thus influence the forearc topography by a drag down effect. However, this drag force should play a minor role in the formation of the forearc subsidence. This is because in our experiments the maxima of the forearc subsidence were observed above the deepest contact point between the subducting and overriding plates at the subduction zone interface (Figures 3b and 3f) rather than within 5–10 cm from the trench where the downward mantle wedge velocity is maximum.
4.2 Comparison With Previous Studies
4.2.1 Analogue Modeling
Several analogue models of subduction have been previously conducted to investigate the topography around subduction zones, in which the overriding plate was either excluded (Husson et al., 2012) or included, but the two involved plates were separated by a thin layer of glucose syrup, with a distance up to 20 mm (scaling to ~120 km for their models) at the surface between the two plates (Guillaume et al., 2013, 2010). In the work of Husson et al. (2012), the model setup consisted of a free or fixed subducting plate but without an overriding plate. Because of the lack of an overriding plate, the surface of the glucose syrup would represent the Earth's surface and the topography of this surface was scanned. Their model results showed that the forearc region is depressed to a depth of up to 0.20 mm. This subsidence corresponds to a dynamic topography induced by vertical motion of the subduction-induced mantle flow. Considering the similar model setup adopted in their models and our models, except for the lack of an overriding plate in their models, and that the maximum subsidence depth in their models (0.2 mm) is far shallower than that in our models (0.65–1.79 mm), these results support our analysis that the force generated by the downward motion of mantle flow in the mantle wedge makes a secondary contribution to the observed forearc topographic subsidence. In the work of Guillaume et al. (2013), a topographic profile parallel to and at ~15 mm away from the trench showed that there is an area of subsidence existing in the middle of the overriding plate. This is comparable with our model result showing a subsidence observed at an average distance of 30 mm away from the trench. Some minor difference between the two studies is probably because in their models there is a window in the middle of the slab but not in our models. Furthermore, in our models the subducting plate is free while in their models the far-field edge of the subducting plate is fixed.
4.2.2 Numerical Modeling
A number of numerical models have previously investigated the parameters that potentially affect the overriding plate topography. These parameters include the coefficient of friction at the subduction zone interface (Bonnardot, Hassani, Tric, Ruellan, et al., 2008; Buiter et al., 2001; Cattin et al., 1997; Hampel & Pfiffner, 2006; Hassani et al., 1997; Zhong & Gurnis, 1992, 1994), the variation of the slab dip angle (Brink, 2005; Hassani et al., 1997), the density contrast between the lithosphere and the asthenosphere (i.e., the age of the lithosphere) (Bonnardot, Hassani, Tric, Ruellan, et al., 2008; Buiter et al., 2001; Hampel & Pfiffner, 2006; Hassani et al., 1997; Zhong & Gurnis, 1992, 1994), the asthenosphere viscosity (Bonnardot, Hassani, & Tric, 2008; He, 2012), the rheology of the subducting plate (Bonnardot, Hassani, & Tric, 2008; Zhong & Gurnis, 1992, 1994), the trenchward velocity of the overriding plate (Buiter et al., 2001; Hampel & Pfiffner, 2006), and the thickness of the overriding plate (He, 2012). Probably, the first numerical subduction modeling studies in a dynamically self-consistent manner to investigate the topography at subduction zones were published more than two decades ago (Zhong & Gurnis, 1992, 1994). In these subduction models, a fault represents the interface between the subducting plate and the overriding plate and a viscous slab subducts into a viscous medium. This setup is similar to that of our models, except that they used a two-dimensional spatial setup (and our models are three-dimensional) and a prescribed interface geometry (which is self-determined in our models). In the work from Zhong and Gurnis (1992, 1994), the backarc region developed a topographic subsidence with a depth of up to ~1.5 km, typically ~150 km away from the trench, when the plates have a Newtonian rheology. They proposed that it was caused by the mantle flow in the backarc extension region. In addition, in their models a subsidence with a ~3 km depth was observed at the trench. This is comparable with and consistent with our model results showing a subsidence depth of 1.0–1.72 mm at the trench, scaling to 2.22–3.82 km in nature. Their model results implied that two main parameters, namely, the age of the lithosphere and the shear stress at the subduction zone interface, have a significant influence on the overriding plate topography. The impact of these two parameters on the overriding plate topography has been tested later (Bonnardot, Hassani, Tric, Ruellan, et al., 2008; Buiter et al., 2001; Cattin et al., 1997; Hampel & Pfiffner, 2006; Hassani et al., 1997). In these models, the age of the lithosphere, represented by the density contrast between the lithosphere and the asthenosphere, and the coefficient of friction along the subduction zone interface have been varied. The observed topographic subsidence in the forearc region was comparable with that in our models. For example, in the work of Hassani et al. (1997), when the density contrast is 100 kg/m3, which is the same as that in our models, an area of topographic subsidence was observed in the forearc region. In their models, the total convergence amount since subduction initiation was only 400 km, indicating that the slab had not yet reached the 670 km discontinuity. Therefore, their models can only be compared with the slab free-sinking phase in our models. They ascribed the forearc topographic subsidence of the overriding plate to the hydrostatic suction (referred to as trench suction in our models) at the subduction zone interface. This explanation for the formation of topographic subsidence in the forearc region is in agreement with our present analysis. Furthermore, in the work of Buiter et al. (2001), a velocity of 4 cm/yr was imposed on a subducting plate with a thickness of 25 km. The model results showed a subsidence with a depth of up to 4 km in the forearc region, which the authors ascribed to slab rollback. This is in general agreement with the findings presented in this study. In addition, Buiter et al. (2001) found that an increase in the density contrast between the lithosphere and the asthenosphere deepens the forearc subsidence to a maximum of 6 km. Such subsidence depths are comparable to the subsidence depths observed in the experiment described in this work (1.44–3.97 km).
In our models, two types of negative topography were observed in the forearc lithosphere: (1) a gradual increase in subsidence toward the trench and (2) superimposed on this first topographic signal a local depression striking parallel to the trench, which was most pronounced during the free sinking phase but gradually decreased during the steady state subduction phase. Through analyzing three forces acting on the forearc lithosphere and comparing our model results with previous investigations of overriding plate topography using geodynamic models, we conclude that the shear force and trench suction force control the first type (gradual, long wavelength) of forearc topographic subsidence, through bending the overriding plate at the trench. Furthermore, the downdip increase of the vertical component of the trench suction force along the subduction zone interface, caused by the progressive slab steepening during the free sinking phase, is the main driver in the formation of the forearc topographic depression. Finally, the vertical component of the subduction-induced mantle flow in the mantle wedge plays a minor role in driving the topographic subsidence observed in the forearc region of the overriding plate.
This project was supported by Discovery grants DP110103387 and DP120102983 from the Australian Research Council (ARC) awarded to W. P. S. Z. C. was supported by APA and IPRS scholarships from the Australian Government and by a JSPS Postdoctoral Fellowship. W. P. S. was supported by a Future Fellowship (FT110100560) from the ARC and a Vici fellowship from the Dutch Science Foundation (NWO). J. C. D. was supported by a DECRA fellowship from the ARC. J. C. D. also acknowledges FCT through project UID/GEO/50019/2013—Instituto Dom Luiz. We would like to thank two anonymous reviewers for their constructive comments that improved the manuscript, and we would like to thank the Editor Jeroen Ritsema for handling this manuscript. All data and methods necessary to understand, evaluate, replicate, and build upon the reported research are presented in the manuscript.
- 2008). Numerical modeling of lithosphere-asthenosphere interaction in a subduction zone. Earth and Planetary Science Letters, 272, 698–708. https://doi.org/10.1016/j.epsl.2008.06.009
- 2008). Effect of margin curvature on plate deformation in a 3-D numerical model of subduction zones. Geophysical Journal International, 173, 1084–1094. https://doi.org/10.1111/j.1365-246X.2008.03752.x
- 2010). Mantle structure and dynamic topography in the Mediterranean Basin. Geophysical Research Letters, 37, L20303. https://doi.org/10.1029/2010GL045001
- 2010). The many surface expressions of mantle dynamics. Nature Geoscience, 3, 825–833. https://doi.org/10.1038/ngeo1020
- 2005). Vertical motions of the Puerto Rico Trench and Puerto Rico and their cause. Journal of Geophysical Research, 110, B06404. https://doi.org/10.1029/2004JB003459
- 2001). A modeling study of vertical surface displacements at convergent plate margins. Geophysical Journal International, 147, 415–427. https://doi.org/10.1046/j.1365-246X.2001.00545.x
- 1997). Quantification of interplate coupling in subduction zones and forearc topography. Geophysical Research Letters, 24, 1563–1566. https://doi.org/10.1029/97GL01550
- 2015). Overriding plate deformation and variability of fore-arc deformation during subduction: Insight from geodynamic models and application to the Calabria subduction zone. Geochemistry, Geophysics, Geosystems, 16, 3697–3715. https://doi.org/10.1002/2015GC005958
- 2016). Does subduction-induced mantle flow drive backarc extension? Earth and Planetary Science Letters, 441, 200–210. https://doi.org/10.1016/j.epsl.2016.02.027
- 1993). Lithospheric buoyancy and collisional orogenesis—Subduction of oceanic plateaus, continental margins, island arcs, spreading ridges, and seamounts. Geological Society of America Bulletin, 105, 715–737.
- 2011). Tectonic geomorphometrics of the western United States: Speculations on the surface expression of upper mantle processes. Geochemistry, Geophysics, Geosystems, 12, Q11002. https://doi.org/10.1029/2011GC003579
- 2013). Three-dimensional dynamic laboratory models of subduction with an overriding plate and variable interplate rheology. Geophysical Journal International, 195, 47–66. https://doi.org/10.1093/Gji/Ggt257
- 2014). Rheology of petrolatum-paraffin oil mixtures: Applications to analogue modelling of geological process. Journal of Structural Geology, 63, 1–11. https://doi.org/10.1016/j.jsg.2014.02.004
- 2015). How weak is the subduction zone interface? Geophysical Research Letters, 42, 2664–2673. https://doi.org/10.1002/2014GL062876
- 2013). A review of observations and models of dynamic topography. Lithosphere, 5, 189–210. https://doi.org/10.1130/L245.1
- 2008). Trench migration, net rotation and slab-mantle coupling. Earth and Planetary Science Letters, 271, 233–240. https://doi.org/10.1016/j.epsl.2008.04.006
- 2010). Mantle flow and dynamic topography associated with slab window opening: Insights from laboratory models. Tectonophysics, 496, 83–98. https://doi.org/10.1016/j.tecto.2010.10.014
- 2013). The dynamics of laterally variable subductions: laboratory models applied to the Hellenides. Solid Earth, 4, 179–200. https://doi.org/10.5194/se-4-179-2013
- 1998). Cretaceous vertical motion of Australia and the Australian-Antarctic discordance. Science, 279, 1499–1504. https://doi.org/10.1126/science.279.5356.1499
- 2001). Residual topography, lithospheric structure and sunken slabs in the central Mediterranean. Earth and Planetary Science Letters, 187, 117–130. https://doi.org/10.1016/S0012-821X(01)00272-2
- 1985). Lower mantle heterogeneity, dynamic topography and the geoid. Nature, 313, 541–545.
- 2000). The influence of plate motions on three-dimensional back arc mantle flow and shear wave splitting. Journal of Geophysical Research: Solid Earth, 105, 28009–28033. https://doi.org/10.1029/2000JB900297
- 2006). Relative importance of trenchward upper plate motion and friction along the plate interface for the topographic evolution of subduction-related mountain belts. Geological Society of London-Special Publication, 253, 105. https://doi.org/10.1144/GSL.SP.2006.253.01.05
- 1997). Study of plate deformation and stress in subduction processes using two-dimensional numerical models. Journal of Geophysical Research, 102, 17,951–17,965.
- 2012). Influence of lithosphere-asthenosphere interaction on the overriding lithosphere in a subduction zone: Numerical modeling. Geochemistry, Geophysics, Geosystems, 13, Q02006. https://doi.org/10.1029/2011GC003909
- 2012). Unraveling topography around subduction zones from laboratory models. Tectonophysics, 526, 5–15. https://doi.org/10.1016/j.tecto.2011.09.001
- 1973). Model experiment of plate movements. Nature Physical Science, 242, 130–134. https://doi.org/10.1038/physci242130a0
- 2007). Trench-parallel flow and seismic anisotropy in the Mariana and Andean subduction systems. Nature, 450, 1222–1226. https://doi.org/10.1038/nature06429
- 2008). Effect of three-dimensional slab geometry on deformation in the mantle wedge: Implications for shear wave anisotropy. Geochemistry, Geophysics, Geosystems, 9, Q01003. https://doi.org/10.1029/2007GC001677
- 1998). Dynamic topography, plate driving forces and the African superswell. Nature, 395, 269–272. https://doi.org/10.1038/26212
- 2010). Analogue modeling of continental collision: Influence of plate coupling on mantle lithosphere subduction, crustal deformation and surface topography. Tectonophysics, 484, 87–102. https://doi.org/10.1016/j.tecto.2009.08.043
- 2014). The impact of slab dip variations, gaps and rollback on mantle wedge flow: Insights from fluids experiments. Geophysical Journal International, 197, 705–730. https://doi.org/10.1093/gji/ggu053
- 2013). Effect of aseismic ridge subduction on slab geometry and overriding plate deformation: Insights from analogue modeling. Tectonophysics, 588, 39–55. https://doi.org/10.1016/j.tecto.2012.12.010
- 2015). How do subduction processes contribute to forearc Andean uplift? Insights from numerical models. Journal of Geodynamics, 96, 6–18. https://doi.org/10.1016/j.jog.2015.04.001
- 2013). Three-dimensional dynamic models of subducting plate-overriding plate-upper mantle interaction. Journal of Geophysical Research: Solid Earth, 118, 775–790. https://doi.org/10.1002/jgrb.50078
- 1998). Mantle flow mechanisms for the large-scale subsidence of continental interiors. Geology, 26, 687–690.
- 2010). Bending mechanics and mode selection in free subduction: A thin-sheet analysis. Geophysical Journal International, 180, 559–576. https://doi.org/10.1111/j.1365-246X.2009.04460.x
- 2004). Topographic response to mantle lithosphere removal in the southern Sierra Nevada region, California. Geology, 32, 245–248. https://doi.org/10.1130/G19958.1
- 2008). Kinematics and flow patterns in deep mantle and upper mantle subduction models: Influence of the mantle depth and slab to mantle viscosity ratio. Geochemistry, Geophysics, Geosystems, 9, Q03014. https://doi.org/10.1130/B26360.1
- 2015). Australian plate motion and topography linked to fossil New Guinea slab below Lake Eyre. Earth and Planetary Science Letters, 421, 107–116. https://doi.org/10.1016/j.epsl.2015.03.036
- 2016). A review of analogue modeling of geodynamic processes: Approaches, scaling, materials and quantification, with an application to subduction experiments. Journal of Geodynamics, 100, 7–32. https://doi.org/10.1016/j.jog.2016.03.009
- 1993). Subduction of the lithosphere and back arc dynamics: Insights from physical modeling. Journal of Geophysical Research: Solid Earth, 98, 16,167–16,185. https://doi.org/10.1029/93JB01094
- 2011). Decoupling during continental collision and intra-plate deformation. Earth and Planetary Science Letters, 305, 435–444. https://doi.org/10.1016/j.epsl.2011.03.028
- 2014). Evolution of 3-D subduction-induced mantle flow around lateral slab edges in analogue models of free subduction analysed by stereoscopic particle image velocimetry technique. Earth and Planetary Science Letters, 403, 368–379. https://doi.org/10.1016/j.epsl.2014.07.007
- 2016). Mantle-induced subsidence and compression in SE Asia since the early Miocene. Geophysical Research Letters, 43, 1901–1909. https://doi.org/10.1002/2016gl068050
- 1992). Viscous flow model of a subduction zone with a faulted lithosphere: Long and short wavelength topography, gravity and geoid. Geophysical Research Letters, 19, 1891–1894.
- 1994). Controls on trench topography from dynamic models of subducted slabs. Journal of Geophysical Research: Solid Earth, 99, 15,683–15,695.