<figure><img src="image_1.png" /></figure>

HAL
open science

# Major Forest Changes in Subtropical China since the Last Ice Age

Qiuchi Wan, Xiao Zhang, Yaze Zhang, Yuanfu Yue, Kangyou Huang, Rachid

Cheddadi, Zhuo Zheng

► To cite this version:

Qiuchi Wan, Xiao Zhang, Yaze Zhang, Yuanfu Yue, Kangyou Huang, et al. Major Forest Changes in Subtropical China since the Last Ice Age. Forests, 2021, 12 (10), pp.1314. 10.3390/f12101314. hal-04934522

HAL Id: hal-04934522

https://hal.science/hal-04934522v1

Submitted on 7 Feb 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Major Forest Changes in Subtropical China since the Last Ice Age

Qiuchi Wan<fnref n="1" />,*, Xiao Zhang<fnref n="1" />,*, Yaze Zhang<fnref n="1" />, Yuanfu Yue<fnref n="2" />, Kangyou Huang<fnref n="1" />, Rachid Cheddadi<fnref n="3" />
and Zhuo Zheng<fnref n="1" />,<fnref n="4" />

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Wan, Q.; Zhang, X.; Zhang, Y.; Yue, Y.; Huang, K.; Cheddadi, R.; Zheng, Z. Major Forest Changes in Subtropical China since the Last Ice Age. *Forests* **2021**, *12*, 1314. https://doi.org/10.3390/f12101314

Academic Editor: Juan A. Blanco

Received: 27 August 2021

Accepted: 22 September 2021

Published: 26 September 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** In the subtropical zone of southern China, there was a considerable conversion of forests from deciduous to evergreen broadleaf in the early Holocene. However, the exact timing of this vegetation change and its relationship to climate are still unclear. We examined a high-resolution pollen record collected in the mid-subtropical zone and then performed a correlation with regional data to reconstruct the history of forest ecosystems since the last deglaciation. Our data show that the expansion of the evergreen plant component already occurred at low elevations during the last deglaciation. The subtropical mountain landscape was not recolonized by evergreen forests until the mid-Holocene at about 8.1 ka BP. Based on fossil pollen reconstruction and climate model simulation, we conclude that the primary increase in evergreen components of subtropical ecosystems was triggered by postglacial temperature increase, and that a complete conversion from deciduous to evergreen forest ecosystems did not occur until Holocene winter temperatures and seasonal temperature contrast reached a threshold suitable for the growth and persistence of evergreen tree species.

**Keywords:** deglaciation; Holocene; pollen analysis; species distribution model; subtropical China; winter temperature

## 1. Introduction

Evergreen broadleaved forests (EBLFs) are now the dominant biome in the subtropical zone in east China with abundant endemic plant species and high biodiversity [1]. The EBLF extends from about 22–23° N to Yangtze River at 31–33° N [2,3]. During the Last Glacial Maximum (LGM), subtropical China was covered by temperate deciduous (TEDE) and cool mixed (COMX) forests in a cooler climate [4,5]. Reconstructed LGM biomes suggest that tropical forests were absent from mainland southern China, or their range was reduced to very limited areas [6]. EBLF and warm mixed forest had retreated to tropical latitudes [7]. Therefore, it is likely that the EBLF was restricted to the southernmost part of China, between 6 and 8°N south of its present range [6–8].

However, exactly when and how the EBLF recolonized subtropical regions since the last deglaciation is poorly known [5,9]. Moreover, the climatic interpretation of forest conversion remains controversial. Although long-term global warming since the last deglaciation is generally considered to be the main cause of forest changes and shifts in zonal vegetation boundaries [10,11], some have claimed that summer rainfall might have played a dominant role [10]. The mechanism of low latitude vegetation response to global

<footnote n="1">¹ Guangdong Key Lab of Geodynamics and Geohazards, School of Earth Sciences and Engineering, Sun Yat-sen University, Zhuhai 519082, China; wanqch@mail3.sysu.edu.cn (Q.W.); zhangyz6@mail2.sysu.edu.cn (Y.Z.); hkangy@mail.sysu.edu.cn (K.H.); eeszzhuo@mail.sysu.edu.cn (Z.Z.)</footnote>

<footnote n="2">² Guangxi Laboratory on the Study of Coral Reef in the South China Sea, Coral Reef Research Center of China, School of Marine Sciences, Guangxi University, Nanning 530004, China; yuanfu.yue@gxu.edu.cn</footnote>

<footnote n="3">³ Institute of Evolutionary Sciences, University of Montpellier, CNRS, IRD, EPHE, 34000 Montpellier, France; rachid.cheddadi@umontpellier.fr</footnote>

<footnote n="4">⁴ Southern Marine Science and Engineering Guangdong Laboratory, Zhuhai 519082, China</footnote>

<footnote>* Correspondence: zhangx636@mail.sysu.edu.cn</footnote>
climate change driven by ice sheet mass and/or insolation is also debated [12,13]. In this study, we examine a high-resolution pollen record of SZY peat bog in the mid-subtropical zone. This pollen record allowed us to reconstruct the timing of the conversion from the LGM deciduous biome to the now dominant EBLF. Climate thresholds for EBLFs were derived from its modern distribution using a species distribution model and past climate variables from transient climate simulations. A reconstruction of biome turnover and climate data helped us identify the key driving climate variables and the effects of post-glacial global warming on EBLF's recolonization of evergreen ecosystems in subtropical China.

## 2. Regional Setting

Study core SZY (119.034° E, 26.777° N; 1007 m) was collected from a mountain peat bog in Jiufeng Mountains, Fujian Province, southeastern China (Figure 1). The summits of nearby mountains are between 1300 and 1600 m. The area is located in the southern part of the mid-subtropical zone, which is strongly influenced by the East Asian Summer Monsoon (EASM). The local mean annual temperature is about 16 °C, and the mean January temperature and daily minimum are 6-7 °C and -4 °C, respectively, at about 1000 m elevation. The average altitude constant is about 0.59 °C/100 m. Annual precipitation varies between 1700 and 1900 mm, with most precipitation falling between late March and July. Evergreen broadleaved forest is the main vegetation type, which can be divided into the following five communities: coniferous/broadleaved mixed forest, warm coniferous forest, evergreen broadleaved forest, bamboo forest and shrubland. In the regions below 1000 m altitude, there are few natural forests left due to heavy human impact. Dominant arboreal plants at the genus level in the local natural or secondary forests of the mountains are *Pinus*, *Castanopsis*, *Cyclobalanopsis*, *Cunninghamia*, *Altingia*, *Cryptomeria*, *Cinnamomum*, *Machilus*, *Schima* and *Eleocarpus*. The common shrub taxa include *Eurya*, *Neolitsea*, *Lindera*, *Ardisia*, *Rhododendron*, *Prunus*, etc. The temperate deciduous and conifer elements, such as *Fagus*, *Quercus*, *Carpinus*, *Acer*, *Tsuga*, *Rhododendron*, *Syplocos*, are generally found at higher altitudes of 1200 m.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Location of the study core SZY.</figcaption></figure>
# 3. Materials and Methods

## 3.1. *Coring, Sampling and Chronology*

The SZY core was taken in 2011 using a Russian Corer. In this study, the upper 273 cm of the peat core (425 cm in total) was investigated. The uppermost sediment (0–32 cm) consists of marsh clay with abundant plant debris and high water content. The main part of the core below 32 cm consists of homogeneous brown-black peat and dark clay. For the upper 273 cm of core, a total of 180 samples (including 132 previously published [5] and 48 newly added) were taken at 1–2 cm intervals, corresponding to a temporal resolution of about 80 years between samples. The age-depth model was constructed using radiocarbon dating with accelerator mass spectrometry (AMS $^{14}$C). Details of radiocarbon samples can be referred to [5]. Based on the calibrated age, the interval extends to 15,000 cal yr BP at 274 cm. The depth-age plot [5] shows a higher sedimentation rate during the Holocene (2.2 cm/100 yr) and a lower rate during the last deglaciation (0.83 cm/100 yr).

## 3.2. *Pollen Analysis*

The pollen preparation method basically followed the traditional laboratory treatment techniques [14] with some improvements such as the use of hydrofluoric acid. Each sample was weighed to 2–3 g of sediment, and one tablet of trace Lycopodium spores (27,637 per tablet) was used for calculation of pollen concentration (grains/g). The main treatment procedure included the hydrochloric acid (10%) to remove carbonates, hydrofluoric acid (40%) to remove silicates and potassium hydroxide to remove organic matter. The palynomorphs were separated from mineral sediment using standard heavy liquid flotation techniques [15]. In general, more than 400 pollen grains were counted, and pollen percentage calculation was based on sum of total pollen excluding aquatic taxa and all spores.

## 3.3. *Biomization*

In order to better illustrate temporal vegetation changes, the pollen taxa were grouped as plant functional types (PFTs). Based on which, the biomes were reconstructed by applying the biomization technique. We basically applied the existing schemes of plant functional types (PFTs) previously used for China [7], except for the exclusion of PFTs from high latitude arctic and boreal zones, and some moderate modifications to fit the regional biomes. In this study, only the following six biomes occurring in the study area were reconstructed: subtropical evergreen broadleaved forest (EBLF), tropical rain forest (TRFO), temperate deciduous broadleaved forest (TEDE), cold mixed forest (CLMX), cool mixed forest (COMX) and alpine shrub and meadow (ALPM). We named the biome of subtropical (warm-temperate) evergreen broadleaved forest as 'EBLF', which is approximately the same with the warm evergreen and mixed forest (WAMF) used in [7] and adopted the biome of ALPM proposed by [4] to identify subtropical mountain vegetation.

## 3.4. *Species Distribution Modelling*

By associating modern presence data and climate data, the species distribution model (SDM) [16] predicts the climatic envelope of the study species or vegetation type and projects it to different time periods to obtain potential distributions. Using this method, the potential distributions of EBLF at the time slices of LGM and the present were predicted. The presence data of EBLF was extracted from a digitalized vegetation map (Supplementary Figure S1). To avoid sampling bias, we divided the distribution range of presence data into several 1° × 1° grids, and randomly picked up only one point from each grid to obtain a filtered subset of presence data (n = 129) using the GridSample function of the R 'dismo' package. Seven key parameters among 19 bioclimatic variables representative of climate conditions from the WorldClim database were used for model calibration (Supplementary Table S1). For the LGM climate conditions, we relied on the results estimated using an Earth System Model, MIROC-ESM [17]. The distribution of EBLF at each time slice was finally predicted using one of the most popular SDM methods, the maximum entropy algorithm
(Maxent) within the R 'dismo' package with default settings based on the filtered presence
data and climate data. The resulting habitat suitability, the value of which ranges from 0
(unsuitable) to 1, was used to indicate the potential distribution of EBLF. By comparing the
predicted distribution model for the present time and the known distribution of EBLF, the
model was validated with a threshold-independent measure of the area under the receiver
operating characteristic curve (AUC). The highest value of 1 (0 being the lowest) represents
a complete match between the modeled distribution and the known distribution.

3.5. *Transient Data Analysis*

TRACE21 is the first state-of-the-art transient simulation of global climate change over the last 21 ka at the T31 spatial resolution [18]. Two climatic variables including mean winter (December–February) temperatures (corresponding to bioclimatic variable 11 in WorldClim) and temperature seasonality (July–January) were used for comparison. We first extracted the temperature values of the grid where SZY is located, and then uniformly correct the simulated values by substituting in the bias caused by the altitudinal lapse-rate of modern climate observation in the study site.

4. Results

The high resolution and well dated record of SZY shows that the vegetation of the last
15,000 years passed through four main phases characterized by four distinct pollen zones
(Figure 2). During the last deglaciation between 15 and 11.5 ka BP, regional vegetation was
dominated by cold-tolerant deciduous forest. The constituent deciduous taxa are *Fagus*,
*Quercus*, *Alnus*, *Betula*, *Corylus*, *Carpinus* and others. Evergreen taxa such as *Cyclobalanopsis*
and *Castanopsis* were also present in this pollen zone, but the percentages are generally
lower than the subsequent Holocene values. Conifers such as *Tsuga* and *Abies* and the
shrub *Rhododendron* (Ericaceae) were more abundant during the time interval earlier than
11.5 ka BP. Herbaceous taxa including *Artemisia* and *Poaceae* fluctuated and sometimes
reached 5% (Figure 2).

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Pollen diagram of selected taxa from the SZY record of the last 15,000 years. All relative abundances were calculated based on total pollen sum.</figcaption></figure>
During the Holocene, at least three zones of pollen assemblage can be identified. The lower zone is from 11.5 to 8.2 ka BP when *Pinus* and Taxodiaceae (possibly Chinese water-cypress *Glyptostrobus*) and wetland sedge increased rapidly. *Alnus* was also relatively abundant, the species of which may belong to the wetland environment taxon in the region (i.e., *A. trabeculosa*). The middle zone dated to 8.1–0.8 ka BP corresponds to the pollen assemblages marked by evergreen broadleaved taxa mainly including *Cyclobalanopsis*, *Castanopsis* and *Altingia*. The wetland sedge pollen of Cyperaceae was extremely abundant in the lower part of this zone, which was followed by an increase in Poaceae and fern spores. The upper most zone corresponding to the last ~0.8 ka BP shows an abrupt increase in Poaceae, *Artemisia*, *Dicranopteris* and *Pinus*, which are considered as the pollen taxa associated with human activities [19]. The total non-arboreal pollen percentage is up to 52% in the uppermost part of the record.

The biome reconstruction result shows a strong dominance of temperate deciduous broadleaved forest (TEDE) between 15 and 11.5 ka BP, and the biome with the second highest score was the cool mixed forest (COMX). An obvious decrease in all biomes except alpine shrub and meadow (ALPM) happened at around 11.5–12.5 ka BP. During the transition from early to mid-Holocene (ca. 11.5–8.2 ka BP), there was a gradual increase in subtropical (warm-temperate) evergreen broadleaved mixed forest (EBLF), although its scores were still lower than those of TEDE.

The onset of the EBLF increase after the obvious biome change at ca. 11.5 ka BP corresponds to the cold Younger Dryas (YD), suggesting that the EBLF's recolonization of the subtropical regions resumed and accelerated after the YD termination. After 8.1 ka BP, EBLF became dominant over TEDE, with several fluctuations of abundance (e.g., at around 3.2 ka BP) until the end of the Holocene (Figure 3).

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Reconstructed biomes of the last 15 ka BP derived from the pollen record SZY in southeastern China. The principal biomes shown are subtropical evergreen broadleaved forest (EBLF), tropical rain forest (TRFO), temperate deciduous broadleaved forest (TEDE), cold mixed forest (CLMX), cool mixed forest (COMX) and alpine shrub and meadow (ALPM).</figcaption></figure>

The species distribution model was used to predict the potential distributions of evergreen broadleaved forest at the LGM and in the present (Figure 4). The predicted present-day distribution of EBLF agrees well with the known distribution (AUC = 0.94), indicating that the model is reliable for drawing further conclusions. The simulated LGM distribution shows that EBLF migrated southward into northern Indochina in the western
part but found potentially suitable conditions (in situ refugia) in the central east (i.e., some large river basins between the Wuyi and Xuefeng mountain ranges).

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Potential distribution of evergreen broadleaved forest (EBLF) in southern China simulated by Maxent. The quantile result of critical climate thresholds for the simulated present distribution is shown below. The suitability value ranges from 0 (unsuitable) to 1. The green to blue distribution ranges represent the predicted most probable climatic-suitable habitats for the EBLF.</figcaption></figure>

# 5. Discussion

## 5.1. *Substantial Forest Transformation during the Holocene*

Our data show that the vegetation in subtropical mountains around 1000 m of elevation during the last deglaciation was dominated by deciduous forest (TEDE) and conifers (COMX). The presence of evergreen components might reflect a kind of mixed forest in which deciduous trees predominated, but there were also a significant number of evergreen taxa. It is also possible that some of the pollen of the evergreen plants was transported from lower altitudes by updraft wind, since the evergreen forest belt in the lower mountains was not very far from the study area during the last deglacial. However, considering the small size of our study site, we suppose that the pollen assemblage of our record largely reflects local vegetation composition [20]. The main transition from a deciduous to an evergreen biome occurred rapidly over a narrow time span at ca. 8.1 ka BP (Figure 3). However, this timing of biome turnover raised controversy as pollen records from other latitudinal positions show different timing. In Taiwan Island, the expansion of subtropical (and/or warm temperate) evergreen forest took place much earlier at about 11.5 ka BP [4]. Fossil pollen data over eastern China show that the timing of vegetation change is not simultaneous, but rather follows a time-transgressive pattern from the south (10.3 ka BP) to the north (6 ka BP) [11].

Effectively, the complete absence of the evergreen broadleaved forest during the LGM in the subtropical zone of the Chinese mainland inferred by previous biome reconstruction [6] is under debate as phylogeographic studies have suggested that species of the EBLF survived the LGM locally in numerous northern refugia [21]. Our pollen record shows that the evergreen forest may have occupied the lower altitude regions even during the glacial periods, taking into consideration the regional altitudinal air temperature lapse rate, which may have steepened during the LGM [22].
Our simulated distribution of evergreen forest shows a shrinkage of evergreen forest distribution, which also persisted in nearby lowland areas as refugia such as the basins between the Wuyi and Xuefeng mountain ranges in southeastern China. Such model simulations imply that the local mountain evergreen forest experienced a complex migration upward from lowland during the early to mid-Holocene, agreeing well with the refugia theory proposed by phylogeographic studies [21,23], although there may have also been some northward shift of constituent evergreen species from the south to the north [21].

## 5.2. Holocene Climatic Changes as a Cause of Subtropical Forest Transformation

The rapid EBLF/TEDE biome change at about 8.1 ka BP suggests that the rate of vegetation change increased as temperature increased and crossed a threshold that delineates the climatic range of the two biomes. Mean minimum temperature and seasonal temperature contrast are commonly considered as effective climate indices for defining the upper limit for the evergreen broadleaved forest [24], as photosynthetic efficiency leads to a reduction in evergreen broadleaved species [25]. From the analysis of Holocene climate changes in the region (Figure 5), we further confirm that the seasonal temperature difference and/or the coldest climate condition must be the most important climatic factors controlling the switch between evergreen/deciduous trees, and that it is much less related to the change of rainfall amount, because both the moisture proxy from stalagmites [26] and the simulated precipitation [27] show a maximum amount centered around 10 ka BP, which contradicts the timing of our pollen records.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Temperature and seasonal difference as the key factors of evergreen forest turnover during the Holocene. Temperature indices including brGDGT-based mean annual temperature reconstruction [27], TRACE21 [18] modeled winter temperature (Twin) and seasonal mean temperature difference between July and January. The values of the two variables were revised from the regional output according to the local air temperature (SZY).</figcaption></figure>

The local air temperature was estimated by the proxy of branched glycerol dialkyl glycerol tetraether (brGDGT) [28], which shows synchronous changes with the reconstructed vegetation changes from the same site of SZY. There was a rapid increase in the
mean annual temperature (MAT) from 13 to 16 °C at about 8 ka BP. The temperature at about 16 °C corresponds to the most appropriate MAT range (14.5–18.5 °C) for the present evergreen forest based on the simulated distribution (Figure 4). This suggests that a local forest change from a deciduous-dominated forest to an evergreen-dominated forest may have been triggered by the temperature change up to the threshold. Note also that the subsequent climate fluctuations after 8 ka BP to the threshold MAT value of about 14.5 °C (e.g., at 3.2 ka BP) repeatedly caused a switch between evergreen/deciduous biomes. This further supports the significant role of temperature threshold in the biome turnover. The vegetation changes in the studied pollen record are closely related to the output of the TRACE21 climate model [18], e.g., the trends in mean winter temperature (Twin) and the seasonal temperature difference between July and January (STdif) from the last deglaciation to the Holocene (Figure 5). The modeled winter temperature in the subtropical climate zone shows a gradual change from 2 to 5 °C in the early Holocene at 11.5–8.1 ka BP (more abrupt since 9 ka BP). Similarly, there was a gradual decrease in the seasonal temperature difference at 11.5–8.1 ka BP from 24.5 to less than 23 °C (Figure 5). Note that these values were corrected to the altitude of the SZY record based on the local average elevation lapse rate. These values at 8.1 ka BP agree very well with the present-day range of the most suitable temperature for evergreen trees (Twin: 5.5–10.5 °C, STdif: 14–23.5 °C). They are also in broad agreement with the meteorological data delineating the modern evergreen broadleaved forest in the subtropical mountains [24,25]. Therefore, there is no doubt that evergreen trees became dominant in the forest when the winter and seasonal thermal conditions reached their thresholds. After the YD at 11.5 ka BP, progressive winter warming and decreasing seasonality led to an increase in evergreen broadleaved elements in the local forest, and the subsequent onset of the Holocene optimum at ~8.1 ka BP forced the general forest turnover.

According to our record, we suggest that the forest change from deciduous to evergreen composition was initiated after the YD in the subtropical region. The greatest extent of EBLF in our study area corresponds to the Holocene optimum, which begins at ~8 ka BP and shows a close relationship with the increase in the winter temperature and the decrease in the seasonal temperature difference to thresholds. Our data are consistent with winter climate dynamics influencing the strength of the Asian summer monsoon [29]. This study provides an improved understanding of the ecologically sensitive region at the interface between the warm temperate and tropical bioclimatic zones. Ice sheet melting exerts an important influence on low latitude climate through associated changes in the global sea level [30], and the sea surface temperature of Indo-Pacific warm pools shows a progressive warming since the beginning of the Holocene [31]. The observed amplitude of subtropical ecosystem changes during global climate warming over the last postglacial period draws our attention to the period when one extended biome with new species replaced another. The biome turnover we identified occurred over a time span comparable to the few decades of ongoing global climate change. Thus, ongoing climate change will likely have a significant effect on ecosystems, especially upland species [32], although the impact is expected to be less significant in upland areas than in lowland areas [33]. Many endemic species, e.g., *Fagus*, *Cryptomeria*, *Tsuga*, *Pseudotsuga*, *Quercus*, in subtropical Asia are near threatened, vulnerable, endangered or even critically endangered according to the IUCN Red List [34]. These species are all fragmented and scattered within the evergreen biome. Ongoing warming is likely to lead to their extinction, resulting in a depletion of biodiversity in Asian ecosystems, which deserve the establishment of new national parks and biosphere reserves.

## 6. Conclusions

This study shows that both ecosystem composition and their range have undergone major changes since the last ice age, especially in the last thousand years. The observed biome transformation in our study site appears to have taken place within a time span comparable to that of ongoing climate warming, providing important knowledge to better
understand ecosystem transformations and manage their future changes. Many endemic evergreen tree species are now threatened with extinction, which also raises many questions regarding the managing of diversity in subtropical ecosystems.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/f12101314/s1, Figure S1: The distribution of the presence points extracted from digitalized EBLF distribution used for predicting the modern distribution of EBLF, Table S1: Estimates of relative contributions of the environmental variables to the Maxent model.

**Author Contributions:** Conceptualization, Q.W., X.Z. and Z.Z.; methodology, Q.W., X.Z., Y.Z., K.H. and Z.Z.; writing—original draft preparation, Q.W., X.Z. and Z.Z.; writing—review and editing, Q.W., X.Z., Y.Z., Y.Y., K.H., R.C. and Z.Z. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the National Natural Science Foundation of China (Grant Nos. 42072205, 41472143, 41902186, 42077414, 41630753, 41472143, 41661144003 and 41702182), Innovation Group Project of Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) (Grant No. 311020002), Project of the influence of Quaternary glacial cycle on biodiversity in Nanling Mountains (Grant No. 2021GJGY001) and the Guangxi scientific projects (No. 2018GXNSFAA281293).

**Data Availability Statement:** Data is contained within the article or supplementary material.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Song, Y.-C.; Chen, X.-Y.; Wang, X.-H. Studies on evergreen broad-leaved forests of china: A retrospect and prospect. *J. East China Norm. Univ.* **2005**, *1*, 1–8.

2. Wu, Z.Y. *Vegetation of China*; Science Press: Beijing, China, 1980.

3. Fang, J.Y.; Wang, Z.H.; Tang, Z.Y. *Atlas of Woody Plants in China: Distribution and Climate*; Higher Education Press: Beijing, China, 2011.

4. Lee, C.-Y.; Liew, P.-M. Late quaternary vegetation and climate changes inferred from a pollen record of Dongyuan Lake in southern Taiwan. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **2010**, *287*, 58–66. [CrossRef]

5. Yue, Y.; Zheng, Z.; Huang, K.; Chevalier, M.; Chase, B.M.; Carré, M.; Ledru, M.-P.; Cheddadi, R. A continuous record of vegetation and climate change over the past 50,000 years in the Fujian province of eastern subtropical China. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **2012**, *365*, 115–123. [CrossRef]

6. Ni, J.; Yu, G.; Harrison, S.P.; Prentice, I.C. Palaeovegetation in China during the late quaternary: Biome reconstructions based on a global scheme of plant functional types. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **2010**, *289*, 44–61. [CrossRef]

7. Yu, G.; Chen, X.; Ni, J.; Cheddadi, R.; Guiot, J.; Han, H.; Harrison, S.P.; Huang, C.; Ke, M.; Kong, Z.; et al. Palaeovegetation of china a pollen data-based synthesis for the mid-holocene and last glacial maximum. *J. Biogeogr.* **2000**, *27*, 635–664. [CrossRef]

8. Zheng, Z.; Yuan, B.; Nicole, P.-M. Paleoenvironments in China during the last glacial maximum and the holocene optimum. *Episodes* **1998**, *21*, 152–158.

9. Ni, J.; Cao, X.; Jeltsch, F.; Herzschuh, U. Biome distribution over the last 22,000yr in China. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **2014**, *409*, 33–47. [CrossRef]

10. Wang, X.; Yao, Y.-F.; Wortley, A.H.; Qiao, H.-J.; Blackmore, S.; Wang, Y.-F.; Li, C.-S. Vegetation responses to the warming at the younger dryas-holocene transition in the hengduan mountains, southwestern china. *Quat. Sci. Rev.* **2018**, *192*, 236–248. [CrossRef]

11. Zhou, X.; Sun, L.; Zhan, T.; Huang, W.; Zhou, X.; Hao, Q.; Wang, Y.; He, X.; Zhao, C.; Zhang, J.; et al. Time-transgressive onset of the Holocene Optimum in the East Asian monsoon region. *Earth Planet. Sci. Lett.* **2016**, *456*, 39–46. [CrossRef]

12. Liu, Z.; Zhu, J.; Rosenthal, Y.; Zhang, X.; Otto-Bliesner, B.L.; Timmermann, A.; Smith, R.S.; Lohmann, G.; Zheng, W.; Elison Timm, O. The holocene temperature conundrum. *Proc. Natl. Acad. Sci. USA* **2014**, *111*, E3501–E3505. [CrossRef]

13. Cheng, Y.; Liu, H.; Wang, H.; Hao, Q. Differentiated climate-driven Holocene biome migration in western and eastern China as mediated by topography. *Earth-Sci. Rev.* **2018**, *182*, 174–185. [CrossRef]

14. Faegri, K.; Iversen, J. *Textbook of Pollen Analysis*; Wiley: New York, NY, USA, 1989.

15. Nakagawa, T. A trial of density pollen fractionation from sediment core samples for the purpose of pollen AMS dating: Evaluation and perspectives. *Summ. Res. Using AMS Nagoya Univ.* **1998**, *9*, 244–252.

16. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum entropy modeling of species geographic distributions. *Ecol. Model.* **2006**, *190*, 231–259. [CrossRef]

17. Hijmans, R.J.; Cameron, S.E.; Parra, J.L.; Jones, P.G.; Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. *Int. J. Climatol.* **2005**, *25*, 1965–1978. [CrossRef]
18. Liu, Z.; Otto-Bliesner, B.L.; He, F.; Brady, E.C.; Tomas, R.; Clark, P.U.; Carlson, A.E.; Lynch-Stieglitz, J.; Curry, W.; Brook, E.; et al. Transient simulation of last deglaciation with a new mechanism for bolling-allerod warming. *Science* **2009**, *325*, 310–314. [CrossRef] [PubMed]

19. Zheng, Z.; Deng, Y.; Zhang, H.; Yu, R.; Chen, Z. Holocene environmental changes in the tropical and subtropical areas of the south China and the relation to human activities. *Quat. Sci.* **2004**, *24*, 387–393.

20. Sugita, S. Pollen representation of vegetation in quaternary sediments: Theory and method in patchy vegetation. *J. Ecol.* **1994**, *82*, 881–897. [CrossRef]

21. Fan, D.; Hu, W.; Li, B.; Morris, A.B.; Zheng, M.; Soltis, D.E.; Soltis, P.S.; Zhang, Z. Idiosyncratic responses of evergreen broad-leaved forest constituents in China to the late Quaternary climate changes. *Sci. Rep.* **2016**, *6*, 31044. [CrossRef]

22. Loomis, S.E.; Russell, J.M.; Verschuren, D.; Morrill, C.; Cort, G.D.; Damsté, J.S.S.; Olago, D.; Eggermont, H.; Street-Perrot, F.A.; Kelly, M.A. The tropical lapse rate steepened during the last glacial maximum. *Sci. Adv.* **2017**, *3*, e1600815. [CrossRef]

23. Wan, Q.; Zheng, Z.; Huang, K.; Guichoux, E.; Petit, R.J. Genetic divergence within the monotypic tree genus platycarya ( Juglandaceae) and its implications for species’ past dynamics in subtropical China. *Tree Genet. Genomes* **2017**, *13*, 73. [CrossRef]

24. Fang, J.; Yoda, K. Climate and vegetation in China v. Effect of climatic factors on the upper limit of distribution of evergreen broadleaf forest. *Ecol. Res.* **1991**, *6*, 113–125. [CrossRef]

25. Song, Y. Perspective of the vegetation zonation of forest region in eastern China. *Acta Bot. Sin.* **1999**, *41*, 541–552.

26. Wang, Y.J.; Cheng, H.; Edwards, R.L.; An, Z.S.; Wu, J.Y.; Shen, C.-C.; Dorale, J.A. A high-resolution absolute-dated late pleistocene monsoon record from hulu cave, China. *Science* **2001**, *294*, 2345–2348. [CrossRef] [PubMed]

27. Cheng, J.; Wu, H.; Liu, Z.; Gu, P.; Wang, J.; Zhao, C.; Li, Q.; Chen, H.; Lu, H.; Hu, H.; et al. Vegetation feedback causes delayed ecosystem response to east Asian summer monsoon rainfall during the holocene. *Nat. Commun.* **2021**, *12*, 1843. [CrossRef]

28. Wang, M.; Zheng, Z.; Man, M.; Hu, J.; Gao, Q. Branched gdgt-based paleotemperature reconstruction of the last 30,000 years in humid monsoon region of southeast China. *Chem. Geol.* **2017**, *463*, 94–102. [CrossRef]

29. Wang, L.; Li, J.; Lu, H.; Gu, Z.; Rioual, P.; Hao, Q.; Mackay, A.W.; Jiang, W.; Cai, B.; Xu, B.; et al. The east Asian winter monsoon over the last 15,000 years: Its links to high-latitudes and tropical climate systems and complex correlation to the summer monsoon. *Quat. Sci. Rev.* **2012**, *32*, 131–142. [CrossRef]

30. DiNezio, P.N.; Tierney, J.E. The effect of sea level on glacial indo-pacific climate. *Nat. Geosci.* **2013**, *6*, 485–491. [CrossRef]

31. Loarie, S.R.; Duffy, P.B.; Hamilton, H.; Asner, G.P.; Field, C.B.; Ackerly, D.D. The velocity of climate change. *Nature* **2009**, *462*, 1052–1055. [CrossRef] [PubMed]

32. Lenoir, J.; Gégout, J.C.; Marquet, P.A.; Ruffray, P.D.; Brisse, H. A significant upward shift in plant species optimum elevation during the 20th century. *Sci. Adv.* **2008**, *320*, 1768–1771. [CrossRef]

33. Bertrand, R.; Lenoir, J.; Piedallu, C.; Riofrio-Dillon, G.; de Ruffray, P.; Vidal, C.; Pierrat, J.C.; Gégout, J.C. Changes in plant community composition lag behind climate warming in lowland forests. *Nature* **2011**, *479*, 517–520. [CrossRef] [PubMed]

34. International Union for Conservation of Nature. The IUCN Red List of Threatened Species; Version 2021-1. 2021. Available online: https://www.iucnredlist.org (accessed on 24 August 2021).