Marine-terminating Outlet Glaciers Research Articles

Spatial variability of marine-terminating ice sheet retreat in the Puget Lowland

Abstract. Understanding drivers of marine-terminating ice sheet behavior is important for constraining ice contributions to global sea level rise. In part, the stability of marine-terminating ice is influenced by solid Earth conditions at the grounded-ice margin. While the Cordilleran Ice Sheet (CIS) contributed significantly to global mean sea level during its final post-Last-Glacial-Maximum (LGM) collapse, the drivers and patterns of retreat are not well constrained. Coastal outcrops in the deglaciated Puget Lowland of Washington State – largely below sea level during glacial maxima, then uplifted above sea level via glacial isostatic adjustment (GIA) – record the late Pleistocene history of the CIS. The preservation of LGM glacial and post-LGM deglacial sediments provides a unique opportunity to assess the variability in marine ice sheet behavior of the southernmost CIS. Based on paired stratigraphic and geochronological work, with a newly developed marine reservoir correction for this region, we identify that the late-stage CIS experienced stepwise retreat into a marine environment between 15 000 and 14 000 years before present, consistent with timing of marine incursion into the region reported in earlier works. Standstill of marine-terminating ice for at least 500 years, paired with rapid vertical landscape evolution, was followed by continued retreat of ice in a subaerial environment. These results suggest rapid rates of solid Earth uplift and topographic support (e.g., grounding zone wedges) stabilized the ice margin, supporting final subaerial ice retreat. This work leads to a better understanding of shallow-marine and coastal-ice-sheet retreat and is relevant to sectors of the contemporary Antarctic and Greenland ice sheets and marine-terminating outlet glaciers.

Sediment discharge from Greenland’s marine-terminating glaciers is linked with surface melt

Sediment discharged from the Greenland Ice Sheet delivers nutrients to marine ecosystems around Greenland and shapes seafloor habitats. Current estimates of the total sediment flux are constrained by observations from land-terminating glaciers only. Addressing this gap, our study presents a budget derived from observations at 30 marine-margin locations. Analyzing sediment cores from nine glaciated fjords, we assess spatial deposition since 1950. A significant correlation is established between mass accumulation rates, normalized by surface runoff, and distance down-fjord. This enables calculating annual sediment flux at any fjord point based on nearby marine-terminating outlet glacier melt data. Findings reveal a total annual sediment flux of 1.324 + /− 0.79 Gt yr-1 over the period 2010-2020 from all marine-terminating glaciers to the fjords. These estimates are valuable for studies aiming to understand the basal ice sheet conditions and for studies predicting ecosystem changes in Greenland’s fjords and offshore areas as the ice sheet melts and sediment discharge increase.

Ice-ocean interactions at the Northeast Greenland Ice stream (NEGIS) over the past 11,000 years

Recent observations have identified increased mass loss from Greenland marine-terminating outlet glaciers (MTOG) with implications for global sea-level rise and wider ocean circulation. The flow of Atlantic-sourced waters to the Greenland margin is thought to be a major control on MTOG behaviour. Investigation of longer-term records of the role of Atlantic-sourced waters on MOTG dynamics are needed to improve understanding of potential future trends in MTOG behaviour. Here we present a multi-proxy study (benthic and planktic foraminifera, dinoflagellate cysts, diatoms, stable isotopes, sea ice biomarkers and sedimentological analyses) from core PS100-198 on the northeast Greenland shelf to investigate the interaction between the Northeast Greenland Ice Stream (NEGIS) and ocean circulation through the Holocene. Proximal glaciomarine conditions at the base of the core indicate deglaciation before 10.9 ka cal BP with the relatively warm Atlantic Water present through advection of the Return Atlantic Current (RAC) across the shelf. The advection of RAC increased through the early Holocene reaching peak subsurface warmth from 8 to 9 ka cal BP. Surface conditions at this time were characterised by heavy sea-ice cover. During the mid-to late Holocene (c. 7–2 ka cal BP) advection of RAC weakened with cooler subsurface waters, but with an amelioration of surface conditions characterised by seasonal sea ice. From c. 2 ka cal BP, during the late Holocene, surface conditions continued to improve with continued seasonal sea-ice cover while subsurface proxies record an increase in RAC advection. The last c. 100 years represent the most ameliorated surface conditions through the Holocene and with subsurface conditions as warm as the early Holocene peak. This coincided with the final break up of ice within 79N fjord and retreat of NEGIS to the Holocene minimum position. Current conditions, therefore, suggest the present-day ice shelf within 79N fjord is most likely susceptible to collapse in the near future. This study highlights the critical influence of Atlantic-sourced waters on the dynamics of major Greenland MTOGs.

Glacier-specific factors drive differing seasonal and interannual dynamics of Nunatakassaap Sermia and Illullip Sermia, Greenland

ABSTRACT Accelerated ice discharge from marine-terminating outlet glaciers accounted for ~48 percent of ice loss from Greenland between 1992 and 2018, and the northwest has been the largest source of dynamic ice loss. Here, we assess the dynamics of two neighboring northwest Greenland glaciers, Nunatakassaap Sermia (NS) and Illullip Sermia (IS), for 2000 to 2020. Retreat rates at NS far exceeded those at IS, and NS accelerated and thinned substantially following the loss of its ice tongue. This was initially driven by loss of buttressing, followed by feedbacks between thinning, surface slope, and effective pressure, as NS retreated into a deeper and wider section of its fjord. At IS, acceleration and thinning were limited due to its location on a bedrock ridge. At NS, net retreat occurred when seasonal retreat persisted through the winter, whereas at IS it resulted from higher summer retreat rates. The timing of seasonal cycles in frontal positions and ice velocities differed markedly between NS and IS: we suggest that ice velocities respond to seasonal meltwater availability at IS and terminus position variations at NS. Overall, our results demonstrate that the dynamic behavior of NS and IS differs markedly at seasonal to decadal timescales and highlight the importance of glacier-specific factors.

Weekly to monthly terminus variability of Greenland's marine-terminating outlet glaciers

Abstract. Seasonal terminus-position variability of Greenland's marine-terminating outlet glaciers is superimposed on multidecadal trends of glacier retreat. To characterize this seasonal variability, we manually digitized terminus positions for 219 marine-terminating glaciers in Greenland from January 2015 through December 2021 using Sentinel-1 synthetic aperture radar (SAR) mosaics. We digitized at a monthly frequency for 199 glaciers and at a 6 d frequency for 20 glaciers. We found that nearly 80 % of glacier termini in Greenland vary significantly on a seasonal basis. For these seasonally varying glaciers, on average, seasonal retreat typically begins in mid-May, and seasonal advance generally commences in early October. The timing of the initiation of the retreat period may be related to the timing of the onset of ice-sheet surface melt. The rate of retreat events peaks in late summer and reaches a minimum in late winter and early spring. The median magnitude of terminus-position seasonality, the difference between glacier length at the dates of peak advance and retreat, is about 220 m. We find a stronger correlation between this magnitude and glacier velocity than between magnitude and glacier width. Terminus-position seasonality can influence longer-term glacier dynamics and, consequently, ice-sheet mass balance. This study contributes to our understanding of terminus-position seasonality for individual glaciers and collectively for glaciers around the entire Greenland Ice Sheet.

TermPicks: a century of Greenland glacier terminus data for use in scientific and machine learning applications

Abstract. Marine-terminating outlet glacier terminus traces, mapped from satellite and aerial imagery, have been used extensively in understanding how outlet glaciers adjust to climate change variability over a range of timescales. Numerous studies have digitized termini manually, but this process is labor intensive, and no consistent approach exists. A lack of coordination leads to duplication of efforts, particularly for Greenland, which is a major scientific research focus. At the same time, machine learning techniques are rapidly making progress in their ability to automate accurate extraction of glacier termini, with promising developments across a number of optical and synthetic aperture radar (SAR) satellite sensors. These techniques rely on high-quality, manually digitized terminus traces to be used as training data for robust automatic traces. Here we present a database of manually digitized terminus traces for machine learning and scientific applications. These data have been collected, cleaned, assigned with appropriate metadata including image scenes, and compiled so they can be easily accessed by scientists. The TermPicks data set includes 39 060 individual terminus traces for 278 glaciers with a mean of 136 ± 190 and median of 93 of traces per glacier. Across all glaciers, 32 567 dates have been digitized, of which 4467 have traces from more than one author, and there is a duplication rate of 17 %. We find a median error of ∼ 100 m among manually traced termini. Most traces are obtained after 1999, when Landsat 7 was launched. We also provide an overview of an updated version of the Google Earth Engine Digitization Tool (GEEDiT), which has been developed specifically for future manual picking of the Greenland Ice Sheet.

A probabilistic framework for quantifying the role of anthropogenic climate change in marine-terminating glacier retreats

Abstract. Many marine-terminating outlet glaciers have retreated rapidly in recent decades, but these changes have not been formally attributed to anthropogenic climate change. A key challenge for such an attribution assessment is that if glacier termini are sufficiently perturbed from bathymetric highs, ice-dynamic feedbacks can cause rapid retreat even without further climate forcing. In the presence of internal climate variability, attribution thus depends on understanding whether (or how frequently) these rapid retreats could be triggered by climatic noise alone. Our simulations with idealized glaciers show that in a noisy climate, rapid retreat is a stochastic phenomenon. We therefore propose a probabilistic approach to attribution and present a framework for analysis that uses ensembles of many simulations with independent realizations of random climate variability. Synthetic experiments show that century-scale climate trends substantially increase the likelihood of rapid glacier retreat. This effect depends on the timescales over which ice dynamics integrate forcing. For a population of synthetic glaciers with different topographies, we find that external trends increase the number of large retreats triggered within the population, offering a metric for regional attribution. Our analyses suggest that formal attribution studies are tractable and should be further pursued to clarify the human role in recent ice-sheet change. We emphasize that early-industrial-era constraints on glacier and climate state are likely to be crucial for such studies.

Late Weichselian and Holocene behavior of the Greenland Ice Sheet in the Kejser Franz Josef Fjord system, NE Greenland

To better understand the past retreat patterns and chronologies of major marine-terminating outlet glaciers, the Late Weichselian and Holocene glacial history of a NE Greenland fjord system is reconstructed using new and previously published swath bathymetry and high-resolution seismic data, supplemented with multi-proxy analysis of sediment cores. The investigated area extends more than 190 km, including Fosters Bugt, Kejser Franz Josef Fjord, Nordfjord and Moskusoksefjord, the drainage routes for Waltershausen Gletscher, one of the largest outlet glaciers in the NE sector of the Greenland Ice Sheet. A complex of moraine ridges identified in Fosters Bugt show former stillstand positions during this period, and, alongside radiocarbon ages, supports the theory of a maximum ice front position in the outer coastal areas during the Younger Dryas-Preboreal period. The retreat moraines distributed within the Kejser Franz Josef Fjord system provide evidence of a stepwise retreat during the last deglaciation, whereby rapid retreat was interrupted by episodes of stillstands and/or readvances of the grounding line. The fjord system was mostly deglaciated prior to 7.8 ka cal BP, with an estimated average retreat rate ranging between 33 and 96 m a−1 since the Younger Dryas. Following this, only minor ice advances in the inner fjords have been identified. The sediment supply from Waltershausen Gletscher dominated the Kejser Franz Josef Fjord, Nordfjord and outer Moskusoksefjord throughout the Holocene period, whereas in the middle and inner Moskusoksefjord the sediment deposits reflect a more local catchment area. The estimated average retreat rate is in the same order of magnitude as the rates from fjords of northern Fennoscandia, indicating the same overall control for the deglaciation, ice melting from increased Northern Hemisphere summer insolation that peaked in Early Holocene.

Brief communication: Preliminary ICESat-2 (Ice, Cloud and land Elevation Satellite-2) measurements of outlet glaciers reveal heterogeneous patterns of seasonal dynamic thickness change

Abstract. Dynamic changes of marine-terminating outlet glaciers are projected to be responsible for about half of future ice loss from the Greenland Ice Sheet. However, we lack a unified, process-based understanding that can explain the observed dynamic changes of all outlet glaciers. Many glaciers undergo seasonal dynamic thickness changes, and classifying the patterns of seasonal thickness change can improve our understanding of the processes that drive glacier behavior. The Ice, Cloud and land Elevation Satellite-2 (ICESat-2) provides space-based, seasonally repeating altimetry measurements of the ice sheets, allowing us to quantify near-termini seasonal dynamic thickness patterns of 37 outlet glaciers around the Greenland Ice Sheet. We classify the glaciers into seven common patterns of seasonal thickness change over a 2-year period from 2019 to 2020. We find small groupings of neighboring glaciers with similar patterns of seasonal thickness change, but, within larger sectors of the ice sheet, patterns of seasonal thickness change are mostly heterogeneous. Future studies can build upon our results by extending these time series, comparing seasonal dynamic thickness changes with external forcings, such as ocean temperature and meltwater runoff, and with other dynamic variables such as seasonal glacier velocity and terminus position changes.

Multi-decadal retreat of marine-terminating outlet glaciers in northwest and central-west Greenland

Abstract. The retreat and acceleration of marine-terminating outlet glaciers in Greenland over the past 2 decades have been widely attributed to climate change. Here we present a comprehensive annual record of glacier terminus positions in northwest and central-west Greenland and compare it against local and regional climatology to assess the regional sensitivity of glacier termini to different climatic factors. This record is derived from optical and radar satellite imagery and spans 87 marine-terminating outlet glaciers from 1972 through 2021. We find that in this region, most glaciers have retreated over the observation period and widespread regional retreat accelerated from around 1996. The acceleration of glacier retreat coincides with the timing of sharp shifts in ocean surface temperatures, the duration of the sea-ice season, ice-sheet surface mass balance, and meltwater and runoff production. Regression analysis indicates that terminus retreat is most sensitive to increases in runoff and ocean temperatures, while the effect of offshore sea ice is weak. Because runoff and ocean temperatures can influence terminus positions through several mechanisms, our findings suggest that a variety of processes – such as ocean-interface melting, mélange presence and rigidity, and hydrofracture-induced calving – may contribute to, but do not conclusively dominate, the observed regional retreat.

Image classification of marine-terminating outlet glaciers in Greenland using deep learning methods

Abstract. A wealth of research has focused on elucidating the key controls on mass loss from the Greenland and Antarctic ice sheets in response to climate forcing, specifically in relation to the drivers of marine-terminating outlet glacier change. The manual methods traditionally used to monitor change in satellite imagery of marine-terminating outlet glaciers are time-consuming and can be subjective, especially where mélange exists at the terminus. Recent advances in deep learning applied to image processing have created a new frontier in the field of automated delineation of glacier calving fronts. However, there remains a paucity of research on the use of deep learning for pixel-level semantic image classification of outlet glacier environments. Here, we apply and test a two-phase deep learning approach based on a well-established convolutional neural network (CNN) for automated classification of Sentinel-2 satellite imagery. The novel workflow, termed CNN-Supervised Classification (CSC) is adapted to produce multi-class outputs for unseen test imagery of glacial environments containing marine-terminating outlet glaciers in Greenland. Different CNN input parameters and training techniques are tested, with overall F1 scores for resulting classifications reaching up to 94 % for in-sample test data (Helheim Glacier) and 96 % for out-of-sample test data (Jakobshavn Isbrae and Store Glacier), establishing a state of the art in classification of marine-terminating glaciers in Greenland. Predicted calving fronts derived using optimal CSC input parameters have a mean deviation of 56.17 m (5.6 px) and median deviation of 24.7 m (2.5 px) from manually digitised fronts. This demonstrates the transferability and robustness of the deep learning workflow despite complex and seasonally variable imagery. Future research could focus on the integration of deep learning classification workflows with free cloud-based platforms, to efficiently classify imagery and produce datasets for a range of glacial applications without the need for substantial prior experience in coding or deep learning.

Rapidly changing glaciers, ocean and coastal environments, and their impact on human society in the Qaanaaq region, northwestern Greenland

Environments along the coast of Greenland are rapidly changing under the influence of a warming climate in the Arctic. To better understand the changes in the coastal environments, we performed researches in the Qaanaaq region in northwestern Greenland as a part of the ArCS (Arctic Challenge for Sustainability) Project. Mass loss of ice caps and marine-terminating outlet glaciers were quantified by field and satellite observations. Measurements and sampling in fjords revealed the important role of glacial meltwater discharge in marine ecosystems. Flooding of a glacial stream in Qaanaaq and landslides in a nearby settlement were investigated to identify the drivers of the incidents. Our study observed rapid changes in the coastal environments, and their critical impact on the society in Qaanaaq. We organized workshops with the residents to absorb local and indigenous knowledge, as well as to share the results and data obtained in the project. Continuous effort towards obtaining long-term observations requiring involvement of local communities is crucial to contribute to a sustainable future in Greenland.

The future sea-level contribution of the Greenland ice sheet: a multi-model ensemble study of ISMIP6

Abstract. The Greenland ice sheet is one of the largest contributors to global mean sea-level rise today and is expected to continue to lose mass as the Arctic continues to warm. The two predominant mass loss mechanisms are increased surface meltwater run-off and mass loss associated with the retreat of marine-terminating outlet glaciers. In this paper we use a large ensemble of Greenland ice sheet models forced by output from a representative subset of the Coupled Model Intercomparison Project (CMIP5) global climate models to project ice sheet changes and sea-level rise contributions over the 21st century. The simulations are part of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). We estimate the sea-level contribution together with uncertainties due to future climate forcing, ice sheet model formulations and ocean forcing for the two greenhouse gas concentration scenarios RCP8.5 and RCP2.6. The results indicate that the Greenland ice sheet will continue to lose mass in both scenarios until 2100, with contributions of 90±50 and 32±17 mm to sea-level rise for RCP8.5 and RCP2.6, respectively. The largest mass loss is expected from the south-west of Greenland, which is governed by surface mass balance changes, continuing what is already observed today. Because the contributions are calculated against an unforced control experiment, these numbers do not include any committed mass loss, i.e. mass loss that would occur over the coming century if the climate forcing remained constant. Under RCP8.5 forcing, ice sheet model uncertainty explains an ensemble spread of 40 mm, while climate model uncertainty and ocean forcing uncertainty account for a spread of 36 and 19 mm, respectively. Apart from those formally derived uncertainty ranges, the largest gap in our knowledge is about the physical understanding and implementation of the calving process, i.e. the interaction of the ice sheet with the ocean.

The contrasting response of outlet glaciers to interior and ocean forcing

Abstract. The dynamics of marine-terminating outlet glaciers are of fundamental interest in glaciology and affect mass loss from ice sheets in a warming climate. In this study, we analyze the response of outlet glaciers to different sources of climate forcing. We find that outlet glaciers have a characteristically different transient response to surface-mass-balance forcing applied over the interior than to oceanic forcing applied at the grounding line. A recently developed reduced model represents outlet-glacier dynamics via two widely separated response timescales: a fast response associated with grounding-zone dynamics and a slow response of interior ice. The reduced model is shown to emulate the behavior of a more complex numerical model of ice flow. Together, these models demonstrate that ocean forcing first engages the fast, local response and then the slow adjustment of interior ice, whereas surface-mass-balance forcing is dominated by the slow interior adjustment. We also demonstrate the importance of the timescales of stochastic forcing for assessing the natural variability in outlet glaciers, highlighting that decadal persistence in ocean variability can affect the behavior of outlet glaciers on centennial and longer timescales. Finally, we show that these transient responses have important implications for attributing observed glacier changes to natural or anthropogenic influences; the future change already committed by past forcing; and the impact of past climate changes on the preindustrial glacier state, against which current and future anthropogenic influences are assessed.

Sediment controls dynamic behavior of a Cordilleran Ice Stream at the Last Glacial Maximum

The uncertain response of marine terminating outlet glaciers to climate change at time scales beyond short-term observation limits models of future sea level rise. At temperate tidewater margins, abundant subglacial meltwater forms morainal banks (marine shoals) or ice-contact deltas that reduce water depth, stabilizing grounding lines and slowing or reversing glacial retreat. Here we present a radiocarbon-dated record from Integrated Ocean Drilling Program (IODP) Site U1421 that tracks the terminus of the largest Alaskan Cordilleran Ice Sheet outlet glacier during Last Glacial Maximum climate transitions. Sedimentation rates, ice-rafted debris, and microfossil and biogeochemical proxies, show repeated abrupt collapses and slow advances typical of the tidewater glacier cycle observed in modern systems. When global sea level rise exceeded the local rate of bank building, the cycle of readvances stopped leading to irreversible retreat. These results support theory that suggests sediment dynamics can control tidewater terminus position on an open shelf under temperate conditions delaying climate-driven retreat.

Calving event size measurements and statistics of Eqip Sermia, Greenland, from terrestrial radar interferometry

Abstract. Calving is a crucial process for the recently observed dynamic mass loss changes of the Greenland ice sheet. Despite its importance for global sea level change, major limitations in understanding the process of calving remain. This study presents high-resolution calving event data and statistics recorded with a terrestrial radar interferometer at the front of Eqip Sermia, a marine-terminating outlet glacier in Greenland. The derived digital elevation models with a spatial resolution of several metres recorded at 1 min intervals were processed to provide source areas and volumes of 906 individual calving events during a 6 d period. The calving front can be divided into sectors ending in shallow and deep water with different calving statistics and styles. For the shallow sector, characterized by an inclined and very high front, calving events are more frequent and larger than for the vertical ice cliff of the deep sector. We suggest that the calving volume deficiency of 90 % relative to the estimated ice flux in our observations of the deep sector is removed by oceanic melt, subaquatic calving, and small aerial calving events. Assuming a similar ice thickness for both sectors implies that subaqueous mass loss must be substantial for this sector with a contribution of up to 65 % to the frontal mass loss. The size distribution of the shallow sector is represented by a log-normal model, while for the deep sector the log-normal and power-law model fit well, but none of them are significantly better. Variations in calving activity and style between the sectors seem to be controlled by the bed topography and the front geometry. Within the short observation period no simple relationship between environmental forcings and calving frequency or event volume could be detected.

Seasonal ice-speed variations in 10 marine-terminating outlet glaciers along the coast of Prudhoe Land, northwestern Greenland

AbstractWe present a 3-year record of seasonal variations in ice speed and frontal ablation of 10 marine-terminating outlet glaciers along the coast of Prudhoe Land in northwestern Greenland. The glaciers showed seasonal speedup initiated between late May and early June, and terminated between late June and early July. Ice speed subsequently decreased from July to September. The timing of the speedup coincided with the onset of the air temperature rise to above freezing, suggesting an influence of meltwater availability on the glacier dynamics. No clear relationship was found between the speedup and the terminus position or the sea-ice/ice-mélange conditions. These results suggest that the meltwater input to the glacier bed triggered the summer speedup. The excess of summer speed (June–August) over the mean for the rest of the year accounted for 0.5–13% of the annual ice motion. Several glaciers showed seasonal frontal variations, i.e. retreat in summer and advance in winter. This was not due to ice-speed variations, but was driven by seasonal variations in frontal ablation. The results demonstrate the dominant effect of glacier surface melting on the seasonal speedup, and the importance of seasonal speed patterns on longer-term ice motion of marine-terminating outlet glaciers in Greenland.

Update of annual calving front lines for 47 marine terminating outlet glaciers in Greenland (1999–2018)

The Greenland ice sheet has been losing mass in response to increased surface melting (Khan et al. 2015; van den Broeke et al. 2017) as well as discharge of ice from marine terminating outlet glaciers (van den Broeke et al. 2009; Box et al. 2018). Marine terminating outlet glaciers flow to the ocean where they lose mass by e.g. iceberg calving. Currently, the mass loss from the Greenland ice sheet is the largest Arctic contributor to global sea-level rise (van den Broeke et al. 2009, 2017; Box et al. 2018). Therefore, monitoring changes in the Greenland ice sheet is essential to provide policy makers with reliable data.
 There is a consensus that most marine terminating outlet glaciers have retreated in recent decades, and that the increased calving rates are a response to recent atmospheric and oceanic warming (e.g. Box et al. 2018; Moon et al. 2018). The rate of dynamic mass loss is determined by changes of the glacier calving front (i.e. its terminus) position, ice thickness and changes in ice flow. Ocean temperature and fjord circulation also influence the calving front stability by melting the glacier below the water line, thinning the ice that is in contact with water (Moon et al. 2014). Change in calving front position is therefore an important indicator for monitoring the dynamic behaviour of the upstream area of the ice sheet, which is further modulated by local topographic features and buttressing effects (Rignot & Kanagaratnam 2006; Nick et al. 2009).
 The Programme for Monitoring of the Greenland Ice Sheet (PROMICE) is dedicated to monitoring changes in the mass budget of the Greenland ice sheet, including monitoring of the calving front lines of marine terminating outlet glaciers. Here, we present an updated collection of annual measurements of end-of-melt-season calving front lines for 47 marine terminating outlet glaciers in Greenland between 1999 and 2018. We also present an example application of the data set, in which we estimate area changes for this group of glaciers since 1999. The Greenland calving front lines were measured from optical satellite imagery obtained from Landsat, Aster, and Sentinel-2 (Table 1). The PROMICE calving front product is freely available for download as ESRI shapefiles.

Semi-automated open water iceberg detection from Landsat applied to Disko Bay, West Greenland.

Changes in Greenland’s marine-terminating outlet glaciers have led to changes in the flux of icebergs into Greenland’s coastal waters, yet icebergs remain a relatively understudied component of the ice-ocean system. We developed a simple iceberg delineation algorithm for Landsat imagery. A machine learning-based cloud mask incorporated into the algorithm enables us to extract iceberg size distributions from open water even in partially cloudy scenes. We applied the algorithm to the Landsat archive covering Disko Bay, West Greenland, to derive a time series of iceberg size distributions from 2000–02 and 2013–15. The time series captures a change in iceberg size distributions, which we interpret as a result of changes in the calving regime of the parent glacier, Sermeq Kujalleq (Jakobshavn Isbræ). The change in calving style associated with the disintegration and disappearance of Sermeq Kujalleq’s floating ice tongue resulted in the production of more small icebergs. The increased number of small icebergs resulted in increasingly negative power law slopes fit to iceberg size distributions in Disko Bay, suggesting that iceberg size distribution time series provide useful insights into changes in calving dynamics.

Response of the flow dynamics of Bowdoin Glacier, northwestern Greenland, to basal lubrication and tidal forcing

ABSTRACTWe use the full-Stokes model Elmer/Ice to investigate the present dynamics of Bowdoin Glacier, a marine-terminating outlet glacier in northwestern Greenland. Short-term speed variations of the glacier were observed, correlating with air temperature and precipitation, and with the semi-diurnal ocean tides. We use a control inverse method to determine the distribution of basal friction. This reveals that most of the glacier area is characterized by near-plug-flow conditions, while some sticky spots are also identified. We then conduct experiments to test the sensitivity of the glacier flow to basal lubrication and tidal forcing at the calving front. Reduction of the basal drag by 10–40% produces speed-ups that agree approximately with the observed range of speed-ups that result from warm weather and precipitation events. In agreement with the observations, tidal forcing and surface speed near the calving front are found to be in anti-phase (high tide corresponds to low speed, and vice versa). However, the amplitude of the semi-diurnal variability is underpredicted by a factor ~ 3, which is likely related to either inaccuracies in the surface and bedrock topographies or mechanical weakening due to crevassing.