Flat-slab subduction, topography, and mantle dynamics in southwestern Mexico

Topography above subduction zones arises from the isostatic contribution of crustal and lithospheric buoyancy, as well as the dynamic contribution from slab-driven mantle flow. We evaluate those effects in southwestern Mexico, where a segment of the Cocos slab subducts horizontally. The eastern part of the volcanic arc - the Trans-Mexican Volcanic Belt - stands at an average elevation of 2.3 km, nearly 1.3 km above the fore-arc. Lateral changes in bulk crustal density are relatively small, and seismic imaging shows that there is little variation in crustal thickness between these two regions. Thus, the elevation difference between the arc and the fore-arc should arise from differences in mantle properties. We present finite element models of flat-slab subduction that provide a simultaneous match to topography, plate velocities, and stress state in the overriding plate. We find that the dynamic effects are primarily controlled by the amount of coupling at the subduction interface and in the mantle wedge, the lack of slab anchoring into the lower mantle, and the absence of continental mantle lithosphere. With a mantle wedge and a subduction interface that are, respectively, 2 and 4 orders of magnitude weaker than the asthenosphere, the flat slab exerts a downward pull that can explain most of the elevation difference between the fore-arc and the arc. We infer that lateral viscosity variations play a significant role in shaping dynamic topography in complex tectonic settings and that sublithospheric dynamics can influence the topography at wavelengths that are significantly shorter than previously recognized.