Progress in numerical modeling of subducting plate dynamics

The core concerns of plate tectonics theory are the dynamics of subducting plates, which can be studied by integrating multidisciplinary fields such as seismology, mineral physics, rock geochemistry, geological formation studies, sedimentology, and numerical simulations. By establishing a theoretical model and solving it with numerical methods, one can replicate the dynamic effects of a subducting plate, quantifying its evolution and the surface response. Simulations can also explain the observations and experimental results of other disciplines. Therefore, numerical models are among the most important tools for studying the dynamics of subducting plates. This paper provides a review on recent advances in the numerical modeling of subducting plate dynamics. It covers various aspects, namely, the origin of plate tectonics, the initiation process and thermal structure of subducting slab, and the main subduction slab dynamics in the upper mantle, mantle transition zone, and lower mantle. The results of numerical models are based on the theoretical equations of mass, momentum, and energy conservation. To better understand the dynamic progress of subducting plates, the simulation results must be verified in comparisons with the results from natural observations by geology, geophysics and geochemistry. With the substantial increase in computing power and continuous improvement of simulation methods, numerical models will become a more accurate and efficient means of studying the frontier issues of Earth sciences, including subducting plate dynamics.