How impurity atoms move through a crystal is a fundamental and renewed issue in condensed matter physics and materials science. Diffusion of oxygen (O) in titanium (Ti) affects the formation of titanium-oxides and the design of Tibased alloys. Moreover, the kinetics of initial growth of titania-nanotubes via anodization of a titanium metal substrate also involves the diffusion of oxygen. Therefore, the understanding of the migration mechanism of oxygen atoms in -Ti is extremely important for controlling oxygen diffusion in Ti alloys. In this work, we show how the diffusion coefficient can be predicted directly from first-principles studies without any empirical fitting parameters. By performing the first-principles calculations based on the density functional theory (DFT) through using the Vienna ab initio Simulation Package (VASP), we obtain three locally stable interstitial oxygen sites in the hexagonal closed-packed (hcp) lattice of titanium. These sites are octahedral center (OC) site, hexahedral center (HE) site, and TiTi bond center crowdion (CR) site with interstitial energies of -2.83, -1.61, and -1.48 eV, respectively. From the interstitial energies it follows that oxygen atom prefers to occupy the octahedral site. From electronic structure analysis, it is found that the TiO bonds possess some covalent characteristics and are strong and stable. Using the three stable O sites from our calculations, we propose seven migration pathways for oxygen diffusion in hcp Ti and quantitatively determine the transition state and diffusion barrier with the saddle point along the minimum energy diffusion path by the climbing image nudged elastic band (CI-NEB) method. The microscopic diffusion barriers (E) from the first-principles calculations are important for quantitatively describing the temperature dependent diffusion coefficients D from Arrhenius formula D = L2v* exp(-((E)/(kBT)), where v* is the jumping frequency and L is the atomic displacement of each jump. The jumping frequency v* is determined from where vi and vj are the vibration frequency of oxygen atom at the initial state and the transition state respectively. This analysis leads to the formula for calculating the temperature dependent diffusion coefficient by using the microscopic parameters (vi and E) from first-principles calculations without any fitting parameters. Using the above formula and the vibration frequencies and diffusion barriers from first-principles calculations, we calculate the diffusion coefficients among different interstitial sites. It is found that the diffusion coefficient from the octahedral center site to the available site nearby is in good agreement with the experimental result, i.e., the diffusion rate D is 1.046510-6 m2s-1 with E of 0.5310 eV. The jump from the crowdion site to the octahedral interstitial site prevails over all the other jumps, as a result of its low energy barrier and thus leading to markedly higher diffusivity values. The diffusion of oxygen atoms is mainly controlled by the jump occurring between OC and CR sites, resulting in high diffusion anisotropy. This finding of oxygen diffusion behavior in Ti provides a useful insight into the kinetics at initial stage of oxidation in Ti which is very relevant to many technological applications of Ti-based materials.