Directional solidification is a common and important process in both scientific research and industrial practice. Owing to the presence of temperature gradients during directional solidification, local remelting and solidification in the mushy zone occurs, resulting in some typical phenomena such as temperature gradient zone melting (TGZM). The TGZM influences the solidifying microstructure and microsegregation significantly. In the present work, a two-dimensional (2D) cellular automaton (CA) model involving the mechanism of both solidification and melting is adopted to investigate the migration phenomena of molten liquid pools in the mushy zone due to the TGZM. The effect of pulling velocity, initial liquid pool position, temperature gradient, and alloy composition on the TGZM kinetics are studied. The simulation results are compared with the analytical predictions, and good agreement between two models is obtained. It is found that under a temperature gradient, the liquid pool always migrates towards the high temperature direction. When the pulling velocity is lower than the critical velocity, the liquid pool migrates through the liquidus into the bulk liquid and the time required for a liquid pool to reach the liquidus increases with pulling velocity increasing. On the other hand, when a pulling velocity higher than the critical value is adopted, the liquid pool moves towards the solidus and the time required for migrating liquid pool to reach the solidus decreases with pulling velocity increasing. For a given pulling velocity, the liquid pools located above the critical position move towards the liquidus, while the others gradually approach to the solidus. When a molten liquid pool migrates towards the liquidus, the migration velocity and liquid pool thickness are found to gradually increase, while the liquid pool composition decreases with time. Inversely, for the molten liquid pool that moves towards the solidus, the migration velocity and liquid pool thickness gradually decrease, while the liquid pool composition increases with time going by. The average migration velocity of liquid pool caused by the TGZM effect increases with temperature gradient increasing and alloy composition decreasing. The CA simulations provide an insight into the complicated interactions among the local temperature, solute distribution and diffusion, and the kinetics of local remelting and solidification in the TGZM process.