Late Holocene Glacial Fluctuations of Schiaparelli Glacier at Monte Sarmiento Massif, Tierra del Fuego (54°24′S)

The Magallanes–Tierra del Fuego region, Southern Patagonia (53–56°S) features a plethora of fjords and remote and isolated islands, and hosts several thousand glaciers. The number of investigated glaciers with respect to the multiple Neoglacial advances is based on a few individual studies and is still fragmentary, which complicates the interpretation of the glacial dynamics in the southernmost part of America. Schiaparelli Glacier (54°24′S, 70°50′W), located at the western side of the Cordillera Darwin, was selected for tree-ring-based and radiocarbon dating of the glacial deposits. One focus of the study was to address to the potential dating uncertainties that arise by the use of Nothofagus spp. as a pioneer species. A robust analysis of the age–height relationship, missing the pith of the tree (pith offset), and site-specific ecesis time revealed a total uncertainty value of ±5-9 years. Three adjacent terminal moraines were identified, which increasingly tapered towards the glacier, with oldest deposition dates of 1749 ± 5 CE, 1789 ± 5 CE, and 1867 ± 5 CE. Radiocarbon dates of trunks incorporated within the terminal moraine system indicate at least three phases of cumulative glacial activity within the last 2300 years that coincide with the Neoglacial phases of the Southern Patagonian Icefield and adjacent mountain glaciers. The sub-recent trunks revealed the first evidence of a Neoglacial advance between ~600 BCE and 100 CE, which so far has not been substantiated in the Magallanes–Tierra del Fuego region.


Introduction
Climatically, the region of Patagonia is situated between the subpolar low pressure trough (~60° S) and the subtropical high-pressure system (~35° S), leading to prevailing westerly winds by this synoptic constellation [1,2]. The mountain range of the Southern Andes forms an obstacle perpendicular to the zonal band of persistent westerlies, resulting in large amounts of precipitation caused by the orographic uplift of humid air masses on the windward side [1,[3][4][5]. The amount of precipitation exceeds several thousand millimeter per year [mma −1 ] on the western side, but decreases rapidly towards the east due to the foehn effect, forming one of the world´s sharpest moisture gradients [6,7]. Though the Pacific Ocean causes moderate temperatures, the vast sum of precipitation nourishes more than 11,000 glaciers south of 45.5° S, covering an area of 22,636 ± 905 km 2 in 2016/17 [8] ( Figure 1A). While the Andes north of the Strait of Magallanes are striking north-south, the mountain range of the Cordillera Darwin (CD) on Tierra del Fuego is predominantly west-east orientated. In addition to the otherwise dominant westerly winds, the CD lies within the sphere of influence of sub-Antarctic air masses, and glaciers may respond differently to climatic changes [9,10]. Although the ~250 km wide CD and the nearby islands comprise over 2100 individual glaciers (2931 ± 155 km 2 [8]), most studies dealing with Holocene glacial variability have focused on the Northern Patagonian Icefield (NPI) and Southern Patagonian Icefield (SPI), while only a few have related to the Magallanes-Tierra del Fuego region [11,12]. Besides the limited accessibility of the fjord-seamed mountain range ( Figure 1B), around 60% of the total glacierized area is either fjord-or laketerminating, resulting in a limited terrestrial evidence of glacial variability [8,10]. In consequence, only a few glaciers have been investigated in detail for their historical variations [13][14][15][16]. Thus, the glaciers of Tierra del Fuego are still considered poorly investigated [17].

Figure 1. A)
Overview of the Patagonian Andes and the Tierra del Fuego with its main ice bodies (black) [8]. On the upwind side of the mountain ranges, blue shades display precipitation amounts exceeding 1000 mma −1 [18]). Light-blue-filled circles and their diameters represent the location and the corresponding annual precipitation of climate stations (data obtained from the Dirección Meteorológica de Chile (DMC) and the Servicio Nacional de Meteorología (SMN)). Grey isotachs represent the mean annual wind speed in m per second (ms −1 ), and grey arrows, the mean annual wind direction for the period between 1980 and 2018 at 1000 mbar pressure level [19]. (B) The fjord-rich west-east striking mountain range of the Cordillera Darwin includes the largest icefield of Tierra del Fuego (Landsat-TM image mosaic, bands 5-4-3, acquisition date: 14-02-2005). (C) According to the west-east elevation profile of the Cordillera Darwin (transect displayed as red line in (B), Monte Sarmiento is the first prominent peak located on the windward, hyper-humid side of the mountain range; glacierized terrain is displayed in blue.
The Little Ice Age (LIA) was a period of global glacial expansion from the 13th century to a maximum between the 17th and 19th centuries, and it is well documented [20,21]. However, the LIA showed a high spatio-temporal variability in a global context [4,20,22], and glacial fluctuations, e.g., for the Southern Hemisphere, do not completely coincide with those from the Northern Hemisphere [23]. In fact, significant differences in spatio-temporal glacial advances have emerged even within the respective hemispheres [20]. As a consequence, the term of a "global phenomenon" is controversial [24]. A considerable variability of the timing of the LIA maximum on decadal and centennial scales has been recorded [11,25]. Hence, a strong need for accurate dating of glacial fluctuations is apparent [11,[26][27][28]. Indeed, the glaciers in the Southern Andes reached their LIA maximum extent at various times across several centuries, culminating between the 17th and 19th centuries [11,29]. Although the number of paleoclimatic studies in Patagonia has increased in recent decades, 85% of the implemented studies are still located on the eastern side of the Andean main crest [30]. Late Holocene glacial fluctuations have scarcely been documented for the glaciers in the Magallanes region; thus, the different phases of the "Neoglaciations" in Southern America are still corroborated by an insufficient number of studies [11,20]. With the exception of the LIA period, almost no information is available on glacial fluctuations within the last 2000 years. The Common Era (CE) is frequently subdivided in paleoclimatic periods based on (historical) observations or reconstructions conversant from the Northern Hemisphere [31]. The LIA is commonly preceded by the Medieval Climate Anomaly (MCA), which is considered to be a warmer period of approximately 400 years in duration (800-1200 CE) [32]. The Dark Age Cold Period (DACP, 300-800 CE) chronologically prior marked a distinctly cooler period which superseded the Roman Warm Period (RWP, 0-300 CE) [32].
Schiaparelli Glacier, within the Monte Sarmiento Massif, is located in one of the most western parts of Isla Grande de Tierra del Fuego ( Figure 1B,C). The glaciers of this pronounced mountain peak descend from over 2000 m asl to almost sea level. West to Mt. Sarmiento, the Magdalena Channel directly connects the Strait of Magellan with the Pacific Ocean. This area has attracted scientific interest since the famous Beagle Expedition accompanied by Robert FitzRoy and Charles Darwin in 1836 [33]. Lithographs from this expedition supposedly show the Schiaparelli Glacier calving into the Magdalena Channel, leading to a dispute over the glacier's extent in the 19th century [34][35][36]. To clarify whether the glacier calved into the sea within recent centuries, Masiokas et al. concluded that "[…] a more detailed analysis of this site, including the tree-ring dating of trees growing on the frontal moraines is needed to resolve this question" [11] (p.257). The formation of moraines can be dated by the germination date of trees colonizing recently exposed glacial deposits [37]. In this study, we performed a tree-ring-based moraine dating of the terminal moraines at Schiaparelli Glacier: maximum tree ages constitute minimum ages for LIA maxima. Since an ecological proxy was used for dating, it is susceptible to uncertainties, e.g., a variable period between glacial recession and floral recolonizations (ecesis time) [38] or uncertainties due to missing the pith of the respective trees (pith offset) [39]. Commonly, dating uncertainty arising from the use of tree-rings on moraines is a quarter of a century or even more [26,40]. To compare phases of culminated glacial activity, the respective influences of the different uncertainty values on the dating accuracy must be disentangled. Consequently, one focus within the study was to provide a robust estimation of the uncertainties that emerge in tree-ring-based dating of glacial deposits. With these results, the plausibility of a fjord-calving glacier in 1836, as is inferred from the historical record, was verified. Furthermore, the study tried to constrain the timing and phases of different Neoglacial advances in Tierra del Fuego. Sub-recent trunks that were incorporated within morainic sediments were radiocarbon dated to interlink local glacial activities with the culminating glacial advances of southernmost South America during the Late Holocene.

Study Site and Regional Climate Conditions
In 2016, the 11 km long northwest facing Schiaparelli Glacier covered an area of 24.3 ± 0.3 km 2 [8] and descended from 2207 m asl (Mt Sarmiento) to 11 m asl. The glacier calves into the moraine-dammed proglacial lake Lago Azul (~1.2 km 2 ). The lake was formed after the glacier receded in the 1940s from its position located close to the terminal moraine system (Figure 2A). Satellite imagery and historical observation reveal that the Schiaparelli Glacier has been subject to a perceptible retreat within the recent decades. This clearly coincides with the widely observed behavior of landand lake-terminating glaciers in the Cordillera Darwin [8]. While attempting to ascend the summit of Mt. Sarmiento in 1913, De Agostini [41] identified the glacier close to the terminal moraines behind a dense forest. During the following decades, the glacial terminus remained in the vicinity of the moraine system [35,36,42]. However, a continuous monitoring of glacial dynamics has only become possible since 1973, when Landsat MSS scenes first became available. As deduced from a TanDEM-X digital elevation model (DEM, acquisition date 13-05-2011, Microwaves and Radar Institute, German Aerospace Center (DLR), Oberpfaffenhofen), the terminal moraine system is composed of at least three major moraine ridges, reaching a maximal height of 35 m. The outermost moraine is located 2.2 km in front of the current glacial terminus. The moraine ridges are cut centrally by an approximately 60 m wide lake discharge, connecting Lago Azul with the Magdalena Channel to the west (Figure 2A, B).
The glacial forefield, especially the depressions between the moraine ridges, was characterized by abundant swampy terrain interspersed with several peat bogs. The hummocky terrain was colonized by a closed tree stand of Nothofagus spp., homogenous in height and diameter, indicating a simultaneous germination, i.e., subsequent to glacial retreat ( Figure 2C). Recently, several tall individual trees had been uprooted by the wind, but no noticeably old and rotten trees were encountered on the ground. This implies that the forest is still at the stage of primary succession. Beyond the morphologically distinguishable moraine system, the forest stand was much more heterogeneous in terms of trunk diameter, species composition, and dead, decayed trees (climax forest).  [8,42] and their corresponding moraine ridges (light grey) and forest cover (green [43]). (B) Cross-section of the terminal moraine system as indicated by the red bar in (A). Three major terminal moraine ridges were identified (TM1-TM3). (C) Peat bog formation between the moraine ridges, colonized with Nothofagus betuloides trees of comparable diameter. (D) Climate diagram calculated as a mean of four CHELSA grids cells covering the end moraine [18]. Schiaparelli Glacier is surrounded by west-east striking lateral moraines at the northern and southern slopes, which reach elevations of 380 meters above sea level (m asl). In the upper hillside of the moraines, the forest stand appeared older and homogenous in its age structure, composed of N. betuloides and N. antarctica. Downslope of the lateral moraines towards the recent glacial margin, the tree stand was more scattered and small in growth, indicating a primary succession after glacial thinning. Above the lateral moraines, the matured forest was replaced by Andean tundra, composed of mosses, gramineous, lichens, and some low shrubs. The overall high annual wind velocities in Tierra del Fuego ( Figure 1A) can intermittently exceed 30 ms −1 on exposed slopes, which can be decisive in forest distribution [44]. In sun-exposed and occasional wind-protected areas, the treeline can reach an elevation up to 550 m asl [15,45].
The late successional, undisturbed forest stand was mainly composed of the endemic evergreen southern beech (Nothofagus betuloides), with scattered individuals of Drymis winteri, Tepualia stipularis, and Nothofagus antarctica [46]. Since Nothofagus spp. saplings are shade intolerant, but can cope well with difficult nutrient recovery or soils containing low nutrient concentrations, they represent the pioneer woody species on recently exposed moraines after deglaciation. In addition, N. betuloides and N. antarctica grow superiorly on waterlogged soils; thus, they frequently occur in the glacial forefield [47][48][49]. Tree germination of Nothofagus spp. is epigeous and occurs in the spring of the year the seeds are produced [16,45]. The species composition represents the ecoclimatic hyper-humid zone of the sub-antarctic broadleaf evergreen rainforest, exceeding precipitation sums of 1000 mma −1 [47,50]. The dominance of N. betuloides is associated not only with a humid oceanic climate, but also with a cold-temperate climate [45].
The permanent cold air emanating from the ice-covered Antarctic continent, combined with the upwelling of cool Antarctic deep water, prevents substantial warming of the polar air south of 50°S and causes overall cool temperatures [51,52]. The daily and seasonal cycle is dampened by the immediate proximity of the ocean and the persistent westerlies. Since there are no climate stations in the vicinity that provide long-term data for a sound climatological classification ( Figure 1A), we had to refer to the CHELSA reanalysis dataset (Climatologies at High Resolution for the Earth's Land Surface Areas), downscaled from ERA-Interim to 30 arc sec [18] ( Figure 2D). Temperatures were relatively homogenous throughout the year (mean temperature 5.9 °C), with an annual amplitude of merely 6.6 °C, and maximum and minimum temperatures during Jan/Feb and Jun/Jul of 9.2/9.0 °C and 2.7/2.6 °C, respectively. Interannual variations can be mostly attributed to the leading modes of atmospheric circulations: the Antarctic Oscillation (AAO), also known as the Southern Annular Mode (SAM) [53,54]. The SAM is defined by pressure anomalies of opposite signs over Antarctica and midaltitudes, influencing the strength and latitudinal position of the westerly wind belt. During positive SAM phases, the southern westerly wind belt strengthens and shifts poleward and vice versa [55,56].

Methods for Dating Glacial Deposits (Based on Tree-Rings and Radicarbon Dates)
As soon as a formerly glacierized terrain becomes exposed due to glacial recession, tree seedlings colonize the bare ground, progressively tapering towards the glacial front [57,58]. The oldest individual trees found on a moraine can be used as a minimum estimation to date the formation of a moraine [38]. Nevertheless, the exact dating of the glacial event is influenced by some uncertainties that must be taken into account [39,59,60]: (i) the pith of the tree was not reached or missed during sampling (pith offset); (ii) the number of tree-rings that are lost when the tree is sampled at a distinct elevation above the ground; (iii) the timespan (ecesis time) between deglaciation (surface exposure), moraine stabilization, soil development, and the colonization by the first trees; and (iv) the probability that the oldest tree was sampled.

Establishment of Tree Ring Site Chronology and Evaluation of Missing Rings
Tree-ring width was measured with a precision of 0.01 mm after preparing the increment cores with a razor blade and chalk for better visibility of light rings. Since the vegetation period, and thus tree growth, extends over the turn of the year, the Schulman convention [61] is used for the Southern Hemisphere, whereby the year in which tree growth starts is used for dating. As recommended by Koch [39], all sampled cores were cross-dated statistically, and the established site chronology was verified against data from the International Tree-Ring Data Bank (ITRDB) using the open source statistical language R [62] and the dplR package [63,64].
If the increment borer missed the pith of the tree, the number of absent rigs must be estimated. The most widely applied geometric model, developed by Duncan [65], uses the curvature and width of the innermost visible ring to estimate the length of the absent radius. The number of missing rings is then calculated by dividing the estimated length of the missing radius with the mean width of the innermost 5, 10, and 15 visible tree rings [60,66]. If multiple radii are available, they are averaged to derive the mean ring width of the respective trees, and to additionally check for differences in eccentric tree growth [67]. The obtained correction factor is then added to the date of the innermost tree-ring to receive the age at coring height. As a prerequisite for this geometrically-based assumption, radial and symmetrical tree growth is required to estimate the length of the missing radius, leading to higher uncertainties with increasing numbers of missing rings [60,68]. So far, age uncertainties induced by the missed pith have only been estimated in a few studies, and must not be generalized across species or continents [60,67]. The effect of the intervals (first 5, 10, and 15 years) for calculating the mean ring width and the inferred pith offset (PO) was evaluated for 24 pairs of increment cores, one of which hit the pith of the respective tree, and another that missed it at a comparable tree height.

Growth Rate Modeling of Nothofagus spp.
Trees cored at a distinct elevation above the root collar do not reflect the date of germination and always underestimate the true tree age [69,70]. Since the root collar is located at or even below the ground level, an exact estimation of the tree age is highly time-consuming and destructive to the individual trees [70]. The study site was located in the Alberto de Agostini National Park as the core zone of the Cape Horn Biosphere Reserve, where it is prohibited to perform destructive sampling, such as taking cross sections at the root crown. Thus, we established a site-specific age-height relationship for Nothofagus spp. by sampling six young trees at different elevations above the ground level and counting the age differences between the distant cores of the respective trees. The dense and already matured forest of the terminal moraine was not suitable for this purpose, because a closed canopy influences the growth rates of young saplings [71,72]. In order to ensure comparable successional growing conditions, we selected the southern lateral moraine, which had been recently exposed after glacial recession, encompassing minimal soil development and scattered individual trees. Growth rates may vary considerably, as the tree stand may include highly stressed individual trees. In such cases, age estimations are prone to underestimating true tree age by more than 35 years [73]. For this reason, we carefully examined the increment cores extracted from the terminal moraine system to identify outliers in radial growth (more than 1.5 standard deviations calculated for the first 15 years in tree growth). To those outliers, an additional age-uncertainty of ±5 years was added, and it was ensured that they were not solely used for the dating of a respective moraine.

Estimation of the Site-Specific Ecesis Time
The ecesis time is spatio-temporally variable and does not follow a superficial dominant factor: it is rather the result of the combination of large-scale and local climatic influences, sufficient seeding availability, e.g., distance of seeding trees, seed dispersal, and the characteristics of the bedding substrate [38,57,74]. Previous estimations of the ecesis time in the Patagonian Andes and Tierra del Fuego have revealed a high temporal variability, ranging from less than one decade [75], to several decades [15], to even a century [34,76]. Even for the same glacier, adjacent slopes are colonized at rates differing by up to 30 years, depending on aspect and inclination [40]. The divergence of the tree colonization prevents a generalized application, and subsequently requires a site-specific determination of the ecesis interval [39,73]. Since a proglacial lake was located between the terminal moraine system and the recent glacial front, it was not possible to assess the ecesis time directly in the glacial forefield. Thus, the only possibility by which to examine tree recolonization after deglaciation was in the areas along the lateral moraines.
To evaluate which of the two lateral moraines was more similar in climatic conditions to the terminal moraine, we calculated the annual incoming solar radiation for the western side of the Sarmiento Massif. The calculation was based on a geometric solar radiation model using a digital elevation model as input, as implemented in the Solar Analyst tool in ArcView GIS (ESRI), developed by Fu and Rich [77,78]. The northern terminal moraine received only one third of the incident solar radiation due to cast shadows and steep slopes, whereas the southern lateral moraine was exposed to similar solar input as the terminal moraine system. Consequently, the ecesis time was estimated on the basis of samples from the southern lateral flank of the glacier.
On the southern lateral moraine at ~200 m asl, a significant change in the tree stand from a more or less dense/matured forest to a scattered tree stand was observed. That ecotone indicated a moraine stage that was exposed recently. Assuming a constant and linear glacial retreat, trees should be progressively more juvenile towards the glacial tongue. The individual trees that appeared the oldest along the moraine were sampled in transects between the recent glacial margin and the dense forest stand ( Figure 3). Along this small lateral moraine, 81 young Nothofagus spp. were sampled at approximately 0.3 m above ground level. At least two cores were taken from each tree using an increment borer, and used for establishing a site chronology. The positions of the sampled trees were mapped using GPS with an accuracy of ±3 m and transferred to a geographical information system. Thereby, the spatial distance between the glacial outlines and sampled trees was precisely determined. Glacial retreat and tree succession was monitored by combining optical satellite imagery and tree-ring sampling for the last 40 years. To ensure that the succession was not disturbed by glacial advances, as they were observed elsewhere in the Cordillera Darwin [8,9], we performed multi-temporal observations of the glacial margins. For mapping glaciers and their variations, multi-spectral satellite imagery of sufficient quality has been available since 1973. However, reliable and sufficient imagery for southernmost South America is rare due to irregular satellite overflights and frequent heavy cloud cover. Therefore, 18 optical satellite scenes acquired from different sensors were used to assess regional glacial extent and variability (Appendix Table A1). Scenes of the Landsat mission (MSS/TM/ETM+ and OLI) and Terra ASTER, acquired from the United States Geological Survey (Level 1GS/TP), were used for a semiautomated mapping of the glacier's extent [8]. Additionally, satellite scenes from the Sentinel-2 MSI program of the European Space Agency (ESA) that has provided data since 2015 were used. All scenes were projected to the local UTM zone (19S) and checked for the orthorectification quality. Ratio images were calculated by dividing the raw digital numbers of the red and shortwave infrared channels [8,79]. We used the equivalent bands 2 and 4 for calculation of the ratio images in the case of the Terra ASTER satellite imagery. Due to the different resolutions of the respective ASTER bands of 15 m and 30 m, Band 4 was resampled to a 15 m resolution, as suggested by Kääb et al. [80]. A suitable threshold value was used after a stepwise increase of 0.1, starting from 1.0, to convert the ratio images into a binary image. By vectorizing the binary image, glacier polygons were generated and compared to false color composites of the mid-infrared, near infrared, and red channels. Afterwards, the glacial outlines at Mt. Sarmiento were subdivided into individual catchments using a void-filled UTM-projected DEM with a resolution of 30 m (SRTM, LP DACC NASA Version 3). The mapping accuracy was influenced by (i) the geolocation and (ii) the scene quality [81,82]. The positional error (i) was reduced by aligning the scenes used to distinct landforms like shorelines and prominent mountain peaks. The Sentinel-2 and most Landsat TM/ETM+ and OLI data matched perfectly, whereas the Terra ASTER and Landsat MSS appeared to have a small offset and were adjusted. Since scene quality (ii) is negatively influenced by seasonal snow cover, satellite scenes between the spring and autumn season, representing the seasonal snow minimum, were mainly used. In the case of partial cloud cover and occultation of the glacier, the next available and suitable satellite scene was used to supplement the cloud-covered parts of the original image.
The change and resulting uncertainty (u) in the position of the glacial margin by comparing two different satellite scenes was dependent on the respective image resolutions (r1 and r2) and the co-registration error (RE) due to the georectification process [83,84].
With a pixel size of 60 m, the Landsat MSS sensor had the lowest spatial resolution compared to the Landsat TM/ETM+/OLI (30 m) or the Terra ASTER and Sentinel-2 scenes (15 and 10 m, respectively). The largest possible uncertainty occurred between Landsat MSS (1973) and Landsat TM (1986) and accounts for 77 m, of which 10 m resulted from the RE. Hence, this uncertainty estimation is the worst-case value; we additionally estimated the precision of our semi-automated derived outlines by applying a positive and negative buffer to the glacial perimeter [85]. The buffer was applied to the glacial outlines before they were subdivided into individual catchments to prevent any uncertainty estimation where glaciers adjoined [81]. A buffer size of ±0.5 pixel of the sensor used (according to ±30, ±15, ±7.5, ±5 m) was well suited for the uncertainty estimation of clean-ice bodies [8,86,87], leading to a smaller maximum uncertainty than assuming that the complete pixel might be faulty.

First-Generation Trees
On older moraines, a very pronounced, old-grown dense forest can have evolved. In such areas, the identification of the oldest tree is not obvious due to its diameter, height, or branch spacing. Through a high number of sampled individual trees (>150 trees) we optimized the likelihood of identifying the oldest tree(s) growing on the moraine ridges [60]. The moraines were colonized by an uniform, dense forest stand which was comparable in diameter and height, with only a few dead trees and an abundance of the Nothofagus spp., typical for postglacial succession and first-germination trees [88]. Trees growing on the distal face of a glacial deposit can be much older than the ones on the proximal side of the moraine, because the latter can only germinate if the glacier completely retreated without any re-advances. Furthermore, the germination can be negatively influenced by katabatic winds [59]. If the glacier is still in the vicinity of the moraines, dumping and sliding processes of debris and meltwater may induce disturbances that additionally impair the natural recolonization. In order to avoid misinterpretations of the glacial position, trees were sampled at the respective moraine crest and at both the distal and proximal moraine faces.

Radiocarbon Ages
Radiocarbon ages ( 14 C BP) received from the logs incorporated into the glacial deposits were calibrated using the CALIB 7.1 software [89] and the SHCal13 calibration curve for the Southern Hemisphere [90]. The median probability of the distribution function was used to report the calibrated radiocarbon age, whereby the calibrated date range is given as 2σ-ranges [91]. If the calibration provided multiple non-overlapping age segments, the period with the highest probability was used. The date notations of the radiocarbon dates are expressed as calibrated years before or after the Common Era (cal BCE/CE) or as calendar years before present (cal ka BP: ka = 10 3 cal. years BP, present = CE 1950).
Since different proxies such as tree-rings or radiocarbon data typically use different date notations, we have reported all dates within the manuscript in a uniform convention, using the Common Era as reference system. The calibrated radiocarbon ages are always denoted with the prefix "cal," and are rounded to the nearest ten years and are expressed as highest probability for better readability. The original, exact calibrated and uncalibrated dates are listed in Table A2 and Table A3.

Dating Uncertainty Introduced by the Estimation of the Pith Offset (PO)
Initially, the influence of the different intervals used to calculate the average ring width of the innermost rings was analyzed. The mean ring width varied insignificantly for the first 5, 10, and 15 years, leading to a maximum deviation that did not exceed 0.8 years if the PO was calculated afterwards. This was mainly attributed to the uniform growing conditions of Nothofagus spp. within the first years and the absence of strong inter-and intraspecies competition, e.g., for light, as the early successional forest does not provide a closed canopy.
The estimated PO and true PO were significantly (r = 0.64, p < 0.01) correlated, as the simple linear regression model revealed ( Figure 4A), but showed higher variance with increasing distance to the pith. This was confirmed by the comparison of the estimated PO with the resultant PO errors (POerror = POestimated− POtrue) in the lower panel of Figure 4A: if the estimated PO was less than 10 years, the PO error resided between −8 and +5 years. This led to an over-or underestimation of the true PO by an average of 2.8 years (σ = 2.0). Due to the low number of observations, we calculated the expanded uncertainty by multiplying the combined standard uncertainty with an appropriate coverage factor (for 0.99 confidence), leading to an average PO deviation of 2.8 ± 1.4 years. For the trees with an estimated PO of less than 10 years, we consequently applied an uncertainty value of ±4 years. Verifiably, the number of incorrectly estimated annual tree-rings increased with growing PO, reaching maximum values of −10 and +12 years, respectively. In accordance with the above procedure, the average PO uncertainty at 0.99 confidence limit accounted to 7.3 ± 2.0 years, leading to an uncertainty value of ±9 years for trees with an PO exceeding 10 years. According to Hochreuther et al. [68], we restricted the reliable PO estimation and disregarded trees with more than 20 calculated missing rings to prevent an over-or underestimation of the true tree age.

Growth Rates of Nothofagus Spp.
The Nothofagus spp. pioneer species located at the lateral moraine were sampled between 0.30 and 0.40 m above their root collars, with an average sampling height of 0.33 m above ground level ( Figure  4B). The age of the lowermost extracted cores was used as the starting reference, but did not represent the true tree age. Subsequently, the temporal difference to the next distant core, which was sampled further up the trunk of the respective tree, was calculated. The cumulative temporal differences (i.e., the effective ageing) depending on the respective height are displayed in Figure 4B. Variations in tree height accumulate as the trees mature, resulting to only minor variations at a young age (less than 0.1 m for the first year of growth) and larger variations (0.3 m) for trees aged over six years. This indicates higher age uncertainties with increasing sampling height above ground level. Finally, annual growth rates were expressed by calculating the differences in height and counted years of the extracted cores ( Figure 4B).
Growth rates are highest as a sapling, reaching values exceeding 25 centimeters per year [cma −1 ]. With increasing maturity, growth rates show a declining trend (r = −0.84, p < 0.01) and reach minimum values of 14 cma −1 . The average growth rate amounts to 20 cma −1 , whereas saplings smaller than 1.2 m show an enhanced annual growth rate of 23 cma −1 , and trees larger than 1.2 m, a distinctly reduced growth rate, corresponding to 18 cma −1 . From a physiological perspective, the sampled trees did not show any growth anomalies such as narrow rings. Additionally, the development of multiple branches at the base of the trunk or compression wood, which would indicate straining and detrimental growth conditions, was not visually apparent. The increase of the tree diameter (~0.3 m above ground level) for the 81 collected trees was quite constant (r = 0.83, p < 0.01) and showed only minor deviations within the first 15 years, but increased with aging of the tree ( Figure 4C). The apparently steady and undisturbed growth patterns at young age on the lateral moraine justify a generalization of the calculated growth rates to the entire set of sampled trees.
Considering the average growth rates for the first three years (24 cma −1 ), the trees that were sampled at the lateral moraine about 0.3 m above the ground must have been at least two years old. Sampling height at the terminal moraine of 1.2 m might vary slightly ±0.1 m (indicated as hatched field in Figure 4B) due the rugged and hummocky terrain or interfering branches. As depicted in Figure 4B, the two year old trees sampled at the lateral moraine reached a coring height of 1.2 m within a further 5 ± 1 years. In consequence, the trees sampled on the terminal moraines would reach, under comparable growing conditions, their coring height within 7 ± 1 years.

Site-Specific Ecesis Time
The average glacial retreat on the southern lateral moraine was calculated as the mean of the orthogonal distance every 30 m to the respective glacial margins. From 1973, Schiaparelli Glacier receded continuously in the lateral direction by 215 m (r = −0.99, p < 0.01) until 2018, which corresponds to an annual loss of approximately 5 ma −1 ( Figure 5A,B). The high error values given in Figure 5 (light grey) describe the largest possible error, caused by the assumption that each of the individual examined pixel pairs of two satellite scenes to be compared is faulty [83], which is, concerning the location of the seedlings, rather unlikely. If, on the other hand, the uncertainty between two satellite scenes were calculated using the applied buffer, which assumes that the accuracy of the glacial boundaries is ±0.5 pixel of the used sensor, the calculated error for the lateral glacial retreat is significantly lower (dark grey).
The germination date as dependent on the distance to the glacial margin in 1973 is displayed on the right ordinate axis in Figure 5B. Trees with a positive distance germinated outside the 1973 glacial extent. All trees with negative values germinated on freshly exposed terrain subsequent to glacial recession. The succession of Nothofagus spp. along the 1.3 km long sample site was rather swift and extensive. The trees that appeared the oldest were sampled, and the simple linear regression analysis (r = −0.93, p < 0.01) revealed that the saplings linearly traced the glacier's retreat (r = −0.99, p < 0.01) ( Figure 5B). According to the two linear models, the germination of the trees took place parallel to the glacial retreat with a time offset of 15 years. Although the colonization of some individual trees on deglaciated deposits can occur very fast (<10 years), the assumption of such a short period for a generalized ecesis period is hardly reasonable. On the one hand, the number of successional trees that germinated shortly after the retreat was relatively small, with the possibility of sampling them later in a dense forest as well. On the other hand, the estimates of the ecesis time recorded on the lateral moraine weretransferred to the frontal moraine system. Due to sliding and dumping processes, the morphological moraine patterns of the lateral moraines were less dissected than those of terminal moraines [92]. Glacial oscillations can lead to complex pushing, thrusting, and dumping processes at the glacial terminus, combined with melt-induced collapses that may result in steadily disturbed moraine formations [93,94]. It is likely that tree colonization will be delayed in such unsettled terrain compared to the relative flat or undisturbed lateral moraine. Therefore, the shortest possible time interval derived from the few first germinated individual trees was not used, but an averaged time delay of 15 years, as revealed by the regression analysis, was used instead. An associated uncertainty value of ±3 years, derived by the standard deviation of the linear model, was applied, resulting in an ecesis interval of 15 ± 3 years.

Tree-Ring-Based Moraine Dates
The root-mean-square error (rmse) for the tree-ring dating of the terminal moraine can be described as the quadratic sum of the independent uncertainty of the PO estimation (±4 or ±9 years), the growth rate model (±1 years), and the ecesis time (±3 years), leading to a minimum and maximum uncertainty of ±5 or ±10 years.
The forest distribution was strongly linked to the presence of glacial deposits which seemed to favor the formation or quality of soils, as the closed forest was exclusively located on glacial debris. At altitudes above 100 m asl, the appearance of trees was solely attributed to the presence of moraines, which simultaneously represented the local tree line ( Figure 6A). The digital elevation model of the TanDEM-X scene indicated at least three terminal moraine ridges north and south of the lake discharge. The field survey revealed an irregular distribution of the forest within the terminal moraines: a vast area north of the lake runoff in which numerous ponds appeared, with trees growing in a scattered pattern. According to the historical information, this area was still glacierized, or had just been recently deglaciated, in the 1940s [36,41] (Figure 2A), which was reflected by the occasional vegetation typical of the successional stage, as was also observed after the subsequent glacial retreat of 1973 ( Figure 6A). This was supported by the comparison of Sentinel-2 scenes acquired from different seasons, which revealed that the generally dominant evergreen N. betuloides was, to a significant extent, substituted by the broadleaf deciduous N. antarctica ( Figure 6A). The latter species is habitually adapted to higher elevations (>100 m asl) but is also abundant in flat and wet terrain, e.g., directly along lakeshores, wetlands, and on recently exposed terrain [16].
The presence of multiple ponds, the simultaneous onset of peat bog formation, and the scattered appearance of N. Antarctica, which often comprises several trunks or forms belts of krummholz [95], appeared unfavorable for tree-ring-based dating of the glacial deposits. In consequence, our moraine dating was concentrated on the area south of the outflow (Figure 6B), as the forest in the area appeared to be older and the terrain was less unsettled. Additionally, the shape of Lago Azul reflects the glacier's principal ice flow direction and is aligned on its elongated axis exactly towards the dated terminal moraine system, supporting the hypothesis that these moraines are representative of former glacial activity.
Along the three terminal moraine ridges (TM1-TM3), the maximum tree age increased with increasing distance to the glacier, whereby the outermost terminal moraine encompassed the oldest trees and the innermost, the youngest ones. This implies a glacial recession without distinct readvances that would superimpose the TM1 moraine since its deposition date. The oldest sampled tree individual, germinated in 1749 CE, was located on the distal side of the TM1 moraine crest in the southeastern part, representing the maximum date of moraine stabilization ( Figure 6B, C). Two additional trees nearby on the proximal slope of the moraine germinated in 1763 CE and 1764 CE, representing the minimum deposition date of glacial sediments. However, the mean germination date of the five oldest individual trees from the proximal moraine slope was about 10 years younger than those from the distal slope ( Figure 6C, lower panel). This suggests that the terrain was either still icecovered for a longer period, or that the glacial front was located nearby. Glacial oscillations, as well as unfavorable site-specific climatic circumstances (e.g., katabatic winds), meltwater streams, or collapse structures introduced by dead ice may impede otherwise rapid colonization. The clearly discernible TM1 moraine ridge diminished in contour in a northwesterly direction and transformed into a diffuse hummocky terrain with multiple small ridges, with the oldest trees in that area dated to 1769 CE.
Tree recruitment on the central moraine TM2 seemed to be most heterogeneous, with the largest age range of germination dates between 1762 CE and 1921 CE ( Figure 6B, C). Similar to the TM1 moraine, trees along the moraine belt towards the lake discharge (in northern direction) were slightly younger. South of the recent lake drainage, the TM2 moraine belt was interrupted by an abrupt terrain flattening and the occurrence of swampy terrain. Simultaneously, tree age was youngest in this area (average germination date 1895 CE), and the individuals appeared stunted and possessed multiple trunks, indicating harsh growing conditions. Some individuals showed clear signs of root exposure, e.g., due to fluvial geomorphic processes, such as the presence of a former lake outflow ( Figure 6B). The oldest trees on the distal slope in the southeastern part of the TM2 moraine dated back to 1760 CE, while the oldest trees towards the north were 20 to 30 years younger, corresponding to 1782 CE and 1786 CE. Tree germination on the distal slope of TM2 coincided with the proximal TM1 moraine side around 1760-1775 CE, whereas tree recruitment on the TM2 crest and proximal side had a relative lag of approximately 15 years ( Figure 6C).
This observation was decisive, since the TM2 moraine belt was almost radially dated before the 19th century, and a calving of the glacier into the Magdalena Channel, as outlined by the Beagle Expedition's lithograph in 1836 CE [33], seems rather unlikely.
Tree recruitment on the terminal moraine closest to the glacier (TM3) was dated to 1850 CE and possessed the smallest age range. Only one older tree was dateable at the easternmost site (1807 CE). The tree position was situated in a transition zone, where the narrow, barely distinguishable lateral moraines fanned out into defined terminal moraines. This individual tree distal to the TM1 moraine was therefore attributed to the proximal slope of the TM2. Around 1900 CE, the glacier receded so far that the proximal TM3 side was recolonized with trees.

Radiocarbon Dates/Glacier Chronology of Southern America
In total, 14 trunks (Table 1) incorporated within the terminal moraine deposits of the Schiaparelli Glacier were discovered, of which five were located directly on the lakeshore south of the lake discharge. The remaining trunks were exposed along the moraine cross section parallel to the river bed. The disordered location of the trunks amidst glacial debris suggested that the trees were killed and deposited during glacial advances. Samples for radiocarbon dating were taken from the trunk exterior to determine the glacial event more precisely.
After tree mortality, sapwood is less resilient to environmental influences than heartwood, which may cause dating uncertainties. However, we observed one trunk still showing morphological signs of an intact sapwood area, leading to the conclusion that the trees were covered rapidly after they were overridden and tilted by the glacial advance. Consequently, we assumed that dating uncertainties were of subordinate importance. With regard to the slope morphology and the low tree line, we considered it unlikely that the trees might be remnants of a forest that collapsed from the lateral slopes after glacial thinning, or was eroded by glacial advance and later transported and accumulated by the glacier. We hypothesize that the trees grew in the apron of the glacier near the terminal moraines, and were buried and transported by glacial advances more or less in situ. Statistically, multiple 14 C dates were available for the same glacial/deposit event. The standard deviations of SC 1-5 could be clustered according to their time overlap and attributed to one single glacial event. The same applied to samples SC 6-9 and SC 10-11, whereas three single glacial oscillations (SC 12, SC 13, SC 14) at 880 CE, 1000 CE, and 1220 CE could be dated.

Late Holocene Glacial Oscillations in Southernmost America
Even though temperatures in the Early Holocene were up to 2 °C higher than in the late 20th century [25], it is highly unlikely that the glaciers disappeared entirely at this time. High precipitation amounts cause a considerable accumulation rate on the hyper-humid side of the Andes. Even under the recent climatic conditions, elevation changes for the glaciers of Tierra del Fuego region have been balanced or even positive above 900m asl [96]. In consequence, the mass and area loss mainly affected the glacial tongues and low elevated glacial areas, leading to a severe glacial retreat. After this prolonged phase of glacial recession in the Early to Middle Holocene (8200-4000 BCE), many glaciers in the world were subject to a glacial resurgence that was initiated between the Middle and Late Holocene (~4200-2500 BCE) [14,97,98]. This period is referred as the "Neoglacial" [97,99,100], and the cooling within the Neoglacial lasted until ~1900 CE. The most recent phase of supra-regional intensified glacial advance is referred as the LIA with regard to glaciation (~1250 CE to ~1950 CE) [20,22]. Irrespective of the definition of the LIA, namely whether a glaciological or a climatological concept is used, the timing and duration is globally and interhemispherically highly variable [22,101].
For the Patagonian Andes the most common glacier chronologies were published in 1965 and 1995/2013 by Mercer and Aniya, respectively, postulating three to five phases of glacial advances, which are named in decreasing age as Neoglacial I-V (NG I-V) [102][103][104]. Nevertheless, lately, studies have provided important evidence for the occurrence of up to seven phases of culminated glacial expansions in the last 10,000 years [14,105,106]. Most of the data on which the glacial chronology of Southern America is based are represented by radiocarbon-dated glacial deposits. The inferred 14 C ages are frequently based on single measurements, and it cannot be ruled out that organic material has been remobilized or overprinted by Late Holocene advances. Therefore, glacial events/chronologies in Southern America are not always well constrained in time [4,30,97], and should be best regarded as broad regional trends [25].
Glacial advances within the last 2500 years are compiled in Figure 7 for the Southern Patagonian Icefield, including adjacent mountain glaciers (~48-51° S) and the Magallanes-Tierra del Fuego region (~52-55° S). Several phases of culminated or enhanced glacial activities can be inferred (displayed as blue vertical bars in Figure 7A,B, Table A2) [4,11,16].
The lower part of Figure 7 compiles paleoclimatic studies from southernmost South America that represent favorable or unfavorable climatic conditions for glacial growth. In fact, no simple interrelationship between glacial oscillations and climatic variables exists, but rather, a complex interweaving of climatic and non-climatic factors. Even glaciers that are in contact in the accumulation zones (located side by side) and are subject to the same local climatic influences can exhibit contrasting frontal variations due to different accumulation-area ratios (AAR) [107,108]. In addition, the calving rate of a respective glacier is decisive for its length variation. Calving fluxes represent inherent dynamic instabilities resulting, for instance, from a complex interaction between subsurface geometry, sub-and englacial water pressure, ambient water temperature, and circulation, tide, and waves [109][110][111], and are frequently decoupled from direct climatic influences [97]. Therefore, the evidence of glacial fluctuations and the imprints of climate signals must be critically assessed, especially in regions with tidewater-and lake-terminating glaciers [97,112].
The timing of the Neoglacial (NG) advance between ~550 BCE to 50 CE in southernmost South America is mostly based on the dating of glacial deposits from Upsala Glacier (~49°40′S) and Occidental Glacier (~48°50′S), descending from the eastern and western side of the SPI, respectively [12,34,102]. There has been so far no clear evidence by morainic deposits that the southern glaciers of the Magellanes-Tierra del Fuego region reacted in phase to that NG event. In contrast, at the Gran Campo Nevado Ice Cap, the lack of evidence for glacial advances between 3000 BCE and 1000 CE, raised the question of a possible shift of the westerlies and resultant drier conditions [113]. However, the new data derived from the Schiaparelli Glacier closes this gap in knowledge, indicating at least one glacial advance during that NG, with the oldest date estimated at 280 cal. BCE. The five radiocarbon dated samples SC 1-SC 5 (Table 1, Figure 7A) suggest several glacial advances before 110 cal. CE. This is corraborated by glacial clay deposition derived from a marine geological fjord core at Marinelli Fjord (54°25′S), suggesting glacial advances between 350 cal. BCE and 150 cal. CE, and a subsequent retreat lasting until ~550 cal. CE ( Figure 7C) [114].
The Neoglacial stage ~300-850 CE is reflected in our study by two advances: four individual trees (SC 6-SC 9, Table 1) and their 14 C ages were indicative for a first advance around 300 cal. CE, followed by a second advance between 670-770 cal. CE (SC 10-SC 11). Although the first glacial advance is so far the only evidence of glacial activity at the beginning of this NG in Tierra del Fuego, the latter one coincides with the glacial advance at Ema Glacier (790 cal. CE) [15]. Ema Glacier descends from the eastern flank of Mt. Sarmiento (54°25′S, Figure 2A). The synchronous advance of the two glaciers reveals a similar glacial behavior to a common climate forcing. The NG event is additionally reaffirmed by several studies in the area between 51°S to 54°S, representing a phase of stable/common climate forcing for the period from ~600 cal. CE to ~900 cal. CE, favorable for glacial advance ( Figure 7C) [10,114,115].
In contrast, recent studies at the fjord of Marinelli Glacier (~54°33′S), northern side of CD) have indicated a massive retreat, with its maximum roughly centered between the neoglacial events at ~300 CE [10,114,116]. At Almirantazgo Fjord, north of Marinelli Glacier (Table A2), a significant meltwater event was detected by the use of a grain-size model between 50 BCE and 750 cal. CE, indicating a rapid retreat of the respective glacier combined with a high meltwater discharge [10]. Vanneste et al. [116] found evidence in elevated glacial dust deposits that the Marinelli Glacier receded so far that it was no longer fjord-terminating, but land-terminating in ~250-350 cal. CE ( Figure 7C, Table A2). These results are in disagreement with the dated glacial advance of Schiaparelli Glacier around 300 CE. Nevertheless, the Marinelli Glacier features some peculiarities, e.g., an exceptional high length reduction since its LIA position. The length reduction along the glacial centerline exceeded 14 km and was around seven times as high as that of Schiaparelli Glacier [8]. Concurrently, and in contrast to the surrounding glaciers, high flow velocity and extreme glacial thinning have been observed within recent decades at Marinelli Glacier [96,117]. It appears that non-climatic factors, e.g., bedrock topography, fjord bathymetry, and calving processes, are considerably influencing the glacier's behavior [9,117]. In consequence, we assume that the glaciers descending from the northern part of the CD may respond differently, or even complementarily, than the glaciers of the hyper-humid western and southern mountain range [9]. This points to the fact that paleoclimatic evidence is not valid when generalized, and therefore must account for differences in terms of glacial exposition/location. This is in line with a review paper dealing with the paleoclimatology of southernmost South America, which highlighted the inconsistent/opposite picture that can be derived from the interpretation of various proxies [30].  Table A2) and tree-ringbased moraine dates from the Tierra del Fuego region [16,118]. Blue vertical bars indicate phases of enhanced glacial activity. (B) Frequency of glacial advances at 25 year intervals, including tree-ringderived advances compiled in Masiokas et al. [11]. (C) Paleoclimate studies representing (un)favorable climatic conditions for glacial expansion (Table A2). (D) High resolution palaeoclimatic reconstruction for South America [119], Southern Annular Mode Index [120], Southern Andes Temperature Anomalies [121] (red lines represent the respective 11 year fast Fourier transform (FFT) filter for each time series), and Temperature Anomalies for the Antarctica [122].
Climate reconstructions for the period of the last two millennia have revealed a globally decreasing temperature trend [101,123]. This can be well observed on the Antarctic continent, where temperature reconstructions for the entire Antarctic show a negative trend until the 20th century ( Figure 7D). This cooling has been reinforced by a progressive shift of the Southern Annular Mode towards positives phases since 1400 CE ( Figure 7D) , which also causes a warming of the Antarctic Peninsula [120].
Between the penultimate and last NG advances, glacial activity is noticeably declining. This is also reflected in the small number of dated glacial advances and climatically unfavorable conditions for glacial growth ( Figure 7B,C). This phase coincides with the period of the Medieval Climate Anomaly in Southern Patagonia [4,115,121]. However, Schiaparelli Glacier still presents glacial activity at the onset of this climate anomaly (SC 13). The high-resolution data series are also rather ambivalent: the warmer period (1150-1350 CE), which can be discerned in the temperature reconstruction of South America (Figure 7D), reveals a temporal offset compared to the culminated glacial advances [119]. During this phase of moderate/higher temperatures, the reconstructed Southern Annular Mode (SAM) [120] has also exhibited less negative values. The core zone of the Southern Hemispheric Westerlies is located between 50 and 55°S. In a positive phase of the SAM, this windbelt contracts towards the pole and weakens the winds between 40 and 50°S, and vice versa [1,124]. After 1350 A.D., a significant decrease in the reconstructed temperature and also a prolonged period of exceptionally negative SAM indices can be considered to mark the (climatological) onset of the LIA. Negative phases of the SAM denote an increasing moisture transport with expanded/enhanced westerly winds north of 50°S, which nourish the large ice masses of the SPI and NPI, and are in consequence favorable for glacial growth.
The latest 14 C-dated glacial event at the Schiaparelli Glacier occurred around 1225 cal. CE (SC 14), and corresponded to the start of the Little Ice Age. Most glacial advances in Southern Patagonia during the LIA are dated between 1200-1350 CE, 1560-1675 CE and between the 18th to 19th century [4,14,118,121] (Figure 7). In the Magallanes region at Pia Glacier (54°46′S), two radiocarbon dates reflect neoglacial maxima of multiple moraine ridges formed between 1140-1340 cal. CE. The last date corresponds to a noticable advance of the Ema Glacier, leading to the conclusion that the glaciers south of the Strait of Magellanes advanced during the LIA in the early 13th to 14th century. Paleoclimatic records support this phase of accentuated glacial growth as a uniformly cool/moist phase ( Figure 7C). A tree-ring-based climate reconstruction performed by Villalba et al. [125] in the southern Andes revealed that the LIA temperatures were about −1.5 °C less than today. Since 1800 CE, temperature have followed a positive trend that was only slightly interrupted at the beginning of the 20th century.
Tree-rings as a proxy for glacial activitiy have provided more accurate dating results in recent centuries compared to the radiocarbon method, which features calibration problems [126]. As Nothofagus spp. are long living species and the moraines are quickly recolonized, the tedious task of seeking sub-recent wood for 14 C dating can be omitted. Therefore, the number of tree-ring-based dated glacial advances has increased since ~1600 CE, as also reflected in the absolute number of glacial advances compiled in Figure 7B. In the Tierra del Fuego region, recently exposed glacial deposits were dated using tree rings at only three locations: Gran Campo Nevado Ice Cap (Lengua Glacier), Santa Inés Island (Alejandro and Beatriz Glaciers), and Mt. Sarmiento Massif (Ema Glacier) ( Figure 1A, 7A) [15,16,118]. The oldest tree-ring-dated glacial advances were substantiated to 1628 CE and 1675 CE at Lenguna Glacier and Alejandro Glacier, respectively. All the subsequent glacial advances were smaller in their glacial demension and did not reach the extent of their precursor. Two reworked stumps embedded until 1577 and 1555 cal. CE represent an advance at Glacier Ema. In the forefield of the Schiaparelli Glacier, no evidence was found for an advance between 1400 CE and 1700 CE. However, the moraine formation ~1750-1850 CE for the glaciers of the Mt. Sarmineto Massif appears to have been synchronous. Additionally, the moraine formation coincides with tree-ring-dated moraines from Santa Inés Island [16] and the Gran Campo Nevado Ice Cap [118]. Centered around ~1750 CE, all studied glaciers of the Magellanes region south of 52°S revealed an advance ( Figure 7A). Due to the southward shift of the Southern Westerlies during the positive phases of the SAM in recent decades [55], or as a prolonged period since the end of the LIA, there has been a slight increase in precipitation at the western side of the Cordillera Darwin of +200 mma −1 , whereas precipitation in Northern Patagonia has been significantly declining (up to -800 mma −1 ) [1]. This could lead to glacial thickening in the accumulation area and compensate for possible warming tendencies in Tierra del Fuego, as observed at the Gran Campo Nevado Ice Cap [127].
Although the TM2 and TM3 moraines are situated close to each other at Schiaparelli Glacier, the moraine at the distal shore of the proglacial lake is about 100 years younger, leading to the conclusion that the glacier stagnated over a longer period in this position or that a readvance took place at the end of the 19th century ( Figure 6C,7A). During a four-year British expedition with the ship "Alert" between 1878-82 CE, the Fuegian-Magellanes region was explored, discovering an advance at Glacier Bay (53°20′ S/72°50′ W) and describing it as "(…) a strange sight, standing in the middle of this terminal moraine, to see, on the one hand, a fresh evergreen forest abounding in the most delicate ferns and mosses; and, on the other, a huge mass of cold blue-veined ice, which was slowly and irresistibly gouging its passage downwards to the sea" [128] (p. 125). This glacial advance, destroying the forest after a period of prolonged glacial inactivity, supports the hypothesis of a late 19th century synchronous advance of the Magallanes glaciers.
In the first half of the 20th century, many glaciers located on the hyper-humid side of the CD were still close or at their Late Holocene maxima [9,15]. Whereas most of the fjord-terminating glaciers are calving in constant positions or have even revealed some readvances in recent decades, the laketerminating glaciers are mostly retreating [8]. This is also valid for Schiaparelli Glacier and Ema Glacier, which were situated in the vicinity of the terminal moraines in the early 1940s [15,35,42]. Since then, these glaciers have been subject to a significant retreat. The glacial recessions along the glacier centerlines of Schiaparelli and Ema Glaciers are displayed in Figure 8A, revealing an almost continous retrat of 2000 m during the last ~75 years. The area of Schiaparelli Glacier decreased by 2.8 ± 0.7 km 2 between 1973 and 2019 ( Figure 8B). The area loss has been accompanied by a rapid thinning of the glacial tongue, exceeding values of −4 ma −1 since 2000 [96]. This indicates that Schiaparelli Glacier has been very sensitive to climatic changes in recent decades, and reacts more rapidly than the surrounding fjord-terminating glaciers.  (Table A1). Uncertainty values were calculated according to Equation (1). For historical information of 1943, an uncertainty value of ±100 m was assumed. B) Area loss of the Schiaparelli Glacier. Uncertainty values were calculated applying a buffer of ±0.5 pixel of the sensor used.

Estimations of the Pith Offset (PO)
Previous studies estimating the error range of PO have observed that for small distances (few missing rings), the true PO was overestimated, and in contrast it was underestimated for larger distances [60,129]. Our results support these observations to the extent that the deviations increased with increasing PO. Nevertheless, our study did not reveal a clear assignment of an over-or underestimation dependent on the distance. The uncertainty value for the PO of the Nothofagus spp. was in good agreement with a study recently performed in High Asia [60] indicating that an estimated PO of more than 20 years may underestimate the true tree age by more than 10 years. If moraine dating is based on tree-rings and glacial fluctuations are to be compared on a decadal scale, we highly recommend disregarding PO estimations that include more than 20 missing rings.

Height-Age Calibration
Our results revealed that high averaged annual growth values of almost 0.25 ma −1 for the young Nothofagus spp. encompassed only minor fluctuations. This is in line with growth rates observed within the hyper-humid zones of the evergreen Magellanic forest on recently exposed glacial deposits even exceeding values of 0.30 ma −1 at the southern part of the CD [118,130]. The single available study in Tierra del Fuego examining the growth rates for southern beech was performed at the small Gran Campo Nevado Ice Cap, and revealed similar but more variable growth rates [118]. However, the sampling site was located 100 m asl higher and received up to three times as much precipitation as our sampling site at Schiaparelli Glacier. This may lead to more unfavorable growing conditions, but highlights the resilience and wide ecological niche of Nothofagus spp. Constraining conditions seemed to be absent within the first years of tree growth, as indicated by a relatively homogenous increase of the tree radii. If tree growth is suppressed, miscalculation or underestimation of true tree age can exceed 35 years when cored at breast height [60,73]. High mountain areas possessing trees at the forest ecotone, crossing the treeline up to the krummholz limit, are particularly affected. Therefore, trees should be carefully checked for tree physiological evidence of aggravated growth conditions before they are incorporated in dating, and must be excluded when appropriate.

Ecesis Time
The ecesis time is the variable highly affecting the accuracy of tree-ring-based moraine dating. At the same time, it is also one of the most difficult parameters to determine, since it is usually estimated on the basis of historical sources and the maximum age of the pioneer trees, which are both susceptible to errors [26,104]. Frequently, the succession of pioneer trees on recently exposed geomorphological terrain, e.g., landslides and floodplains subsequent to lake outburst floods, is used as an evaluation of the ecesis time [26,66]. The germination of seedlings is strongly dependent on the bedding substrate, which does not necessarily have to be the same for, for example, floodplains and moraines, and can possibly lead to further estimation mistakes. This high variability in the ecesis period is compiled in Figure 9 for the southern Patagonian Andes. Figure 9. Overview of different ecesis intervals for trees colonizing recently exposed glacial deposits in the Patagonian Andes. Numbers highlighted in red represent the intervals; numbers in square brackets represent the respective reference study (NPI, SPI: Northern and Southern Patagonian Icefield; GCN: Gran Campo Nevado; ISI: Santa Inés Island). [9,15,16,26,27,34,40,75,76,118,[130][131][132][133][134][135][136][137][138][139][140] In the Magellanes-Tierra del Fuego region, the ecesis time appears to be swift, at least for the hyper-humid sites. Bahía Pía, including Pía Glacier, is located 90 km southeast of Mt. Sarmiento, on the southern slope of the Cordillera Darwin (Figure 9). A tree-ring-based ecesis time of 20 years was estimated in the glacier's forefield by accounting the maximum tree age [9]. However, the authors did not provide any information about the exact coring height, PO, and growth rate model. Therefore, the ecesis time could be even shorter. Indeed, a second study carried out 15 years later at Bahía Pía estimated a time offset of only 4 years between glacial recession and tree re-colonization [130]. By comparing aerial photographs and dendrochronology, an ecesis interval of 15-20 years was estimated at the Gran Campo Nevado Ice Cap (Figure 9) [118]. The same method was applied to reconstruct the glacial history at the most western part of the Strait of Magellanes at Santa Inés Island, where a minimum period of 32 years was used ( Figure 1A, 9) [16]. The only significantly longer time span was estimated for the eastern slope of the Sarmiento Massif (Glacier Ema), namely <50 years [15]. Since no exact site-specific estimation was performed, this maximum ecesis period could be slightly shorter, and would rather match the regional values. Nevertheless, those differences appear reasonable when compared with other locations, e.g., Mt. Tronador (~41° S), north of Puerto Montt (Figure 9). Despite the short distance between the glaciers descending from this northern Patagonian stratovolcano, completely different intervals of up to 70 years have been measured [27,75,121]. The colonization of freshly exposed moraines after glacial recession is a complex interaction of different climatic and nonclimatic factors that should not be considered invariant in space and time. Major periods of tree establishment have been associated with warmer phases during the last century, which complicates the calculation of a "typical" ecesis time [141]. Additionally, the proximity to the glacial tongue and the katabatic winds may negatively affect tree recruitment and colonization. Therefore, recolonization periods that were estimated in the 20th and 21st centuries should not be automatically transferred to the period of the Little Ice Age. Our sampling site was located at an elevation of 200-250 m asl, whereas the terminal moraine system to which the ecesis time was applied is located between 1-40 m asl. Consequently, we assume that even despite a possible warming trend, the ecesis time should be in the same order. This presumption is corroborated by observations at Lengua Glacier (Gran Campo Nevado Ice Cap), where killed and tilted trees mark a glacial advance in 1872. The subsequent recolonization with trees was revealed to be in the same magnitude as recent estimations under current climate conditions [118]. Nevertheless, a changing climate is likely to have long-term effects on tree recruitment. Particularly in the transition of some forest ecotones, e.g., in the lee of the Andes, variations in temperature and/or precipitation can have serious effects on forest resilience and composition. This Schiaparelli Glacier study site was located within the ecoclimatic optimum for Nothofagus betuloides, within the cool and hyper-humid climate of Tierra del Fuego [142]. As a consequence, minor climate changes should only affect the production and survival rate of these seedlings in a subordinate way.

Conclusions
The use of trees and their corresponding germination dates represents a very reliable proxy for the exact dating of glacial deposits in the past centuries. The time span subsequent to the recession of Schiaparelli Glacier, moraine stabilization, soil formation, and individual trees that colonize the exposed material (ecesis time) was within 15 ± 3 years: very short.
Our study demonstrated the high potential and precision of determining local ecesis time by combining optical satellite imagery with dendrochronological methods. We recommend establishing linear models for the geolocation and age of the pioneer species and for the glacial retreat. In the case of constant glacial retreat, the trees will taper towards the glacial margin and both slopes of the linear models will have comparable values. These results were more reliable than usage of the shortest possible ecesis time that was generalized on the basis of the maximum age of a few individual trees.
The Nothofagus spp. at our study site was revealed to be a fast and homogenous growing pioneer species, indicating only minor uncertainties when cored at breast height (±1 year). However, under constrained growth conditions, growth rates might vary strongly, leading to higher uncertainties. Therefore, the trees, unless a dendro-climatological analysis is required, should be cored as close to the root collar as possible [16,60].
The largest uncertainty values resulted from missing the pith of the tree and the subsequent errors due to estimating the missing rings. The uncertainty value, introduced for estimation of a small number of missing rings (<10) was comparatively small, but increased with the number of missing rings (>10) from ±4 to ±10 years, respectively. For this reason, trees with an estimated number of more than 20 missing rings should not be included in the interpretation, which is in line with recent results for Himalayan tree species [60]. In our study, it was possible to perform moraine dating with an uncertainty value of less than 10 years by precisely determining the influencing and constraining variables resulting to formation of three dated moraines in ~1755 CE, ~1770 CE, and ~1880 CE, which is in accordance with the so far southernmost tree-ring-based moraine dates [16]. Furthermore, it was refuted that the Schiaparelli Glacier calved into the fjord, as was supposed based on a lithograph of the Beagle expedition, accompanied by C. Darwin and R. FitzRoy in 1836 CE [33]. The dating of sub-recent trunks imbedded into the glacial deposits revealed several Neoglacial advances within the last 2500 years that coincide with at least three periods of enhanced glacial activity at the Southern Patagonian Icefield and adjacent mountain glaciers, suggesting a fairly synchronous pattern. However, the number and spatial coverage of studies in southernmost South America is still very fragmented.
The presence of the Neoglacial advance (~550 BCE to 50 CE) was radiocarbon dated in glacial deposits within the Magellanes-Tierra del Fuego region for the first time. The following phases of culminated glacial advances were in phase with previously dated glacial events. Nevertheless, Schiaparelli Glacier was revealed to be advancing with an (small) offset during phases of a supraregional glacial recession. Additionally, the Schiaparelli Glacier even showed a contrary pattern: while a prominent glacier north of the CD was decreasing, our studied glacier was revealed to be advancing. This leads to the conclusion that further case studies are needed to assess the convergence and homogeneity of glacial activity across the region. Acknowledgments: The authors are grateful to Dr. Ricardo Jaña (Instituto Antarctica Chileno) and Inti Gonzalez (Centro Quarternario, CEQUA, Punta Arenas) for the logistical and administrative effort that made the field campaigns possible. Stephanie Weidemann kindly helped with sampling of buried wood for 14 C dating. We are grateful to all members of the field party to Mt. Sarmiento in April 2016. Special thanks to M.L. Fanti and W.K. Wauzi for their inspiration and stimulating discourses. The Landsat and ASTER data were kindly provided via USGS Earth Explorer. The authors would like to thank the German Aerospace Center and European Space Agency for providing the TanDEM-X data and Sentinel-2 data, respectively. The data analysis and interpolation were carried out using the open source statistical language R (www.r-project.org) with additional packages. We thank the reviewers for their careful reading of the manuscript and their constructive remarks.

Conflicts of Interest:
The authors declare no conflict of interest. The funding agency had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.