[1] Numerical simulations of mantle convection with semidynamic plate subduction in 3-D geometry is performed to investigate the role of harzburgite layers in the morphology of subducting plates and the behavior of oceanic crustal layers. The results show that a "buckled" crustal layer is observed under the "heel" of the stagnant slab that begins to penetrate into the lower mantle, regardless of the magnitude of the viscosity of the harzburgite layer, Δη hrz , when the factor of viscosity increase at the boundary of the upper and lower mantle, Δη lm , is larger than 60-100. As Δη hrz increases, the curvature of buckling is larger. When Δη lm and Δη hrz are larger, the volumes of crustal and harzburgite materials trapped in the mantle transition zone (MTZ) are also larger, although almost all of the materials penetrate into the lower mantle. These materials are trapped in the MTZ for over tens of millions of years. Citation: Yoshida, M. (2013), The role of harzburgite layers in the morphology of subducting plates and the behavior of oceanic crustal layers, Geophys. Res. Lett., 40,[5387][5388][5389][5390][5391][5392]