Climate fluctuations in the past and in the future are likely to result in population expansions,shifts,or the contraction of the ecological niche of many species,and potentially leading to the changes in their geographical distributions.Prediction of suitable habitats has been developed as a useful tool for the assessment of habitat suitability and resource conservation to protect wildlife.Here,we model the ancestral demographic history of the extant modern Chinese Muntjac Muntiacus reevesi populations using approximate Bayesian computation (ABC) and used the maximum entropy model to simulate the past and predict the future spatial dynamics of the species under climate oscillations.Our results indicated that the suitable habitats for the M.reevesi shifted to the Southeast and contracted during the Last Glacial Maximum,whereas they covered a broader and more northern position in the Middle Holocene.The ABC analyses revealed that the modern M.reevesi populations diverged in the Middle Holocene coinciding with the significant contraction of the highly suitable habitat areas.Furthermore,our predictions suggest that the potentially suitable environment distribution for the species will expand under all future climate scenarios.These results indicated that the M.reevesi diverged in the recent time after the glacial period and simultaneously as its habitat\'s expanded in the Middle Holocene.Furthermore,the past and future climate fluctuation triggered the change of Chinese muntjac spatial distribution,which has great influence on the Chinese muntjac\'s population demographic history.