Remotely sensed time series of rapid terrace formation: Laguna del Viedma valley (Patagonia)

. The Patagonian Andes were affected by a range of geophysical drivers of landscape incision during the Last Glacial Interglacial Transition and Early Holocene. Deciphering drivers of river system response during this period is complex, and magnitudes and timescales of landscape change are poorly constrained. Herein, a remotely sensed time series of modern terrace formation is investigated from the Laguna del Viedma valley as a modern analogue of Late Quaternary landscape evolution in Patagonia. The aim of the research was to constrain the timing of terrace formation following lake level fall of the Laguna del Viedma over a 35 year period from 1985-2019. The objectives were to: 1) use satellite imagery from 1985-2019 to document glacier and lake changes in the study area; 2) map landforms of the Laguna del Viedma valley; and 3) analyse terrace elevations. In total seven terrace surfaces were distinguished, with the oldest four pre-dating the ALOS PALSAR DEM (February 2000) used. Landform evidence shows the highest, and vegetated, T1 terrace surface (+40-75 m) grades to the highest lake level and was likely the elevation of the valley floor during Holocene neoglacials. Viedma glacier recession then led to a phase of lake regressions/transgressions with an overall trend of lake level fall. The DEM shows ~20 m incision from the 1985 floodplain level (T3) to the T4 level floodplain by 2000. This constrains a minimum rate of incision of 1.33 m/yr, however, the satellite time series demonstrates rapid T3 terrace formation, with evidence for mass movements contributing to lateral terrace erosion by 1986. The implications of the data are discussed within the context of the Late Quaternary palaeohydrology of Patagonia where lake level falls of 10s to 100s of metres occurred within many large river systems of the Patagonian Andes from 42-52⁰ S. The data herein demonstrate that base level falls from sudden lake drainage events were likely a major driver of rapid landscape change in Patagonia during deglaciation.


Introduction
Base level is a key concept in fluvial geomorphology (e.g., Davis, 1902;Schumm, 1993). As defined by Schumm (1993) base level is an imaginary level surface that limits fluvial incision. In inland settings, Faulkner et al. (2016) recognised three local to regional-scale types of base level: enclosed depressions, lake levels, or rivers downstream of tributary valleys. Base level, therefore, can be controlled by rise/fall in eustatic sea level (Schumm, 1993;Harvey, 2002), rise/fall of lake level altitude (Harvey, 2002;Bowman et al., 2010) or incision/aggradation in downstream river valleys (Faulkner et al., 2016).
Changes to lake levels and downstream riverbed elevations can be particularly dynamic during deglaciation and paraglacial response of Late Quaternary terrestrial landscapes. During deglaciation the magnitude of ice dammed lake falls in Patagonia can be in the order of 10s to 100s of metres (e.g., Bell, 2008;García et al., 2014;Solari et al., 2012;Turner et al., 2005;Thorndycraft et al., 2019). Rates of base level change can occur over abrupt (10 1 -10 2 yr) timescales, for example in response to changes in sediment supply (Faulkner et al., 2016), while glacial lake outburst flood (GLOF) events can cause rapid (<1 yr), geologically instantaneous, lake level falls (Jacquet et al., 2017). Furthermore, high magnitude GLOFs can cause major erosion and channel incision along the flood pathway (Baker, 2009;O'Connor et al., 2020;Benito and Thorndycraft, 2020), lowering base level at confluences with tributaries otherwise unaffected by the flood event. The interpretation of rates of Late Quaternary landscape change can, therefore, be challenging in formerly glaciated valleys, not least given constraints on dating precision. This is especially so where there are multiple drivers, such as the Patagonian Andes where isostatic uplift (Dietrich et al., 2010), tectonics and dynamic topography (Guillaume et al., 2013;Dávila et al., 2019;Ávila and Dávila, 2019) can lead to uplift or subsidence.
In this paper, the rate of terrace formation in response to base level falls/transgressions resulting from ice dammed lake drainage/refilling is investigated using satellite remote sensing as a modern analogue for Late Quaternary landscape evolution in Patagonia. The study area is the lower Viedma glacier, an outlet glacier of the Southern Patagonia Icefield, and the ice-dammed Laguna del Viedma valley (Fig. 1). The aim of the paper is to test the hypothesis that incision up to 10s of metres could occur over rapid (sub-decadal) timescales in response to lake level fall. The objectives are to: 1) use satellite imagery from the period 1985-2019 to document landform, glacier and lake changes in the study area; 2) use remotely sensed imagery to map landforms of the lower Viedma glacier and Laguna del Viedma valley; and 3) analyse terrace elevations from the ALOS PALSAR DEM, which uses SRTM datasets collected in February 2000. The implications of the Laguna del Viedma data are discussed with respect to the interpretation of Late Quaternary landforms resulting from deglaciation of the former Patagonian Ice Sheet.

Study area
The Viedma Glacier (Fig. 1) is one of the major outlet glaciers draining the eastern ice-shed of the Southern Patagonia Icefield. The glacier terminates in a calving ice-margin at the western end of Lago Viedma at ~250 m a.s.l. At the southern ice margin of the lower Viedma, the ice blocks a SW-NE aligned drainage basin (Fig. 2) forming the ice-dammed Laguna del Viedma. During the local Late Glacial Maximum the Viedma ice limits extended 10s of km east of its present position (Ponce et al., 2019) so the catchment of the present day Laguna del Viedma would have been beneath thick ice at this time.
The emergence of ice-free conditions in the Laguna del Viedma valley is undated but the presence of trimlines, ice-scoured bedrock, tributary valley moraine complexes (Fig. 2), and dense valley forest indicate the likely extent of Holocene neoglacial advances based on moraine records elsewhere. For example, to the north of the Viedma glacier, the Torres glacier, which shows a similar valley glacier landsytem, formed moraines at ~9.7, 6.9, 6.1, 4.5, and 0.53 ka (Reynhout et al., 2019). The presence of an arête on the eastern catchment watershed shows the modern glaciers are now mainly disconnected from the icefield (Fig. 1 The drainage area of the Laguna del Viedma valley is 149 km 2 . The contemporary planform of the main valley floor is predominantly braided, and the three main rivers that feed the Laguna del Viedma are fed from glacial lakes (Fig. 2).

Geomorphological mapping
Landform mapping was carried out at 1:5,000 scale using the high resolution ESRI base layer imagery (cf. Bendle et al., 2017a) in ArcMap 10.3 and ArcGIS Pro using UTM projection (Zone 18S). A selection of images are shown in figure 3. Recourse was made to the 3D imagery available in ArcGIS Earth to cross-check breaks of slope. Landform identification focused on glacial, fluvial and slope landforms where lack of forest cover allowed, and followed standard nomenclature used in geomorphological mapping in Patagonia (e.g., Bendle et al., 2017a, Leger et al., 2020. Standard GIS mapping protocols were used (Smith et al., 2014;Chandler et al., 2018).

Satellite change detection
To map the position of the Viedma Glacier ice margin and spatial extent of Laguna del Viedma, remote sensing imagery was used, with the timeframe spanning 1985 to 2019. Ice margins and lake extent were mapped manually, at 1:25,000 scale using ArcMap v.10.3, following standard GIS mapping protocols (e.g., Smith et al., 2014;Chandler et al., 2018). In total, 26 images were sourced from Landsat and Sentinel missions (Table 1), with a greater frequency of images in the latter part of the record. Eight of the 25 images had partial cloud cover (Table 1). Sometimes cloud cover precluded mapping the Viedma glacier snout, but in all the selected images the spatial extent of Laguna del Viedma could be inferred. All the LANDSAT imagery was collected by the LANDSAT 5 TIR scanner, except for the January 1985 image, which used the MSS scanner.

Terrace elevations
Landform elevations were determined using the ALOS PALSAR Radiometric Terrain Correction dataset, a digital elevation model (DEM) product with a spatial resolution of 12.5 m, considered more robust than the 30 m ASTER GDEM because more pixels covered the terrace remnants. A study evaluating vertical accuracy of open source DEMs for gully erosion concluded the root mean square error (RMSE) for the ALOS PALSAR data was 1.7 m, compared to 1.9 m for the ASTER GDEM (Khal et al., 2020). The source DEM for the ALOS PALSAR dataset was the 30 m resolution SRTMGL1 DEM, collected in February 2000. To determine the relative heights of mapped terrace levels relative to the valley floor, elevations were extracted from mapped terrace surfaces and the median value for each was obtained. This was achieved by firstly mapping the terrace landform, then clipping the

Glacial and lake geomorphology of the study area
Prominent lateral moraines and a terminal moraine were mapped beyond the 1985 ice extent of the Viedma Glacier (Fig. 4). The terminal moraine, with arcuate ridges and kettle holes, forms an isthmus extending out into lago Viedma. This marks a position 2.5 km down-valley from the 1985 snout position and likely dates to the last Holocene readvance (neoglacial) dated to ~0.53 ka in the Fitz Roy valley to the north (Reynhout et al., 2019). At the time of moraine formation, the right margin of the Viedma glacier would have extended further into the Laguna del Viedma valley than at present, as indicated by the eroded bedrock beneath the treeline (Fig. 1), and a moraine to the west of the current ice margin. The satellite time series shows that this moraine was overtopped by lake waters during high stands of Laguna del Viedma and multiple, closely spaced, shorelines (formed by lake level falls) are imprinted upon it (Fig. 3D). This suggests the landform could have been reworked, but the orientation of the ridge and presence of large boulders indicate a moraine rather than fluvioglacial origin.
Up-valley from the Laguna del Viedma, tributary valleys contain prominent moraines (Fig. 2), with the youngest being unvegetated, while older moraine sets are covered in trees. This suggests multiple Holocene periods of glacier readvance (neoglacials). The geochronology of the Torres glacier moraines, preserved in the Fitz Roy valley to the north of the Viedma glacier, indicates Holocene glacier readvances at 9.7±0.4, 6.9±0.2, 6.1±0.3, 4.5±0.2 and 0.53±0.06 ka (Reynhout et al., 2019). is unvegetated. Moraines of the Upsala glacier (Lago Argentino basin) to the south of Laguna del Viedma record other neoglacials at 2.5-2.0 ka and 1.5-1.1 ka . Kaplan et al. (2016), based on Holocene glacier reconstructions from the Lago Argentino basin, conclude that maximum Holocene glacier advances occurred 6.1 to 4.5 ka. The undated moraines (Fig. 2) are likely better correlated to those of the Torres rather than Upsala glacier because the narrower valley glacier landsystem of the Torres glacier is a more appropriate analogue than the large outlet Upsala glacier. Therefore, the youngest of the forested moraines in the Laguna del Viedma valley is correlated herein to the 4.5 ka moraine (Kaplan et al., 2016;Reynhout et al., 2019), while it is highly likely the unvegetated moraines correlate to the final (0.53 ka) neoglacial (Reynhout et al., 2019). Interpretation of the study area mapping and regional Holocene glacier chronology therefore indicates that a 7.6 km length of the Laguna del Viedma valley was likely ice-free between up-valley ice limits and the Viedma glacier moraine since at least 4.5-6.0 ka, allowing formation of Laguna del Viedma. The highest shorelines visible on satellite imagery indicate a maximum lake elevation of ~530 m a.s.l.

Change detection
Glacier outlines for the lower Viedma glacier are shown in figure 4 for five snapshots from the satellite time series . The general pattern of change is for glacier snout recession from its 1985 extent to 2019, but there is an acceleration in rate of recession over the last decade. The glacier receded 2.7 km from the inferred 0.53 ka to 1860 terminal moraine position (Davies and Glasser, 2012;Reynhout et al., 2019)  The spatial extent of Laguna del Viedma is shown in figue 5 for selected years. In general, there is a decrease in the area of the lake from 1985 through to 2019 but there is evidence for cycles of emptying and refilling (lake regression/ transgression), as is usual for ice-dammed lakes in the Patagonian Andes (Dussaillant et al., 2010;Jacquet et al., 2017). Evident from the lake outlines is the effect of head-cutting erosion following lake drainage events (Jacquet et al., 2017), for example, an arm of the 1998 lake extends further up-valley than the 1985 shoreline position.

Terraces of the Laguna del Viedma valley
The geomorphological map of the Laguna del Viedma terraces is presented in figure 6. The main river occupies a braidplain inset beneath a staircase of at least seven, distinguishable terrace surfaces. The highest T1 terrace is well-vegetated with trees and shrubs (Fig. 3A). The largest preserved area of T1 occurs in the west of the study area (Fig. 6) and occupies a position in between two rivers entering the lower Laguna del Viedma valley. Other preserved T1 surfaces occur on the southern valley side, where they are dissected by tributary streams entering the main valley. Where vegetation cover is sparser a braided channel pattern is visible on the terrace surface (Fig. 3A). Horseshoe-shaped scarp headwalls occur along the edges of the T1 terrace surfaces, indicating head-cutting erosion by mass movement events. The most easterly terrace remnant has a steeply dipping scarp slope, facing eastwards down-valley, characteristic of a delta.
The spatial extent of the unvegetated T2 terrace surface is less extensive than T1 and is represented by three preserved remnants all on the north side of the valley. The furthest eastwards terrace occupies a broadly similar position to the furthest eastwards T1 terrace remnant. Horseshoe-shaped scarp headwalls also occur along the edges of the T2 terrace surfaces (e.g., Fig. 3A). The T3 terrace extends further eastwards than T1 and T2, with two large terraces preserved on both the northern and southern sides of the valley (Fig. 3C), the easterly scarp slope of the latter is imprinted with shorelines. At the confluence with a southerly tributary a delta surface at the T3 level also preserves eroded shorelines (Fig. 3B). These terraces mark the position of the 1985 lake level (Fig. 5) at the beginning of the satellite time series. Smaller remnants of both T3 and T4 terrace levels occur up-valley, for example in the centre of the valley upstream from the confluence located at 73°9' W (Fig. 6). The T5-T7 terraces are all preserved at the eastern end of the study area (Fig. 6)  floor, with the greatest height difference downvalley to the east as illustrated by cross-sections XS-1 to XS-3 (Fig. 7) and the terrace long-profile (Fig. 8). The average dimensionless slope of the terrace is calculated as 0.02 but appears to decrease down-valley where it grades to ~520-530 m a.s.l., coinciding with the highest elevation shorelines in the valley. Terrace T2 is located ~10-20 m below T1, and the T3 surface, which grades to the 1985 lake level, is ~30 m below T2. There is then a fall of 10 m (upstream) to ~30 m (downstream) from terrace T3 to surface T4, the key incision event captured by the year 2000 elevation data and satellite record. The formation of terraces T5-T7 (Fig. 6) post-date the ALOS PALSAR DEM and so height differences were not calculated.

Rates of landscape change in the Laguna del Viedma valley
The geomorphology of the Laguna del Viedma valley demonstrates incision of unconsolidated sediments on the valley floor with at least seven terrace surfaces identified through remote sensing, and a peak total incision of ~70 m in a reach upstream of the 1985 delta position. Establishing the timeframe of this total incision, prior to the satellite time series is uncertain but given the geomorphological evidence -the position of the ice dammed lake and terraces below mature forest ( Fig. 1) developed on post-glacial soils, the T1 floodplain would clearly have been active during Holocene neoglacials. One interpretation of the geomorphology is that the vegetated T1 terrace coincides with the vegetated moraines up-valley, while the unvegetated T2 terrace is associated with the youngest unvegetated moraines. As discussed in Section 4.1, based on the Torres glacier chronology the age of the youngest vegetated moraines are correlated to the 4.5 ka neoglacial, while it is highly likely the unvegetated moraines date to the most recent 0.53 ka neoglacial (Reynhout et al., 2019). Based on an age of 0.53 ka for the unvegetated moraines and terrace T2, an average incision rate of ~0.08-0.12 m/yr from T2 to the valley floor can be constrained.
The satellite imagery time series provides robust data on more recent rates of landscape change. In particular, the fall of ~20 m from T3 (an active floodplain in 1985) to T4 prior to the source SRTM DEM data collection (2000) provides a minimum rate of incision of 1.33 m/yr. However, the satellite imagery ( Fig. 9) indicates much of the incision occurred quickly. Fig. 9A shows the active floodplain at T3 feeding into Laguna del Viedma in January 1985. By September 1986 significant scarp slopes had already formed on the T3 terrace (Fig. 9B) and slumping of the scarp slope can be seen in the October 1986 image (Fig. 9C). Over the rest of the time series we can observe lateral erosion, in particular to the north of the river, likely due to lateral fluvial erosion coupled with mass movement processes. In terms of vertical incision, therefore, we can infer from the time series rapid incision in the order of metres to tens of metres immediately following lake level fall from its 1985 level. For comparison, Tobal et al. (2021) calculate an average incision rate for the río Deseado (46-47ᴼ S) outwash terraces of ~0.135 m/ka since 1.2 Ma in response to tectonic and climate forcing. The impact of rapid base level falls on rates of incision has been demonstrated elsewhere through field, laboratory flume and remote sensing studies. The latter was employed by Jacquet et al. (2017) to quantify bed elevation changes following drainage of lago Cachet II, an ice dammed lake formed in a tributary valley of the Colonia glacier (Northern Patagonia Icefield). The study demonstrated that during a phase of 21 episodic GLOF drainage events since April 2008 there was a net loss of ~25x10 6 m 3 of sediment, with >40 m of incision into the unconsolidated valley floor sediments.
Dam decommissioning projects also provide opportunities for investigating geomorphic response to rapid base level fall. Wilcox et al. (2014) monitored the removal of Condit Dam on the White Salmon River in northwest USA. The 38 m high dam, which impounded 1.6 million m 3 of water, was breached by blasting a 5 m hole in its base. Following the breach, the reservoir drained in 90 minutes and ~10% of its 1.8 million m 3 of sediment was eroded during this initial phase of the event. Sediment mobilisation included slumping of the reservoir sediments, which helped contribute to hyperconcentrated flows downstream. Erosion of the reservoir sediments proceeded over weeks to months, and by week 15 it was calculated that >60% of reservoir sediment had been eroded (Wilcox et al., 2014).
Laboratory flume experiments have been used to investigate the geomorphic processes of rapid base level falls. Experiments were carried out by Cantelli et al. (2004) to investigate the geomorphological impact of dam removal. During an experimental run with complete dam removal, Cantelli et al. (2004) observed erosional narrowing, which resulted in initially high sediment discharges. Due to maximum boundary shear stresses in the middle of the channel, reducing to zero at the channel margins, incision rates were greatest in the centre promoting vertical downcutting and narrowing of the channel. This created a positive feedback where streamlines were focused towards the centre of the channel increasing erosion rates. Channel narrowing ceased when lateral sediment transport, from the increasingly higher banks, slowed channel degradation. Once degradation was stopped the channel widened (Cantelli et al., 2004).
In summary, therefore, the available DEM and satellite time series from Laguna del Viedma demonstrates rapid vertical incision from T3 to T4 with subsequent valley floor widening, consistent with observations from other ice dammed lake drainage events (e.g., Jacquet et al., 2017), and both field and flume dam removal students (Cantelli et al., 2004;Wilcox et al., 2014).

Implications for Late Quaternary landscape evolution
The Laguna del Viedma and Cachet valleys, dammed by outlet glaciers of the Southern and Northern Patagonia icefields respectively, can be considered modern analogues for landscapes of deglaciation in the Patagonian Andes. Late Quaternary deglaciation of Patagonia featured large outlet glaciers, which had the potential to dam tributary valleys. Glacier recession led to the formation of large proglacial lakes (Caldenius, 1932;Turner et al., 2005;Hein et al., 2010;Solari et al., 2012;García et al., 2014, FIG. 9. Selected LANDSAT imagery from the time series, showing a sector focused on the up-valley limit of the Laguna del Viedma in 1985. Key landforms and geomorphic processes are indicated. 2019; Thorndycraft et al., 2019), most of which either totally drained e.g., the Cisnes valley (García et al., 2019), or partially drained by punctuated drainage events e.g., lago General Carrera/Buenos Aires (Bell, 2008). These lake falls were in the order of 10s to 100s of metres and the largest could have caused catastrophic flooding, as demonstrated by the 110,000 m 3 /s flood reconstructed in the Baker valley (Benito and Thorndycraft, 2020). Across Patagonia, therefore, there was a spatially and temporally complex pattern of lake drainage events (Thorndycraft et al., 2019), which resulted in rapid base level lowering leading to the style of landscape evolution recorded in the Laguna del Viedma valley. Examples of Late Quaternary paraglacial landscapes (Ballantyne, 2002) from 46-47° S in Patagonia (Fig. 10), with significant valley incision, are presented in figures 11 and 12. The terrace sequence of the Río Jeinimeni valley (Fig. 11A) exhibits up to ~100 m of total incision in the reach photographed. At the 25 and 20 ka timeframes of the PATICE reconstruction of the former Patagonian Ice Sheet, Davies et al. (2020) show the valley to be ice free and therefore the valley has a similar geomorphic setting to the Laguna del Viedma, with a large outlet glacier blocking an ice-free river draining glaciated valleys (Fig. 10). Drainage from the Jeinimeni valley would have been blocked (creating an ice dammed lake) by the Buenos Aires ice lobe, the snout of which was located ~40 km to the east of the valley mouth until at least 17 ka (Bendle et al., 2017b). The landscape history of the palaeolake is yet to be detailed but the terrace sequence (Fig. 9a), based on the observed record from Laguna del Viedma, likely records a series of rapid incisions in response to lake level falls as the Buenos Aires ice-lobe thinned and receded.  (Thorndycraft et al., 2019). The valley floor at the Chacabuco site ( Fig. 12) was influenced by two base level falls: 1) ~11.7-12.6 ka: the drainage of lago Chelenko (Thorndycraft et al., 2019), a large palaeolake that joined lago General Carrera/Buenos Aires with lago Cochrane/ Puerreydon through the Baker valley (grey-blue shading); and 2) ~9.6±0.8 ka: the Bertrand-Baker catastrophic flood (Benito and Thorndycraft, 2020;Benito et al., 2021), which eroded bedrock along the indicated flow pathways.
The río Fénix Chico provides an example of a valley formed in a proglacial geomorphic setting (Fig. 10). The development of a proglacial lake at ~400 m a.s.l. at this locality occurred from ~18.0 ka (Bendle et al., 2017b). The geomorphological evidence of deltas and palaeoshorelines around present day lago General Carrera/Bueno Aires (Fig. 10) demonstrate a lake level fall of ~100 m (e.g., Bell, 2008), triggered by a drainage reversal from Atlantic to Pacific drainage (Turner et al., 2005), with timing constrained to 14.2-12.6 ka by Thorndycraft et al. (2019). Here, vertical incision following the lake level fall likely created instability in the saturated lacustrine sediments so the Fénix Chico landscape is dominated by mass movement processes (Fig. 11b), rather than fluvial terrace formation, as also documented at lago Puerreydón to the South (Pánek et al., 2020). The drainage reversal event led to the rapid formation of the Fénix Chico valley, which contrasts with the Laguna del Viedma setting, where catchment drainage has fed into the lake up-valley for millennia since ice-free conditions were established.
The final example, presented in figure 12, is from the Chacabuco valley in Chile, where there is evidence for multiple drivers of incision. The altitude of raised deltas (Fig. 12), ~40-50 m higher than equivalent landforms further east, demonstrate isostatic uplift (see Fig. 10 in Thorndycraft et al., 2019). The valley floor delta flowed into Palaeolake Chelenko, which drained ~12.6-11.7 ka (Thorndycraft et al., 2019). Furthermore, the catastrophic flood of the río Baker (see Fig. 10 for the flood pathway), which occurred ~9.6±0.8 ka (Benito et al., 2021), triggered geologically instantaneous bedrock incision of ~20 m downstream at the Baker-Chacabuco confluence (Benito and Thorndycraft, 2020). Maximum total incision at the site reaches ~90 m, consistent with the evidence from the modern Laguna del Viedma and Cachet II (Jacquet et al., 2017)  at this locality. The vertical incision downstream of the terraced tributary valley (Fig. 12) is consistent with rapid vertical incision of the Cantelli et al. (2004) flume experiments, however the terraces in the lower reach of the tributary valley could be a complex response to both rapid base level fall and isostatic uplift. The Chacabuco valley, therefore, has experienced a range of spatially variable rates of isostatic uplift as well as rapid base level falls. It is challenging to date Late Quaternary rapid incision events as the centennial, at best, uncertainty ranges on typical dates from fluvial terraces are larger than the decadal/sub-decadal timescales it took to form the Laguna del Viedma terraces. However, careful spatial and temporal analyses of landform records in the Patagonian Andes may yield insights to decipher likely controls (e.g., isostatic uplift versus base level fall) and relative rates of landscape incision (e.g., rapid versus centennial/millennial scales).

Conclusions
Satellite time series  and geomorphological mapping of the Laguna del Viedma valley (49° S), dammed by the Viedma glacier, demonstrates rapid terrace formation in response to lake level falls in the order of tens of metres. The available DEM shows that the formation of the ~20 m T3 terrace level occurred at a mean rate of 1.33 m/yr over a ~15 year time period, but the satellite evidence shows that most of the incision likely occurred more rapidly based on the apparent size of the scarp slope evident within months of the drainage event on satellite imagery.
The data have implications for the interpretation of rates of landscape evolution based on Late Quaternary landform records of the Patagonian Andes during the Late Glacial Interglacial Transition, given the widespread occurrence of palaeolake drainage events across the region. Examples of landscapes from the río Baker catchment (45-46° S), which we know to have been affected by rapid base level falls (e.g., Thorndycraft et al., 2019;Benito and Thorndycraft, 2020) were presented to demonstrate that sudden base level falls likely played an important role in landscape incision processes in addition to geophysical drivers of uplift such as isostatic uplift or neotectonics. FIG. 12. ALOS PALSAR DEM with transparent greyscale slope DTM overlay from the Chacabuco valley, Chile (see Fig. 1 for locations), illustrating a delta which fed into Palaeolake Chelenko (Thorndycraft et al., 2019) and river terrace sequence. Some key spot heights (m a.s.l.) are added to help demonstrate the amount of incision.
sensing and GIS class, so I thank students who took the course and contributed to discussions of the data. My thoughts on Late Quaternary landscape evolution in Patagonia have benefitted from discussions in the field with numerous colleagues and postgraduate students including: J. Bendle, G. Benito, E. Cooper, B. Davies, J.L. García, J. Martin, A. Nuevo-Delaunay, C. Méndez, A. Palmer, J. Pike and C. Sancho. I thank two anonymous reviewers whose comments improved the original manuscript.