Fifty Years of Tidewater Glacier Surface Elevation and Retreat Dynamics along the South-East Coast of Spitsbergen (Svalbard Archipelago)

Tidewater glaciers on the east coast of Svalbard were examined for surface elevation changes and retreat rate. An archival digital elevation model (DEM) from 1970 (generated from aerial images by the Norwegian Polar Institute) in combination with recent ArcticDEM were used to compare the surface elevation changes of eleven glaciers. This approach was complemented by a retreat rate estimation based on the analysis of Landsat and Sentinel-2 images. In total, four of the 11 tidewater glaciers became land-based due to the retreat of their termini. The remaining tidewater glaciers retreated at an average annual retreat rate of 48 m year−1, and with range between 10–150 m year−1. All the glaciers studied experienced thinning in their frontal zones with maximum surface elevation loss exceeding 100 m in the ablation areas of three glaciers. In contrast to the massive retreat and thinning of the frontal zones, a minor increase in ice thickness was recorded in some accumulation areas of the glaciers, exceeding 10 m on three glaciers. The change in glacier geometry suggests an important shift in glacier dynamics over the last 50 years, which very likely reflects the overall trend of increasing air temperatures. Such changes in glacier geometry are common at surging glaciers in their quiescent phase. Surging was detected on two glaciers studied, and was documented by the glacier front readvance and massive surface thinning in high elevated areas.

Aerodynamic roughness length of crevassed tidewater glaciers from UAV mapping

The aerodynamic roughness length (z0) is an important parameter in the bulk approach for calculating turbulent fluxes and their contribution to ice melt. However, z0 estimates for heavily crevassed tidewater glaciers are rare or only generalised. This study used uncrewed aerial vehicles (UAVs) to map inaccessible tidewater glacier front areas. The high-resolution images were utilised in a structure-from-motion photogrammetry approach to build digital elevation models (DEMs). These DEMs were applied to five models (split across transect and raster methods) to estimate z0 values of the mapped area. The results point out that the range of z0 values across a crevassed glacier is large, by up to 3 orders of magnitude. The division of the mapped area into sub-grids (50 m × 50 m), each producing one z0 value, accounts for the high spatial variability in z0 across the glacier. The z0 estimates from the transect method are in general greater (up to 1 order of magnitude) than the raster method estimates. Furthermore, wind direction (values parallel to the ice flow direction are greater than perpendicular values) and the chosen sub-grid size turned out to have a large impact on the z0 values, again presenting a range of up to 1 order of magnitude each. On average, z0 values between 0.08 and 0.88 m for a down-glacier wind direction were found. The UAV approach proved to be an ideal tool to provide distributed z0 estimates of crevassed glaciers, which can be incorporated by models to improve the prediction of turbulent heat fluxes and ice melt rates.

A fully-coupled 3D model of a large Greenlandic outlet glacier with evolving subglacial hydrology, frontal plume melting and calving

We present the first fully coupled 3D full-Stokes model of a tidewater glacier, incorporating ice flow, subglacial hydrology, plume-induced frontal melting and calving. We apply the model to Store Glacier (Sermeq Kujalleq) in west Greenland to simulate a year of high melt (2012) and one of low melt (2017). In terms of modelled hydrology, we find perennial channels extending 5 km inland from the terminus and up to 41 and 29 km inland in summer 2012 and 2017, respectively. We also report a hydrodynamic feedback that suppresses channel growth under thicker ice inland and allows water to be stored in the distributed system. At the terminus, we find hydrodynamic feedbacks exert a major control on calving through their impact on velocity. We show that 2012 marked a year in which Store Glacier developed a fully channelised drainage system, unlike 2017, where it remained only partially developed. This contrast in modelled behaviour indicates that tidewater glaciers can experience a strong hydrological, as well as oceanic, control, which is consistent with observations showing glaciers switching between types of behaviour. The fully coupled nature of the model allows us to demonstrate the likely lack of any hydrological or ice-dynamic memory at Store Glacier.

A fully-coupled 3D model of a large Greenlandic outlet glacier with evolving subglacial hydrology, frontal plume melting and calving

We present the first fully coupled 3D full-Stokes model of a tidewater glacier, incorporating ice flow, subglacial hydrology, plume-induced frontal melting and calving. We apply the model to Store Glacier (Sermeq Kujalleq) in west Greenland to simulate a year of high melt (2012) and one of low melt (2017). In terms of modelled hydrology, we find perennial channels extending 5 km inland fromthe terminus and up to 41 km and 29 km inland in summer 2012 and 2017, respectively. We also report a hydrodynamic feedback that suppresses channel growth under thicker ice inland and allows water to be stored in the distributed system. At the terminus, we find hydrodynamic feedbacks exert a major control on calving through their impact on velocity. We show that 2012 marked a year inwhich Store Glacier developed a fully channelised drainage system, unlike 2017, where it remained only partially developed. This contrast in modelled behaviour indicates that tidewater glaciers can experience a strong hydrological, as well as oceanic, control, which is consistent with observations showing glaciers switching between seemingly dominant types. The fully coupled nature of the model allows us to demonstrate the likely lack of any hydrological or ice-dynamic memory at Store.

Sensitivity of subglacial drainage to water supply distribution at the Kongsfjord basin, Svalbard

By regulating the amount, the timing, and the location of meltwater supply to the glacier bed, supraglacial hydrology potentially exerts a major control on the evolution of the subglacial drainage system, which in turn modulates ice velocity. Yet the configuration of the supraglacial hydrological system has received only little attention in numerical models of subglacial hydrology so far. Here we apply the two-dimensional subglacial hydrology model GlaDS (Glacier Drainage System model) to a Svalbard glacier basin with the aim of investigating how the spatial distribution of meltwater recharge affects the characteristics of the basal drainage system. We design four experiments with various degrees of complexity in the way that meltwater is delivered to the subglacial drainage model. Our results show significant differences between experiments in the early summer transition from distributed to channelized drainage, with discrete recharge at moulins favouring channelization at higher elevations and driving overall lower water pressures. Otherwise, we find that water input configuration only poorly influences subglacial hydrology, which instead is controlled primarily by subglacial topography. All experiments fail to develop channels of sufficient efficiency to substantially reduce summertime water pressures, which we attribute to small surface gradients and short melt seasons. The findings of our study are potentially applicable to most Svalbard tidewater glaciers with similar topography and low meltwater recharge. The absence of efficient channelization implies that the dynamics of tidewater glaciers in the Svalbard archipelago may be sensitive to future long-term trends in meltwater supply.

Aerodynamic roughness length of crevassed tidewater glaciers from UAV mapping

The aerodynamic roughness length (z0) is an important parameter in the bulk approach for calculating turbulent fluxes and their contribution to ice melt. However, for heavily crevassed tidewater glaciers z0 estimations are rare or only generalized. This study used unmanned aerial vehicles (UAVs) to map inaccessible tidewater glacier front areas. The high-resolution images were used in a structure-from-motion photogrammetry approach to build digital elevation models (DEMs). These DEMs were applied to five different models (split across transect and raster methods) to estimate z0 values of the mapped area. The results point out that the range of z0 values across a glacier is large, with up to three (locally even four) orders of magnitude. The division of the mapped area into sub-grids (50 m x 50 m), each producing one z0 value, best accounts for the high spatial variability of z0 across the glacier. The z0 estimations from the transect method are in general higher (up to one order of magnitude) than the raster method estimations. Furthermore, wind direction (values parallel to the ice flow direction are larger than perpendicular) and the chosen sub-grid size turned out to have a large impact on the z0 values, again presenting a range of up to one order of magnitude each. On average, z0 values between 0.08 m and 0.88 m for a down-glacier wind direction were found. The UAV approach proved to be an ideal tool to provide distributed z0 estimations of crevassed glaciers, which can be incorporated by models to improve the prediction of turbulent heat fluxes and ice melt rates.

Harbor Seals as Sentinels of Ice Dynamics in Tidewater Glacier Fjords

Tidewater glaciers calve icebergs into the marine environment which serve as pupping, molting, and resting habitat for some of the largest seasonal aggregations of harbor seals (Phoca vitulina richardii) in the world. Although they are naturally dynamic, advancing and retreating in response to local climatic and fjord conditions, most tidewater glaciers around the world are thinning and retreating. Climate change models predict continued loss of land-based ice with unknown impacts to organisms such as harbor seals that rely on glacier ice as habitat for critical life history events. To understand the impacts of changing ice availability on harbor seals, we quantified seasonal and annual changes in ice habitat in Johns Hopkins Inlet, a tidewater glacier fjord in Glacier Bay National Park in southeastern Alaska. We conducted systematic aerial photographic surveys (n = 55) of seals and ice during the pupping (June; n = 30) and molting (August; n = 25) periods from 2007 to 2014. Object-based image analysis was used to quantify the availability and spatial distribution of floating ice in the fjord. Multivariate spatial models were developed for jointly modeling stage-structured seal location data and ice habitat. Across all years, there was consistently more ice in the fjord during the pupping season in June than during the molting season in August, which was likely driven by seasonal variation in physical processes that influence the calving dynamics of tidewater glaciers. Non-pup harbor seals and ice were correlated during the pupping season, but this correlation was reduced during the molting season suggesting that harbor seals may respond to changes in habitat differently depending upon trade-offs associated with life history events, such as pupping and molting, and energetic costs and constraints associated with the events.

Complex multi-decadal ice dynamical change inland of marine-terminating glaciers on the Greenland Ice Sheet

Greenland's future contribution to sea-level rise is strongly dependent on the extent to which dynamic perturbations, originating at the margin, can drive increased ice flow within the ice-sheet interior. However, reported observations of ice dynamical change at distances >~50 km from the margin have a very low spatial and temporal resolution. Consequently, the likely response of the ice-sheet's interior to future oceanic and atmospheric warming is poorly constrained. Through combining GPS and satellite-image-derived ice velocity measurements, we measure multi-decadal (1993–1997 to 2014–2018) velocity change at 45 inland sites, encompassing all regions of the ice sheet. We observe an almost ubiquitous acceleration inland of tidewater glaciers in west Greenland, consistent with acceleration and retreat at glacier termini, suggesting that terminus perturbations have propagated considerable distances (>100 km) inland. In contrast, outside of Kangerlussuaq, we observe no acceleration inland of tidewater glaciers in east Greenland despite terminus retreat and near-terminus acceleration, and suggest propagation may be limited by the influence of basal topography and ice geometry. This pattern of inland dynamical change indicates that Greenland's future contribution to sea-level will be spatially complex and will depend on the capacity for dynamic changes at individual outlet glacier termini to propagate inland.

High-concentration sediment plumes, Horseshoe Island, western Antarctic Peninsula

The variability in sediment concentration and spatial distribution of meltwater discharges from tidewater glaciers can be used to elucidate climatic evolution and glacier behaviour due to the association between sediment yield and glacier retreat (e.g. Domack & McClennen 1996). In an accelerated deglaciation environment, higher sediment concentrations in the water column can change the glacimarine costal dynamics and affect productivity and sea floor ecosystems (e.g. Marín et al. 2013). In the Antarctic Peninsula Region, meltwater or turbid plumes were previously believed to be rare or without an important role in the sedimentary glacimarine environment (e.g. Griffith & Anderson 1989), but recent studies have shown that this is a common phenomenon in subpolar and transition polar climates (Yoo et al. 2015, Rodrigo et al. 2016). In the current climate change scenario, accelerated glacier retreats and mass losses can produce an increasing input of glacial meltwater into the fjord regions, a situation that is not yet well evaluated in the Antarctic Peninsula. In this short note, after in situ observation of an unusual waterfall from the southern side of the main western tidewater glacier (Shoesmith Glacier) of Horseshoe Island (Lystad Bay), Marguerite Bay (Fig. 1), we report high turbidity values associated with plumes from the glacier, whose values were higher than reported data from subpolar/transition polar Antarctic climates.

Using sub-daily timelapse imagery to investigate the behaviour of Narsap Sermia, SW Greenland.

<p>Greenland’s tidewater glaciers (TWG) have been retreating since the mid-1990s, contributing to mass loss from the Greenland Ice Sheet and sea level rise. Satellite imagery has been widely used to investigate TWG behaviour and determine the response of TWGs to climate. However, multi-day revisit times make it difficult to determine short-term processes such as calving and shorter-term velocity changes that may condition this. </p><p>Here we present velocity, calving and proglacial plume data derived from hourly time-lapse images of Narsap Sermia, SW Greenland for the period July 2017 to June 2020 (n=13,513). Raw images were orthorectified using the <em>Image GeoRectification And Feature Tracking toolbox</em> (ImGRAFT; Messerli & Grinsted, 2015) using a smoothed ArcticDEM tile from 2016 (RMSE=44.4px). TWG flow velocities were determined using ImGRAFT feature tracking, with post-processing adjusting for varying time intervals between image acquisitions (if >1 hour) and removing outliers (>x2 mean). The high temporal resolution of the imagery also enabled the manual mapping of proglacial plume sizes from the orthorectified images and the recording of individual calving events by visually comparing images.</p><p>Results show a total retreat of approximately 700 m, with a general velocity increase from ~15 m/d to ~20 m/d over the investigated time period and highly variable hourly velocities (±12m/d). The number of calving events and plume sizes remain relatively stable from year to year throughout the observation period. However, later in the record plumes appear earlier in the year and the size of calved icebergs increases significantly, which suggests a change in calving behaviour. </p>