Introduction
Situated in the north-eastern Antarctic Peninsula, James Ross Island contains more than 100 glaciers, which covered around 1780 km2 in 2009 (Davies and others, Reference Davies, Carrivick, Glasser, Hambrey and Smellie2012). Mount Haddington Ice Cap, with its outlets, accounts for 80% of the glacierized area of the island coalescing with Dobson Dome in the southern part on the Ulu Peninsula (Fig. 1a). The northern part of Ulu Peninsula, the mostly deglaciated part of James Ross Island, contains small alpine glaciers, which account for <2% of the total glacierized area on island. These glaciers lie 10–30 km from the north-eastern coast of the Antarctic Peninsula and ~73–91 km from its northern tip. Their location on the eastern side of the Antarctic Peninsula, dominated by the ice-bound Weddell Sea, provides a polar continental climate with cold and dry winters and above-freezing summer conditions (Turner and others, Reference Turner2016; Ambrožová and others, Reference Ambrožová, Láska and Kavan2020). The near-zero air temperature during the melt season in this region, together with low-elevated accumulation areas and small ice volume renders land-terminating glaciers on the Ulu Peninsula vulnerable to rising temperatures (Davies and others, Reference Davies, Carrivick, Glasser, Hambrey and Smellie2012; Engel and others, Reference Engel, Láska, Nývlt and Stachoň2018). As a result, these glaciers have a short response time and are highly sensitive to air temperature (Davies and others, Reference Davies2014; Engel and others, Reference Engel, Láska, Kavan and Smolíková2023), and may provide insights into the recent mass changes of similar small glaciers in the Antarctic Peninsula region. In addition, these glaciers contribute significantly to the runoff of local streams and preserve adjacent fresh water ecosystems, which represent the most important terrestrial habitats of Antarctic life (e.g. Hawes and Brazier, Reference Hawes and Brazier1991; Kavan and others, Reference Kavan, Nedbalová, Nývlt, Čejka and Lirio2021).
Due to the remote location of glaciers in the northern part of James Ross Island, little information exists regarding recent changes in their topography and mass balance. However, remote-sensing studies indicate that ice caps, valley and cirque glaciers on the Ulu Peninsula experienced large retreat rates between 1979 and 2006 (Engel and others, Reference Engel, Láska, Kavan and Smolíková2023 and references therein). These data also indicate that land-terminating glaciers on James Ross Island retreated faster from 1988 to 2001 than between 2001 and 2009 (Davies and others, Reference Davies, Carrivick, Glasser, Hambrey and Smellie2012; Silva and others, Reference Silva2020). A further decrease in the retreat rates was reported for small glaciers on the Ulu Peninsula (Engel and others, Reference Engel, Kropáček and Smolíková2019) and for outlet glaciers around James Ross Island (Lippl and others, Reference Lippl2019). A period of reduced retreat rates terminated in the mid-2010s when retreat rates started to increase in the northern part on the Ulu Peninsula (Engel and others, Reference Engel, Kropáček and Smolíková2019, Reference Engel, Láska, Kavan and Smolíková2023). The rapid rate of glacier retreat at the end of 20th century is in accordance with a long-term increase in the near surface air temperature (e.g. Vaughan and others, Reference Vaughan2003) while the subsequent decelerated retreat reflects the temporary cooling at the turn of the millennium (Turner and others, Reference Turner2016, Reference Turner2020; Oliva and others, Reference Oliva2017). The return to a rapid retreat rates after 2015 coincides with a prominent increase in both summer and mean annual air temperature (Engel and others, Reference Engel, Láska, Kavan and Smolíková2023). Although annual precipitation values were highly variable, winter precipitation increased during the 2010s in the northern part of the Antarctic Peninsula (Carrasco and others, Reference Carrasco, Bozkurt and Cordero2021).
Just two case studies have been published on the surface mass balance (SMB) of the glaciers on James Ross Island. The initial study describes changes in SMB of Davies Dome (Fig. 1b) and Lookalike Glacier (referred to as Whisky Glacier or IJR-45 in previous publications) during the period 2009–15 (Engel and others, Reference Engel, Láska, Nývlt and Stachoň2018), while the second focuses on the effect of air temperatures and snow height changes on the SMB of Triangular Glacier (Engel and others, Reference Engel, Láska, Kavan and Smolíková2023). The data indicate a predominantly positive SMB until the mass-balance year 2014/15 and persistent glacier mass loss afterwards (Engel and others, Reference Engel, Láska, Nývlt and Stachoň2018, Reference Engel, Láska, Matějka and Nedělčev2022). The change in mass-balance trends around 2015 was attributed to the prominent warming along the eastern coast of the Antarctic Peninsula and the pronounced response of small glaciers to melt-season temperatures (Engel and others, Reference Engel, Láska, Kavan and Smolíková2023). The sensitivity of these glaciers to air temperature has also been observed elsewhere on the north-eastern side of the Antarctic Peninsula (Abram and others, Reference Abram2013; Davies and others, Reference Davies2014; Marinsek and Ermolin, Reference Marinsek and Ermolin2015) but their response to precipitation (sensu Navarro and others, Reference Navarro, Jonsell, Corcuera and Martín-Español2013) remains unknown. Meteorological and glaciological observations on Triangular Glacier have shown that glacier surface decreases mainly by snow melt, but wind transport of snow can be also effective (Engel and others, Reference Engel, Láska, Kavan and Smolíková2023).
This study aims to reveal the response of small glaciers on James Ross Island to recent atmospheric warming in the northern Antarctic Peninsula region. SMB data collected on Davies Dome and Lookalike Glacier since the mass-balance year 2015/16 are presented and compared with the available SMB records from the northern Antarctic Peninsula region. Recent changes in geometry of these glaciers are assessed using satellite imageries, digital elevation models (DEMs) and ice thickness data derived from a ground-penetrating radar (GPR) survey. Finally, the air temperature records from Johann Gregor Mendel Station (63°48’S, 57°53’W, 10 m a.s.l.; hereafter referred to as Mendel) and Esperanza Base (63°24’S, 56°59’W, 13 m a.s.l.) weather stations are analysed to connect the observed changes in glacier mass to the mean summer air temperature and number of positive degree-days (PDD).
Study sites
Davies Dome (GLIMS Glacier ID: G301945E63889S), Lookalike Glacier (G302049E63932S) and Triangular Glacier (G302151E63856S) are located in the predominantly ice-free northern part of James Ross Island, which is situated in the precipitation shadow of the high-elevation (1500–2100 m a.s.l.) plateau regions of the northern Antarctic Peninsula (Fig. 1a). The northern part of James Ross Island, the Ulu Peninsula, is composed of lava-fed delta sequences, which overlie Cretaceous marine sediments (Mlčoch and others, Reference Mlčoch, Nývlt and Mixa2019). The sedimentary strata form a low-relief landscape, while the volcanic rocks build prominent cliff-bounded mesa landforms with surfaces at c. 300–600 m a.s.l. (Jennings and others, Reference Jennings2021). Small cold-based ice caps cover the highest parts of these mesas, while headwall cliffs bound the upper part of cirque glaciers (Fig. 1b). The mean annual air temperature is −6.7°C at sea level, but is <−9.0°C at higher-elevation (Ambrožová and others, Reference Ambrožová, Láska, Hrbáček, Kavan and Ondruch2019; Kaplan-Pastíriková and others, Reference Kaplan-Pastíriková, Hrbáček, Uxa and Láska2023). Precipitation is modelled to be between 300 and 700 mm in water equivalent (w.e.) per year (van Wessem and others, Reference van Wessem2016; Palerme and others, Reference Palerme2017) but the prevailing south to south-westerly winds strongly affects the distribution of snow cover, transporting snow from flat upper portions of glaciers to their lee sections (Engel and others, Reference Engel, Láska, Nývlt and Stachoň2018; Kavan and others, Reference Kavan, Nývlt, Láska, Engel and Kňažková2020).
Davies Dome is an ice dome with a flat-topped surface at c. 500 m a.s.l. (Fig. 2a). The dome is located on a volcanic plateau with an elevation of ~400–450 m a.s.l. and its southern part is drained by a 2.5 km long marine-terminating outlet that terminates in Whisky Bay. In 2010, the maximum ice thickness of the dome determined from GPR data was 83 ± 2 m (Engel and others, Reference Engel, Nývlt and Láska2012). A cumulative mass gain (0.11 ± 0.37 m w.e.) was observed over the period 2009–15 on the dome but the mass change on its sea-terminating outlet remained unknown (Engel and others, Reference Engel, Láska, Nývlt and Stachoň2018).
Lookalike glacier is a valley glacier that is gently inclined (mean slope of 6°) to the NNE. Headwall cliffs and steep slopes bound the upper part of the glacier and its eastern side, respectively. Prominent ice-cored moraines rim the glacier tongue from the north-west (Fig. 2b). The ice thickness of Lookalike Glacier increases towards its central part, where it was measured to be nearly 160 m thick in 2010 (Engel and others, Reference Engel, Nývlt and Láska2012). Over the 2009–15 period, the glacier had a positive cumulative SMB of 0.57 ± 0.67 m w.e. (Engel and others, Reference Engel, Láska, Nývlt and Stachoň2018).
Triangular Glacier is located in the western slope of Lachman Crags plateau at the foot of headwall cliffs that are 120–200 m high (Fig. 2c). This cirque glacier is inclined to the south-west and its surface has a mean slope of 15°. Massive ice-cored moraines delimit the north-western margin of the glacier, and more gentle slopes prevail along the glacier's southern boundary. The glacier, with a maximum measured ice thickness of 103 ± 5 m, experienced a cumulative mass loss of −1.66 ± 0.83 m w.e. over the period 2015–2020 (Engel and others, Reference Engel, Láska, Kavan and Smolíková2023).
Data and methods
Annual mass-balance was measured using the glaciological method between late January and mid-February over the period 2015/16–2020/21 (whenever we refer to the annual mass balance we are specifically referring to the mass-balance year). Observations included the measurement of the height and the position of ablation stakes drilled into the glaciers. The ablation stakes are regularly distributed across the surface of the glaciers except for the crevassed sea-terminating outlet of Davies Dome. Since the mass loss due to frontal ablation of this outlet are not encompassed by the glaciological method, the total mass balance of the entire glacier is less positive or more negative than the calculated SMB estimate. The inter-annual height differences of the stakes were converted to meters of water-equivalent based on variable density for snow above the equilibrium line and glacier ice in the ablation zone (cf. Huss, Reference Huss2013). We converted the values in the ablation and accumulation zones using densities of 900 and 500 ± 90 kg m–3, respectively (Engel and others, Reference Engel, Láska, Nývlt and Stachoň2018). As glaciological measurements are made only once a year during the summer, we use this approach assuming that the net elevation difference between subsequent years corresponds to ice and snow in the ablation and accumulation areas, respectively.
Subsequently, we extrapolated point balances over the entire glacier area using the natural neighbour technique (de Smith and others, Reference de Smith, Goodchild and Longley2015). This technique was selected as the most appropriate interpolation algorithm for our dataset based on cross-validation results (see Engel and others, Reference Engel, Láska, Nývlt and Stachoň2018 for details). Finally, we quantified uncertainties of the calculated balance values following the recommendations of Huss and others (Reference Huss, Bauder and Funk2009). We applied an uncertainty of ±90 kg m–3 to the snowpack density in the accumulation area to address the uncertainty in the SMB determination at the individual stakes (Engel and others, Reference Engel, Láska, Nývlt and Stachoň2018). The uncertainty that results from spatial interpolation was estimated by the std dev. of the cross-validation residuals obtained from the interpolation of the annual mass-balance grids. The uncertainty in the cumulative SMB was calculated as the square root of the sum of squares of the annual errors for the period 2015/16–2020/21.
Based on the annual mass-balance grids, equilibrium-line altitude (ELA) and accumulation area ratio (AAR) values were determined for the individual mass-balance years within the period 2015/16–2020/21. The ice surface elevation was derived from the Reference Elevation Model of Antarctica (Howat and others, Reference Howat, Porter, Smith, Noh and Morin2019) converted to orthometric heights based on the Earth Gravitational Model 2008 1 min geoid height file (Pavlis and others, Reference Pavlis, Holmes, Kenyon and Factor2012). The orthometric (geoid) heights allowed for the comparison of the determined ELAs with the values reported by Engel and others (Reference Engel, Láska, Nývlt and Stachoň2018) for the period 2009/10–2014/15. The balanced-budget AAR was estimated based on the relation between the annual SMB and AAR values (Cogley and others, Reference Cogley2011). Mean values of ELA, SMB and AAR are reported with the standard error of the mean.
To evaluate the volume of Davies Dome including its marine-terminating terminus, we analysed the spatial changes of the complete glacier (volume has only been reported for the ice dome by Engel and others, Reference Engel, Nývlt and Láska2012). The position of the glacier front was delineated using WorldView-2/3 and Pléiades imageries from 2015/16 and 2020, respectively, and the ice surface elevation was recorded from the Reference Elevation Model of Antarctica and Pléiades DEM (see the details in Table 1). GPR was used to determine the ice thickness and the volume of both the dome and outlet. The GPR data were collected using an unshielded 50 MHz antenna and RAMAC CU-II control unit (MALÅ GeoScience, 2005). The ice-thickness data retrieved from two-way travel times assuming a radio-wave velocity of 168 m μs–1 (Narod and Clarke, Reference Narod and Clarke1994) were interpolated over the glacier surface using a minimum curvature method (Engel and others, Reference Engel, Nývlt and Láska2012). The uncertainties associated with ice thickness and glacier volume were determined following the recommendations presented by Lapazaran and others (Reference Lapazaran, Otero, Martín-Español and Navarro2016a, Reference Lapazaran, Otero, Martín-Español and Navarro2016b) and Martín-Español and others (Reference Martín-Español, Lapazaran, Otero and Navarro2016).
22 ground control points were applied to improve the registration of the Reference Elevation Model of Antarctica strip and Pléiades DEMs.
Recent changes in the extent and the ice surface elevation of Davies Dome and Lookalike Glacier were assessed based on the 2006 DEM (Czech Geological Survey, 2009) photogrammetrically derived from British Antarctic Survey (BAS) aerial photographs, the Reference Elevation Model of Antarctica and Pléiades DEM (see details in Table 1). The differences between the ice surface DEMs summed with the volume estimates determined from the GPR survey at Lookalike Glacier (Engel and others, Reference Engel, Nývlt and Láska2012) and Davies Dome (Fig. 3) yield the ice volume changes since 2006. The ice-thickness data derived from the GPR survey at Davies Dome were used to model the glacier bed topography following the procedures described by Engel and others (Reference Engel, Láska, Kavan and Smolíková2023). The bed topography was calculated by subtracting the GPR-derived ice-thickness from the Reference Elevation Model of Antarctica surface elevation.
We used the hourly air temperatures from the Mendel station record, which extends back to the beginning of SMB measurements on James Ross Island, as reference climate data. The air temperature was measured at 2 m above the ground using an EMS33H probe and Minikin TH datalogger (EMS Brno, Czech Republic) with an accuracy of ±0.15°C. The hourly values of air temperature were averaged to produce daily, summer (December-January-February) and annual means for the period 2009/10–2020/21. Furthermore, a PDD approach was applied for the surface melt parameterization (Hock, Reference Hock2003), using hourly air temperature observations. The number of PDD is derived as the integral of the positive air temperature (T) in degrees Celsius (°C) over a selected time interval (A):
Using Eqn (1), we computed the annual sum of PDD for each mass-balance year and compared it with the annual SMB data.
The long-term climate variability in the northern Antarctic Peninsula was assessed based on the monthly mean air temperatures from Esperanza station (location shown in Fig. 1a) reported in the SCAR MET READER database (Turner and others, Reference Turner2004). This dataset was used to calculate the mean annual air temperatures and to determine the linear trends for individual decades in the period 1980–2021.
Results
Climate conditions
The mean annual air temperature record from Esperanza station shows a warming through the period 1980–2021 (Fig. 4a). The highest warming trend occurred during the decades 1980–89 (1.33°C decade−1) and 1990–99 (2.60°C decade−1). The temporary cooling followed during the decade 2000–09 (−1.51°C decade−1) that was the coldest period in the northern Antarctic Peninsula region since the 1980s (Fig. 4b). This cooling culminated at the turn of the 2000s when series of low annual temperatures was recorded. During the subsequent decade 2010–19, air temperature started to increase again (1.23°C decade −1).
The mean annual air temperature at Mendel station in the northern Ulu Peninsula was −6.6°C in the period 2009/10–2020/21. Air temperature has gradually increased since 2009 (Fig. 5), which can be seen by the increase in mean temperatures between the periods 2009/10–2014/15 (−7.4°C) and 2015/16–2020/21 (−6.1°C). The warmest mass-balance years were 2010/11 and 2016/17, with the mean annual air temperature spanning between −5.7 and −3.9°C. The coldest years were 2009/10 and 2012/13, with the annual means of −8.5 and −8.2°C, respectively.
The mean summer (DJF) air temperature varied throughout the entire period, with a minimum of −1.2°C in summer 2013/14 and a maximum of 1.6°C in summer 2020/21. An increase in the mean summer temperatures, which is most relevant to glacier surface melting, was observed between the periods 2009/10–2014/15 (−0.5°C) and 2015/16–2020/21 (0.4°C). There is considerable variability in the annual sum of PDDs, which fluctuates between a minimum of 126.6 K d−1 in 2012/13 and a maximum of 354.1 K d−1 in 2016/17. The annual PDD sum was almost twice higher in this year compared with the coldest years 2009/10 and 2012/13.
The negative correlation between annual SMB and PDD (Fig. 6) indicates the sensitivity of glaciers to air temperature fluctuations. The glaciers respond more sensitively to summer temperatures as indicated by the linear regression analysis between annual SMB and the summer sum of PDD (Fig. 6b). The higher slope of the regression line obtained for summer PDD implies pronounced changes in glacier mass that occur within the warm period of the year and in the low-altitude zones of the glaciers. The highest correlation was found for Davies Dome (r = − 0.965), while a slightly lower correlation was calculated for Triangular Glacier (r = − 0.868). These findings coincide well with the previous studies reporting the high sensitivity of small glaciers on James Ross Island (Davies and others, Reference Davies2014; Engel and others, Reference Engel, Láska, Kavan and Smolíková2023) and in adjacent regions of the northern Antarctic Peninsula (Braun and Hock, Reference Braun and Hock2004; Skvarca and others, Reference Skvarca, De Angelis and Ermolin2004; Jonsell and others, Reference Jonsell, Navarro, Bañón, Lapazaran and Otero2012; Navarro and others, Reference Navarro2023).
Glacier area and volume changes
Lookalike Glacier extended over an area of 2.31 ± 0.02 km2 in 2015, storing 0.2413 ± 0.0136 km3 of ice. The glacier extent in 2015 was 3.8% smaller compared with the ice area in 2006 (Engel and others, Reference Engel, Nývlt and Láska2012), which represents a mean area loss of 0.010 ± 0.008 km2 a−1 (−0.39% a−1). From October 2015 to January 2020, the glacier extent shrank to 2.29 ± 0.02 km2 experiencing an area reduction of 1.0% at a mean rate of −0.005 ± 0.013 km2 a−1 (−0.23% a−1). Over this period, Lookalike Glacier lost 0.096 km3 (0.4%) of ice volume at a mean rate of −0.12% a−1. Its land-terminating glacier front retreated <30 m in total over a period of 4.3 years at a mean retreat rate of 6.8 ± 1.3 m a−1 and its mean ice surface elevation slightly decreased (Table 2).
Davies Dome shrank from 6.49 ± 0.03 km2 in 2006 to 6.27 ± 0.05 km2 in 2016 (−3.4%) and then to 6.24 ± 0.02 km2 in 2020 (−0.6%). Its total area reductions were larger than those of Lookalike Glacier but the mean loss rates for the periods 2006–2016 (−0.32% a−1) and 2016–2020 (−0.18% a−1) were comparable to the values obtained for the Lookalike Glacier. The maximum measured ice thickness of 92.3 ± 8.3 m occurs in the middle part of the marine-terminating outlet (at 231 m a.s.l.) of Davies Dome, while the ice dome on a volcanic plateau has the maximum thickness of 88.0 ± 7.4 m (Fig. 7a). The ice volume of Davies Dome decreased by 0.00404 km3 (1.8%) from November 2016 to January 2020 (−0.56% a−1) mainly due to surface melt. The volume loss due to marine-terminating front retreat accounts for 4.4 ± 0.4 105 m3, or 11% of the total ice volume loss. The calving front retreated by 61.0 ± 2.8 m (19.0 ± 1.6 m a−1) and its mean height decreased by 4.1 ± 1.9 m over a period of 3.2 years. The most prominent calving event since the summer 2017/18 occurred on 12 February 2020 (Fig. 8).
Surface mass balance of glaciers
The observed glaciers experienced a negative annual SMB during most of the 6-year period. The land-terminating Lookalike Glacier had a negative balance in 2015/16, near-zero balance in 2016/17, positive balances in 2017/18–2018/19 and negative balances in the last 2 years (Table 3). The surface mass loss culminated in 2020/21 when the most negative annual balance (−0.70 ± 0.32 m w.e.) since 2009/10 was recorded. This negative value contributed substantially to the cumulative SMB of −0.89 ± 0.86 m w.e. over the years 2015/16–2020/21. The highly negative SMB (−0.89 ± 0.62 m w.e.) in 2020/21 was also determined for the neighbouring Triangular Glacier, where the ELA was located above the highest point of the glacier (329 m). This loss added to the past annual balance values (Engel and others, Reference Engel, Láska, Kavan and Smolíková2023) yields a cumulative SMB of −2.59 ± 0.96 m w.e. over the same 6-year period (the mean ELA was determined to be >255 ± 23 m).
The Davies Dome plateau glacier, which includes a marine-terminating outlet, experienced a net surface mass loss in all years between 2015/16 and 2020/21, except in 2018/19. In the latter year, a slightly positive SMB was determined for both the separate ice dome (0.09 ± 0.16 m w.e.) and the whole glacier including its marine-terminating outlet (0.01 ± 0.30 m w.e.). As the rest of the investigated land-based glaciers on James Ross Island, Davies Dome experienced the largest surface mass loss in 2020/21 for both the separate ice dome (−0.77 ± 0.15 m w.e.) and for the complete glacier with outlet (−0.83 ± 0.19 m w.e.). A surface lowering was recorded at all the ablation stakes in this particular year, which is the first year since the beginning of the mass-balance measurements in 2009 with zero AAR and ELA above the upper glacier limit (Table 3). A cumulative SMB of −1.55 ± 0.60 m w.e. was calculated for Davies Dome (without its outlet) over the period 2015/16–2020/21.
The mean surface mass loss of 0.26 ± 0.11 m w.e. a−1 measured for Davies Dome represents the ice dome only. The marine-terminating outlet decreases the annual SMB values by 0.10–0.06 m w.e. (Table 2) and further losses arise from iceberg calving and submarine melting at the glacier front. Frontal ablation remains unknown but the mean annual volume loss due to ice front retreat from 2015/16 to 2019/20 was 0.14 ± 0.01 million m3, or 0.06 m w.e. Considering this value for the years 2015/16–2020/21, the mass loss of the complete glacier exceeds 0.4 ± 0.2 m w.e. a−1. The volume loss due to frontal retreat accounts for 13–29% of the surface mass loss of Davies Dome and implies the minimum contribution of frontal ablation assuming low ice surface velocities and small number of calving events observed at the glacier terminus. If we consider the similar contribution of both the surface and frontal ablation in the northern Antarctic Peninsula region (Navarro and others, Reference Navarro, Jonsell, Corcuera and Martín-Español2013; Osmanoğlu and others, Reference Osmanoğlu, Navarro, Hock, Braun and Corcuera2014), the total mass loss of Davies Dome with its outlet increases to 0.7 ± 0.2 m w.e. a−1.
A west–east gradient in SMB dominated the annual mass-balance spatial distribution on Lookalike Glacier over most of the investigated period (Fig. 9). A less pronounced asymmetric pattern was observed in 2016/17 and 2019/20, when SMBs were close to zero and the SMB distribution was largely controlled by elevation (Figs 9b, e). The SMB distribution of Davies Dome was more complicated, with increasing mass loss towards the periphery of the dome (i.e. elevation controlled SMB) only observed in the markedly negative mass-balance year of 2020/21 (Fig. 10f). In all the previous years between 2015/16 and 2019/20, the SMB increased along the SW–NE transect attaining positive values on the north-eastern slope of the dome (Figs 10a–e).
Equilibrium-line altitude
The annual ELA determined for Lookalike Glacier ranged between 310 and 378 m except for 2020/21, when enhanced ablation raised it to 449 m (Table 3). The equilibrium line was nearly parallel to contours in 2016/17 only (Fig. 9b). In the rest of the period, it had a complex spatial pattern and crossed several contour lines. The annual values yielded a mean ELA of 362 ± 18 m over the 6-year period. The AAR values ranged from 0.71 to 0.02, with an average of 0.42 ± 0.10. At Davies Dome, the annual ELA ranged from 335 to 438 m before the last year, when the ablation zone expanded across the entire glacier (Fig. 10f). Including the very high ELA of 2020/21, which was above the maximum altitude (519 m) of the glacier, the arithmetic mean of the six annual values yields a value of 427 ± 22 m. AAR was often quite low, with a value of zero in 2020/21, and nearly-zero in 2015/16 and 2019/20. Combined with the remaining annual values, AAR was, on average, 0.24 ± 0.11. The relation between the SMB and AAR, observed over 12 years, indicates the balanced-budget AAR of 0.58 and 0.48 for Lookalike Glacier and Davies Dome, respectively.
The mean ELA over the period 2015/16–2020/21 was 51 and >34 m higher compared with the previous 6-year period at Lookalike Glacier and Davies Dome, respectively (Table 4). The observed shift in ELA follows an increase in the mean summer air temperature of 0.9°C on the glaciers between the 6-year periods. Despite this shift, the current ELA probably remains lower compared with the period 2005–09 when the mean summer air temperature on the Ulu Peninsula was 0.3°C higher (Ambrožová and others, Reference Ambrožová, Láska, Hrbáček, Kavan and Ondruch2019). Furthermore, it equals the regional ELA at the turn of the millennium when the temporary cooling phase initiated (Turner and others, Reference Turner2016). The mean AAR decreased between the periods 2009/10–2014/15 and 2015/16–2020/21 by 0.26 and 0.2 at Lookalike Glacier and Davies Dome, respectively. As a result, the mean AAR values for the recent 6-year period are lower than the balanced-budget AAR indicating that glaciers are out of balance with the current climate conditions. This implies further glacier retreat and volume loss if the current climate conditions in the region persist.
Data for the period 2009/10–2014/15 after Engel and others (Reference Engel, Láska, Nývlt and Stachoň2018).
An extreme surface ablation occurred on the investigated glaciers during the summer of 2020/21. Enhanced melt expanded to the upper part of Lookalike Glacier causing slush swamps and sheet flows below the snow-covered ice divide. Open cryoconite holes formed at an exceptional elevation of 440 m a.s.l., and supraglacial channels (hundreds of meters long and up to 1 m deep) emerged along both margins of the glacier. At the same time, snow cover melted over 100% of Triangular Glacier exposing glacier ice to melt conditions and prompting surface meltwater runoff across the glacier. Finally, surface melt extended across the entire flat top surface of Davies Dome (~500 m a.s.l.), which is normally a dry snow zone. The meltwater was partly retained within the snowpack, but it eventually flowed from the percolation zone forming supraglacial stream networks on the north-eastern slope of the dome and on the marine-terminating outlet.
Discussion
Accelerated glacier mass loss
The retreat rates determined for Lookalike Glacier and Davies Dome indicate regional changes in glacier evolution at the beginning of the 21st century. In 2015/16, the glaciers remained close to their 2006 extent experiencing reductions of 3.8% (−0.4% a−1) and 3.4% (−0.3% a−1) over the period 2006–2015/16. These values coincide with that determined for Triangular Glacier (3.2% and −0.3% a−1), being an order of magnitude lower compared with the reductions reported for Lookalike Glacier (10.6%), Davies Dome (20.7%) and Triangular Glacier (30.8%) over the period 1979–2006 (Engel and others, Reference Engel, Nývlt and Láska2012, Reference Engel, Láska, Kavan and Smolíková2023). The ongoing decrease of area loss (−0.2% a−1) derived for both Lookalike Glacier and Davies Dome over the period 2015/16–2020, and the accelerated retreat of Triangular Glacier (−0.9% a−1), indicate a shift in glacier evolution around the mid-2010s (Fig. 11). This shift is also demonstrated by changes in the rates of glacier volume loss, that nearly doubled between the periods 2006–2015/16 (0.1 to −0.3% a−1) and 2015/16–2020 (−0.1 to −0.6% a−1). The timing of the shift is also constrained with the glaciological measurements, which yield predominantly negative annual SMBs of the observed glaciers since 2015/16. The surface mass loss over the period 2015/16–2020/21 is higher than the surface mass gain over the previous 8-year interval (Table 4).
The observed shift towards mass loss indicates a change in the mass balance of the investigated glaciers. However, uncertainty in the annual SMB values and the short observation period do not allow statistically significant trend to be determined. Nevertheless, the recent onset of enhanced glacier mass loss, directly following a sustained period of increased air temperatures, implies that the previous period of reduced warming and mass losses (Turner and others, Reference Turner2016; Oliva and others, Reference Oliva2017) has terminated. Increased air temperatures since the mid-2010s (Carrasco and others, Reference Carrasco, Bozkurt and Cordero2021) and enhanced temperature-controlled processes such as surface snow melt (Zheng and others, Reference Zheng, Zhou and Liang2019; Zhou and others, Reference Zhou, Liu and Zheng2021) and glacier ablation (Engel and others, Reference Engel, Láska, Kavan and Smolíková2023; WGMS, Reference Zemp, Gärtner-Roer, Nussbaumer, Welty, Dussaillant and Bannwart2023) in the north-eastern Antarctic Peninsula indicate a return to a long-term positive trends in air temperature (e.g. Turner and others, Reference Turner2020), snow melt (Abram and others, Reference Abram2013) and ice mass loss (Carrivick and others, Reference Carrivick, Davies, Glasser, Nývlt and Hambrey2012). However, mass loss of the glaciers on the Ulu Peninsula is driven by summer temperatures (Engel and others, Reference Engel, Láska, Kavan and Smolíková2023) that are related to the highly variable phase of the Southern Annular Mode (Marshall and Thompson, Reference Marshall and Thompson2016) and, therefore, short-term periods of less negative or positive mass balances are still possible in the upcoming years. Ongoing atmospheric warming, and enhanced glacier retreat and thinning observed on James Ross Island since 2020/21 corroborate our hypothesis of a return to the long-term negative mass-balance trend.
The largest losses of ice volume determined for Davies Dome result from surface melt and frontal ablation at the marine-terminating outlet. Higher volume losses of marine-terminating glaciers compared to land-terminating ice masses were observed around the northern Antarctic Peninsula, but ocean-calving glaciers have been considered much more dynamic (Skvarca and others, Reference Skvarca, Rott and Nagler1995; Carrivick and others, Reference Carrivick, Davies, Glasser, Nývlt and Hambrey2012; Cook and others, Reference Cook, Vaughan, Luckman and Murray2014). The limited changes of Davies Dome may be attributed to the location of its accumulation zone at a relatively high elevation (>400 m a.s.l.), short floating section of its outlet (200–700 m) and narrow terminus (325 m or >2% of the glacier perimeter) affected by calving (Fig. 7b). By contrast, marine-terminating outlets on James Ross Island are mostly kilometres long and wide (Lippl and others, Reference Lippl2019) being more vulnerable to reduced sea-ice extent in the Weddell Sea (Parkinson, Reference Parkinson2019) and related changes in air temperature (Turner and others, Reference Turner2020; Carrasco and others, Reference Carrasco, Bozkurt and Cordero2021).
The increased mass loss observed on the Ulu Peninsula since 2015/16 corresponds with the onset of negative annual SMBs determined for Bahía del Diablo Glacier on Vega Island, the north-eastern Antarctic Peninsula (Fig. 12). One year later, the shift to accelerated mass loss was observed on glaciers on the South Shetland Islands north from the Peninsula (WGMS, 2021). The timing of the change in mass-balance trends around the northern Antarctic Peninsula is consistent with the acceleration of glaciers on the eastern and western Antarctic Peninsula (Tuckett and others, Reference Tuckett2019; Wallis and others, Reference Wallis, Hogg, van Wessem, Davison and van den Broeke2023). Increased ice loss since 2017 has been reported for Pine Island and Thwaites glaciers in the Amundsen Sea sector of West Antarctica (Joughin and others, Reference Joughin, Shapero, Smith, Dutrieux and Barham2021; Surawy-Stepney and others, Reference Surawy-Stepney, Hogg, Cornford and Davison2023). Irrespective of the unique climate conditions in Antarctica, recently accelerated global glacier mass loss has been detected in most of the glacierized regions around the world (Hugonnet and others, Reference Hugonnet2021).
Glacier–climate interaction
The observed changes in the glacier extent are consistent with the air temperature record from Esperanza station (Figs 4b, 11). The largest retreat rates reported by Engel and others (Reference Engel, Nývlt and Láska2012, Reference Engel, Láska, Kavan and Smolíková2023) for the period 1979–2006 coincide with the most prominent temperature increase in the Antarctic Peninsula region prior 1999 (e.g. Vaughan and others, Reference Vaughan2003; Turner and others, Reference Turner2005; Jones and others, Reference Jones2019). The subsequent decrease in retreat rates observed over the period 2006–2015/16 resulted from cold conditions during the 2000s. This temporary cooling and its consequences for glacier melting in the northern Antarctic Peninsula were reported by Navarro and others (Reference Navarro, Jonsell, Corcuera and Martín-Español2013) and Oliva and others (Reference Oliva2017).
The cooling phase on Ulu Peninsula culminated between 2012 and 2014 when cold summer conditions, with the mean summer air temperature of <−1.0°C, reduced seasonal ablation and annual surface mass gains. The subsequent temperature increase resulted in considerably warmer conditions during the period 2015/16–2020/21 as indicated by predominantly positive mean summer temperatures and 1.3°C higher mean annual temperature compared with the mean for the previous 6-year period. The observed increase in the annual and summer air temperatures coincides with the acceleration of warming (Carrasco and others, Reference Carrasco, Bozkurt and Cordero2021) and more frequent heat waves across the Antarctic Peninsula (González-Herrero and others, Reference González-Herrero, Barriopedro, Trigo, López-Bustins and Oliva2022). In addition to the warming trend, an increase in rainfall events in the northern part of the Peninsula since 2014 (Carrasco and others, Reference Carrasco and Cordero2020) contributed to glacier surface melting.
Whereas glacier surface melting varies from year to year with summer temperature, snow accumulation (a major positive component of SMB), is a function of precipitation. In situ measurements of precipitation using a disdrometer have been initiated at Mendel station in the summer 2020/21 and numerical modelling remains the only source of information about local precipitation (Matějka and Láska, Reference Matějka and Láska2022). The model simulates small precipitation increase in the northern part of James Ross Island over the period 1979–2014 (1–3 mm w.e. a−2), but without a significant trend (van Wessem and others, Reference van Wessem2016). At regional scale, precipitation records are available from O'Higgins station in the northern part of Antarctic Peninsula and from Frei and Prat stations located on the South Shetland Islands. These station records have revealed a large year to year variability of precipitation amounts, an increase in precipitation around the northern Antarctic Peninsula since 2000, and the largest positive trend of winter precipitation over the last 10 years (Carrasco and Cordero, Reference Carrasco and Cordero2020). The strong variability in precipitation is consistent with modelled results by van Wessem and others (Reference van Wessem2016), but biases remain in the identified trends over the northern Antarctic Peninsula.
The distribution of cumulative SMB on the glaciers (Fig. 13) confirms the significance of snowdrift and redistribution by wind on the glacier mass changes. The wind-driven changes in the distribution of snow were suggested as a reason for the differences in the surface mass gain in the accumulation zones of Davies Dome and Lookalike Glacier (Engel and others, Reference Engel, Láska, Nývlt and Stachoň2018). The removal of snow and its effect on the mass-balance pattern was recently reported for Triangular Glacier (Engel and others, Reference Engel, Láska, Matějka and Nedělčev2022). A shift of the accumulation area on Davies Dome to its north-eastern slope (Fig. 13b) confirms the significance of snowdrift for SMB of small glaciers in the Antarctic Peninsula region (Navarro and others, Reference Navarro, Jonsell, Corcuera and Martín-Español2013). The accumulation zone covered almost the entire dome over the cold period 2009/10–2014/15 but moved to its lee-side surface during the subsequent years when negative SMB prevailed. Despite the north-eastern aspect of this surface being directly exposed to solar radiation, the accumulation of snow is enhanced by south-westerly winds (Kavan and others, Reference Kavan, Nývlt, Láska, Engel and Kňažková2020). As a result, the accumulation zone of Davies Dome was not located on the coldest flat top of the dome but at surfaces with the highest accumulation of snow. The lee-side factor prevails in the SMB distribution on the glacier over the complete observation period 2009/10–2020/21, though the effect of the altitude remains visible as a minor accumulation zone at the highest part of the dome.
Implications of the observed changes
Significant relationships between the annual values of the SMB, ELA and AAR (Fig. 14) indicate that all these variables can be used to analyse the impacts of climate conditions on the investigated glaciers. The plots of SMB versus ELA and AAR reflect higher climate sensitivity of Triangular Glacier that has narrow accumulation zone in its upper section and more simple geometry than Lookalike Glacier and Davies Dome. The lower ELA determined for Triangular Glacier (Table 4) and the higher slope of the regression line obtained for its AAR (Fig. 14b) imply pronounced changes in glacier mass under less extreme climate perturbations. The location of this cirque glacier at a low elevation range suggests that its ELA will rise more frequently in the future above the upper limit of the glacier, implying its earlier disappearance compared with Lookalike Glacier and Davies Dome (Engel and others, Reference Engel, Láska, Kavan and Smolíková2023). Triangular Glacier is the smallest ice mass in the Antarctic periphery subject to annual glaciological measurements and represents a category of glaciers that is most vulnerable to regional atmospheric warming (e.g. Simões and others, Reference Simões, Dani, Bremer, Aquino and Arigony-Neto2004; de Woul and Hock, Reference de Woul and Hock2005; Huss and Fischer, Reference Huss and Fischer2016).
Despite the larger ice volume of Davies Dome and its relatively low rate of mass loss, its gently sloped accumulation area renders this glacier particularly susceptible to possible future rises in ELA. It has been suggested that glaciers with low-lying ablation areas and large flat accumulation areas could lose their accumulation area rapidly under rise in ELA (Davies and others, Reference Davies, Carrivick, Glasser, Hambrey and Smellie2012). The SMB measurements on Davies Dome document the validity of this view as an increase of ~35 m in mean ELA between the periods 2009/10−2014/15 and 2015/16−2020/21 resulted in 90% loss of accumulation area (Figs 13, 15). Further rises in ELA would probably speed-up surface melt considerably through the ice dome thinning and related enhanced reduction of accumulation area. The mass-balance modelling suggests that extreme negative mass-balance rates can be expected on ice cap-like glaciers which cannot retreat to higher elevations (Bolibar and others, Reference Bolibar, Rabatel, Gouttevin, Zekollari and Galiez2022). The lack of topographical feedback from glacier retreat and non-linear response of ice melt and mass balance to increasing air temperatures might be largely relevant on James Ross Archipelago and South Shetland Islands where flat-top glaciers prevail.
Besides the strong vulnerability of ice cap-like glaciers with marine-terminating outlet to rise in ELA, another implication of these glaciers is that mass loss due to frontal ablation is highly variable and sensitive to changes at the ice front (Davies and others, Reference Davies, Carrivick, Glasser, Hambrey and Smellie2012; Dryak and Enderlin, Reference Dryak and Enderlin2020). The rapid retreat of the Davies Dome ice front provides an opportunity to study glacier mass balance during the transition of its calving front to land-terminating terminus. A scarcity of direct observations makes it difficult to predict the timing of this transition. However, persistent retreat of ice front since the 2010s suggests that it could take place within years. Once the terminus recedes out of the sea, a decrease in retreat rate and mass loss can be expected. This change, however, remains poorly known making the future mass balance of ice caps with marine-terminating outlets uncertain. These glaciers prevail in Antarctic periphery and frequently occur in Alaska, Greenland and other Arctic regions.
The observed change in SMB trends in the mid-2010s and subsequent increases in the rate of ice loss confirms that small glaciers are highly sensitive to changes in air temperature. This sensitivity has been demonstrated by recently accelerated retreat of small ice bodies across the glacierized regions including those with similar climate settings (Hugonnet and others, Reference Hugonnet2021). The accelerated rate of ice loss of these glaciers increases mass deficit exceeding the predicted mass changes and the associated sea-level rise (e.g. Edwards and others, Reference Edwards2021; Rounce and others, Reference Rounce2023). The potential contribution of small glaciers to sea-level rise is limited (Bahr and Radić, Reference Bahr and Radić2012; Marzeion and others, Reference Marzeion2020), but a further increase in their mass loss would accelerate their dynamic mass loss. Consequently, this would decrease their longevity and result in a higher than expected rate of sea-level rise (Marzeion and others, Reference Marzeion2020; Edwards and others, Reference Edwards2021). Irrespective of their contribution to sea-level change, small glaciers cover most Antarctic islands and represent unique ecosystem (Anesio and Laybourn-Parry, Reference Anesio and Laybourn-Parry2012). Moreover, melt-water discharge from these glaciers has a significant impact on the development of terrestrial ecosystems (Jung and others, Reference Jung, Jeen, Lee and Lee2022) and deglaciated landscape (Ruiz-Fernández and others, Reference Ruiz-Fernández2019).
Conclusions
Lookalike Glacier and Davies Dome in the northern part of James Ross Island have lost their ice volume at an accelerated rate since 2015. After a period of reduced glacier retreat between 2006 and 2015, the land-terminating Lookalike Glacier lost −0.1% a−1 of its volume from October 2015 to January 2020 and Davies Dome, together with its marine-terminating outlet, shrank in volume by −0.6% a−1 between November 2016 and January 2020. The mean annual volume loss due to calving front retreat over a period of 4.3 years was estimated to be 1.4 ± 0.1 105 m3, or 11% of the total volume loss from Davies Dome.
Lookalike Glacier and Davies Dome experienced a surface mass loss of −0.89 ± 0.86 and −1.55 ± 0.60 m w.e., respectively, losing mass at the mean rate of −0.15 ± 0.13 and −0.26 ± 0.11 m w.e. a−1. Surface mass loss for the period 2015/16–2020/21 exceeds the surface mass gain reported for the previous 6-year period resulting in negative cumulative SMBs over the years 2009/10–2020/21 for both Lookalike Glacier (−0.33 ± 1.09 m w.e.) and Davies Dome (−1.44 ± 0.70 m w.e.). The mean annual SMBs over the period 2015/19–2020/21 were −0.15 ± 0.13 m w.e. for Lookalike Glacier and −0.26 ± 0.11 m w.e. for Davies Dome.
The mean ELAs at Lookalike Glacier (362 ± 18 m a.s.l.) and Davies Dome (>427 ± 22 m a.s.l.) during the period 2015/16–2020/21 were 51 and >34 m higher compared with the period 2009/10–2014/15. Between these 6-year periods, the mean AAR for Lookalike Glacier and Davies Dome decreased from 0.68 ± 0.09 to 0.42 ± 0.10 and from 0.44 ± 0.09 to 0.24 ± 0.11, respectively. The mean AARs are lower than the balanced-budget ratio values, which indicate that the glaciers are out of balance with the current climate conditions.
The transition from the surface mass gain to mass loss around 2014/15 is attributed to the changes in the near surface air temperature trends in the region. The mean annual temperature increased from −7.4°C in the period 2009/10–2014/15 to −6.1°C in the period 2015/16–2020/21. This enhanced rate of mass loss reflects warmer melt-season conditions in the latter period, which had both mean summer air temperature ( + 0.9°C) and PDD sum ( + 33.1 K) higher compared with those of the period 2009/10–2014/15. The high sensitivity of the investigated glaciers to air temperature changes in the Antarctic Peninsula region confirms their vulnerability to the acceleration of regional warming associated with the occurrence of high-temperature extremes and heat wave events over the last few years.
Data availability
The Reference Elevation Model of Antarctica created by the Polar Geospatial Center (PGC) is available from https://doi.org/10.7910/DVN/SAIK8B and mass-balance records for glaciers on the South Shetland Islands and Vega Island from https://doi.org/10.5904/wgms-fog-2021-05. The Pléiades stereo-pair used in this study was provided by the Pléiades Glacier Observatory initiative of the French Space Agency (Centre National d'Études Spatiales, CNES). Other data are available from the author upon request.
Acknowledgements
The research was funded by the Czech Science Foundation (project no. GC20-20240S) and the Ministry of Education, Youth and Sports (projects VAN 2022 and LM2015078). Kind support in the field was highly appreciated from Klára Ambrožová, Filip Hrbáček, Daniel Nývlt, Jakub Ondruch, Peter Váczi and the crew at the J.G. Mendel Station. Chris Stringer is thanked for his suggestions and detailed language check of our manuscript. Finally, we thank Francisco Navarro and an anonymous referee for their thorough reviews and scientific editor Ali Graham for his constructive comments.