Numerical modeling of subduction zones and implications for mantle convection