Overriding Plate Velocity Control on Surface Topography in 2‐D Models of Subduction Zones

• Published in: Geochemistry, geophysics, geosystems : G3 geophysics, geosystems : G3, 2020-04, Vol.21 (4), p.n/a
• Main Authors: Cerpa, Nestor G., Arcay, Diane
• Format: Article
• Language: English
• Subjects: dynamic topography; Earth Sciences; Free surfaces; Geologic depressions; Geophysics; Kinematics; Ocean surface topography; overriding plate kinematics; Physiographic features; Sciences of the Universe; Subduction; subduction dynamics; Subduction zones; Tectonics; Topography; Upper mantle; Velocity
• ISSN: 1525-2027; 1525-2027
• DOI: 10.1029/2019GC008900

We study the dynamically induced flexural topography in subduction numerical mechanical models. We focus on the topographic changes at the overriding plate (OP) surface induced by variations in OP kinematics, particularly when the subducting plate (SP) has a stationary motion after having reached the rigid base of the upper mantle. Our models consist of two viscoelastic plates with free surfaces and an isoviscous mantle. Friction is imposed along the planar subduction interface. We first characterize the main topographic features at a constant OP velocity, using spatial definitions based on geometrical estimations of the volcanic arc position. The models exhibit the formation of a bulge in the forearc area followed landwards by a depression and a smaller second bulge, both bracketing the arc region. The steady‐state distance to the trench of these features increases with OP velocity. Their amplitude is affected by the far‐field OP tectonic regime that depends on kinematics, and plates and subduction interface strength. We next test the effect of sudden changes in OP velocity. An OP acceleration yields a transient topographic tilt, during which the outer forearc quickly subsides whereas the arc region uplifts, and that is followed by reverse slower motions. An OP slowdown induces opposite motions. The rates of elevation change during the tilt are approximately proportional to velocity variations and mainly sensitive to the SP strength. The rates are higher than 0.1 mm/yr for velocity changes higher than 1 cm/yr. We suggest that topographic accommodations of OP velocity changes should be considered when quantifying nonisostatic topography. Key Points The topography amplitude and position on top of the overriding plate are controlled by its kinematics and the plate and interface strengths Changes in overriding plate velocity yield transient strong (rates ∼0.1 mm/yr) surface vertical motions from the trench to the arc region The strength of the subducting plate has the strongest impact on the topographic response after a change in overriding plate kinematics