This study, using species distribution modeling (involving a new approach that allows for uncertainty), predicts the distribution of climatically suitable areas prevailing during the mid-Holocene, the Last Glacial Maximum (LGM), and at present, and estimates the potential formation of new habitats in 2070 of the endangered and rare Tertiary relict tree Davidia involucrata Baill. The results regarding the mid-Holocene and the LGM demonstrate that south-central and southwestern China have been long-term stable refugia, and that the current distribution is limited to the prehistoric refugia. Given future distribution under six possible climate scenarios, only some parts of the current range of D. involucrata in the mid-high mountains of south-central and southwestern China would be maintained, while some shift west into higher mountains would occur. Our results show that the predicted suitable area offering high probability (0.5‒1) accounts for an average of only 29.2% among the models predicted for the future (2070), making D. involucrata highly vulnerable. We assess and propose priority protected areas in light of climate change. The information provided will also be relevant in planning conservation of other paleoendemic species having ecological traits and distribution ranges comparable to those of D. involucrata.