Molecular dynamic simulations were conducted to study the self-interstitial migration in zirconium. By defining the crystal lattice point, at which more than one atom fall in the Wigner-Seitz cell of the lattice point, for the location of interstitial atoms (LSIA), three types of events were identified for LSIA migration: the jump remaining in one 1120 direction (ILJ), the jump from one 1120 to another 1120 direction in the same basal plane (OLJ) and the jump from one basal plane to an adjacent basal plane (OPJ). The occurrence frequencies of the three types were calculated. ILJ was found to be the dominant event in the temperature range (300K to 1200K), but the occurrence frequencies of OLJ and OPJ increased with increasing temperature. Although the three types of jumps may not follow Brownian and Arrhenius behavior, on the whole, the migration of the LSIAs tend to be Brownian-like. Moreover, the migration trajectories of LSAs in the hcp basal-plane are not what are observed if only conventional one-or two-dimensional migrations exist; rather, they exhibit the feature we call fraction-dimensional. Namely, the trajectories are composed of line segments in 1120 directions with the average length of the line segments varying with the temperature. Using Monte Carlo simulations, the potential kinetic impacts of the fraction-dimensional migration, which is measured by the average number of lattice sites visited per jump event (denoted by nSP E ), was analyzed. The significant differences between the nSP E value of the fraction-dimensional migration and those of the conventional one-and two-dimensional migrations suggest that the conventional diffusion coefficient, which cannot reflect the feature of fraction-dimensional migration, cannot give an accurate description of the underlying kinetics of SIAs in Zr. This conclusion may not be limited to the SIA migration in Zr and could be more generally meaningful for situations in which the low dimensional migration of defects has been observed.