Factors Contributing to Slab Locations and Geometries in Reconstructions of Past Mantle Flow

Individual sinking slabs present markedly different geometries between 410 and 660 km depths, from vertical slabs penetrating the lower mantle to slabs stagnating above the lower mantle. The proposed factors determining these contrasted geometries include mantle viscosity and the magnitude and evolution of trench retreat. Here, we assess the success of paleo-geographically driven global mantle flow models in matching slabs in tomographic models between 400 km and 1,000 km depth. We quantify the spatial match between predicted present-day mantle temperature anomalies and vote maps of tomographic models. We investigate the sensitivity of the spatial match to input parameters of the mantle flow model: imposed tectonic reconstruction, model start age, and viscosity contrast between the upper and lower mantle. We evaluate the visual match between model slabs and tomographic vote maps for three circum-Pacific regions with contrasted slab dip angles between 400 km and 1,000 km depth. Predicted model slabs better match slabs inferred from tomography when there is an increase in viscosity at 660 km depth. The temporal evolution of the models and the global match at present day suggest that the subduction history could be refined in the global tectonic reconstructions that we considered. For example, we suggest that the subduction to the east of Japan should be offset by approximately 100 km to the west at ∼80 Ma to match the anchoring of a continuous slab into the lower mantle suggested by tomography.