The recrystallization behavior of cold-rolled, commercial-purity titanium was studied experimentally and with Monte-Carlo (MC) modeling. Utilization of EBSD-OIM as input for MC modeling resulted in realistic predictions of recrystallization kinetics, microstructure and texture, which were in good agreement with experimental results. MC modeling of recrystallization kinetics predicted that the non-uniform stored energy distribution, heterogeneous nucleation of recrystallization and recovery in combination leads to a negative deviation from linear JMAK kinetics. It was found that concurrent recovery that takes place during recrystallization is an important process that controls both the overall recrystallization kinetics and the deviation of linear JMAK kinetics. On the other hand, the non-uniformly distributed stored energy itself has little effect on the negative deviation from JMAK kinetics but intensifies the deviation when heterogeneous nucleation is combined. Modeling results also revealed that heterogeneous nucleation of recrystallized grains and their early impingement in local areas of high deformation are essential for producing a log-normal distribution of grain size and a typical recrystallization texture of rolled titanium.