Quaternary glaciolacustrine deposits around a Triple Junction site: Paleolakes at the foot of the Northern Patagonian Ice field (Argentina and Chile)

The area involved by the triple junction between the South American, Nazca and Antarctic plates activity was affected by Quaternary glaciations. Before 12,800 yrs BP an extended ice field occupied the top of the Patagonian Andes, irradiating glaciers towards the east and the west dominantly. Towards the east, the ice melted in piedmont lakes; towards the west, fjords melted into the Pacific Ocean. The Upper-Pleistocene climate amelioration caused the recession of those glaciers. Some piedmont lakes reversed their Atlantic outflow towards to the Pacific Ocean. The glaciers retreat caused the fluvial reactivations along crustal former faults that were located below the ice. The Patagonian ice field became therefore split into present Northern and Southern fields. At the second largest lake of South America, the Buenos Aires-General Carrera Lake, the water level dropped from about 500 m over present mean sea level to 230 m. Several glaciolacustrine deposits from this area are indicating significant variations caused by climatic changes, volcanism and tectonics, differing in spatial and temporal magnitudes. The triple junction activity involved subduction of the Chile Ridge below the continental South American plate, volcanic activity and faulting. During the glacier melting the Baker River captured three eastern-moving glacial systems towards the southwest, towards the Pacific Ocean. This rapid event is thought to occur 12,800 yrs BP. The lowering of these glaciolacustrine systems should be also interpreted in terms of the tectonic activity in the region and considering other processes operating in the lakes and within the watersheds.

Significant changes in the watersheds of the lakes of Patagonia and Tierra del Fuego have been reported (Ton-That et al., 1999;Kaplan et al., 2004;McCulloch et al., 2000McCulloch et al., , 2005Del Valle et al., 2007;Solari et al., 2012: Horta et al., 2017Bendle et al., 2017;Thorndycraft et al., 2019;Bourgois et al., 2019). Concerning the Patagonian Ice Field, it became divided about 12,800 years ago (Turner et al., 2005). This caused a reversal in the direction of the

Introduction
During the Quaternary, climatic changes affected both hemispheres although they were not simultaneous and with the same magnitudes (Kaplan et al., 2004). The question remains whether these changes occurred gradually or abruptly (Broecker, 1998), or if the Holocene glacial advances occurred worldwide and simultaneously (Douglass et al., 2005). The behavior of glaciers is today matter of significant concern in relation to their effects on the present sea-level rise (Zemp et al., 2019). In the last years several reports have focused on the evolution of Quaternary lakes in relation to the Patagonian ice field (Ton-That et al., 1999;Kaplan et al., 2004;McCulloch et al., 2005;Hein et al., 2009Hein et al., , 2010Horta et al., 2011Horta et al., , 2017Van Daele et al., 2016;Bendle et al., 2017;Thorndycraft et al., 2019;Bourgois et al., 2019;Clague et al., 2020), presently divided into the Northern Ice Field (attached to the San Valentín Ice Field) and the Southern Ice drainage of two different lake complexes: the Buenos Aires-General Carrera, and the Pueyrredon-Posadas-Cochrane, with deactivations and activations of some river valleys (Isla et al., 2015;Glasser et al., 2016). Jökulhlaups events would have taken place during the evolution of some these lakes (Turner et al., 2005). Two glacier readvances were proposed to occur at this region ca. 8,200 and 6,200 yrs BP and assigned to a northward migration of the Southern Westerlies (Douglass et al., 2005). Many authors interpreted the geomorphological variations in particular areas of the Patagonian Andes in terms of climatic changes (Kaplan et al., 2004;McCulloch et al., 2005;Hein et al., 2009Hein et al., , 2010Horta et al., 2011Horta et al., , 2017Bendle et al., 2017;Thorndycraft et al., 2019;Bourgois et al., 2019). However, many changes could be also triggered by tectonic and volcanic activities associated to the dynamics of the triple junction site at the pacific coast in Patagonia Georgieva et al., 2016), that migrated significantly from south to north since the Miocene (ca. 15 Ma). This work comprises the description of glaciolacustrine deposits corresponding to Upper Pleistocene lakes that include the Buenos Aires-General Carrera, Bertrand, Plomo, Posadas, Pueyrredon and Cochrane lakes (Fig. 1). It has been assumed that both glacial valleys were occupied by a single and extensive lake named Chalenko Lake (Thorndycraft et al., 2019;Bourgois et al., 2019), although ice-walled lakes at different water levels could have existed (Martinod et al., 2016). The altitudes of these deposits were measured and analyzed in order to estimate the maximum altitude of these ancient paleolakes, and their lowering to present lakes. The particular evolution of these lakes in relation to the Triple Junction dynamics is compared to other piedmont lakes that had also significant changes in their watersheds at Patagonia and Tierra del Fuego.

Geological setting
The Buenos Aires-General Carrera Lake is the second largest lake in South America. It is occupying a glacial valley oriented W to E across the Patagonian Andes Cordillera where several Quaternary glaciations have occurred (Ton-That et al., 1999;Hein et al., 2009;Lagabrielle et al., 2010). The orientation of these lakes parallel to the dominant westerly winds causes seiches, mainly at the Argentine sector where the regional topography is not an obstacle to the wind action (Isla et al., 2010;Richter et al., 2010). At the southern shore of the General Carrera-Buenos Aires Lake, the Las Dunas River has built seven Gilberttype deltas (Bell, 2009). The oldest is at 650 m over present sea level (MSL) and assumed to be constructed ca. 13,000 years ago. The youngest non-operative delta (second delta) was probably constructed 8,500 years ago. The present delta is assumed to be stable at an altitude of 201 m for the last 6,200 years (Bell, 2009). The altitudes of these delta terraces above MSL along the General Carrera-Buenos Aires Lake shoreline were analyzed to discern its evolution (Bourgois et al., 2016) although their interpretations are not so simple (Martinod et al., 2016). An ice cover persisted in General Carrera-Buenos Aires Lake area until at least 16,934±116 cal yrs BP (Bendle et al., 2017). Moraine successions corresponding to the maximum glacier extensions (23,000 yrs BP; Fénix V Moraine) were established either for the General Carrera-Buenos Aires (Kaplan et al., 2004) and Pueyrredon-Posadas-Salitroso-Ghio lakes (Hein et al., 2010). Between 16,000 and 15,600 years BP, the Chacabuco valley, located between General Carrera-Buenos Aires and Pueyrredon-Posadas valleys, was reported free of ice and colonized by herbs, shrubs and evergreen forests (Villa-Martínez et al., 2012).
This lake complex is surrounded by basaltic rocks of Miocene-Pliocene ages (Meseta del Lago Buenos Aires Plateau; Panza and Nullo, 1994;Lagabrielle et al., 2010). The volcanic activity in this area remained until Pleistocene (Mazzoni and Rabassa, 2010;Ghiglione et al., 2016), with ash layers interfingered with glacial moraines (Fidalgo and Riggi, 1965;Ton-That et al., 1999;Kaplan et al., 2004), and with several intrusions of different basaltic dykes, as lopolites or laccolites  that could have caused local variations of the topographic altitudes. The recent alkaline volcanic activity (<14 Ma) is related to Chile Ridge subduction at the Taitao Peninsula (Guivel et al., 1999;Lagabrielle et al., 2000Lagabrielle et al., , 2004. This particular case is a Ridge-Trench-Trench (RTT) junction dominated by horizontal displacements and produces tilting of crustal blocks along a network of strike-slip faults on the upper plate . The northward migration of the Triple Junction along this active margin has caused a bimodal composition of volcanic rocks: mafic to intermediate lavas (48-56% SiO 2 ) and acidic lavas ranging from dacites to rhyolites (66-73% SiO 2 ; Lagabrielle et al., 2000). After 3 Ma a tectonic crisis led to the general uplift of the piedmont areas in the eastern side of the Andes and the rapid subsidence of narrow transverse (E-W) valleys that were later occupied by the piedmont lakes (General Carrera-Buenos Aires, and Cochrane-Pueyrredon-Posadas; Lagabrielle et al., 2010). The oblique subduction of the Chile Ridge caused the migration of the triple junction site from south to north up to the present location at latitude 46º30' S (Folguera and Ramos, 2002;Lagabrielle et al., 2010;Russo et al., 2010). The opening of an astenospheric window produced the dynamic uplift of the retroarc migrating from the piedmonts to the Extra Andean Patagonia Georgieva et al., 2016) causing a tectonic crisis for the last 3 Ma . This uplifting trend migrated northwards as the Chile Ridge subducted below South American Plate (Ghiglione et al., 2016). The relative motion between South American and Nazca plates was calculated about 8.5 cm/yr, and 2 cm/yr between South American and the Antarctic plates (Boutonnet et al., 2010). The northward movement is caused by the asymmetric growing of this active ridge, even the architecture of the continental margin is different to the north and south of the present triple junction. To the north, the continental margin is characterized by a small, Plio-Quaternary, accretionary prism with much of the older material removed by tectonic erosion .
To the south, the accretionary prism widens abruptly as the deformation front migrated seawards and a broader zone of trench turbidites is involved in the deformation (Ranero et al., 2006). Considering the magnitudes of these tectonic effects in the region, the rol of the glacio-isostatic rebounds have been overestimated.

Methods
Glaciolacustrine deposits identified from field mapping were located on a GIS database (Arc View). Heights of the tops of these deposits were established with a handled GPS in order to recognize the deposits corresponding to the highest levels of the former lakes (Table 1). Radiocarbon dates are reported in the original units, either calibrated or without calibration ( Table 2). The thicknesses of these sedimentary sequences were measured in each location. Seven samples from these glaciolacustrine deposits were analyzed for grain sizes distribution, combining sieving and the pipette methods. Sieving was performed every 0.5 phi units, from coarse to very fine sand. Results are reported in a triangle plot with the sand, silt and clay percentages. Altitudes were referred to mean sea level (MSL) avoiding estimations to the present level of the lakes, as these large lakes, oriented W-E, are subject to differences of several meters due to barometric-triggered seiches and wind effects (Richter et al., 2010).

Glaciolacustrine sequences
Several glaciolacustrine sequences were identified and analyzed in relation to their altitudes, varve spacing and geographic location. Four were described in the Argentine territory while other four sequences are reported from the Chilean territory (Fig. 2, Table 1). Argentine sequences are described from E to W; the Chilean from N to S. Grain sizes of the sampled varves vary significantly from 0 to 35% of sand. These variations were interpreted according to the lake dimensions, dynamics and distance from the glacial sources.
The samples were taken both in the Argentine (Fenix Chico River, Route 41 and Jeinimeni River) and in Chilean (Fachinal, Chacabuco Bridge) territories. Their compositions are dominated by clay and silt fractions.

Fénix River terrace (top about 460 m above MSL)
On the way from Perito Moreno to Los Antiguos (Fig. 3), the route 43 runs close to the headlands of the Fenix River. This terrace is about 460 m over MSL and corresponds to the former headlands of the Deseado River watershed. Presently, it is the continental divide between the Atlantic and Pacific watersheds. The outcrop is very small; a till is  outcropping at the base of the profile (Fig. 3A-B). On top of this sequence convolute bedding denotes physical perturbations (gravity-dominated phenomena are assumed) of the ancient lake during deposition.

Lake Buenos Aires terrace (top about 354 m above MSL)
This is the original sequence described by Caldenius (1932) where he counted and measured 530 varves. The deposit is 22.6 m thick with maximum varve thicknesses of 0.57 m. The sample, collected from the middle of the sequence, is composed of silt and clay with small amounts of sand (Fig. 4); muscovite grains are dominant. The sequence was revisited and described in detail by Bendle et al. (2017). An ash layer intercalated in the sequence is assumed to belong to the Hudson Volcano (170 km to the NW of this site) and dated in 17,378±118 cal yrs BP (Bendle et al., 2017). However, this tephra layer was not accepted as belonging to the Hudson Volcano (Bourgois et al., 2019). Five different sites from the Eastern Buenos Aires Lake sum 994±36 varve-years (Bendle et al., 2017;Fig. 3C, D). Based on a 10-year moving average of the varve thickness, an average ice retreat of 5-10 m/yr occurred until 17,322±115 cal yrs BP that accelerated to 15-18 m/yr to the end of the sequence (Bendle et al., 2017). Some of the assumptions to estimate these retreats were also questioned (Bourgois et al., 2019).

Route 41 (South from Los Antiguos; top
terrace about 514 m above MSL) About 13 km to the south of Los Antiguos village, route 41 is on top of a former lacustrine terrace. Further south, perched glacial features have been described in spatial relation to Monte Zeballos basalts . This glaciolacustrine terrace are composed of gravel (coarser than 3 mm) alternating with sand and suggesting the coastal facies of the former lake (Fig. 3E, F). In a 1.5 m thick deposit 16 varves were counted (about 0.10 m each). In the sample, from the middle of the sequence, the silt fraction is dominant (56%) with less proportions of clay (29%) and sand (15%; Fig. 4).

Jeinimeni River (top terrace about 262 m above MSL)
Between the Argentine Customs control and the Argentina-Chile boundary, the route runs close to a 20 m tall cliff along the Jeinimeni River valley. This valley was affected by glaciers and landslides flowing from the basalts of the Meseta del Lago Buenos Aires . Several outcrops of a former lake are accessible from the route FIG. 4. Grain-size composition of the samples collected from the glaciolacustrine facies. Red squares are related to the Buenos Aires Lake, blue triangles to the General Carrera Lake, and black circles to the Chacabuco Bridge. (Fig. 3G, H). The sample collected was composed of silt (62%), clay (30%) and only 7% of sand.

Fachinal terrace (top about 304 m above MSL)
The Fachinal terrace is attached to the southern coast of the General Carrera Lake, Chile (Fig. 5A). Along the route at the southern shore of the lake, there are several outcrops. At the top of one of these sequences, there is a small quarry with varves of 1-2 cm thick (Fig. 5B). The glaciolacustrine deposit is close to the mountain and composed of coarser material: 73% silt, 16% sand and 9% of clay. Thirty km to the SW, the Las Dunas River delta has been subject to a particular analysis of its poorly sorted gravel spanning from well-rounded to subangular grain sizes, from cobbles (up to 200 mm) to coarse sand (Bell, 2009). Two moraines are related to the modern delta; they were dated on 8.5±0.7 kyrs BP (outer moraine) and 6.2±0.8 kyrs BP (inner moraine; Douglass et al., 2005). At the southern coast of the General Carrera Lake two ancient shorelines were analyzed in terms of their altitude variations, indicating a tilting rate of 0.6 m/km for the higher shoreline, and 0.45 m/km at the lower coast (Turner et al., 2005).

Bertrand Lake (about 236 m above MSL)
Along the route from Bertrand Lake towards the south, there are saltuary outcrops with glaciolacustrine facies. Moraines close to Bertand Lake have been dated between 14,100 and 15,900 yrs BP (Thorndycraft et al., 2019). Varves are about 1 cm thick and composed of fine sediments (Fig. 5C, D).

Borde Baker (about 230 m above MSL)
At the western border of the Baker River, a huge landslide discovered a glaciolacustrine sequence close to the touristic resort Borde Baker. The landslide area, difficult to access, is about 200 m long and 30 m height. A glaciofluvial deposit is on top of the cliff (Fig. 5E, F).

Chacabuco Bridge (about172 m above MSL)
On the route from El Furioso to Cochrane, several glaciolacustrine sequences are outcropping. One of the characteristic of this area is the right angles of the Baker River in its flow towards the south, indicating a tectonic control. The section was described in detail by Turner et al. (2005) and recently reinterpreted by Thorndycraft et al. (2019). An ancient Chacabuco Paleolake was reported as older than the General Carrera-Buenos Aires and Pueyrredon-Posadas systems, and at an altitude of 620 m above MSL (Glasser et al., 2016). At the place called Puente Chacabuco (Chacabuco Bridge), the outcrop is more than 400 m long and 8 m height. The thicknesses of the varves are of approximately 1 m each (Fig.  5G, H). Three layers were sampled at this site; they are composed of 32-47% of silt; 32-42% of clay and 9-34% of sand (Fig. 4). The samples from this outcrop are the coarser sampled with thicker varves.

Altitude distribution
At the Argentine side, the maximum altitude of the glaciolacustrine sequences diminishes towards the west. At the Chilean side, the altitudes of the former lake levels diminish northwards (Fig. 2). At the Argentine territory the altitudes are higher as the distances to the present lakes increase. The highest levels are at about 500 m over MSL suggesting the maximum extension of the largest paleolake.

Landscape evolution
The decrease in the water level of this great lake caused the formation of several Gilbert-type deltas in relation to each level (Bell, 2009). Beach plains also formed attached to the most modern delta (Bell, 2009). The tilting of the shoreline deposits of the General Carrera Lake was analyzed exclusively in relation to the distance from the glacier center (Turner et al., 2005).
According to the ash layer recognized close to the top of the sequences at lake Buenos Aires terrace, dated in 17,378±118 cal yrs BP, and assumed to belong to Hudson volcano (known as Ho), the receding rates of the glaciers were moderate (5-10 m/yr) before the eruption and accelerated to 15-18 m/yr after that age (Bendle et al., 2017). These watersheds were also subject to intense faulting derived from the interactions of the subduction of the active ridge since 3 Ma to Present (Lagabrielle et al., , 2004.

Discussion
Variations in lake levels can be assigned to climatic changes, but also to variations in the watersheds triggered by tectonic movements, fluvial captures or volcanic activity (Isla et al., 2015;Martinod et al., 2016). Palinological studies in some of the piedmont lakes in the Patagonian Andes indicate significant variations of the westerly winds that could led to dessication intervals and increments of their strength between 11,200 and 4,500 years BP suggesting a broadening of the wind belt (Van Daele et al., 2016). The more complete varve sequence considered in this study (Eastern Lake Buenos Aires terrace), however, is only reporting approximately 994 years. For the Eastern Buenos Aires Lake, several lava deposits of Upper Pleistocene age were reported leading to the deactivation of the Deseado River (Ton-That et al., 1999;Kaplan et al., 2004). The retreat of the ice lobes about 10 kyrs BP was proposed for the activation of the Baker River and the reversal of the drainage of the lacustrine system towards the Pacific Ocean (Isla et al., 2015). During the interval of this drainage reversal, several paleolakes (Cisnes-Ñirehuao, Deseado, Tranquilo, Colonia, Cochrane-Pueyrredon-Chacabuco) were defined and mapped (Davies et al., 2020). A temporal reversal via the Bayo River was proposed (Martinod et al., 2016) operating during the interval between 15,000 and 12,600 years BP (Bourgois et al., 2019). Catastrophic GLOF events have also occurred at this watershed (Ulloa et al., 2018;Thorndycraft et al., 2019;Bañales-Seguel et al., 2020). The reactivation of tectonic faults induced by the subduction of the Chile Ridge, mostly related to the Cachet and Exploradores fault zones, can explain these recent geomorphological changes (Agurto-Detzel et al., 2014;Georgieva et al., 2016). An isostatic rebound was proposed to explain the altitudes of the Las Dunas River terraces (Bourgois et al., 2016;Martinod et al., 2016), disregarding the tectonic tilting (Turner et al., 2005). These tectonic processes were reported for the Meseta del Lago Buenos Aires where an inversion of a paleo-piedmont occurred about 3 Ma . The basalts of Meseta del Lago Buenos Aires and Meseta de Chile Chico are genetically related to a subducted segment of the Chile Ridge, enclosed between the Tres Montes and Esmeralda fracture zones (Lagabrielle FIG. 6. Sketch map of the evolution of the Northern Patagonian Ice Field and the associated glaciolacustrine lakes. A. The extension of former lakes during the Pleistocene was located at the 500 m over MSL. B. Present conditions for the approximately 230 m height of the Buenos Aires-General Carrera Lake and the Pueyrredon-Posadas lacustrine complex (modified after Glasser et al., 2004Glasser et al., , 2005Bendle et al., 2017Bendle et al., ). et al., 2004. Tilting effects were also reported for the former shorelines of ice-dammed lakes at the Magellan Strait region . The isostatic rebound of the Northern Patagonian Ice Field (NPIF) should be minimized considering that it was neglected in Tierra del Fuego, where glaciations were assumed to be more extended than at the NPIF (Mörner, 1991). The basaltic extrusions and intrusions occurred during the last 3 Ma conditioned the morphology where the glaciers were limited to flow . Lava deposits of Upper Pleistocene age were reported to the east of Buenos Aires Lake (Ton-That et al., 1999). These deposits caused that some Argentine rivers (Pinturas, Deseado) diminished significantly their watersheds. As consequence the Baker River incorporated stepwise several lakes in its watershed about 10 kyrs, and led to the progradation of its delta system close to Caleta Tortel (Isla et al., 2015). An evolution model of these steps in the retreat of the ice fronts and lake level drops has been proposed based on digital elevation models (Fig. 7). A fault oriented N-S was mapped at the Eastern border of the present NPIF (Cachet Fault). This block is assumed to move towards the north (Georgieva et al., 2016). The final stages of the drainage of the General Carrera- FIG. 7. Evolution model of the ice-front retreatments of the original Patagonian Ice Field (PIF), lake-level variations, blocking of the flow towards the East (Fenix, Pinturas, Deseado rivers) and the glaciolacustrine-deposits at the Buenos Aires Lake (modified from Ton-That et al., 1999, andGeorgieva et al., 2016). NPIF: Northern Patagonian Ice Field. Buenos Aires Lake system has been reported with the catastrophic flooding of the Baker River about 10,000 yrs BP (Thorndycraft et al., 2019;Davies et al., 2020). Other glacial lake outburst floods (GLOFs) have occurred recently along the Baker (Ulloa et al., 2018) and Exploradores rivers (Bañales-Seguel et al., 2020). Recent GLOFs are of special concern in relation to the historic mountain-glacier contributions to the sea level (Zemp et al., 2019) Several glaciolacustrine sequences have been reported along the Southern Andes of Patagonia and Tierra del Fuego from 40 to 54° S, although the altitudes and ages of these paleolakes are not similar. Some of these sequences are suggesting reversals in the drainage patterns: a. The Aluminé and Moquehue lakes (38°55' S; Fig. 1) were though to drain towards the northeast to the Pacific Ocean by the watershed of the Biobio River. After the eruption of the Sollipulli Volcano about 2,900±60 yrs BP, the Alumine Lake is draining towards the SE (Isla et al., 2019). The Aluminé River is flowing to the Collón Cura, Limay and Negro rivers, towards the Atlantic Ocean. b. The great Elpalafquen paleolake (41° S; Fig. 1) comprises several modern piedmont lakes (Nahuel Huapi, Gutiérrez, Mascardi, Fonck and Steffen, and the small shallow lakes Cari Lafquen Grande and Cari Lafquen Chica). This great lake originally drained towards the Atlantic Ocean via the Limay River. Afterwards, it began to flow to the Pacific Ocean via the Manso River about 13,200 yrs BP (Del Valle et al., 2007;Davies et al., 2020). At the area between Pilcaniyeu and Maquinchao, east of the present piedmont lakes, two dacitic tephras indicate that this ancient lake remained between 13,800 and 11,200 yrs BP (Tatur et al., 2002). c. Another paleolake, the Caldenius (47°50' S; Fig. 1) comprised the Belgrano, Azara, Mogote, Nansen, Volcán and Burmeister lakes, of former Atlantic destination that opened towards the Pacific Ocean about 10,000 yrs BP (Del Valle et al., 2007). Two lakes, related to present Belgrano and Burmeister lakes, of Upper Pleistocene age (about 31,200 yrs BP) were positioned 880 m above MSL. Between 11,371 and 10,819 yrs BP, a lake close to Northern Belgrano Lake, reached levels of 900-920 m above MSL. About 6,900 yrs BP the level dropped to split into two-lake settings separated by moraines (Horta et al., 2017).
d. The Great Tehuelche paleolake (51°50' S; Fig. 1) occupied the area of Torres del Paine. It comprised several of present lakes: Nordenskjold, Pehoe, Sarmiento, del Toro and Porteño. Between 24,291 and 17,500 yrs BP the paleolake level was between 260 and 240 m above MSL (Solari et al., 2012). Between 14,800 and 12,600 there was a glacial readvance while the lake was at an altitude of 165-150 m above MSL (Solari et al., 2012). During that interval several temporal lakes were recognized (Blanca, Skyring, Otway; Davies et al., 2020). 7,113 years ago the area of the former Great Tehuelche Lake began to drain towards the Pacific Ocean, via the Seno Ultima Esperanza (Solari et al., 2012). Two glaciolacustrine terraces were discriminated along the Almirante Montt Gulf: the oldest (150-125 m height) attached to the Arauco Morainic Complex (Arauco paleolake), and the modernmost (30 m altitude), attached to the present oceanic coast (Sagredo et al., 2011). Between 1,200 and 200 years ago the del Toro Lake was at altitudes about 80 m above present MSL (Solari et al., 2010). e. Another extended piedmont lake has occupied the central portion of the Magellan Strait and Bahía Inútil (38°55' S) during the end of the Pleistocene (MacCulloch et al., 2005). This paleolake Magallanes-Bahía Inútil, occupying both glacial valleys, persisted about 15,000 yrs BP (Davies et al., 2020). Two glacial advances were discriminated in relation to the Mount Reclus eruption event. f. In the Grande Island of Tierra del Fuego, changes in the glaciolacustrine watersheds have also taken place. The ancient paleolake called Fueguino Lake (54°04' S; Fig. 1), comprising the Yehuin, Chepelnut, Yakush and Penny lakes, was draining towards the Atlantic Ocean during the Upper Pleistocene (Del Valle et al., 2007). The Fagnano Lake captured this watershed after 7,800 yrs BP, and therefore today it is draining towards the Pacific Ocean, via Seno Almirantazgo Channel (Del Valle et al., 2007). The comparisons between the altitudes and ages of these Upper-Pleistocene glaciolacustrine sequences should be analyzed considering also the differential tectonic trends and the distance to volcanic features (Ton-That et al., 1999;Lagabrielle et al., 2010;Georgieva et al., 2016). Glacio-isostatic rebounds should be disregarded as a major effect on topographic variations in comparison to the tectonic and volcanic-triggered effects (lopolites, laccolites). However, these comparisons should consider other issues as watershed changes, precipitation variations, local climate and wind-induced lake dynamics.
6. Conclusions 1. The great Buenos Aires-General Carrera paleolake reversed its flow at the same time that its level dropped from approximately 500 to 230 m over mean sea level about 12, 800 yrs BP. This altitudinal difference should be interpreted considering climatic changes, local tectonic effects (faults, including tilting) probably induced by the subduction of the Chile Ridge. 2. The geomorphological changes that affected this former glacial valley transformed it into a great lacustrine basin and led to the capture of another glacial valley discharge occupied today by the Cochrane, Pueyrredon and Posadas lakes. 3. As the drainage of the glacial valleys reversed their discharges direction towards the Pacific Ocean (via the Baker River), the rivers flowing towards the Atlantic Ocean (Deseado, Pinturas) had experienced significant diminution on their discharges since 12,800 yrs BP. 4. Much of these geomorphologic changes should be also assigned to the reactivation of faults caused by the subduction of the active Chile ridge at the Triple Junction region, and the volcanic activity, occurred since 3 Ma to Present. Therefore, it is not recommended to extrapolate the timing of the geomorphological variations of a region assigned to climatic changes without considering other factors (tectonic, volcanic activity, lake dynamics, watershed changes).