Vegetation Response to Holocene Climate Change in the Qinling Mountains in the Temperate–Subtropical Transition Zone of Central–East China

Global warming is having a profound influence on vegetation and biodiversity patterns, especially in alpine areas and high latitudes. The Qinling Mountain range is located in the transition zone between the temperate and subtropical ecosystems of central–east China and thus the vegetation of the area is diverse. Understanding the long-term interactions between plant diversity and climate change can potentially provide a reference for future landscape management and biodiversity conservation strategies in the Qinling Mountains region. Here, we use a pollen record from the Holocene sediments of Daye Lake, on Mount Taibai in the Qingling Mountains, to study regional vegetation changes based on biomes reconstruction and diversity analysis. Temperature and precipitation records from sites close to Daye Lake are used to provide environmental background to help determine the vegetation response to climate change. The results indicate that climate change was the main factor influencing vegetation and palynological diversity in the Qinling Mountains during the Holocene. The cold and dry climate at the beginning of the early Holocene (11,700–10,700 cal yr BP) resulted in a low abundance and uneven distribution of regional vegetation types, with the dominance of coniferous forest. During the early Holocene (10,700–7,000 cal yr BP), temperate deciduous broadleaf forest expanded, palynological diversity and evenness increased, indicating that the warm and humid climate promoted vegetation growth. In the middle Holocene (7,000–3,000 cal yr BP), the climate became slightly drier but a relatively warm environment supported the continued increase in palynological diversity. After ∼3,000 cal yr BP, palynological diversity and the evenness index commenced a decreasing trend, in agreement with the decreased temperature and precipitation in the Qinling Mountains. It’s noteworthy that human activity at this time had a potential influence on the vegetation. During the past few centuries, however, palynological diversity has increased along with the global temperature, and therefore it is possible that in the short-term ongoing climatic warming will promote vegetation development and palynological diversity in the area without human interference.


INTRODUCTION
Climate change has a major impact on the global environment and biota, including humans (Shaver, 2000;Root et al., 2003;Peng et al., 2004;Parmesan, 2007). Global temperatures have risen substantially since the last glacial maximum (Xiao et al., 2020), especially during the past few decades, due to increasing atmospheric greenhouse gas concentrations. As an important component of global biodiversity, plant diversity plays an essential role in maintaining ecosystem stability and environmental protection (Ehrlich and Wilson, 1991;Thomas et al., 2004). However, plant diversity is currently undergoing rapid changes due to climatic warming and the increased frequency of extreme climatic events, as well as the reduction of living habitats (Walther et al., 2002;Wake and Vredenburg, 2008;Barnosky et al., 2011). However, several studies have shown that increased temperatures are conducive to plant growth at high latitudes (Peuelas and Filella, 2001), while at the same time, the vegetation at high latitudes and altitudes is more sensitive to climate change (Pepin et al., 2015;Ellis et al., 2021). Human activities also have important direct impacts on terrestrial ecosystems (Vitousek, 1994;Crutzen and Stoermer, 2000;Crutzen, 2002;Mottl et al., 2021). Humans altered the land surface, ecosystem dynamics and biodiversity by clearing the native vegetation for agriculture and living space, which has led to increased species extinctions, with rates accelerating during the past few centuries and with potentially irreversible effects (Braatz et al., 1992;Pimm et al., 2014;Lewis and Maslin, 2015). Therefore, the long-term paleovegetation and paleoclimatic records is of great significance to explore regional patterns of vegetation and climate change to distinguish the effects of climate change and human activity on vegetation. Such records may also provide a potential analogue of the interaction of climate change and vegetation in the future.
Pollen analysis has been widely used to study the history of vegetation and plant diversity on different timescales (Jackson and Blois, 2015;Birks et al., 2016a,b;Felde et al., 2020), and fossil pollen records can provide important information about the response of vegetation to past changes in climate and land use (Willis et al., 2004(Willis et al., , 2010Rudaya et al., 2020). The Qinling Mountain range is the natural geographic and climatic boundary between North and South China. The area is influenced mainly by the East Asian summer monsoon and the vegetation is highly sensitive to climate change (Liu et al., 2008;Zhang et al., 2021). Mount Taibai, which is the highest peak in the Qinling Mountain range, has numerous alpine lakes, and the pollen assemblages in the sediments of these lakes are a valuable archive of information on climate change and vegetation succession. Previous studies of the sediments have included analyses of biological indicators such as pollen and diatoms to reconstruct changes in vegetation cover and limnological conditions in the Mount Taibai area (Zhao et al., 1999;Zhang et al., 2001;Cui et al., 2003;Li et al., 2008Li et al., , 2017Cheng et al., 2017Cheng et al., , 2020Huang et al., 2020). Additionally, organic geochemical indicators have been used to characterize climate change on different timescales (Liu et al., 2005;Wang et al., 2016Wang et al., , 2019. However, these studies were not conducted on open-water lake sediments and were focused on the vegetation of the middle-late Holocene and the modern vegetation; moreover, they were low temporal resolution. Detailed investigations of the response of plant diversity to climate change based on highresolution continuous sediments is important for biodiversity conservation in the area. Lake sediments have the advantages of high stratigraphic continuity and resolution, which make them well suited for studies of regional climate and vegetation changes Wagner et al., 2000;Rioual et al., 2001). Daye Lake is an enclosed alpine lake above the forest line on Mount Taibai, and here we report the results of a high-resolution biomes and palynological diversity reconstruction from the Holocene sediments of the site. The principal aims of the study were as follows: (1) To investigate the Holocene evolution of vegetation and climate based on the pollen records from Daye Lake in the Qinling Mountains region; (2) to reveal changes in the palynological diversity and infer the plant diversity changes; and (3) to explore the relationship between vegetation and climate change, to provide a basis for predicting possible future changes in the vegetation and plant diversity of the region.

Study Area
Mount Taibai, with a maximum altitude of 3,767 m a.s.l., is the highest part of the Qinling Mountain range and the highest mountain in central and eastern China. The lithology is granite formed by the large-scale incursion of acidic magma during the Yanshan Movement in the Cretaceous (Tian and Huang, 1990). There are numerous Quaternary relict glacial landforms in the Mount Taibai area, with periglacial landforms generally developed in areas below the snow line of ∼2,600 m a.s.l. (Pang, 1998). However, the snow line has risen above Baxiantai, which is the highest peak of Mount Taibai (3,767 m a.s.l.).
There are pronounced differences in precipitation, temperature, topography, and vegetation distribution on the southern and northern sides of the Mount Taibai. Moreover, diverse temperate and subtropical plant species, as well as numerous medicinal herbs, survive in the region (Pang, 1998). The combination of thermal and hydrological conditions results in a vertical distribution of vegetation belts in the area. The area below 780 m a.s.l. is extensively cultivated and the natural vegetation has been substantially altered by human activities. Above this cultivated zone, the following vegetation zones on the northern slope of Mount Taibai Figure 1) is a small, closed glacier lake with a surface area of 0.5 ha and maximum water depth of 18 m, formed after the retreat of the Quaternary glacier activity, located above the tree line on the northern slope of Mount Taibai. It is mainly supplied FIGURE 1 | (A) Location of the study site (DY: Daye Lake) and other sites (S1-S4) in the surrounding region that are mentioned in the text (S1: Daihai Lake; S2: Gonghai Lake; S3: Qinghai Lake; S4: Dajiuhu wetland). The vegetation zones in China (https://geodata.pku.edu.cn) are also shown; (B) digital elevation map of the study area; (C) the picture of the Daye Lake.
by precipitation which is concentrated in summer. The regional climate is influenced by the East Asian summer monsoon, and the annual mean air temperature and precipitation are 8.1 • C and 720 mm, respectively. The coldest month (January) and warmest month (July) have respective average temperatures of −4 and 19.3 • C, based on a 30-year meteorological dataset (CE 1981(CE -2010 from Taibai county meteorological station (1,545 m a.s.l) (see Figure 2 in Zhang et al., 2021). Based on the lapse rates for the northern slopes of the Qinling Mountains proposed by Ma et al. (2017), the mean annual, January and July temperatures are 2.8, −7.3, and 13.6 • C, respectively. Daye Lake is surrounded by alpine shrubland and meadow. The shrub taxa are dominated by Salix cupularis, Rhododendron capitatum, and Spiraea, and the local herb taxa are mainly Ranunculaceae, Brassicaceae, and Poaceae (Zhang et al., 2021).

Sediments and Dating
A 2.35 m-long sediment core (DYH16A) was retrieved from Daye Lake in April 2016 using a UWITEC piston corer. The undisturbed core was transported to the laboratory and subsampled at 1 cm-intervals. Seventeen samples of bulk organic sediments were analyzed by accelerator mass spectrometry (AMS) at Arizona University, United States, and Peking University, China ( Table 1). The ages of bulk sediment samples may be influenced by old carbon due to the complexity of carbon sources. Therefore, we used the intercept method after fitting a 2nd-order polynomial curve to get the carbon reservoir age which is comparable to the age difference between pollen concentrates and bulk organic from the surface sample (see the Supplementary Material). Finally, 1540 years was selected as the reservoir age of Daye Lake. A chronological model (Figure 2) was developed using the R package rbacon, version 2.5.0 (Blaauw and Christen, 2011) based on calendar ages before present (BP = 1950 CE) and calibrated with the IntCal20 dataset (Northern Hemisphere) (Reimer et al., 2020). Note that only the top 1.26 m of the core were discussed in this study.

Pollen Analysis
The top 1.26 m of core DYH16A at 1 cm-intervals were analyzed in this study. A known number of Lycopodium spores were added to 0.5 g of dried sediment to calculate pollen concentrations (grains g −1 ). Samples were then treated sequentially with HCl, KOH, and HF, following the procedures described in Faegri and Iversen (1989) and Moore et al. (1991). Pollen identification was carried out under a light microscope at 400× magnification with reference keys (e.g., Wang et al., 1995;Tang et al., 2016). At least 500 terrestrial pollen grains were counted each sample. Pollen percentages were calculated based on the sum of terrestrial taxa, excluding aquatic taxa and spores. Pollen diagrams were draw and constrained incremental sum of squares (CONISS) analysis was conducted using Tilia software of version 2.0 (Grimm, 1987).
A lower-resolution pollen record from the top 1.03 m of the core DYH16A has been published by Zhang et al. (2021) to reconstruct the regional plant abundances based on the REVEALS model. This record spans the past ∼7000 years which FIGURE 2 | AMS 14 C dates for core DYH16A from Daye Lake based on the R package rbacon, version 2.5.0 (Blaauw and Christen, 2011) and the resulting age-depth model. The red line represents mean calendar years. The gray shading and the blue line indicate calendar ages with the 95% confidence range and the age errors, respectively. the sediment accumulation rate was relatively uniform. The higher-resolution pollen record presented here spans the entire Holocene and extends to the base of the core where the sediment accumulation rate was lower. But the entire core is continuously deposited which can reflect the Holocene changes of vegetation and climate in the region.

Biome Reconstruction
The pollen record from Daye Lake was used for biomes reconstruction in order to explore the vegetation response to Holocene climate change. The procedures used for biomization analysis follow those described in Prentice et al. (1996). First, pollen types (taxa) were assigned to one or more plant functional types (PFTs) ( Table 2), and the PFTs were then classified into biomes (Table 3). Finally, the affiliated scores for each biome were calculated based on the pollen percentages and the PFT sets and were then used to represent the abundance of a particular biome. Assignment of the identified pollen taxa from Daye Lake to PFTs and biomes follows Yu et al. (2000) and Ni et al. (2014) for China, together with reference to the modern regional vegetation distribution. Biome definition depends on the influence of the dominant PFTs with the highest scores. Biome affinity scores provide an objective evaluation of the vegetation composition and have been applied to several surface and fossil pollen datasets (Tarasov et al., 1998;Yu et al., 2000;Rudaya et al., 2008;Ni et al., 2014). Biome modeling was conducted using PPPBase software (Guiot and Goeury, 1996), based on three matrixes (PFT × taxon, biome × PFT and biome × taxon matrix) which comprised 51 pollen types, 13 PFTs and 5 biomes in this study (Tables 2, 3).

Pollen-Based Plant Diversity
The reconstructions of palynological diversity can be used to directly interpret the composition and structure of the vegetation. In order to reconstruct the palynological diversity history over the Holocene in Qingling Mountains region from different aspects, Hill's indices (Hill, 1973) and rarefaction analysis (Hurlbert, 1971) are used to estimate the diversity richness including palynological richness, Shannon diversity, and Simpson diversity. Palynological richness is the simple count of the number of taxa which represents the total number of taxa in a community to describe the vegetation community; the number of common taxa is known as the Shannon diversity index; and the number of dominant taxa is a measure of the Simpson diversity (Rudaya et al., 2020;Kulkarni et al., 2021). We also calculated evenness index to reflect the community structure and to assess palynological diversity at the regional scale and variations in vegetation composition (Pielou, 1966). The value of evenness is between 0 (only one species in the pollen assemblage) and 1 (each species with the same abundance in the pollen assemblage), the low value indicates the greater difference in the abundance of the plant types in pollen assemblages (Li and Zhao, 2018). All of these indices were calculated using the "vegan" package, version 2.5-7, in R (Oksanen et al., 2020).

Pollen Assemblages in the Holocene Sediments of Daye Lake
The pollen assemblages from Daye Lake (Figure 3) were conducted with CONISS analysis and divided into five significant zones, i.e., three pollen zones and two subzones, suggested by broken stick model. The zones are described below.
Zone I (11,700-10,700 cal yr BP). This zone is characterized by the highest percentages of coniferous trees (average of 42.18%). Among them, Pinus, Abies, and Picea have high percentages, with averages of 27.60, 9.18, and 3.23%, respectively, while Tsuga has the lowest (average of 2.17%). Deciduous broadleaf trees are relatively poorly represented (30. 65%). Among the shrub taxa, Ephedra is well represented (average of 1.26%), while Rosaceae, Fabaceae, and Euphorbiaceae are poorly represented. Herbs are dominated by Artemisia which has an average of 21.45%. The pollen concentration is the lowest during the entire Holocene.
Zone II (10,700-3,800 cal yr BP). This zone is characterized by the lowest proportion of coniferous trees (average of 22.66%) among the three zones and the significant increase of deciduous broadleaf trees (from 30.65 to 48.62%). The obvious increase of Corylus/Carpinus (with average values of 18.51%) and evergreen oak (Quercus_evergreen) are at the expense of the decrease of Pinus (from 27.60 to 9.88%). However, the herb pollen proportion  decrease, which is largely due to the decreases in Artemisia (from 9.31 to 5.30%) and Amaranthaceae (from 2.18 to 1.32%), whereas shrub taxa such as Rosaceae (from 2.16 to 5.66%) and Fabaceae (from 0.20 to 1.33%) increase. The pollen concentration is relatively high in this zone. Owing to the different trends of Pinus, Abies, Corylus/Carpinus, Fagus, and Artemisia within this zone, two subzones, i.e., Zone IIa (10,700-6,000 cal yr BP) and Zone IIb (6,000-3,800 cal yr BP) is further divided. Zone IIa shows the lowest abundance of coniferous trees (20.49%) and the highest value of deciduous broadleaf trees (51.48%) in the Holocene. In Zone IIb, Artemisia reach the lowest percentages in the entire record (averages of 5.04%). Zone III (3,800-0 cal yr BP). There is a minor increase in the percentages of coniferous trees and herbs (average of 26.87 and 23.51%, respectively). The abundance of deciduous broadleaf trees and shrubs shows a decreasing trend (average of 40.80 and 6.10%, respectively). This zone can also be further divided into Zone IIIa (3,800-1,200 cal yr BP) and Zone IIIb (1,200-0 cal yr BP). Zone IIIa is characterized by continuous decrease in deciduous broadleaf trees, increase in coniferous trees and herbaceous plant (which is mainly due to increases in Artemisia and Amaranthaceae). However, Zone IIIb is characterized by the increase of coniferous trees, decrease of deciduous broadleaf trees, and a substantial increase in herbs, especially in Artemisia (from 7.99 to 13.62%, the maximum value over the Holocene) and Poaceae (from 2.04 to 3.09%). The pollen concentration in Zone IIIa is high whereas the one in Zone IIIb is low.

Biome Reconstruction From the Holocene Sediments of Daye Lake
The biome reconstruction from Daye Lake (Figure 4) shows that the coniferous forests (CO) vegetation type displays the highest affinity scores during the first stage of the early Holocene (11,700-10,700 cal yr BP). This vegetation reconstruction suggests drought-and cold-tolerant needles communities played a major role in the regional vegetation cover. However, the temperate deciduous broadleaf forests (TEDE) become dominated during the Holocene interval after which suggests a warm and moist environment during the Holocene. There is a pronounced inverse relationship between TEDE and CO throughout the Holocene.
The Holocene period experiences drastic changes in the biome composition, characterized by a marked increase in TEDE, herbaceous vegetation (STEP) and warm mixed forest (WAMF) scores, which is accompanied by a decrease in the CO scores. The drastic decrease of CO (which is corresponding to the apparent reduce of Pinus, see Figure 3), and the visible increase of TEDE suggests a warm and moist environment during the period of∼10,700 and 6,000 cal yr BP. The period of 6,000-3,800 cal yr BP is characterized by the minimum scores of CO, decreased scores of TEDE, and the lowest scores of STEP, but a visible increased scores of WAMF. CO increase slightly compared with the previous stage from 3,800 to 1,200 cal yr BP. The increasing trends are observed in herbs and WAMF. During the past ∼1200 years, WAMF decrease and STEP increase from low to high values. There is an obvious low value of TEDE between 1000 and 600 years and then the scores increase significantly, whereas conifers show the opposite trend.

Changes in Palynological Diversity Reconstructed From the Holocene Sediments of Daye Lake
Diversity analysis based on pollen records from lake sediments can accurately reflect the richness and composition of the regional vegetation (Stirling and Wilsey, 2001;Birks et al., 2016a). The Shannon and Simpson diversity indexes show similar temporal trends and are positively correlated with evenness ( Figure 5). The palynological richness shows an increasing trend, with pronounced fluctuations, throughout the Holocene. Therefore, Shannon diversity and evenness index are employed in the following discussion. During 11,700-10,700 cal yr BP, the palynological diversity increases substantially, corresponding to a trend of increasing evenness. Palynological richness, diversity and evenness index are in increasing trends until 6,000 cal yr BP; however, there is an interval of pronounced low palynological diversity and evenness index during 9,000-8,200 and 7,800-7,000 cal yr BP based on the raw data (gray line in Figure 5). Palynological diversity increases substantially from 6,000 to 3,800 cal yr BP, and then decreases slightly between 3,800 and 1,200 cal yr BP. The decreasing trend is observed from ∼1,200 cal yr BP, which is followed by another increase during the past ∼600 years.

Vegetation Response to Holocene Climate Change in the Qinling Mountains
Climate change is a major factor influencing the distribution and abundance of vegetation (Bartlein et al., 1986;Ma et al., 2013); however, the response of vegetation and plant diversity to climate change shows strong regional variations. For example, temperature is the major influence on diversity in high-altitude mountains and high-latitude regions (Hou, 1983;Zhang et al., 2020), while precipitation is more important in arid and semiarid regions such as the Loess Plateau of China (Shang and Li, 2010;Sun et al., 2013). Elsewhere, such as in southwest China, temperature and precipitation both have major influences on the vegetation Trivedi et al., 2020). In the case of east-central China, the Qingling Mountains region, the relationship between climate change and vegetation over the Holocene need to be investigated.
The Qinling Mountain range is a typical region influenced by the East Asian summer monsoon, where plant growth is influenced predominantly by temperature and precipitation (Gansert, 2004;Zhang et al., 2006). The response of the vegetation at Daye Lake to climate change is manifested by changes in community composition and in the relative abundance of the dominant taxa (Figures 3, 4). Coniferous trees such as Abies and Picea can tolerate low temperatures and drought, and generally grow at high altitudes (Achuff and Roi, 1977;Liu Z. L. et al., 2002). By contrast, deciduous broadleaf forest (e.g., Quercus_deciduous, Carpinus_Corylus, and Ulmaceae) is favored by a warm and moist environment Xu et al., 2019;Trivedi et al., 2020). Evergreen trees such as Quercus_evergreen and Cyclobalanopsis are widely distributed in subtropical and tropical regions, especially in locations with high temperatures and precipitation (Srivastava et al., 2021). Herbaceous plants such as Artemisia and Amaranthaceae are drought-tolerant (Kolb and Sperry, 1999;Uotila et al., 2021).
In order to determine the vegetation response at Daye Lake to Holocene climate change, we compared the pollen assemblages, biome scores, Shannon diversity and evenness index with paleoclimatic records from other sites in the region. The cited records in the study comprise the following: monsoon precipitation and temperature reconstructions based on pollen records from Gonghai and Daihai Lakes, both to the north of the Qinling Mountains (Xu et al., 2010;; temperature and precipitation reconstructions based on records of branched fatty alcohol ratio BNA 15 and pollen assemblages from the Dajiuhu wetland to the south of the Qinling Mountains, respectively (Huang et al., 2013;Sun et al., 2019); alkenonebased temperature and pollen-based precipitation records from Qinghai Lake, to the northwest of the Qinling Mountains (Herzschuh et al., 2010;Hou et al., 2016; Figure 6).
In the first stage of the early Holocene (11,700-10,700 cal yr BP), coniferous forests with Pinus, Abies and Picea were abundant in the study area, as well as steppe habitats dominated by Artemisia. Deciduous broadleaf taxa were poorly represented in the vegetation (Figures 3, 4). The dominance of coniferous forest indicates that the climate during this period was cold and dry, which is confirmed by the paleoclimatic record from sites elsewhere in the region (Figure 6).
During the early Holocene (10,700-7,000 cal yr BP), there were significant decreases in Pinus and Artemisia, and evergreen trees such as Quercus_evergreen began to appear (Figure 3), and there was a substantial increase in deciduous broadleaf and warm mixed forests. The pollen percentages and biome results indicate that the region was occupied by deciduous broadleaf forests, suggesting that the climate in the early Holocene was warm and humid. At the same time, temperature and precipitation commenced an increasing trend and reached their highest levels in the region of the southern and northern China, and the northeastern Qinghai-Tibet Plateau (Figure 6). However, the climate change at this stage is not stable, and there are fluctuations on a centennial time scale. For example, the abundance of deciduous broadleaf taxa decreased at 9,000-8,200 cal yr BP and 7,800-7,000 cal yr BP in Daye Lake, revealing deteriorating climate, while it corresponded to the declining of temperature and precipitation in Qinghai Lake and Dajiuhu wetland, respectively. And the Gonghai Lake in the northern China also has a low precipitation.
In the middle Holocene (7,000-3,000 cal yr BP), coniferous forest expanded slightly, as indicated by increases in Pinus and  (Xu et al., 2010). (B) Pollen-based precipitation reconstruction from Gonghai Lake . (C) Temperate deciduous broadleaf forest scores based on the biome reconstruction from Daye Lake. (D) Shannon index of palynological diversity from Daye Lake. (E) Evenness index from Daye Lake. (F) Temperature reconstruction based on the branched fatty alcohol ratio BNA 15 from the Dajiuhu wetland (Huang et al., 2013). (G) Pollen-based precipitation reconstruction from the Dajiuhu wetland (Sun et al., 2019). (H) U K 37 -inferred temperature record from Qinghai Lake (Hou et al., 2016). (I) Pollen-based precipitation reconstruction from Qinghai Lake (Herzschuh et al., 2010).
Picea pollen percentages, while percentages of thermophilous deciduous broadleaf taxa gradually decreased. The representation of Artemisia reached the lowest level in the profile. The hydrothermal combination at this time were probably worser than during the previous stage; nevertheless, the climate was still conducive to the growth of warm-temperate mixed forest (Figure 4). The substantial climatic differences existed between the northern, western, and southern parts of China at this time. At Qinghai Lake in the northeastern Tibetan Plateau, to the northwest of the Qinling Mountains, there was strong decrease in temperature which reached its lowest level during the Holocene, probably because of intensification of the Westerlies (Hou et al., 2016); while precipitation decreased substantially. However, temperatures at sites to the north of the Qinling Mountains were close to those sites to the south. At Daihai Lake, temperatures reached their highest level during the Holocene, while in the Dajiuhu wetland the temperature remained at a high level. This dipolar pattern of precipitation is a feature of the climate of China on various timescales in north and south of China and is discussed in detail in . Overall, our results of pollen percentages and biome reconstruction may indicate that the changes of precipitation pattern at Daye Lake during the Holocene were between northern and southern China.
During the late Holocene, especially after 3,000 cal yr BP, Pinus and Artemisia commenced a pronounced increasing trend and there was the increased representation of herbaceous vegetation (Figure 3). The biome reconstruction shows an increase in coniferous forests and herb-dominated vegetation, accompanied by decreases in deciduous broadleaf and warm mixed taxa (Figure 4), indicating a climatic deterioration. It is noteworthy that decreases in temperature and precipitation occurred at other sites adjacent to Qingling Mountains region (Figure 6).

Relationship Between Palynological Diversity and Holocene Climate Change in the Qinling Mountains
Pollen is an important indicator for the reconstruction of past plant diversity. There have been many studies on plant diversity based on the palynological diversity (Xiao et al., 2008;Liang et al., 2019;Rudaya et al., 2020;Mottl et al., 2021). Even though there are several uncertainties involved in reconstructing plant diversity from fossil pollen records. For example, palynological diversity cannot represent the full range of plant diversity due to the limitations of pollen identification, sample size, and the complexity of depositional processes (Odgaard, 2001;Weng et al., 2006;Meltsov et al., 2011;Birks et al., 2016a;Liang et al., 2019). Several taxa can only be identified to family level, and poorly represented entomophilous taxa can lead to the underestimation of plant diversity. However, the vegetation of the Qinling Mountains is dominated by woody plants which are identified mainly at the genus level (4 families + 16 genera of trees and shrubs in Figure 3). And at least 500 pollen grains were counted for each sample in this study to ensure a statistically adequate and representative sample of observations. Moreover, Zhang et al. (2021) suggested that the pollen loading in Daye Lake is predominantly wind-transported, indicating that the palynological diversity reconstructed in this study is representative of the vegetation of the Qinling Mountains region. Nevertheless, several studies have shown that palynological diversity is positively correlated with plant diversity (Odgaard, 1999;Meltsov et al., 2011;Birks et al., 2016a;Felde et al., 2016). Overall, despite these potential limitations of palynological diversity, we assume that the pollen records from Daye Lake can provide a reliable record of changes in plant diversity in the Qingling Mountains region during the Holocene.
The record of palynological diversity at Daye Lake is in good agreement with records of regional climate change (Figure 6), indicating that palynological diversity was highly dependent on the evolution of the regional climate. The low values of palynological diversity and evenness index at the beginning of the early Holocene are consistent with the relatively low temperature and precipitation indicated by the paleoclimate reconstructions. The increased palynological diversity and evenness index during the early Holocene is also consistent with the inferred warm and humid climate, and numerous studies have demonstrated that intervals of warm and humid climate are characterized by high vegetation cover and high palynological diversity in northwest, southwest and Loess Plateau regions of China (Fang et al., 2003;Zhang et al., 2005;Xiao et al., 2008;Zhao and Ding, 2014). Notably, palynological diversity and evenness index at Daye Lake continued to increase during the middle Holocene (7,000-3,000 cal yr BP), even though the precipitation decreased; temperatures remained at a high level (Figures 6A, F), the environmental conditions in the Qingling Mountains continued to favor a high degree of palynological diversity, indicating that the vegetation in the Qinling Mountains was more sensitive to temperature than to precipitation. The palynological diversity also commenced a decreasing trend after ∼3,000 cal yr BP, and it's clear that the environmental conditions in the area became less favorable for certain plant species, leading to their decreased representation or disappearance. Cheng et al. (2017) also suggested that the alpine timberline on Mount Taibai moved upward to a higher altitude during ∼4,200-2,000 cal yr BP and started to move downward from ∼2,000 cal yr BP. The visible decrease in the palynological diversity of the Daye Lake between 1,200 and 600 cal yr BP may also be resulted from the cold and dry environment which was consist with the sudden decrease in temperature and precipitation in the northern and southern of the Qinling Mountains within the recognizable range of the age error. At the same time, archeological evidence indicated that this period corresponded to the cultural and economic prosperity of the Song Dynasty (Archaeological research institute of Sui-Tang Dynasties of Shanxi Academy of Archaeology, 2018) which revealed the possible influence of human activities on the vegetation. Interestingly, during the past few centuries, palynological diversity increased, accompanied by an increase in temperature and a decrease in precipitation (Figure 6). This seems to support evidence that temperature had been a major factor affedting vegetation evolution in Qinling Mountains region, and that the rising temperatures resulting from global warming had a positive effect on palynological diversity in the region. However, the possible roles played by human activities should also be considered. For example, during the past several hundred years, the increases in Artemisia, Amaranthaceae, and Poaceae pollen percentages in the sediments can be attributed to the development of agriculture in low-altitude areas and open land may be one of the reasons for the increase in palynological diversity (Woodbridge et al., 2021).
Overall, the results of this study indicate that plant diversity in the Qinling Mountains during the Holocene was driven primarily by climate change. The relatively warm and humid environment to the south of the Qinling Mountains has resulted in a higher plant diversity in southern China than in northern China (Sun et al., 2009;Shang and Xing, 2016). Sun et al. (2013) concluded during the warm period of the Holocene, subtropical vegetation in the south expanded and migrated northward into the Qinling Mountains region. The impact of climate on vegetation is also exerted via extreme climatic events such as droughts, floods, and abrupt cooling, which lead to species extinctions and the reduction of biodiversity (Srivastava et al., 2021). Park et al. (2019) showed that in coastal East Asia successive droughts during the ∼8.2 ka caused forest damage and a major reduction of plant diversity. Additionally, Tao et al. (2013) demonstrated that in northwest China the cold and dry climatic event at 4.2 ka resulted in the conversion of grassland to desert, with a greatly reduced plant diversity and vegetation structure. However, the pollen record from Daye Lake seems not to contain clear evidence for the influence of abrupt climatic events on palynological diversity (Figure 6). Notably, climate change will affect the vegetation phenology which will subsequently impact the plant diversity. For example, changes in seasonal temperatures will alter the flowering time of plants, with increases in annual average temperature overall resulting in an earlier flowering time (Hart et al., 2014;Yu et al., 2016).
There is a long history of human impacts on the composition and structure of vegetation, associated with the development of agriculture and the use of fire (Ruddiman, 2003;Rudaya et al., 2020;Kulkarni et al., 2021). In southern China significant agricultural activity began ∼6000 years ago (Cheng et al., 2018), while the large-scale alteration of the vegetation composition and ecosystems has occurred during the past few millennia (Lewis and Maslin, 2015;McMichael, 2021). In the Qinling Mountains, agricultural activity began as early as ∼8000 cal yr BP, associated with cultures such as the Peiligang and Yangshao (Zhang et al., 2018). However, even though a long history of human activities was documented in this area, there is no ambiguous evidence for human disturbance of the vegetation at Daye Lake, possibly due to its high altitude (3,615 m a.s.l.), location above the forest line (Zhang et al., 2021). Liu H. et al. (2002) suggested that the manifest influence of human activities on vegetation is concentrated on the low latitudes in Mount Taibai.
In summary, the temperature of the Qinling Mountains has increased in recent decades due to global warming, which will impact the regional vegetation composition and palynological diversity Huang et al., 2020). Additionally, the palynological diversity of the Qinling Mountains region may increase with the ongoing climate warming in the future when the direct human interference is absent.

CONCLUSION
Our results show that the climate was the main factor affecting regional vegetation changes and palynological diversity. At the beginning of the early Holocene (11,700-10,700 cal yr BP), the cold and dry climate of the region resulted in predominantly coniferous forest vegetation, with a low palynological diversity and evenness index. During the early Holocene (10,700-7,000 cal yr BP), the warm and humid climate with fluctuations on a centennial time scale promoted the growth and expansion of deciduous broadleaf forest, and there was a corresponding increase in palynological diversity. During the middle Holocene (7,000-3,000 cal yr BP), although the climate became drier, the relatively warm environmental conditions continued to promote vegetation growth and high palynological diversity. During the past ∼3000 years, the temperature and precipitation in the Qinling Mountains decreased substantially, which led to reduction in palynological diversity, and human activities may play a potential role. However, during the past few centuries palynological diversity has increased along with rising temperatures, which may indicate that ongoing warming in the region may have a beneficial effect on vegetation growth and palynological diversity and the natural plant diversity may be in an increasing trend without human interference in the future in the region.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS
YZ: conceptualization, methodology, analysis, and writing -original manuscript. QC: conceptualization, methodology, reviewing, and editing. YH: reviewing and editing. AZ: supervision, project administration, providing revisions and comments, and funding acquisition. DW: providing revisions and comments. All authors contributed to the article and approved the submitted version.