Submarine Melting Research Articles

Ocean Modification and Seasonality in a Northern Ellesmere Island Glacial Fjord Prior to Ice Shelf Breakup: Milne Fiord

AbstractUnderstanding the impact of the break‐up of northern Ellesmere Island ice shelves on fjord dynamics is limited by a lack of ocean observations prior to ice loss. Based on profiling and mooring data collected between 2011 and 2019, we describe the oceanography of Milne Fiord prior to the 2020 breakup of the Milne Ice Shelf, including, sources of ambient water, under‐ice hydrography, glacial modification, and seasonality. Ambient waters originate in the Canada Basin but are modified by interaction with the 100 m thick ice shelf at the fjord entrance, and extensive glacier tongue at the head. Properties within the 436 m deep fjord are depth‐dependent, with freshwater surface runoff and subglacial discharge accumulating in the upper 50 m of the fjord each summer behind the ice shelf. Freshwater export is restricted to a basal channel in the ice shelf, resulting in the fjord being more stratified and warmer than waters offshore year‐round. Below the ice shelf and above a 260 m deep sill unrestricted exchange allows warm Atlantic Water to penetrate to the 150 m deep glacier grounding line where submarine melting occurs. Basal meltwater spreads down‐fjord close to its depth of origin due to the strong stratification. Interannual warming of fjord deepwater is likely driven by shoaling of the Arctic Ocean thermocline and spillover at the sill, highlighting the link between fjord properties and regional oceanography. Further breakup of the ice shelf is predicted to substantially alter fjord dynamics, with consequences for ocean forcing of the Milne Glacier.

Glacier–plume or glacier–fjord circulation models? A 2-D comparison for Hansbreen–Hansbukta system, Svalbard

AbstractUp to 30% of the current tidewater mass loss in Svalbard corresponds to frontal ablation through submarine melting and calving. We developed two-dimensional (2-D) glacier–line–plume and glacier–fjord circulation coupled models, both including subglacial discharge, submarine melting and iceberg calving, to simulate Hansbreen–Hansbukta system, SW Svalbard. We ran both models for 20 weeks, throughout April–August 2010, using different scenarios of subglacial discharge and crevasse water depth. Both models showed large seasonal variations of submarine melting in response to transient fjord temperatures and subglacial discharges. Subglacial discharge intensity and crevasse water depth influenced calving rates. Using the best-fit configuration for both parameters our two coupled models predicted observed front positions reasonably well (±10 m). Although the two models showed different melt-undercutting front shapes, which affected the net-stress fields near the glacier front, no significant effects on the simulated glacier front positions were found. Cumulative calving (91 and 94 m) and submarine melting (108 and 118 m) along the simulated period showed in both models (glacier–plume and glacier–fjord) a 1:1.2 ratio of linear frontal ablation between the two mechanisms. Overall, both models performed well on predicting observed front positions when best-fit subglacial discharges were imposed, the glacier–plume model being 50 times computationally faster.

The Impacts of a Subglacial Discharge Plume on Calving, Submarine Melting, and Mélange Mass Loss at Helheim Glacier, South East Greenland

AbstractAlmost half of the Greenland ice sheet's mass loss occurs through iceberg calving at marine terminating glaciers. The presence of buoyant subglacial discharge plumes at these marine termini are thought to increase mass loss both through submarine melting and by undercutting that consequently increases calving rates. Plume models are used to predict submarine melting and undercutting. However, there are few observations that allow these relationships to be tested. Here, we use airborne lidar from the terminus of Helheim Glacier, SE Greenland to measure the bulge induced at the surface by the upwelling plume. We use these measurements to estimate plume discharge rates using a high‐resolution, three‐dimensional plume model. Multiyear observations of the plume are compared to a record of calving from camera and icequake data. We find no evidence to suggest that the presence of a plume, determined by its visibility at the surface, increases the frequency of major calving events and also show that mass loss at the terminus driven directly by plume discharge is significantly less than mass loss from major calving events. The results suggest that the contribution of direct plume‐driven mass loss at deep marine‐terminating glaciers may be less than at shallower termini.

Factors Controlling Terminus Position of Hansbreen, a Tidewater Glacier in Svalbard

AbstractThe mechanism of glacier recession and its climatic controls are complex processes that differ across the Arctic region. Here, we investigate factors influencing front variations of Hansbreen, a glacier terminated in Hornsund fjord (SW Svalbard). We apply remote sensing data to observe glacier front fluctuations between 1992 and 2015 and compare them to atmospheric and oceanographic data, sea water depth at the terminus and surface velocity. Rate of subglacial meltwater discharge approximated by the seasonal positive degree‐day index (PDD) together with sea thermal conditions appear to be the main factors responsible for the fluctuations of the front of Hansbreen, while water depth at the front plays a secondary role. Taking into account ocean and air thermal conditions, the studied period has been divided into warm, cold and moderate years. The glacier retreated considerably throughout a bedrock overdeepening in the very warm period 2012–2014. This recession coincided with a slower ice flow due to intense subglacial runoff and increased submarine melting. The long‐term retreat was interrupted by glacier advances in colder years, regardless of water depth at the front. The slower recession rate was the combined effect of decreased subglacial melting and increased glacier movement associated with lower subglacial runoff. Although the seasonal PDD is a good indicator of the front fluctuations, the duration of the retreat and advance periods are strongly correlated with the sea surface temperature. Expected climate warming and an increase of water temperature in the West Spitsbergen Current will stimulate further recession of Hansbreen in future.

Influence of glacier runoff and near-terminus subglacial hydrology on frontal ablation at a large Greenlandic tidewater glacier

AbstractFrontal ablation from tidewater glaciers is a major component of the total mass loss from the Greenland ice sheet. It remains unclear, however, how changes in atmospheric and oceanic temperatures translate into changes in frontal ablation, in part due to sparse observations at sufficiently high spatial and temporal resolution. We present high-frequency time-lapse imagery (photos every 30 min) of iceberg calving and meltwater plumes at Kangiata Nunaata Sermia (KNS), southwest Greenland, during June–October 2017, alongside satellite-derived ice velocities and modelled subglacial discharge. Early in the melt season, we infer a subglacial hydrological network with multiple outlets that would theoretically distribute discharge and enhance undercutting by submarine melt, an inference supported by our observations of terminus-wide calving during this period. During the melt season, we infer hydraulic evolution to a relatively more channelised subglacial drainage configuration, based on meltwater plume visibility indicating focused emergence of subglacial water; these observations coincide with a reduction in terminus-wide calving and transition to an incised planform terminus geometry. We suggest that temporal variations in subglacial discharge and near-terminus subglacial hydraulic efficiency exert considerable influence on calving and frontal ablation at KNS.

Vertical directionality and spatial coherence of the sound field in glacial bays in Hornsund Fjord.

Arctic glacial bays are among the loudest natural environments in the ocean, owing to heavy submarine melting, calving, freshwater discharge, and ice-wave interactions. Understanding the coherence and vertical directionality of the ambient sound there can provide insights about the mechanisms behind the ice loss in these regions. It can also provide key information for operating technologies such as sonar, communication, and navigation systems. To study the unexplored sound coherence and vertical directionality in glacial bays, a vertical hydrophone array was deployed, and acoustic measurements were made at four glacier termini in Hornsund Fjord, Spitsbergen, in June and July 2019. The measurements show that the sound generated by melting glacier ice is more dominant in the upper portion of the water column near the glacier terminus. The melt water from the submarine melting and the freshwater discharge from the glacier create a glacially modified water duct near the sea surface. This disrupts the inter-sensor vertical coherence in the channel. However, some coherence across the duct is preserved for sound arising from spatially localized events at low frequencies. Overall, the observations in this study can help improve the understanding of the submarine melting phenomenon in glacial bays.

Vertical line array measurements of the sound radiated by melting glaciers in Hornsund Fjord

Marine-terminating glaciers worldwide are melting rapidly in response to climate shifts, resulting in the delivery of freshwater into the oceans. Submarine melting at the glacier-water interface accounts for a significant component of the freshwater delivery from the glacier. This melting produces a distinct acoustic signature, providing a potentially viable modality to monitor glacial ice melting on a large scale via acoustic sensing. In order to evaluate the utility of ambient noise oceanography as a tool to quantify glacial ice-melt, in June 2019 we deployed a vertical hydrophone array and made acoustic measurements at some glacier termini in Hornsund Fjord, Spitsbergen. Quantification via array processing proves to be challenging due to (1) the space- and time-varying sound-speed profile in the underwater channel, and the way it refracts sound in an unknown manner, (2) limited vertical resolution of the array due to its aperture, and (3) interference from other noise sources such as melting icebergs, for example, contaminating the melt sound recordings. We present preliminary results from the processing which reveal different acoustic levels arising from submarine melting at different glaciers. The sound from the melt seems to be more dominant in the upper layers of the water at the glacier-sea interface.

Depth dependence of acoustic signals produced by bubble release events in melting glacier ice

Melting glaciers and ice sheets are important contributors to sea level rise, but the rates at which they are losing mass are difficult to measure precisely. Using the acoustic signals produced by the rapid release of pressurized air from bubbles contained in the ice is a promising avenue for obtaining quantitative estimates of the rate of submarine melting at the termini of tidewater glaciers. The amplitude and character of the observed average signal generated by the bubble release events have been found to depend strongly upon the hydrostatic pressure of the water the bubbles are released into, and therefore on the depth at which said events occur; this dependence must be better understood before the goal of measuring melting rates can be achieved. Data from field experiments performed at the Hornsund fjord in Svalbard, Norway, in which the acoustic melting signal from blocks of glacier ice was recorded at varying depths, will be presented, and compared to the output of an idealized model for the acoustic signal produced by the release of bubbles from melting ice.

Numerical Modeling Shows Increased Fracturing Due to Melt-Undercutting Prior to Major Calving at Bowdoin Glacier

Projections of future ice sheet mass loss and thus sea-level rise rely on the parametrization of iceberg calving in ice sheet models. The interconnection between submarine melt-induced undercutting and calving is still poorly understood, which makes predicted contributions of tidewater glaciers to sea level rise uncertain. Here, we compare detailed 3-D simulations of fracture initiation obtained with the Helsinki Discrete Element Model (HiDEM) to observations, prior to a major calving event at Bowdoin Glacier, Northwest Greenland. Observations of a plume surfacing at the calving location suggest that local melt-undercutting influenced the size of the major calving event. Therefore, several experiments are conducted with various local and distributed (front-wide) undercut geometries. Although the number of undercut experiments is limited by computational requirements, one of the conjectured undercut geometries reproduces the crevasse leading to the observed major calving event in great detail. Our simulations show that undercutting leads to initiation of wider fractures more than 100 m upstream of the terminus, well beyond the directly undercut region. When combining a moderate distributed undercut with local amplified undercuts at the two observed plumes, fracture initiation also increases in between the local undercuts. Thus, our results agree with previous studies suggesting the existence of a 'calving amplifier' effect by submarine melt, both upglacier and across-glacier. Consequently, the simulations show the potentially large impact of submarine melt-induced undercutting on iceberg size.

Surface emergence of glacial plumes determined by fjord stratification

Abstract. Meltwater and sediment-laden plumes at tidewater glaciers, resulting from the localized subglacial discharge of surface melt, influence submarine melting of the glacier and the delivery of nutrients to the fjord's surface waters. It is usually assumed that increased subglacial discharge will promote the surfacing of these plumes. Here, at a western Greenland tidewater glacier, we investigate the counterintuitive observation of a non-surfacing plume in July 2012 (a year of record surface melting) compared to the surfacing of the plume in July 2013 (an average melt year). We combine oceanographic observations, subglacial discharge estimates and an idealized plume model to explain the observed plumes' behavior and evaluate the relative impact of fjord stratification and subglacial discharge on plume dynamics. We find that increased fjord stratification prevented the plume from surfacing in 2012, show that the fjord was more stratified in 2012 due to increased freshwater content and speculate that this arose from an accumulation of ice sheet surface meltwater in the fjord in this record melt year. By developing theoretical scalings, we show that fjord stratification in general exerts a dominant control on plume vertical extent (and thus surface expression), so that studies using plume surface expression as a means of diagnosing variability in glacial processes should account for possible changes in stratification. We introduce the idea that, despite projections of increased surface melting over Greenland, the appearance of plumes at the fjord surface could in the future become less common if the increased freshwater acts to stratify fjords around the Greenland ice sheet. We discuss the implications of our findings for nutrient fluxes, trapping of atmospheric CO2 and the properties of water exported from Greenland's fjords.

How Oceanic Melt Controls Tidewater Glacier Evolution

AbstractThe recent rapid retreat of many Arctic outlet glaciers has been attributed to increased oceanic melt, but the relationship between oceanic melt and iceberg calving remains poorly understood. Here, we employ a transient finite element model that simulates oceanic melt and ice break‐off at the terminus. The response of an idealized tidewater glacier to various submarine melt rates and seasonal variations is investigated. Our modeling shows that for zero to low oceanic melt, the rate of volume loss at the front is similar or higher than that for intermediate oceanic melt rates. Only very high melt rates lead to increasing volume losses. These results highlight the complex interplay between oceanic melt and calving and question the general assumption that increased submarine melt leads to higher calving fluxes and enhanced retreat. Models for tidewater glacier evolution should therefore consider calving and oceanic melt as tightly coupled processes rather than as simple, additive parameterizations.

Analysis of Antarctic Peninsula glacier frontal ablation rates with respect to iceberg melt-inferred variability in ocean conditions

AbstractMarine-terminating glaciers on the Antarctic Peninsula (AP) have retreated, accelerated and thinned in response to climate change in recent decades. Ocean warming has been implicated as a trigger for these changes in glacier dynamics, yet little data exist near glacier termini to assess the role of ocean warming here. We use remotely-sensed iceberg melt rates seaward of two glaciers on the eastern and six glaciers on the western AP from 2013 to 2019 to explore connections between variations in ocean conditions and glacier frontal ablation. We find iceberg melt rates follow regional ocean temperature variations, with the highest melt rates (mean ≈ 10 cm d−1) at Cadman and Widdowson glaciers in the west and the lowest melt rates (mean ≈ 0.5 cm d−1) at Crane Glacier in the east. Near-coincident glacier frontal ablation rates from 2014 to 2018 vary from ~450 m a−1at Edgeworth and Blanchard glaciers to ~3000 m a−1at Seller Glacier, former Wordie Ice Shelf tributary. Variations in iceberg melt rates and glacier frontal ablation rates are significantly positively correlated around the AP (Spearman'sρ= 0.71,p-value = 0.003). We interpret this correlation as support for previous research suggesting submarine melting of glacier termini exerts control on glacier frontal dynamics around the AP.

Twenty-first century ocean forcing of the Greenland ice sheet for modelling of sea level contribution

Abstract. Changes in ocean temperature and salinity are expected to be an important determinant of the Greenland ice sheet's future sea level contribution. Yet, simulating the impact of these changes in continental-scale ice sheet models remains challenging due to the small scale of key physics, such as fjord circulation and plume dynamics, and poor understanding of critical processes, such as calving and submarine melting. Here we present the ocean forcing strategy for Greenland ice sheet models taking part in the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), the primary community effort to provide 21st century sea level projections for the Intergovernmental Panel on Climate Change Sixth Assessment Report. Beginning from global atmosphere–ocean general circulation models, we describe two complementary approaches to provide ocean boundary conditions for Greenland ice sheet models, termed the “retreat” and “submarine melt” implementations. The retreat implementation parameterises glacier retreat as a function of projected subglacial discharge and ocean thermal forcing, is designed to be implementable by all ice sheet models and results in retreat of around 1 and 15 km by 2100 in RCP2.6 and 8.5 scenarios, respectively. The submarine melt implementation provides estimated submarine melting only, leaving the ice sheet model to solve for the resulting calving and glacier retreat and suggests submarine melt rates will change little under RCP2.6 but will approximately triple by 2100 under RCP8.5. Both implementations have necessarily made use of simplifying assumptions and poorly constrained parameterisations and, as such, further research on submarine melting, calving and fjord–shelf exchange should remain a priority. Nevertheless, the presented framework will allow an ensemble of Greenland ice sheet models to be systematically and consistently forced by the ocean for the first time and should result in a significant improvement in projections of the Greenland ice sheet's contribution to future sea level change.

A Decade of Variability on Jakobshavn Isbrae: Ocean Temperatures Pace Speed Through Influence on Mélange Rigidity.

The speed of Greenland's fastest glacier, Jakobshavn Isbrae, has varied substantially since its speedup in the late 1990s. Here we present observations of surface velocity, mélange rigidity, and surface elevation to examine its behaviour over the last decade. Consistent with earlier results, we find a pronounced cycle of summer speedup and thinning followed by winter slowdown and thickening. There were extended periods of rigid mélange in the winters of 2016-17 and 2017-18, concurrent with terminus advances ~6 km farther than in the several winters prior. These terminus advances to shallower depths caused slowdowns, leading to substantial thickening, as has been noted elsewhere. The extended periods of rigid mélange coincide well with a period of cooler waters in Disko Bay. Thus, along with the relative timing of the seasonal slowdown, our results suggest that the ocean's dominant influence on Jakobshavn Isbrae is through its effect on winter mélange rigidity, rather than summer submarine melting. The elevation time series also reveals that in summers when the area upstream of the terminus approaches flotation, large surface depressions can form, which eventually become the detachment points for major calving events. It appears that as elevations near flotation, basal crevasses can form, which initiates a necking process that forms the depressions. The elevation data also show that steep cliffs often evolve into short floating extensions, rather than collapsing catastrophically due to brittle failure. Finally, summer 2019 speeds were slightly faster than the prior two summers, leaving it unclear whether the slowdown is ending.

Meltwater Intrusions Reveal Mechanisms for Rapid Submarine Melt at a Tidewater Glacier

AbstractSubmarine melting has been implicated as a driver of glacier retreat and sea level rise, but to date melting has been difficult to observe and quantify. As a result, melt rates have been estimated from parameterizations that are largely unconstrained by observations, particularly at the near‐vertical termini of tidewater glaciers. With standard coefficients, these melt parameterizations predict that ambient melting (the melt away from subglacial discharge outlets) is negligible compared to discharge‐driven melting for typical tidewater glaciers. Here, we present new data from LeConte Glacier, Alaska, that challenges this paradigm. Using autonomous kayaks, we observe ambient meltwater intrusions that are ubiquitous within 400 m of the terminus, and we provide the first characterization of their properties, structure, and distribution. Our results suggest that ambient melt rates are substantially higher (×100) than standard theory predicts and that ambient melting is a significant part of the total submarine melt flux. We explore modifications to the prevalent melt parameterization to provide a path forward for improved modeling of ocean‐glacier interactions.

Subglacial Plumes

Buoyant plumes form when glacial ice melts directly into the ocean or when subglacial meltwater is discharged to the ocean at depth. They play a key role in regulating heat transport from the ocean to the ice front, and in exporting glacial meltwater to the open ocean. This review summarizes current understanding of the dynamics of these plumes, focusing on theoretical developments and their predictions for submarine melt rates. These predictions are sensitive to ocean temperature, the magnitude and spatial distribution of subglacial discharge, the ambient stratification, and, in the case of sub–ice shelf plumes, the geometry of the ice shelf. However, current understanding relies heavily on parameterizations of melting and entrainment, for which there is little in the way of validation. New observational and experimental constraints are needed to elucidate the structure of the plumes and lend greater confidence to the models.

Tides modulate crevasse opening prior to a major calving event at Bowdoin Glacier, Northwest Greenland

AbstractRetreat of calving glaciers worldwide has contributed substantially to sea-level rise in recent decades. Mass loss by calving contributes significantly to the uncertainty of sea-level rise projections. At Bowdoin Glacier, Northwest Greenland, most calving occurs by a few large events resulting from kilometre-scale fractures forming parallel to the calving front. High-resolution terrestrial radar interferometry data of such an event reveal that crevasse opening is fastest at low tide and accelerates during the final 36 h before calving. Using the ice flow model Elmer/Ice, we identify the crevasse water level as a key driver of modelled opening rates. Sea water-level variations in the range of local tidal amplitude (1 m) can reproduce observed opening rate fluctuations, provided crevasse water level is at least 4 m above the low-tide sea level. The accelerated opening rates within the final 36 h before calving can be modelled by additional meltwater input into the crevasse, enhanced ice cliff undercutting by submarine melt, ice damage increase due to tidal cyclic fatigue, crevasse deepening or a combination of these processes. Our results highlight the influence of surface meltwater and tides on crevasse opening leading to major calving events at grounded tidewater glaciers such as Bowdoin.

Distinct Frontal Ablation Processes Drive Heterogeneous Submarine Terminus Morphology

AbstractCalving and submarine melt drive frontal ablation and sculpt the ice face of marine‐terminating glaciers. However, there are sparse observations of submarine termini, which limit estimates of spatially varying submarine melt. Here we present a detailed survey of a west Greenland glacier to reveal heterogeneity in submarine terminus morphology. We find that the majority of the terminus (~77%) is undercut, driven by calving in the upper water column and submarine melting at depth. The remaining ~23% of the terminus is overcut, driven by calving alone. We use observations of six subglacial discharge outlets, combined with a plume model, to estimate spatially varying discharge fluxes. While small discharge fluxes (<43 m3/s) feed numerous, deeply undercut outlets with subsurface plumes, ~70% of the net subglacial flux emerges at the terminus center, producing a vigorous, surface‐reaching plume. This primary outlet drives large, localized seasonal retreat that exceeds calving rates at secondary outlets.

Sensitivity of Tidewater Glaciers to Submarine Melting Governed by Plume Locations

AbstractThe response of tidewater glaciers to ocean warming remains a key uncertainty in sea level rise predictions. Here we use a 3‐D numerical model to examine the response of an idealized tidewater glacier to spatial variations in submarine melt rate. While melting toward the center of the terminus causes only a localized increase in mass loss, melting near the lateral margins triggers increased calving across the width of the glacier, causing the terminus to retreat at several times the width‐averaged melt rate. This occurs because melting near the margins has a greater disruptive impact on the compressive stress arch that transfers resistance from the side walls to the body of the glacier. We suggest that the rate of terminus advance or retreat may thus be governed by the difference between ice velocity and submarine melting in the slow‐flowing zones away from the glacier center.

Estimating Greenland tidewater glacier retreat driven by submarine melting

Abstract. The effect of the North Atlantic Ocean on the Greenland Ice Sheet through submarine melting of Greenland's tidewater glacier calving fronts is thought to be a key driver of widespread glacier retreat, dynamic mass loss and sea level contribution from the ice sheet. Despite its critical importance, problems of process complexity and scale hinder efforts to represent the influence of submarine melting in ice-sheet-scale models. Here we propose parameterizing tidewater glacier terminus position as a simple linear function of submarine melting, with submarine melting in turn estimated as a function of subglacial discharge and ocean temperature. The relationship is tested, calibrated and validated using datasets of terminus position, subglacial discharge and ocean temperature covering the full ice sheet and surrounding ocean from the period 1960–2018. We demonstrate a statistically significant link between multi-decadal tidewater glacier terminus position change and submarine melting and show that the proposed parameterization has predictive power when considering a population of glaciers. An illustrative 21st century projection is considered, suggesting that tidewater glaciers in Greenland will undergo little further retreat in a low-emission RCP2.6 scenario. In contrast, a high-emission RCP8.5 scenario results in a median retreat of 4.2 km, with a quarter of tidewater glaciers experiencing retreat exceeding 10 km. Our study provides a long-term and ice-sheet-wide assessment of the sensitivity of tidewater glaciers to submarine melting and proposes a practical and empirically validated means of incorporating ocean forcing into models of the Greenland ice sheet.