Revisit the Medieval Warm Period and Little Ice Age in Proxy Records from Zemu Glacier Sediments, Eastern Himalaya: Vegetation and Climate Reconstruction

: The Late Holocene fossil pollen records from the Zemu glacier, located in Yabuk, North Sikkim, in the eastern Himalayas, effectively generated quantitative climate reconstructions based on the transfer function model. The transfer function model was developed by establishing a modern pollen–climate calibration set from the temperate alpine belt of North Sikkim. A redundancy analysis was carried out to detect the pattern of variation of climatic variables in the modern pollen datasets. The mean annual precipitation (MAP) and mean temperature of the warming month (MTWA) had the strongest inﬂuence on the composition of the modern pollen samples among the climatic variables considered in the analysis. Proxy data in the form of fossil pollen records were analyzed for reconstructing past climate based upon the relationships between modern pollen vegetation assemblages and climatic patterns. Transfer functions for MAP and MTWA were developed with the partial least squares (PLS) approach, and model performance was assessed using leave-one-out cross-validation. The validated model was used to reconstruct MAP and MTWA for the last 2992 cal years BP (1042 BC) in North Sikkim. The variability observed in the reconstructions was analyzed for past global climatic events. It was further compared with the available regional and hemispheric proxy-based climate reconstructions. The reconstructions captured comparable Medieval Warm Period (MWP) and Little Ice Age (LIA)-like events from the Zemu glacier region. The fossil pollen data and climate reconstructions were further compared with the mineral magnetism data of the subsurface sediment proﬁle.


Introduction
The decadal-to-millennial-scale climate variability records of the Quaternary period are of significant interest to scientists; their chronology and mechanisms have been widely discussed, prompting extensive analysis and recommendations [1]. However, further detailed analysis is needed in order to obtain a deeper understanding of the authenticity of past climate change studies. An introspective look at the dynamics and functioning of the past climatic system will allow us to understand its behavior and patterns in the future. This will improve greatly if long-term climate data are available beyond the existing instrumental records. Thus, palaeoclimatic data are essential to obtain comprehensive information about the past climate system, without which reliable forecasting of future conditions will not be possible [2]. High-elevation mountainous regions also represent unique areas for the detection of climatic change and the assessment of climate-related impacts [3]. One reason for this is that, as the climate changes rapidly with altitude over relatively short horizontal distances, so do vegetation and hydrology [4]. Past climate records from such high-elevation mountainous regions and remote areas can be used to evaluate the regional climate and its association with global climate dynamics. These high-altitude environments comprise glaciers, snow, permafrost, water, and the uppermost  The vegetation of Sikkim has three discrete botanical zones: tropical, temperate and alpine, distributed according to the elevation and characteristics of the vegetation [112]. The vegetation of the state is classified as low, middle, and upper hill forests, followed by rhododendron-conifer zones and Alpine scrubs and grassland [112]. The present study area falls under Rhododendron-conifer (cold temperate or sub-alpine zone, 2700-3600 m) and alpine scrub and grasslands (3600-4300 m and above). The dominant tree taxa along the sampling sites are mostly conifers: Abies densa, Juniperus recurva, Larix griffithiana, J. squamata, Picea spinulosa, and Tsuga dumosa, and broad-leaved taxa: Betula utilis and Rhododendron spp. At the termination of the tree line above 3600 m, J. pseudosabina and J. recurva grow as scrubs on hill slopes. The stunted height of Rhododendron spp was observed growing on exposed rock crevices at altitudes of 3600 m or above. Ephedra gerardiana were found growing only on the hilltops around Thangu. A few species belonging to genera such as Ranunculus, Anemone, Delphinium, Rhus, Potentilla, Primula, Fragaria, Cassiope and Allium were seen on the gentle slopes of open meadows. The detailed vegetation of Zemu Valley in Sikkim Himalaya has been thoroughly enumerated by Smith and Cave [113].

Modern Pollen Sites and Assemblages
The modern pollen sampling sites considered in the present study are well represented in the north-west transect (Lachen to Zemu Glacier) but the densities are relatively low in the northern transect of Lachen towards Thangu (Figure 1). These sites cover vegetation zones ranging from sub-tropical to temperate to alpine, with both coniferous and broadleaved elements. The climatic variation of these sites ranges with gradients in MAT (1.4-10.6 • C), MTCO (−7.4-2.6 • C), MTWA (9-16 • C), MAP (61-145.5 mm), and elevation (2600-4311 m) based on gridded climate data sets ( Figure 3). The details of the gridded climate data have been described in Supplementary Information.

Ordination
The unconstrained ordination using PCA showed that the first two axes capture 49.1% (axis 1, 31.1%; axis 2, 18.0%) of the total variance. This suggests that only a few variables dominate the data structure, but they are considered high in such high-elevation and climatically contrasting regions. The constrained ordination using RDA with pollen and climate data (MAP and MTWA) showed the eigen values of the first and second axes are 0.1444 and 0.0421, explaining 77.43% of the total variance in the pollen assemblage data (Figures 4 and 5). The RDA was also carried out with MAP and MTWA separately. The first RDA axis captures 12.23% variance in the pollen datasets and shows significant correlation with MAP (r = 0.696, p < 0.02), while with MTWA, the first axis captures 9.93% variance with significant correlation (r = 0.667, p < 0.046). The RDA with MAP and MTWA variables shows a similar pattern to the PCA, suggesting that both the variables have a strong influence on the underlying structure of the modern pollen data.

Modern Pollen Sites and Assemblages
The modern pollen sampling sites considered in the present study are well represented in the north-west transect (Lachen to Zemu Glacier) but the densities are relatively low in the northern transect of Lachen towards Thangu (Figure 1). These sites cover vegetation zones ranging from sub-tropical to temperate to alpine, with both coniferous and broad-leaved elements. The climatic variation of these sites ranges with gradients in MAT (1.4-10.6 °C), MTCO (−7.4-2.6 °C), MTWA (9-16 °C), MAP (61-145.5 mm), and elevation (2600-4311 m) based on gridded climate data sets ( Figure 3). The details of the gridded climate data have been described in Supplementary Information.  and climate data (MAP and MTWA) showed the eigen values of the first and second axes are 0.1444 and 0.0421, explaining 77.43% of the total variance in the pollen assemblage data (Figures 4 and 5). The RDA was also carried out with MAP and MTWA separately. The first RDA axis captures 12.23% variance in the pollen datasets and shows significant correlation with MAP (r = 0.696, p < 0.02), while with MTWA, the first axis captures 9.93% variance with significant correlation (r = 0.667, p < 0.046). The RDA with MAP and MTWA variables shows a similar pattern to the PCA, suggesting that both the variables have a strong influence on the underlying structure of the modern pollen data. In RDA, some pollen taxa are grouped together and some are separated ( Figure 5). The RDA axis 1 segregated coniferous taxa (Juniper, Larix, Abies, Picea and Tsuga) and broadleaved taxa (Juglandaceae, Quercus, Salix, Fabaceae and Oleaceae) on the left side. The coniferous taxa (Pinus) and broadleaved taxa (Betula, Alnus, Acer, Corylus, Rhododendron and Magnoliaceae) were separated on the right side of axis 1. It was shown that apart from Pinus, all other conifer taxa were grouped together by axis 1. The RDA axis 1 also differentiated in the pollen assemblage of wild grass (Poaceae < 50 µ) and cultivated grass (Poaceae > 50 µ) on two sides of the axis. Most of the shrub and herb elements were concentrated in the center of the ordination and did not show clear segregations. The modern pollen sites were separated in RDA according to climatic influence (Figure 5). The first RDA axis segregates sites with MAP (60.9-79.0 mm) and MTWA (9.0-10.3 °C) from sites with MAP (116.5-145.5 mm) and MTWA (13.2-16.4 °C). The sites separated by RDA axis 1 are comparable with the altitudinal gradient. The sites (SKM081, SKM087, SKM090, SKM101, SKM105 and SKM0108) segregated on the right side of axis 1 ( Figure  5) have low MAP and MTWA, are from a higher elevation (4016-4311 m), and are mostly closer to the Zemu glacier and within its vicinity, which is in an alpine climate. The sites (SKM022, SKM031, SKM034, SKM039, SKM046, SKM049 and SKM054) with high MAP and MTWA segregated on the left side of axis 1 ( Figure 5) are mostly from a lower elevation (2600-3303 m) and are influenced by a cool, temperate climate. The sites (SKM058, SKM060, SKM065 and SKM116), which are laid in intermediate and close to the center of both axes, and sites (SKM069, SKM074, SKM077), isolated on the left side of axis 1, are mostly influenced by both the higher and lower elevation climatic conditions of the region, with moderate MAP (95.9-109.7 mm) and MTWA (10.5-12.5 °C). The altitudinal gradient of these sites ranges from 3501 to 3915 m, with a sub-alpine climate, and seems segregated under the influence of RDA axis 2 ( Figure 5). Similar segregation was observed when the modern pollen data were grouped in zonation ( Figure 4) using unconstrained chord-distance clustering (CONISS) on a square-root-transformed modern pollen percentage assemblage.

Chronology and Palaeo-Vegetation History
Based on the derived chronology using linear regression, the Yabuk sediment profile covers the last 2992 cal years BP (1042 BC), confirming that these sediments belong to the Late Holocene. So, the present fossil records represent the vegetation and climate records since the Late Holocene. The chronology detail of the fossil sediment profile is provided In RDA, some pollen taxa are grouped together and some are separated ( Figure 5). The RDA axis 1 segregated coniferous taxa (Juniper, Larix, Abies, Picea and Tsuga) and broadleaved taxa (Juglandaceae, Quercus, Salix, Fabaceae and Oleaceae) on the left side. The coniferous taxa (Pinus) and broadleaved taxa (Betula, Alnus, Acer, Corylus, Rhododendron and Magnoliaceae) were separated on the right side of axis 1. It was shown that apart from Pinus, all other conifer taxa were grouped together by axis 1. The RDA axis 1 also differentiated in the pollen assemblage of wild grass (Poaceae < 50 µ) and cultivated grass (Poaceae > 50 µ) on two sides of the axis. Most of the shrub and herb elements were concentrated in the center of the ordination and did not show clear segregations.
The modern pollen sites were separated in RDA according to climatic influence ( Figure 5). The first RDA axis segregates sites with MAP (60.9-79.0 mm) and MTWA (9.0-10.3 • C) from sites with MAP (116.5-145.5 mm) and MTWA (13.2-16.4 • C). The sites separated by RDA axis 1 are comparable with the altitudinal gradient. The sites (SKM081, SKM087, SKM090, SKM101, SKM105 and SKM0108) segregated on the right side of axis 1 ( Figure 5) have low MAP and MTWA, are from a higher elevation (4016-4311 m), and are mostly closer to the Zemu glacier and within its vicinity, which is in an alpine climate. The sites (SKM022, SKM031, SKM034, SKM039, SKM046, SKM049 and SKM054) with high MAP and MTWA segregated on the left side of axis 1 ( Figure 5) are mostly from a lower elevation (2600-3303 m) and are influenced by a cool, temperate climate. The sites (SKM058, SKM060, SKM065 and SKM116), which are laid in intermediate and close to the center of both axes, and sites (SKM069, SKM074, SKM077), isolated on the left side of axis 1, are mostly influenced by both the higher and lower elevation climatic conditions of the region, with moderate MAP (95.9-109.7 mm) and MTWA (10.5-12.5 • C). The altitudinal gradient of these sites ranges from 3501 to 3915 m, with a sub-alpine climate, and seems segregated under the influence of RDA axis 2 ( Figure 5). Similar segregation was observed when the modern pollen data were grouped in zonation ( Figure 4) using unconstrained chord-distance clustering (CONISS) on a square-root-transformed modern pollen percentage assemblage.

Chronology and Palaeo-Vegetation History
Based on the derived chronology using linear regression, the Yabuk sediment profile covers the last 2992 cal years BP (1042 BC), confirming that these sediments belong to the Late Holocene. So, the present fossil records represent the vegetation and climate records since the Late Holocene. The chronology detail of the fossil sediment profile is provided in Table 1. The rate of sedimentation calculated was mostly similar from top to bottom in depth, and its average rate was 0.08 cm/year. Based on the CONISS of the percentage of pollen data, three zones, YAB-I, YAB-II and YAB-III, were established to represent the major phase of vegetation change. Two zones, YAB-II and YAB-III, are further divided into two sub-zones to represent the minor vegetation changes within the broad phase ( Figure 6). The vegetation changes in each zone are described in Supplementary Information.  Based on the CONISS of the percentage of pollen data, three zones, YAB-I, YAB-II and YAB-III, were established to represent the major phase of vegetation change. Two zones, YAB-II and YAB-III, are further divided into two sub-zones to represent the minor vegetation changes within the broad phase ( Figure 6). The vegetation changes in each zone are described in Supplementary Information.

Pollen-Climate Calibration and Transfer Function
The RDA shows that there is a statistically significant relationship between modern pollen assemblages and the climatic variables MAP and MTWA ( Figure 5). A transfer function was developed for MAP and MTWA as these were the significant climatic variables entered in the ordination analysis, RDA. Prior to final model development, we found four samples (SKM022, SKM069, SKM074 and SKM101) producing extreme values in the analyses. These sites were removed from the training set prior to performing PLS regression and calibration. The improved pollen-climate PLS transfer function was finally developed based on the screened dataset ( Table 2). The performance of transfer function models for MAP and MTWA assessed with components retained in the PLS regression approach and assessed by leave-one-out cross-validation showed the sixth component for MAP and the fifth component for MTWA to be the best performing models (Table 2, Figure  7). The sixth component of MAP showed the lowest value (17 mm) of RMSEP and the highest value (0.585) of r 2 between observed and predicted MAP ( Table 2). A strong relationship (r 2 = 0.993, p < 0.01) was observed between observed MAP and the values pre-

Pollen-Climate Calibration and Transfer Function
The RDA shows that there is a statistically significant relationship between modern pollen assemblages and the climatic variables MAP and MTWA ( Figure 5). A transfer function was developed for MAP and MTWA as these were the significant climatic variables entered in the ordination analysis, RDA. Prior to final model development, we found four samples (SKM022, SKM069, SKM074 and SKM101) producing extreme values in the analyses. These sites were removed from the training set prior to performing PLS regression and calibration. The improved pollen-climate PLS transfer function was finally developed based on the screened dataset ( Table 2). The performance of transfer function models for MAP and MTWA assessed with components retained in the PLS regression approach and assessed by leave-one-out cross-validation showed the sixth component for MAP and the fifth component for MTWA to be the best performing models (Table 2, Figure 7). The sixth component of MAP showed the lowest value (17 mm) of RMSEP and the highest value (0.585) of r 2 between observed and predicted MAP ( Table 2). A strong relationship (r 2 = 0.993, p < 0.01) was observed between observed MAP and the values predicted for the study sites based on observed MAP (Figure 7). In the case of MTWA, the fifth component showed the lowest value (1.1537 • C) of RMSEP and the highest value (0.696) of r 2 between observed and predicted MTWA ( Table 2). A strong relationship (r 2 = 0.982, p < 0.01) was observed between the observed MTWA and the values predicted for the study sites based on the observed MTWA ( Figure 7). A comparison of observed MAP and MTWA and residuals (predicted minus those observed for MAP and MTWA) showed that the calibration model predicts values relatively well within the range of observed MAP and MTWA ( Figure 7). Table 2. Model performance statistics as assessed by leave-one-out cross-validation for the using first x (x = 1-6) components of the partial least square regression (PLS). Component marked with * are used for further analysis and for climate (MAP and MTWA) reconstruction.

Late Holocene Quantitative Palaeoclimate Reconstructions
The successful validation of the transfer function models based on standard statistical criteria [38] shows prospects for the quantitative reconstruction of both precipitation (MAP) and temperature (MTWA) variables in the present study. This has led to the reconstruction of MAP and MTWA over 2992 cal years BP (1042 BC) for the Zemu glacier region of Yabuk, North Sikkim, eastern Himalaya (

Late Holocene Quantitative Palaeoclimate Reconstructions
The successful validation of the transfer function models based on standard statistical criteria [38] shows prospects for the quantitative reconstruction of both precipitation (MAP) and temperature (MTWA) variables in the present study. This has led to the recon-  The reconstructed MAP showed variability during several periods as higher (wet) and lower (dry) episodes ( Figure 8). The beginning of the reconstructed MAP, i.e., during 2992-2188 cal years BP (BC 1042-238), shows a lowering of precipitation, with the lowest (90 mm) in 2870 cal years BP (BC 920). The precipitation was higher in 1551 cal years BP and suddenly decreased until 1157 cal years BP (AD 793). Then, high precipitation was

Magnetic Parameters, Fossil Pollen and Quantitative Palaeoclimate
The reconstructed quantitative precipitation and temperature records were compared with the pollen zones developed in the fossil data and supplemented with magnetic parameters, which are described in detail in Supplementary Information (Figure 9).

Magnetic Parameters, Fossil Pollen and Quantitative Palaeoclimate
The reconstructed quantitative precipitation and temperature records were compared with the pollen zones developed in the fossil data and supplemented with magnetic parameters, which are described in detail in Supplementary Information (Figure 9).

Modern Pollen Distribution
The samples collected for modern pollen analysis have a distinct pattern of dominance of species that is also segregated also in terms of the variation in the altitude of sampling. The instances of grouping and/or clustering in both CONISS and ordination statistics following the same trends are evident. The first group (Group I) of samples or clusters included higher altitude samples (SKM90, SKM81, SKM87, SKM77, SKM101,

Modern Pollen Distribution
The samples collected for modern pollen analysis have a distinct pattern of dominance of species that is also segregated also in terms of the variation in the altitude of sampling. The instances of grouping and/or clustering in both CONISS and ordination statistics following the same trends are evident. The first group (Group I) of samples or clusters included higher altitude samples (SKM90, SKM81, SKM87, SKM77, SKM101, SKM116, SKM108), occupying an altitudinal range of (3608-4300 m) and samples (SKM22) of lower altitude (2600 m). These samples showed a high dominance of alpine taxa. However, taxa such as Tsuga, Quercus pollen, and a few instances of Larix pollen were also present, probably arriving from adjoining regions.
In the second group (Group II), the samples (SKM54, SKM74, SKM65, SKM58, SKM34, SKM69, SKM39, SKM46, SKM49, SKM31, and SKM60) collected at various sites in the valley had an altitude range of 2800-3700 m. These sets of samples recorded the dominance of subalpine tree taxa, a few shrubs and herbs, and larger grass pollen and aquatic taxa. This variation in taxa can again be attributed to the lower temperatures at this altitude and increased soil moisture due to snow melt and a higher amount of precipitation. This sampling area is exposed to variable climatic conditions, allowing mixed species to appear in the sediments. However, a slight shift towards higher altitude taxa cannot be overruled. There is one sample (SKM 105), which was collected from the distant region of Chopta Valley, that has similar taxa but high levels of Rhododendron in the sample, consequently forming a separate node in the second group or cluster of samples. However, the ordination statistics and RDA imply considerations for species grouping, including climatic parameters. Hence, it is imperative that all the data and sample groups are separated not only based on the species and climate variables (MAP and MTWA), but also according to the altitude of sampling. To further strengthen the results of RDA, the CONISS results are also comparable, giving a definitive pattern to the species distribution in the region.

Pollen-Climate Relationships
The relationship between pollen assemblages and the environmental variables established using the ordination method showed the significant influence of MAP and MTWA on vegetation. The climatic variable MAP with axis 1 showed a high correlation compared to MTWA, which is slightly less correlated with axis 1, but both are statistically significant (p < 0.05). This suggests that the control of precipitation over the distribution of surface pollen could be the primary environmental driver controlling vegetation and species variations. It has been observed that sites with arid or semiarid climatic conditions, with vegetation cover and composition, are mostly controlled by moisture availability [62]. The analysis of the satellite-based normalized difference vegetation index (NDVI) compared with rainfall data from states in India showed a strong relationship between them. The study showed that NDVI is affected positively by the amount and intensity of rainfall [88]. This relationship is stronger during the monsoon season, the period that contributes the maximum percentage of water availability in annual rainfall. The NDVI precipitation result also supports the fact that precipitation has a stronger effect on vegetation distribution in the present high-elevation study area of the eastern Himalayas. In North Sikkim, precipitation in the form of rain or snow occurs throughout the year. Thus, in the present study region, precipitation has a strong influence throughout the year, which supports the present findings that highlight the significant influence of precipitation on the regional vegetation. However, in Sikkim, the rainfall/precipitation decrease sharply at higher elevations (>2000 m) in all seasons except during winter [111]. Hence, the availability of soil moisture for most taxa is high at all times, either due to rainfall or snow. However, due to the variation in altitude, the precipitation and temperature distribution at all sampling sites vary within a wide range. In the glacier region (>4000 m), the temperature decreases with snow cover, and precipitation also loses its intensity. In the cool, temperate region (2600-3300 m), both precipitation and temperature are usually at their peak during the south-west monsoon months (June to September). The rise in temperature could lead to high evapotranspiration and become the limiting factor for vegetation. The significant relationship observed with temperature could be possible, because the growth of the plants varies greatly along both latitudinal and altitudinal gradients. In the present study, the influence of both temperature and precipitation might be due to the fact that both temperature and precipitation vary greatly with altitudinal gradients. Due to this altitudinal-climatic distribution, both temperature and precipitation become the important limiting factors for the growth of plants in this region. The altitudinal influence on climate cannot be ruled out in mountainous and high-altitude regions. As the elevation increases, the temperature decreases in most of the mountainous region [4,132,133], in addition to precipitation, as also observed in the present study sites (Table S1) and reported for the entire region of Sikkim by [111]. Similar inferences on pollen-temperature-precipitation have been observed from adjoining regions such as the Tibetan Plateau, China, and Mongolia [57,59,62,[134][135][136][137] on modern pollen-climate calibration, along with the statistical parameters observed in the transfer function model. The influence of precipitation as a primary environmental factor can also be explained by the fact that modern pollen sites are located in high-elevation areas with arid or semiarid climatic conditions, where vegetation cover and composition are mostly controlled by moisture availability.

Environmental Geomagnetism and Climate Reconstruction
The reconstruction of mean annual precipitation (MAP) has been based on magnetic minerals within a soil in various studies using different quantitative techniques [96,[138][139][140][141][142]. Magnetic minerals such as magnetite/maghemite, hematite, and goethite abundances in soils have been identified to have a quantitative relationship with precipitation [96,142,143]. Magnetic palaeo-precipitation proxies have been studied widely in the Quaternary to recent Loess-paleosol systems in China [138,144], New Zealand [145], Alaska [146], North America [147][148][149], and Russia [150,151]. These records exhibit that the relationship between magnetic minerals and MAP differs spatio-temporally [96]. The Loess-Paleosol sequences of China have shown that variation and aeolian dust fluxes are influenced by regional monsoonal climate patterns and increases in pedogenic processes due to high precipitation and temperature [96,138,[152][153][154]. Thus, the strong relationship between iron oxide mineralogy in soil sequences and climatic variations has been well established globally [96]. However, the magnetic palaeo-precipitation proxies used to estimate MAP have been said to have certain uncertainties [96,155,156].
In this study, we do not attempt a direct palaeo-precipitation reconstruction using mineral magnetic data. We have attempted to make a comparison of the reconstructed MAP using palynological data with magnetic mineral variations. Mineral magnetism data are an additional tool to identify the variations in the palaeoclimate of the Zemu glacier. Similar variations comparable to the reconstructed climate (MAP, MTWA) at Zemu glacier were found. There were many controlling factors, such as soil moisture, glacial melt, detrital deposition, podogenic processes, and regional sediment influx at the Zemu glacier. All these factors, along with the climatic parameters, control the magnetic mineral variations in the sediment. Here, the fossil pollen data also correspond to vegetation changes along with climatic variations at the Zemu glacier. The fossil pollen data in zone YAB-I is less, but an increase in the ferromagnetic mineral content is concurrent with the high MAP and MTWA and high pollen counts. In the other two zones, the vegetation gives a pattern similar to the changes in MAP and MTWA, along with the mineral magnetic. The glacier also demonstrated changes regarding the type of sediment influx at the sampling site during high precipitation periods in zones YAB-II and YAB-III, which are composed of coarse gravel with sand and silty clays. Such changes are visible in the lithology profile of the samples collected near the Zemu glacier. It is evident that material coming down from the snout of the glacier and the surrounding topography influenced the magnetic mineralogy. The vegetation, climate, and sediment flux in the valley during the Late Holocene are clearly linked together in our record. Both palynological and environmental geomagnetic data provide strong evidence of the influence of MAP and MTWA as climatic factors working in the Zemu glacier valley.
The environmental magnetism parameters are significantly low in terms of their numerical values throughout the sediment profile. They successfully capture the signals and oscillations in the environment and climate at the Zemu glacier. These parameters were analyzed for all the samples, which included samples not analyzed for palynomorphs due to the high sand content. Though we managed to obtain the entire profile data, it is possible that including all the sample data reduced the possibility of easily identifying the magnetic behavior of the minerals in the profile. However, it was important that we supported our climatic interpretation and reconstruction based on other supportive evidence.

Evaluation of Climate Reconstructions
The climate of the Late Holocene, especially the last 1000-1500 years, sets the background for 20 th -century climate warming. The climate records during this period are required to perceive the natural variability in the climate system. The proxy-based long climate records emphasize the interest in quantitative climate data beyond the range of instrumental records. Long-term climate reconstructions derived from various well-dated proxy data indicate that the past millennium is the best documented interval with both historical and climate data [157]. The present pollen-based quantitative climate (MAP and MTWA) reconstructions successfully captured well-known global climate events of the recent millennium. These climatic signals helped us assess the background for recent century-long climate warming in the largest glaciated region of the eastern Himalayas.
In the present reconstructed MTWA, the high temperature covering the time span of 929-472 cal years BP (AD 1021-1478) is comparable with the MWP. After attaining a gradual rise in temperature during 509-472 cal years BP (AD 1441-1478), the significant lowering of temperature until 421 cal years BP (AD 1529) could be related to the cooling period near the beginning of the LIA. This cooling period observed in the present MTWA reconstruction is comparable with a widely known event, the Spörer Minimum, observed during AD 1460-1550 [158]. Jiang and Xu [158] clearly demonstrate that if the Spörer Minimum extended from AD 1460 to 1550, it could be a specious result; if its extent was AD 1400-1510, it could be determined as an actual feature of solar variability during that time. The slight decrease in temperature after attaining the increasing trends in our MTWA reconstruction during 296-233 cal years BP (AD 1654-1717) is comparable with a cold event widely known as the Maunder Minimum (AD 1645-1715) in the middle of the LIA [159]. The Maunder Minimum refers to a period of low solar activity [160]. The physical link between low sunspot activity and cooling temperatures has not been established, but the coincidence between the Maunder Minimum and the deepest trough of the LIA is suggestive of such a connection [161]. The present reconstructed MTWA has not captured this event significantly compared to the Spörer Minimum. Thus, the overall lowering of temperature covering a time span of 421-170 cal years BP (AD 1529-1780) suggests that the LIA event was present around the Zemu Glacier Valley. The LIA episodes observed in the reconstructed MTWA showed a sudden decrease in temperature in the early part of the initiation of the cold period, and gradually the temperature increased, with a slight decrease during the middle of the LIA. The reconstructed temperature around 107 cal years BP (AD 1843) shows the arrival of warming after the cold episodes of LIA in the Zemu glacier valley and is similar to the Current Warming Period (CWP).
The reconstructions from our records captured comparable well-known climatic events (MWP and LIA), but the timing identified is indeed not the same as that recorded in various studies (described in Sections 4.5 and 4.6). The MWP was a time of warm climate during AD 1100-1300, which peaked in the 11th century and ended abruptly around 1300 AD [8]. Similarly, the LIA event is a known widespread climatic cooling event that lasted approximately from the 14th to the 18th centuries AD [12]. Later, Mann et al. [10] argued that MWP might be related to other climate events around the world and occur from about AD 950 to 1250, and that the cooling event LIA lasted from AD 1400 to 1700. However, Graham et al. [11] reported the duration of MCA/MWP as ca. 900-1350 AD and the LIA as ca. 1500-1850 AD. The present climate reconstruction shows cool (dry) and warm (wet) variability within the LIA and MWP, respectively. We observed cold (dry) climatic records punctuated within the 'warm (wet)' period of MWP. Similarly, warm (wet) records are observed during the 'cool (dry)' period of LIA. This illustrates the fact that 'Medieval Warm Period' and 'Little Ice Age' are moveable terms depending on region. However, the present finding supports that the worldwide climatic variability of MWP with a warm (wet) climate and LIA with cold (dry) periods also exists in this remote region of the Zemu glacier valley in the eastern Himalayas. Wanner et al. [162], based on a synthesis review of mid-tolate Holocene global climate change, concluded the existence of warm and cold periods, especially hemispheric-scale warming leading to the MWP and subsequent cooling into the LIA. The possible role of this climatic variability could be due to various natural factors. These could be fluctuations in solar activity, large tropical volcanic eruptions, internal variability in oceanic oscillations, including the El Niño Southern Oscillation (ENSO) and the North Atlantic Oscillation (NAO), changes in the thermohaline circulation, and complex feedback mechanisms between the ocean, atmosphere, sea ice, and vegetation [162]. The strong coupling between the biogenic proxies and reconstructed total solar irradiance during the last millennium from the northeastern Arabian Sea also suggested solar control on the intensity of the Indian monsoon and likely the MWP and LIA as well. Moreover, the large eruptions in the tropics and high latitudes were the primary drivers of interannual to decadal temperature variability in the Northern Hemisphere during the past 2500 years [163]. Based on the model simulations, Polanski et al. [157] stated that the eastern Himalayas are affected by summer variations in the temperature gradient between the Bay of Bengal and Indian subcontinent and by a zonal band of intensified Indian-East Asian monsoon link north of 25 • N.
Apart from these broad global events of MWP and LIA, our climate reconstruction is also comparable with other climatic events that occurred elsewhere in the globe [164] but were not widely explained compared to MWP and LIA. These events were recorded in the reconstruction before the MWP period. The cooling of temperature at the beginning of the reconstructed MTWA (2992-2870 years BP, BC 1042-920) might be linked with the Sub-Atlantic Cold Period (SACP) that occurred during BC 975-250 [164]. The rise in temperature in the present MTWA reconstruction during 2188-1551 cal years BP (BC 238-399) is similar to the Roman Warm Period (RWP) observed during BC 250-AD 450 [164]. Another event known as the Dark Age Cold Period (DACP), which took place during AD 450-950 [164], is also comparable with the cooling events observed in our temperature reconstruction during the later part of 1157-1021 cal years BP (AD 793-929). However, we do not support the claim that these events were strongly captured in our reconstruction due to the fact that the palynological samples are smaller at the beginning of the chronology. However, at the same time, we can argue that the present reconstruction has the ability to capture climatic fluctuations that occurred during the Late Holocene elsewhere. In future studies, the observation can be strengthened with an increase in palynological sample numbers and a higher-resolution chronology. Thus, greater emphasis was given to MWP and LIA events, which were compared with other records in detail.

Climate Reconstructions Comparison
In the present study, we provide the first quantitative climate (precipitation and temperature) reconstruction of the Zemu glacier valley in the eastern Himalayas using the transfer function model. In order to look at our climate reconstructions of the Zemu glacier from a larger geographical perspective, we compared them with proxies-based climate reconstructions carried out in India and adjoining areas, including south Asian continents and the northern hemisphere. For this comparison, two approaches were considered: (i) qualitative climate inference made from different proxy-based studies and (ii) quantitative climate reconstruction from various proxies available from NOAA data contributors (www.ncdc.noaa.gov, accessed on 10 January 2022) and obtained from personal contribution on request. For both approaches, comparisons were made for past climatic events and variability recorded during the Late Holocene.

Climate Reconstructions Comparison-Qualitative Approach
The climatic variability during the MWP and LIA time periods recorded in the present pollen-based climate reconstructions is comparable with other pollen and proxy-based studies carried out in different regions of India, the Himalayas, and the adjoining oceanic region. The time span and climatic inference made for the MWP and LIA periods in these studies [16,84,85,87,158, are presented in Figures 10 and 11, respectively, and are described in the following section. We have given an extensive explanation and comparison of these records in Supplementary Information. Quaternary 2023, 6, x FOR PEER REVIEW

Climate Reconstructions Comparison-Quantitative Approach
The quantitative climate (precipitation and temperature) reconstructions developed in the present study were compared graphically with available Late Holocene reconstructions ( Figure 12). For this comparison, reconstructed temperature and precipitation data using other proxy data were procured from NOAA data contributors and also provided by the author of the research paper.
The MAP reconstruction is compared with available precipitation, drought, and ice core records. The MAP is compared with 3500 years of annual precipitation reconstruction from the Northeastern Tibetan Plateau based on tree-ring data [205] and the north central China 1800-year decadal precipitation index reconstruction using stalagmite and historical document records [206] shows similar trends in the later period of MWP and LIA. The reconstructed Palmer drought severity index (PDSI) from the Monsoon Asia Drought Atlas (MADA) [207] from the nearest grid point (88.75/28.75) in the present study site also shows consistency in the trends during LIA. The MAP has been also compared with ice core records of Rongbuk Glacier [173], Dasuopu [208,209], Guliya [210], Dunde [211], and Puruogangri core I and II [212]. Apart from Puruogangri core II, all other ice core records compared with the present reconstructed MAP shows similar variability in the reconstruction.
Similarly, the present reconstructed temperature (MAT) has been compared with temperature reconstructions from adjoining areas and regional, Arctic, hemispheric and global records ( Figure 13). The tree-ring-based June-August temperature reconstruction from Bhutan [16] shows cold conditions during the early 15th, 16th, and late 17th to early 18th centuries. Another tree-ring-based temperature reconstruction from Nepal [200] also shows a lowering of temperature during AD 1670-1770, which falls in the period of LIA. In addition, temperature reconstruction for Asia [213] using multi-proxy records, temperature reconstruction using stalagmite records from Beijing, China [214], temperature reconstruction in China using proxy records with methods of principal component regression and partial least squares regression [215], and Arctic-wide temperature reconstruction for the past 2000 years [216] also show comparable variability with the present reconstructed temperature. The temperature reconstructions available from various hemispheric and global levels include: multi-proxy reconstructions of the extra-tropical northern hemisphere  • N) mean temperature [217]; reconstructions of northern hemisphere temperatures and climate forcing over the past 1000 years using multi-proxy records [218]; low-frequency temperature signals reconstructed using Regional Curve Standardized tree-ring chronologies from 14 sites in the northern hemisphere extra-tropics [219]; northern hemisphere high-elevation tree-ring-based temperature reconstruction using the regional curve standardization technique [220]; multi-proxy northern hemisphere temperature reconstruction using inverse regression-truncated EOF climate field reconstruction [221]; northern hemisphere temperature reconstruction calculated by combining low-resolution proxies with tree-ring data [222]; multi-proxy northern hemispheric surface temperature reconstruction [223]; and multi proxy global surface temperature [224] also show similar temperature variability to the present temperature reconstruction in different time periods, especially during MWP and LIA ( Figure 13). However, the differences in variability in both the temperature and precipitation comparative analysis with the present climate reconstruction might be due to the differences in the proxy data used, the reconstructed seasons, and the resolution of the proxies.

Climate Reconstructions and Fluctuation of Sikkimese Glaciers
The present findings of LIA in the climate reconstructions might have linkages with the advances and retreats of the Zemu glacier and other glaciers in Sikkim, Himalaya. In 1995, a study team at the Khangchendzonga Biosphere Reserve reported that since the LIA, the Zemu glacier, Onglokthang glacier, and Rathong Chu glacier in North Sikkim had retreated by 3-4 km,~500 m, and~600 m, respectively [225]. A similar response of the cold and dry period of climate on glacier advancement has been recorded in the glacier region of Sithikher Bog, western Himalayas [179], and the Sonapani glacier advance stage II in the Lahaul valley, western Himalayas [226]. The climate reconstruction from this largest glacier of the eastern Himalayas also shows consistency with the climatic observations of the Late Holocene made from the other largest glacier of the western Himalayas [177]. Chaujar [186] also reported that the northern and southern hemispheric glaciers started their retreat in the mid-18th century, which indicates the end of the LIA climatic period. This temporal linkage suggests the possible existence of a common trend in the mountain regions of both hemispheres and the Himalayas, and, thus, its occurrence is a global phenomenon.
Observations made during the early 1940s and 1970s [106,227,228] also suggested the expansion and then receding of the Zemu glacier. The valley has its maximum depth at the center, where the thickness of the glacier ice exceeds 300 m. This suggests that the glacier once occupied the entire valley, extending from the northern to the southern rock walls, and that it has undergone considerable lateral shrinkage, giving rise to the northern and southern trenches. Subsequent movement of the glacier, confined by the lateral moraine ridges, has eroded the bedrock to a greater extent in the central zone than that underlying the trenches from which the glacier has receded to its present position [106].
The Zemu Glacier, located in in Sikkim, the Himalayas, which was visited by the 2nd German Himalaya Expedition in 1931, has been receding since 1891 [227,228], which correlates with evidence from our records of warming trends in temperature. These observations support the evidence of LIA-like events and recent warming phases recorded in our quantitative temperature reconstruction.
The climate oscillation recorded in the present climate reconstruction, especially during the well-known global climate periods of MWP (warm-humid climate) and LIA (cold-dry climate), could be related to glacier behavior and tree line fluctuations. During the warm and humid period, the glaciers might have retreated in this region, and, as a consequence, tree line might ascend to higher elevations. During the cold-dry period, the glacier of this region as well as the tree line might have advanced towards the lower elevations. Thus, the pollen-based climate reconstruction from the glacier region enabled a chronology of glacier activity for the LIA and more recent periods to be established.
The model simulations support the suggestion that the LIA may have been brought about by low northern hemisphere orbital forcing during the Late Holocene, with coincidentally unusually low solar activity and a high number of major volcanic events [162]. So, we can contemplate that this region of the Himalaya was also influenced by solar fluctuations and volcanic activities that occurred elsewhere during the LIA period. The tree-ring studies carried out in the adjacent Bhutan Himalaya [16] also reported that the cooling of the climate during LIA was influenced by solar and volcanic eruptions.

Conclusions
The calibration set for modern pollen-climate (temperature and precipitation) provided the development of quantitative climate reconstructions from sub-surface sediment pollen spectra from the Zemu glacier region, eastern Himalayas. The development of a transfer function model using a linear model based on Partial Least Square regression led to the reconstruction of the mean annual precipitation (MAP) and the mean temperature of the warmest month (MTWA). Precipitation has the strongest influence on the composition of the modern pollen samples among the other climatic variables considered in the analysis. The climate reconstruction based on pollen data, supplemented with environmental magnetism analysis, cumulatively indicates the changing climatic scenario around the Zemu Glacier during the later part of the Holocene. In the absence of long-term instrumental records from this largest glacier region of the eastern Himalayas, proxy-based quantitative climate reconstructions provide an understanding of the relative contribution of changes in temperature and precipitation in glacier environments. The present study provided valuable insight into past precipitation and temperature changes, which helped us identify shifts in climatic conditions as demonstrated by various research outputs, from regional to hemispheric. The climatic variability reconstructed in this study clearly demonstrates close agreement between oscillatory climatic phenomena of the past in the region of the North Atlantic Ocean and the distant Himalayas. Our climate reconstructions capture well-documented world-wide climatic events, such as the MWP and LIA, as well as those of the even earlier Dark Ages Cold Period and Roman Warm Period. Since the present climate reconstructions have been carried out at a remote location in the largest glaciated region of the eastern Himalayas, the observed warming during the late 19 th century and then the cooling trend during the early 20th century are most likely of natural origin and not due to any anthropogenic influences. The overall assessment of quantitative climate reconstruction in the Zemu glacier area and its temporal trends, which are comparable with regional, continental, hemispheric perspectives and from a global perspective, should look forward to a more comprehensive analysis with both spatial and temporal coverage.

Data Availability Statement:
The data that support the findings of this study are available from the corresponding author upon request.

Acknowledgments:
The authors (NM and SKS) are thankful to the Director, Birbal Sahni Institute of Palaeosciences, Lucknow for providing the necessary support and facilities to them for the compilation and to publish this work (BSIP No. 11/2015-16). Thanks, are also due to Late R. Ramesh, Physical Research Laboratory, Ahmadabad for providing us C-14 dates of the sediments. The author (NM) is grateful to the Indian Institute of Geomagnetism, Navi Mumbai, India for providing her permission and facilities to work on environmental geomagnetism. Author (SKS) is thankful to Pradeep Kumar, CCF and Usha Lachungpa, Principal Research Officer (WL) from Department of Forest, Government of Sikkim for providing all the support and permission to carry out the field work in North Sikkim. Authors are also thankful to NOAA for providing contributed climate reconstruction data. One of the authors NM is highly grateful to the C.S.I.R., Government of India for providing her financial support during the period of her PhD dissertation. Author NM is also grateful to Department of Science and Technology (DST), New Delhi for the financial support provided to her during the compilation of this manuscript, under the WOS(A) scheme for the project No. SR/WOS-A/ES-18/2014(G). Author SKS also wish to thanks porters especially Tshering Lachenpa and Pema Lachenpa who helped him during the field expedition to this high-altitude glacier terrain of the North Sikkim. The authors are indebted to three anonymous reviewers for their critical comments and valuable suggestions.