An inventory of glacial lakes in the South Shetland Islands (Antarctica): temporal variation and environmental patterns

: Monitoring and inventorying proglacial lakes in the Maritime Antarctica region is essential for understanding the effects of climate change on these environments. This study uses Landsat images to create a map of lakes in ice-free areas of the South Shetlands Islands (SSI) for 1986/89, 2000/03 and 2020, and verifi c ation of patterns of change in lake areas and numbers. Normalized water difference index (NDWI) products, image segmentation, fi eld records, and cartographic products from other studies were used to validate the results. Results show a 60% increase in the number of lakes from 1986/89 to 2000/03; and a 55% increase from 2000/03 to 2020. There was a 52% increase in lake areas from 1986/89 to 2000/03; a 79% increase from 2000/03 to 2020; and a 173% increase from 1986 to 2020. From 1986 to 2020, the most signifi c ant changes were a decrease in the average elevation and distance from glaciers and an increase in distance from the sea. In 2020, SSI lakes were predominantly coastal and ice-marginal, with an E and S orientations, fl at surfaces, and a low declivity.


INTRODUCTION
Glacial lakes are expanding rapidly, and their evolution is connected to glaciers behavior, which, in general, have been retreating due to increasing mean air temperature (Carrivick & Tweed 2013, Wang et al. 2016, Shugar et al. 2020. Global warming exposes over deepened terrain and causes glacial melting, increasing the amount of water available to lakes and raising the number and size of glacial lakes (Che et al. 2004, Chen et al. 2007). Shugar et al. (2020) reported an increase in glacial lake number and total area by 53% and 51%, respectively, in an inventory mapping of glacial lakes worldwide. Hence, inventories of glacial lakes are critical for assessing the infl uence of climate change on the glacial regime (Raj & Kumar 2016).
The rapid development of remote sensing technologies has driven the monitoring of changes in glacial lakes, mainly due to the vastness of these areas and logistical diffi culty in accessing them (Quincey et al. 2007, Wang et al. 2013, Raj & Kumar 2016, Wang et al. 2020. Landsat series is one of the satellites with the most extensive temporal and spatial record for observing Earth's surface, and it provides products suitable for mapping glacial lakes (Nie et al. 2017, Wang et al. 2020. Several authors have successfully used Landsat scenes for this purpose in various parts of the Cryosphere, focusing on the Himalayas in particular. Zhang and Wang et al. (2020) for the mountainous regions of Asia.
Some glacial lake characteristics are outlined and considered in studies attempting to organize inventories using remote sensing techniques. Some lakes are seasonal -disappearing or freezing at certain times of the year, making a transient and dynamic environment (Carrivick & Tweed 2013, Jawak 2015 -and turbidity modifies reflectance in the visible and near-infrared wavelengths (Gardelle et al. 2011). There is also interference from clouds, mountain shadows, and other sources of error (Sheng et al. 2016). Furthermore, using different data sources and spatial resolutions can lead to inconsistencies when analyzing changes in glacial lakes (Salerno et al. 2012, Nie et al. 2017. Aside from the inventory, gathering data on lake location, area, and number requires considering the environmental importance of these water bodies in the Antarctic system. Environmental configuration, according to Carriwick & Tweed (2013), has a significant influence on the growth, contraction, filling, emptying, and persistence of the proglacial lakes, which is critical for understanding past environmental conditions and predicting future environmental changes. For the Antarctic Peninsula, Bozkurt et al. (2021) expect an increase in annual-mean near-surface temperatures of about 0.5-1.5 °C for the  period. This temperature increase may cause changes in icefree areas, leading to new glacial lakes.
Other authors have developed glacial lake inventories in ice-free areas of Maritime Antarctica. In the SSUm, Rosa et al. (2021) and Oliveira et al. (2021) created glacial lake inventories for Nelson and King George Islands (KGI) from 1986 to 2020. However, these studies did not assess the distance between lakes and glaciers or between lakes and the coastline.
These data are critical for understanding environmental evolution and the increase in the number and size of lakes in proglacial areas (i.e., near glaciers) or other parts of the analyzed areas. This study uses Landsat images to create a map of lakes in ice-free areas of the SSI from 1986 to 2020, with subsequent verification of patterns of change in the lake area, number, and location. It investigates the evolution of lakes on each island and between the SSI islands. It addresses the contrasts in the studied period by observing variations in morphometric and spatial parameters. Furthermore, conducting a lake inventory for all ice-free areas of SSI, which has not previously been developed, and presenting a semi-automatic method for lake mapping -allowing researchers to monitor lakes in other sectors with ice-free areas.

Study area
The SSI is an archipelago located between latitudes 61°S and 63°S, and longitudes 63°W and 53°W. SSI's largest islands are King George, Nelson, Robert, Greenwich, Snow, Low, Deception, Livingston, Smith, Elephant, and Clarence (Ochyra, 1998). This study did not include the northern SSI (e.g., Elephant and Clarence). The Drake Passage separates the archipelago from South America, and the Bransfield Strait separates it from the Antarctic Peninsula (AP) (Figure 1).
The South Shetland Islands are in the Maritime Antarctic region where land is covered by ice caps and permanent snowfields -although some peninsulas, such as Fildes Peninsula in King George Island (KGI) and Byers in Livingstone Island, are free of permanent ice (Watcham et al. 2011). According to Bremer (2008, ice-free areas of SSI are drained by intermittent channels, and morphodynamics is associated with physical and chemical processes caused by snow and ice melting. According to Lopez-Martinez et al. (2012), the Smith (Mt. Foster, 2012 m) and Livingston (Mt. Friesland, 1770 m) islands have mountainous relief, whereas Snow (305 m), Low (180 m), Robert (385 m), and Nelson (332 m) islands have low and flat relief. The authors also note that only 10% of the 4700 km² of raised platforms of marine origin are free of ice.
Several studies have shown glacier retreat on the SSI (Rückamp et al. 2011, Szio & Bialik 2018, Pudelko et al. 2018, Perondi et al. 2020, which has been linked to an increase in air temperature in the AP region from 1948 to the first decades of the 21st century , Turner et al. 2005). The SSI have calving glaciers (Silva et al. 2019). Pudełko et al. (2018) show water-terminating (lake or lagoon) and land-terminating glaciers in King George Island. According to Jonsell et al. (2012), these glaciers are more sensitive to climate change as the ice is near the melting point under pressure.
The annual mean air temperature is -2.8°C, and the austral summer (December to March) mean air temperature is slightly higher than 0°C, with strong winds and high weather variability (Ferron et al. 2004

Imagery
The summer season is when there is less snow cover and the water supply is influenced by glacier melting and liquid precipitationaccording to Zhang et al. (2015), this configures the best period to make lake inventories -this study uses images acquired on different dates during the austral summers of 1986/89, 2000/03 and 2020, as shown in Table I.

Digital Elevation Model (DEM)
The DEM used in this study is from the Antarctic Elevation Reference Model (REMA) project (Howat et al. 2019), which is a time-stamped Digital Surface Model (DSM) of Antarctica that includes data from WorldView-1, WorldView-2, and WorldView-3, as well as a small number from GeoEye-1. Data were acquired between 2009 and 2017, with the majority collected during the austral summer seasons of 2015 and 2016. Each REMA DEM is vertically registered to satellite altimetry measurements (Cryosat-2 and ICESat), resulting in absolute uncertainties of less than 1 meter over most areas (Howat et al. 2019). Due to REMA-8 gaps, there are some missing-data pixels in the DEM data. Hence, neighborhood-analysis and interpolation were performed to create a multi-source DEM-derived product (REMA-2 and TanDEM-X DEM) to fill in all no-data values. TanDEM-X, REMA-8, and REMA-2 have high resolution (12 m, 8.8 m, and 2.2 m, respectively) and show high vertical accuracy (RMSE 2.9 m, 0.6 m, and 3.5 m, respectively) (DLR 2019, Howat et al. 2019). The final product is a DEM with an 8-meter spatial resolution.

Field data and supporting materials
During the austral summers of 2013, 2014, 2015, and 2020, four fieldwork campaigns were conducted on King George Island (Admiralty Bay and Fildes Peninsula) and Nelson Island (Harmony Point) to collect lake photographs and GPS coordinates ( Figure 2).
Supporting materials are used to verify lake pre-existence and validate polygons whose characteristics do not determine whether the vector generated refers to a lake. SPOT images from 1988 and 2000 (at the end of the summer season) are used to assess lake boundaries and a Sentinel-2 MSI image from 2020. Google Earth's high resolution satellite imagery base map, available in QGIS 4.3 (Zhang et al. 2015, Khadka Hemisphere. An error in image registration of ± 0.5 or ± 1 pixel is admitted, with these values being also found by other authors (Gardelle et al. 2011, Nie et al. 2013, Li et al. 2020. Due to clouds in some of the ice-free areas, mosaics of images from different dates are created for the images obtained in the 1980s and 2000 (Nie et al. 2015. In order to detect glacial lakes, all shaded areas in images are removed. A steepest slope > 30% mask is created using DEM data. This masking procedure is essential for eliminating shadow pixels since the spectral response of water is very similar to shadow values recorded by sensors, as demonstrated by Mostafa and Abdelhafiz (2017). The normalized water difference index (NDWI), given by the [(Green -NIR) / (Green + NIR)] formula (McFeeters 1996), is calculated after mask application. Green and NIR wavelengths correspond to B2 and B4 in Landsat TM and ETM+ sensors, and B3 and B5 in the OLI sensor, respectively. NDWI aims to distinguish bodies of water from other land targets (Li et al. 2011), with values ranging from 0.60 to 0.95, usually corresponding to the surface of lakes (Huggel et al. 2002).
Cropped NDWI images are segmented using the region growing method, with a maximum of five iterations and an area of 8 neighboring pixels considered (Rosa et al. 2021). Lakes are then manually vectorized using a combination of segmentation data and color compositions of radiometrically corrected images, yielding a method based on an integrated approach (Huggel et al. 2002, Li et al. 2011, Rai et al. 2017, Khadka et al. 2018, Wang et al. 2020. It should be noted that glacial lakes are identified and edited by one single researcher. Due to ice and snow-covered lakes, the methodological assumptions of Carrivick & Quincey (2014) are applied in some images, and vector files from the first mapped period (1986/89) are copied to the following period (2020).
Lake vectors are validated using photographs, GPS points, SPOT, Sentinel-2, and Google Earth imagery (section 2.4). The control points obtained in the field (Figure 2a-e) are critical so that vectorization can be performed in other environments due to the similarity of NDWI values and visual interpretation.

Lake spatial interpretation
Following lake mapping, each polygon is labeled with a number (Gardelle et al. 2011) and a set of morphometric and spatial attributes such as elevation, slope, aspect, distance from the glacier, and distance from the ocean are calculated. These parameters are selected because they are the primary determinants of the geomorphological and sedimentary characteristics of the lakes (Rubensdotter & Rosqvist 2009, Carrivick & Tweed 2013), indicating potential causes for changes in the attributes and location of the lakes, and facilitating data interpretation over time. Elevation, slope, and aspect parameters are derived from DEM using QGIS 3.4 software -a 100-meter buffer is created around lakes, and the "basic statistics for numeric fields'' tool is used to determine the mean value of each parameter in a given class.
Spatial parameters -distance from the glacier and distance from the ocean -are divided into distance ranges. The distance to the sea is divided into five intervals: 0-500 m, 500-1000m, 1000-1500 m, 1500-2000 m, and 2000-2500 m. The distance to the glacier is divided into eleven intervals of 500 m due to its greater amplitude, which ranges from 500 to up to 5000 m.
Collected data is subjected to descriptive statistical analyses (frequency tables, measures of central tendency, and variability -mean and standard deviation) to assess the quantity and characteristics of lakes in the studied years, as well as inferential statistical analyses: Kolmogorov-Smirnov Normality test, to support the use of parametric tests; and Student's t-test for independent samples, to assess statistical differences between parameter averages for the 1980s and 2020s (increase or decrease in parameters over the years). All analyses are performed in Microsoft Excel and SPSS v.21 (Statistical Package for Social Sciences, version 21), and for all tests, a significance level of 5% is established. Figure 3 summarizes the methodology used in this study.

Lake identification
Total results for the SSI For the years 1986/89, 155 lakes are identified by the method, with 84 (54%) obtained by the combination of shadow masking, NDWI, and segmentation (semi-automatic method); and 71 (46%) identified through validation using the supporting materials listed in section 2.4. For the years 2000/03, 248 lakes are identified, with 161 obtained using the semi-automatic method and 87 using validation. For the 2020 image, 385 lakes are identified, with 328 (85%) using the semi-automatic method and 57 using validation. The number of lakes increased by 60% from 1986/89 to 2000/03 and 55% from 2000/03 to 2020. From 1986 to 2020 (34 years of monitoring), the number of lakes increased by 148%.

Results by islands
As for increases in area, the three islands with the greatest increases in total area from 1986/89 to 2000/03 are Livingston (0.32 km²), King George (0.31 km²), and Greenwich (0.12 km²). In percentage, the three islands with the highest increase are Greenwich (300%), Nelson (167%), and Livingston (114%). From 2000/03 to 2020, the three islands with the greatest increase in absolute area value are Livingston (0.79 km²), Deception (0.47 km²), and King George (0.36 km²). In percentage, the three islands with the highest increase are Snow (800%), Livingston (132%), and Nelson (94%). Smith Island has no area values because no lakes are identified in the assessment period. Figure 6 depicts a graph showing the absolute increase in area per island over the two periods analyzed.
Significant differences are observed (Table  III) for the elevation data (p-value = 0.045, Student's t-test, 5% significance level) with the average being higher in 1989 than in 2020, for distance from the ocean data (p = 0.033) with the average being higher in 2020 than in 1989, and for distance from glaciers data (p-value = 0.043) with the average being higher in 1989 than in 2020.

Lake location patterns in islands
Between 1989 and 2020, lakes were mostly located at elevations lower than 197.6 m. In 2020, lakes began to occur at lower elevations on most islands compared to 1989 (Table IV). Exceptions are Robert and Livingston islands, with more lakes located at higher mean elevation ranges (with little variation) in 2020 than in 1989. When compared to other islands, Deception Island shows the greatest variation in lake location by mean elevation between 1989 and 2000 ( Figure  8). Greenwich, Nelson and Robert Islands show the smallest variation in lake location by mean elevation.
As for slope, there is a predominance of lakes located in slopes percentage ranging from 10% to 15% on islands (Figure 8). There are two exceptions: King George and Deception islands, which have lakes in areas with a higher slope than the other islands in 1989 and 2020. Greenwich Island has the greatest variation in lake location by slope class, while Nelson Island has the smallest variation.
Considering the distance from the ocean parameter (Figure 8), the islands of Nelson, Robert, Livingston, and Low are highlighted as they present new lakes located farther away from the coast . When years and locations of lakes are compared, KGI, Deception, and Greenwich islands present lakes with the closest proximity to the coast in 2020. Livingston Island shows the greatest variation in the environmental context in this parameter, followed by Deception Island. Robert Island shows the smallest variation in terms of the distance from the ocean parameter.
When comparing years, the average distance between lakes and glaciers on the islands KGI, Nelson, Deception, Greenwich, Low, and Livingstone has decreased over time. Only Robert Island shows new lakes (2020) that are located further away from glaciers. Lake location on KGI Island varies the most concerning the average distance from the glacier margin ( Figure  8).   Results also show that the percentage of lakes identified using the semi-automatic method has increased over time, reducing the need for validation, reducing errors, and promoting effective monitoring. In a dynamic environment with a large spectral mixture of targets, mapping gains from using Landsat 8 OLI images are substantial. Landsat images, however, are limited in their ability to identify small lakes due to their spatial resolution. Thus, given that the vast majority of recent lakes are small, it is noteworthy that the use of Sentinel-2 images may expand the possibility of monitoring lakes in Maritime Antarctica, as several authors have demonstrated the importance of understanding the behavior of proglacial lakes in conjunction with glacier retreat, e.g., Himalayas (King et al. 2019, Watson et al. 2020. Following the assumptions of Gardelle et al. (2011), the smallest glacial lake detected in this study is 4 pixels (0.0036 km²) in size, indicating that our study may have underestimated the number of lakes in the icefree areas of the South Shetlands. Wetlands, which are found around lakes and have NDWI values similar to water bodies, add to the complexities of identifying lakes using segmentation and NDWI methods. Therefore, it is recommended that wetlands areas be mapped -an inventory that has yet to be completed for the ISS -and that they are linked to the presence of lakes. This is most evident in portions with the most ice-free areas of SSI as a result of deglaciation processes (John 1972), such as Harmony Point and Stansbury Peninsula (Nelson Island); Fildes, Barton, and Potter Peninsulas (King George Island); and Byers Peninsula (Livingston Island). Lake mapping in these areas is also hampered by shadow interference since their terrain relief is linked to subglacial erosive action. As observed in fieldwork, there are shadows in sectors of cirque valleys, arches, rocky promontories, and U-shaped valleys.
Turbidity may also have hampered lake identification, particularly in areas adjacent to the ice, which have high sedimentary production due to paraglacial activity (Ballantyne 2002) and highlight the rapid evolution of landscapes undergoing glacier shrinkage (Williams & Koppes 2019). In addition to coping with the high cloud cover in the region, which limits satellite image acquisition, proglacial lakes with high sediment input are limiting for semi-automatic vector identification, characterizing this environment as one of the most complex and difficult to map in Maritime Antarctica. As a result, the inventory presented in this research required manual verification, with vectorization of glacial lakes with the aid of high-resolution images available for the study area.
Due to the transience and dynamics of proglacial environments, a satellite image represents only one moment during the ablation period, and the limit of a lake is only an instantaneous state on a given day (Wang et al. 2016, Rosa et al. 2021. Almost all paraglacial landforms and landscapes are transitory (Slaymaker 2009). Therefore, proglacial lakes can expand, coalesce, decrease in size, disappear, or become detached from an ice margin (Carrivick & Quincey 2011).

Increase in the number and area of lakes
Over the 34-year analysis period, the gain in lake number (230 new lakes) and area (3.13 km²) observed represents considerable environmental changes for the ice-free areas of the SSI. Variations in lake area show that the period 2003-2020 was the most significant for the percentage increase in the lake area, while the period 1989-2000 was the most significant for the percentage increase in lake number. It's important to note that there is evidence of an increase in mean air temperature for both periods and evidence of glacier retreat -the formation of new ice-free areas results in environments that are more prone to lake formation.
According to some studies in KGI, for example, a high retreating trend was observed for some glaciers in both periods, and these changes may be related to both the increasing trend of regional average air temperature and the variability of average winter air temperature since the 1980s (Rosa et al. 2021). In recent decades, Rafa et al. (2018) linked glacier area loss to melting degree days in the twentieth and twenty-first centuries in KGI. The year 2020 is highlighted in particular because it sets the record for the highest air temperature in the history of the Antarctic. The WMO recorded a temperature of 18.3°C for the 12 hours (LST) of February 6, 2020, at Base Esperanza in the northern tip of the Antarctic Peninsula, according to Rocha et al. (2021). These changes are most likely associated with an upward trend in annual Positive Degree-Days (Pudełko et al. 2018), as well as an increase in mean surface air temperature , Turner et al. 2005, and in sea surface temperature -SST (Vaughan 2003) for the Bellingshausen Sea.
It is emphasized that the climatic factor is critical in causing glacier retreat, exposing new ice-free areas, and causing melting to form lakes. Robert, Nelson, and KGI are three specific cases that are examined. Robert Island, despite having an increase of 4 lakes from 2003 to 2020, these are small lakes near the edge of the glacier and in only two ice-free areas that already had lakes in 1989 -demonstrating no abrupt environmental change. Nelson Island, which had lakes in 6 new ice-free portions, and KGI, which has developed 10 new ice-free areas in the period of analysis, all with the presence of newly mapped lakes.
Other studies for the KGI have shown an increase in ice-free areas and marginal and coastal lakes due to glacier retreat. Several studies have shown glacier retreat on the KGI (Rosa et al. 2015Perondi et al. 2020). KGI has ice-marginal environments with predominantly flat to gently undulating terrain -in some parts -, which favors the enlargement of lakes over time. This can also be seen in ice-free areas on other islands. Carrivick and Quincey (2014) showed for areas of Greenland that the expansion of the lakes occurs due to the melting of water storage in lakes. One of the consequences of this process, according to the authors, is a reduction in water and sediment supply to the ocean.
KGI, which concentrated most of the lakes in 1986/89 (57.8%), experienced a decrease in value in 2020 (43.9%), indicating the occurrence of important processes such as glacier retreat and the expansion of ice-free areas in the other islands of the South Shetlands. Despite this decrease in lake number proportion, KGI continues to be the island with the most lakes and the greatest increase in lake number (79). Fildes Peninsula, the largest ice-free peninsula in SSI, is located in the KGI and has an ice-free area of approximately 28.8 km² (Andrade et al. 2018). In 2020, Fildes Peninsula alone had 62 lakes, more than any other island in the study area except Livingston (73 lakes).
Livingston Island had the greatest absolute and relative increase in the number of lakes and the greatest absolute increase in area between 1989 and 2000. Between 2000 and 2020, Livingston Island ranked third in total lake growth and first in total area growth. Livingston Island nearly equaled the KGI in terms of area value in 2020. High values of increase in area and number of lakes for Livingston Island are due to the decrease in thickness and shrinkage in glaciers on this island, according to Molina et al. (2007). Furthermore, Byers Peninsula, located in Livingston Island, has been exposed for a long time and is responsible for the variations in distances observed between the lakes and coastline and between lakes and glacier fronts when the 1989-2020 periods are compared.
Results show that Deception, Snow, Greenwich, and Nelson islands also experienced substantial gains in the lake area and/or number. Ice-free areas on these islands may be in topographic areas that allow for the formation of lake basins. Therefore, mapping the evolution of glacier retreat and geomorphological mapping of ice-free areas can aid in understanding the processes that promote lake growth. According to Lopéz-Martinez et al. (2012), the SSI has 33 different types of periglacial relief, demonstrating the great diversity of relief forms and processes in Maritime Antarctica compared to continental Antarctica.
Greenwich Island differs from the others in that it is the only one to show the period 1989-2000 as the greatest gain in absolute lake area -rather than 2000-2020; additionally, for 1989-2000, Greenwich is the island with the highest percentage of increase in area and the second-highest percentage of gain in number, which is consistent with Stastna (2010). According to Stastna (2010), the Chilean Prat station recorded a rate of increase in mean air temperature of + 0.38 °C per decade for Greenwich Island -which showed the greatest increase in area between 1981and 2000. Fatras et al. (2020 report 20.4% of the decrease in the area for four glaciers in the north of Greenwich Island between 1956 and 2019, notably the portions that recorded the development of lakes. Glaciological factors also explain the occurrence of lakes. Islands with no lakes (such as Smith) or few lakes (such as Robert, Low, and Snow) are characterized by ice caps that distinguish them from other islands in the study area. However, Snow Island is highlighted in this lake inventory because it has experienced the most dramatic changes -going from zero to 22 lakes in 34 years. The presence of lakes on Snow Island can alter the sedimentary dynamics of this environment by interfering with the passage of melting water from a glacier and reducing the speed of water flow, resulting in sedimentation of eroded material, according to Carrivick & Tweed (2013).

Morphometric parameter analysis
There is a pattern in lake location in the SSI from 1989 to 2020, with lakes not occupying levels higher than 250 m. Lakes have an average declivity of 13.1%, and they do not occur more than 2500 meters from the coast or more than 10,000 meters from glaciers. For 2020, there are primarily coastal lakes, marginal to the ice, eastoriented, on flat surfaces, and a low percentage of slope in the vicinity of the lake.
The majority of new lakes are proglacial, verified by a decrease in the average distance from glaciers over the years for all islands except Robert. For 2020, most lakes (58.3%) are in recent ice-marginal environments (class 0-500m) formed by changes in terrestrial glaciers. Lakes with the greatest average distances from glaciers are found on KGI, Nelson, and Livingston islands, all of which have significant ice-free areas, as explained in the previous section. On the island of Nelson, the appearance of proglacial lakes in 9 new ice-free areas corroborates the significant decrease in the average distance from the glacier -from 2388.9 m to 918.4 m between 1989 and 2020. However, the predominance of lakes in Harmony Point and Stanbury point has resulted in an increase in the average distance from the glaciers due to the retreat of the ice masses. Therefore, a future analysis by ice-free sectors is recommended for those parameters.
The comparative analysis identified significant changes in distance from the ocean (p = 0.033). Data show two trends in terms of the distance between lakes and the coast. There is a rate of increase in the average distance from the coastline, related to glaciers with terrestrial fronts that were retreating (e.g., Nelson, Robert, Low, and Livingston). Compared in time, the lakes of Deception, Greenwich, and KGI were closer to the coast. Oliveira (2021) observed for KGI that tidal-ending glaciers receded faster, causing lakes to form near the shoreline. Such changes can lead to glacier melting and shrinking processes (Jonsell et al. 2012) and the transformation of marine glaciers to non-marine glaciers and new ice-marginal lakes (Perondi et al. 2020, Oliveira et al. 2021.
The increase in lakes in classes above 1000 meters from the ocean may be related to deglaciation processes and landscape evolution of ice-free areas. With the continuation of glacier retreat, the formation of latero-frontal moraines, and changes in the proglacial landforms, lakes can form in areas further away from the ocean and at higher altitudes. Furthermore, advances in the radiometric resolution of Landsat 8 OLI images may have increased the identification of periglacial lakes, raising the average of lakes farther away from the coastline.
In general, there is a decrease in the average hypsometric (81.8 m to 75.6 m) of lake locations during the 1989-2020 period. Only on the islands of Livingstone, and Robert, the hypsometric mean shows an increase in this period -these islands have lakes further away from sea level (i.e., at lower elevations). Robert Island is also the only one to show an increase in the average distance between lakes and glaciers.
When the slope of the land surrounding the lakes is considered, most of the lakes have slopes ranging from 2 to 15% -as also found by Oliveira et al. (2021). From 1989 to 2020, there were no changes in the average slope value of the areas surrounding the lakes, which remained at 13.1%. Deception Island has the highest slope values, which is explained by its characteristics. The morphology of Deception Island is influenced by two main geomorphological processesvolcanic activity and glacial action, with glaciers covering approximately 57% of the island surface. An almost complete ring of steep hills encircles the sunken interior of Port Foster bay and the outer coast. Thus, melting ice and snow, periglacial activity, slope processes, and the presence of permafrost have also contributed to the formation and evolution of the island relief, with the development of ice streams and lakes that occupy several craters (Smellie et al. 2002).
As for slope orientation, there are changes primarily to north and south for lakes identified in 1986 and to east and west in 2020. According to Andrade et al. (2014), slopes oriented north have higher solar radiation values than slopes facing south. The issue of solar radiation influences the melting of snow and glaciers, resulting in the formation of lakes, but not in the same way that meltwater channels flow, for example, which may explain the formation of lakes in other aspects. Local morphometry issues interfere with meltwater accumulation and lake formation in periglacial environments due to subglacial reliefs, whereas paraglacial processes predominate in proglacial zones, with rapid post-depositional change. Therefore, a DEM from 2016 may not accurately represent either the past or the current proglacial environment.
Given the diverse scenarios presented and the fact that SSI is subjected to similar climatic conditions, the variations may be related to glaciological and environmental factors. Livingston, Greenwich, and Deception have more diverse contexts of location in the parameters assessed. Less variation can be found on islands with more glacial coverage and flatter exposed areas (e.g., Robert & Nelson). In the context of SSI, Perondi et al. (2020), Rosa et al. (2021) have already shown that lake formation is influenced by its location as well as other factors such as glaciers slope, aspect, altitude, and coastline distance, which was also proven in this research.
Understanding the geomorphological scenario in which the lakes are inserted allows for making inferences about glacier retreat patterns and understanding how the ISS proglacial system has evolved and behaved. As ice-free areas are highly vulnerable to climate change, as demonstrated by López-Martinez et al. (2016), their mapping and monitoring should be developed. Furthermore, as the lake inventory proposed in this research may be used to support other studies, knowing the morphometric parameters is essential to make inferences about climate variations (Rubensdotter & Rosqvist 2009); provide the interpretation of processes and reconstructions of past environments and processes (Carrivick & Tweed 2013); and identify 'marginal ice' or 'contact ice' lakes and distal proglacial lakes as they are distinct from each other in terms of geomorphology (Carrivick & Tweed 2013).

CONCLUSION
This first inventory of all lakes located in icefree areas of the South Shetland Island and their morphometric and spatial parameters, is important because this information can be used in other studies, such as paleoclimatic paleoenvironmental studies. Future research may reveal the impact of the observed changes in the modification of local hydrological flow regimes.
Temporal lake monitoring in SSI, as well as glacier monitoring (Cook et al. 2014), is a challenge due to some factors observed while surveying available satellite images, such as partial or total cloud coverage in the images, lakes with dimensions smaller than those possible to be detected in Landsat images, and the phenology of lakes, which implies the freezing of the liquid surface. Given the limited availability of satellite images, the importance of the inventory and monitoring presented in this study is highlighted since there is a limitation for the inclusion of other periods for temporal monitoring.
Given the difficulty of obtaining optical images as previously stated, the integrative approach proposed in this study, which included satellite image processing products and fieldwork in conjunction with previously published mappings, was vital for the Antarctic environment. Furthermore, the points obtained in the field were critical in allowing validation in some free-of-ice areas and allowing extrapolation to other ISS areas.
Data on lake and area percentage gain shows a considerable amplitude; however, it is noteworthy that all islands show an increase in data. The greater percentage gain and number of lakes in a specific ice-free area represents different hydrosedimentological conditions for the proglacial portions in the face of glacier retreat. Because of large amounts of sediments and meltwater present in paraglacial processes, lake formation usually results in sediment damming and disruption of the transfer of sedimentary material and meltwater to the sea, representing important environmental changes for the area.
As previously documented by other SSI studies, increases in lake number and area are linked to glacier retreat. The decrease in lake predominance in KGI indicates that other islands, such as Deception, Greenwich, Livingston, and Snow, are undergoing pronounced environmental changes, suggesting locations for future research and monitoring.
In general, it is concluded that glacial melting water and the formation of new ice-free areas are the primary sources of the majority of lake expansion and formation from 1986 to 2020. This is because the distance between lakes and glaciers decreased in the year 2020. Furthermore, increases in average distance from the ocean demonstrate that many ice-free areas of the SSI are in different stages of proglacial landscape evolution. These glaciers changed their fronts from marine to land since 1986/89, resulting in the formation of new lakes.
Ultimately, it is emphasized that lake type classification, mapping of fluvial drainage networks in ice-free areas, and investigation of glacial retreat patterns can help to understand better contrasts presented by the values of lake variation per island evidenced in this research. However, such analysis requires the use of higher resolution satellite images and the implementation of additional fieldwork. Changes in slope and orientation of the lake's area suggest that future research should concentrate on geomorphological analyses of the local relief, associating it with paleo reliefs or paraglacial environments.