[1] Changes of plate motion may have induced subduction initiation (SI), but the tectonic history of SI is different from one subduction zone to another. Izu-Bonin-Mariana (IBM) SI, accompanied by strong backarc spreading and voluminous eruption of Boninites, contrasts with the Aleutians which shows neither. Using finite element models, we explore visco-elasto-plastic parameters and driving boundary conditions for SI evolution. With an imposed velocity, we find three different evolutionary modes of SI: continuous without backarc spreading, continuous with backarc spreading and a segmented mode. With an increase in the coefficient of friction and a decrease in the rate of plastic weakening, the amount of convergence needed for SI increases from $20 to $220 km, while the mode changes from segmented to continuous with backarc spreading and eventually to continuous without backarc spreading. If the imposed velocity boundary condition is replaced with an imposed stress, the amount of convergence needed for SI is reduced and backarc spreading does not occur. These geodynamic models provide a basis for understanding the divergent geological pathways of SI. First, IBM evolution is consistent with subduction of an old strong plate with an imposed velocity which founders causing intense backarc spreading and Boninitic volcanism. Second, the New Hebrides SI is in the segmented mode due to its weak plate strength. Third, the Puysegur SI is in the continuous without backarc spreading mode with no associated volcanic activities. Fourth, the Aleutians SI has neither trench rollback nor backarc spreading because the slab is regulated by constant ridge-push forces.