The effect of variable thermal diffusivity on kinematic models of subduction