Background: Subalpine ecosystems at high altitudes and latitudes are particularly sensitive to climate change. In South Korea, the prediction of the species richness of subalpine plant species under future climate change is not well studied. Thus, this study aims to assess the potential impact of climate change on species richness of subalpine plant species (14 species) in the 17 mountain national parks (MNPs) of South Korea under climate change scenarios' representative concentration pathways (RCP) 4.5 and RCP 8.5 using maximum entropy (MaxEnt) and Migclim for the years 2050 and 2070. Results: Altogether, 723 species occurrence points of 14 species and six selected variables were used in modeling. The models developed for all species showed excellent performance (AUC > 0.89 and TSS > 0.70). The results predicted a significant loss of species richness in all MNPs. Under RCP 4.5, the range of reduction was predicted to be 15.38-94.02% by 2050 and 21.42-96.64% by 2070. Similarly, under RCP 8.5, it will decline 15.38-97.9% by 2050 and 23.07-100% by 2070. The reduction was relatively high in the MNPs located in the central regions (Songnisan and Gyeryongsan), eastern region (Juwangsan), and southern regions (Mudeungsan, Wolchulsan, Hallasan, and Jirisan) compared to the northern and northeastern regions (Odaesan, Seoraksan, Chiaksan, and Taebaeksan). Conclusions: This result indicates that the MNPs at low altitudes and latitudes have a large effect on the climate change in subalpine plant species. This study suggested that subalpine species are highly threatened due to climate change and that immediate actions are required to conserve subalpine species and to minimize the effect of climate change.