Submarine Melt Rates Research Articles

Improved Parameterizations of Vertical Ice‐Ocean Boundary Layers and Melt Rates

AbstractBuoyancy fluxes and submarine melt rates at vertical ice‐ocean interfaces are commonly parameterized using theories derived for unbounded free plumes. A Large Eddy Simulation is used to analyze the disparate dynamics of free plumes and wall‐bounded plumes; the distinctions between the two are supported by recent theoretical and experimental results. Modifications to parameterizations consistent with these simulations are tested and compared to results from numerical and laboratory experiments of meltwater plumes. These modifications include 50% weaker entrainment and a distinct plume‐driven friction velocity in the shear boundary layer up to 8 times greater than the externally‐driven friction velocity. Using these updated plume parameter modifications leads to 40 times the ambient melt rate predicted by commonly used parameterizations at vertical glacier faces, which is consistent with observed melt rates at LeConte Glacier, Alaska.

Modelled frontal ablation and velocities at Kronebreen, Svalbard, are sensitive to the choice of submarine melt rate scenario

Abstract Both submarine melt and calving are important for the overall mass balance of marine-terminating glaciers, but uncertainty is rife with regards to the magnitude of the processes. Modelling allows for these processes to be investigated without the need to visit inaccessible ice marginal zones. This study looks at the impact of different submarine melt and sea-ice back pressure scenarios on modelled calving activity and dynamics at Kronebreen, Svalbard, by running separate summer and winter simulations with various submarine melt parameterisations and sea-ice characteristics. It is found that submarine melt is an important driver of seasonal variation in modelled glacier dynamics and calving activity, with the choice of sliding law also exerting a significant influence on results.

‘Stable’ and ‘unstable’ are not useful descriptions of marine ice sheets in the Earth's climate system

Investigations of the time-dependent behavior of marine ice sheets and their sensitivity to basal conditions require numerical models because existing theoretical analyses focus only on steady-state configurations primarily with a power-law basal shear stress. Numerical results indicate that the choice of the sliding law strongly affects ice-sheet dynamic behavior. Although observed or simulated grounding-line retreat is typically interpreted as an indication of marine ice sheet instability introduced by Weertman (1974), this (in)stability is a characteristic of the ice sheet's steady states – not time-variant behavior. To bridge the gap between theoretical and numerical results, we develop a framework to investigate grounding line dynamics with generalized basal and lateral stresses (i.e. the functional dependencies are not specified). Motivated by observations of internal variability of the Southern Ocean conditions we explore the grounding-line response to stochastic variability. We find that adding stochastic variability to submarine melt rates that produced stable steady-state configurations leads to intermittently advancing and retreating grounding lines. They can also retreat in an unstoppable manner on time-scales significantly longer than the stochastic correlation time-scales. These results suggest that at any given time of their evolution, the transient behavior of marine ice sheets cannot be described in terms of ‘stable’ or ‘unstable’.

Spatial variation in acoustic field due to submarine melting in glacial bays

Climate-change induced melting is leading to accelerating ice loss at tidewater glaciers worldwide. A significant component of the freshwater flux from these glaciers arises from submarine melting at the glacier-ocean interface. This melting causes a distinct acoustic signature due to the release of pressurized bubbles underwater, opening up the possibility of monitoring this phenomenon on a large scale using passive acoustic systems. To evaluate the use of sound in monitoring submarine glacier melting, we made acoustic measurements using vertical hydrophone arrays in four glacial bays in Svalbard in 2019. As the recording array was moved away from the glacier, the variation in the recorded acoustic level due to melting ice did not follow a uniformly decreasing trend as one might expect. Moreover, the acoustic intensities at different glaciers were clustered at different levels. These observations indicate that the geometry of the glacier-ocean interface, thermohaline structure of the underwater channel and presence of floating ice in the bay played a role in determining the acoustic field. Most of this variation can be explained through propagation modeling. Moving forward, this model-based interpretation of the field will play an integral part in inverting the sound to estimate the submarine melt rate.

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.

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.

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.

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.

Direct observations of submarine melt and subsurface geometry at a tidewater glacier.

Ice loss from the world's glaciers and ice sheets contributes to sea level rise, influences ocean circulation, and affects ecosystem productivity. Ongoing changes in glaciers and ice sheets are driven by submarine melting and iceberg calving from tidewater glacier margins. However, predictions of glacier change largely rest on unconstrained theory for submarine melting. Here, we use repeat multibeam sonar surveys to image a subsurface tidewater glacier face and document a time-variable, three-dimensional geometry linked to melting and calving patterns. Submarine melt rates are high across the entire ice face over both seasons surveyed and increase from spring to summer. The observed melt rates are up to two orders of magnitude greater than predicted by theory, challenging current simulations of ice loss from tidewater glaciers.

Sensitivity of a calving glacier to ice–ocean interactions under climate change: new insights from a 3-D full-Stokes model

Abstract. Iceberg calving accounts for between 30 % and 60 % of net mass loss from the Greenland Ice Sheet, which has intensified and is now the single largest contributor to global sea level rise in the cryosphere. Changes to calving rates and the dynamics of calving glaciers represent a significant uncertainty in projections of future sea level rise. A growing body of observational evidence suggests that calving glaciers respond rapidly to regional environmental change, but predictive capacity is limited by the lack of suitable models capable of simulating calving mechanisms realistically. Here, we use a 3-D full-Stokes calving model to investigate the environmental sensitivity of Store Glacier, a large outlet glacier in West Greenland. We focus on two environmental processes: undercutting by submarine melting and buttressing by ice mélange, and our results indicate that Store Glacier is likely to be able to withstand moderate warming perturbations in which the former is increased by 50 % and the latter reduced by 50 %. However, severe perturbation with a doubling of submarine melt rates or a complete loss of ice mélange destabilises the calving front in our model runs. Furthermore, our analysis reveals that stress and fracture patterns at Store's terminus are complex and varied, primarily due to the influence of basal topography. Calving style and environmental sensitivity vary greatly, with propagation of surface crevasses significantly influencing iceberg production in the northern side, whereas basal crevasses dominate in the south. Any future retreat is likely to be initiated in the southern side by a combination of increased submarine melt rates in summer and reduced mélange strength in winter. The lateral variability, as well as the importance of rotational and bending forces at the terminus, underlines the importance of using the 3-D full-Stokes stress solution when modelling Greenland's calving glaciers.

Spatio-temporal variations in seasonal ice tongue submarine melt rate at a tidewater glacier in southwest Greenland

ABSTRACTSubmarine melting of tidewater glaciers is proposed as a trigger for their recent thinning, acceleration and retreat. We estimate spring submarine melt rates (SMRs) of Kangiata Nunaata Sermia in southwest Greenland, from 2012 to 2014, by examining changes in along-fjord freeboard and velocity of the seasonal floating ice tongue. Estimated SMRs vary spatially and temporally near the grounding line, with mean rates of 1.3 ± 0.6, 0.8 ± 0.3 and 1.0 ± 0.4 m d−1across the tongue in 2012, 2013 and 2014, respectively. Higher melt rates correspond with locations of emerging subglacial plumes and terminus calving activity observed during the melt season using time-lapse camera imagery. Modelling of subglacial flow paths suggests a dynamic system capable of rapid re-routing of subglacial discharge both within and between melt seasons. Our results provide an empirically-derived link between the presence of subglacial discharge plumes and areas of high spring submarine melting and calving along glacier termini.

Calving controlled by melt-under-cutting: detailed calving styles revealed through time-lapse observations

AbstractWe present a highly detailed study of calving dynamics at Tunabreen, a tidewater glacier in Svalbard. A time-lapse camera was trained on the terminus and programmed to capture images every 3 seconds over a 28-hour period in August 2015, producing a highly detailed record of 34 117 images from which 358 individual calving events were distinguished. Calving activity is characterised by frequent events (12.8 events h−1) that are small relative to the spectrum of calving events observed, demonstrating the prevalence of small-scale calving mechanisms. Five calving styles were observed, with a high proportion of calving events (82%) originating at, or above, the waterline. The tidal cycle plays a key role in the timing of calving events, with 68% occurring on the falling limb of the tide. Calving activity is concentrated where meltwater plumes surface at the glacier front, and a ~ 5 m undercut at the base of the glacier suggests that meltwater plumes encourage melt-under-cutting. We conclude that frontal ablation at Tunabreen may be paced by submarine melt rates, as suggested from similar observations at glaciers in Svalbard and Alaska. Using submarine melt rate to calculate frontal ablation would greatly simplify estimations of tidewater glacier losses in prognostic models.

A two-dimensional glacier–fjord coupled model applied to estimate submarine melt rates and front position changes of Hansbreen, Svalbard

ABSTRACTWe have developed a two-dimensional coupled glacier–fjord model, which runs automatically using Elmer/Ice and MITgcm software packages, to investigate the magnitude of submarine melting along a vertical glacier front and its potential influence on glacier calving and front position changes. We apply this model to simulate the Hansbreen glacier–Hansbukta proglacial–fjord system, Southwestern Svalbard, during the summer of 2010. The limited size of this system allows us to resolve some of the small-scale processes occurring at the ice–ocean interface in the fjord model, using a 0.5 s time step and a 1 m grid resolution near the glacier front. We use a rich set of field data spanning the period April–August 2010 to constrain, calibrate and validate the model. We adjust circulation patterns in the fjord by tuning subglacial discharge inputs that best match observed temperature while maintaining a compromise with observed salinity, suggesting a convectively driven circulation in Hansbukta. The results of our model simulations suggest that both submarine melting and crevasse hydrofracturing exert important controls on seasonal frontal ablation, with submarine melting alone not being sufficient for reproducing the observed patterns of seasonal retreat. Both submarine melt and calving rates accumulated along the entire simulation period are of the same order of magnitude, ~100 m. The model results also indicate that changes in submarine melting lag meltwater production by 4–5 weeks, which suggests that it may take up to a month for meltwater to traverse the englacial and subglacial drainage network.

Greenland iceberg melt variability from high-resolution satellite observations

Abstract. Iceberg discharge from the Greenland Ice Sheet accounts for up to half of the freshwater flux to surrounding fjords and ocean basins, yet the spatial distribution of iceberg meltwater fluxes is poorly understood. One of the primary limitations for mapping iceberg meltwater fluxes, and changes over time, is the dearth of iceberg submarine melt rate estimates. Here we use a remote sensing approach to estimate submarine melt rates during 2011–2016 for 637 icebergs discharged from seven marine-terminating glaciers fringing the Greenland Ice Sheet. We find that spatial variations in iceberg melt rates generally follow expected patterns based on hydrographic observations, including a decrease in melt rate with latitude and an increase in melt rate with iceberg draft. However, we find no longitudinal variations in melt rates within individual fjords. We do not resolve coherent seasonal to interannual patterns in melt rates across all study sites, though we attribute a 4-fold melt rate increase from March to April 2011 near Jakobshavn Isbræ to fjord circulation changes induced by the seasonal onset of iceberg calving. Overall, our results suggest that remotely sensed iceberg melt rates can be used to characterize spatial and temporal variations in oceanic forcing near often inaccessible marine-terminating glaciers.

Estimating Spring Terminus Submarine Melt Rates at a Greenlandic Tidewater Glacier Using Satellite Imagery

Oceanic forcing of the Greenland Ice Sheet is believed to promote widespread thinning at tidewater glaciers, with submarine melting proposed as a potential trigger of increased glacier calving, retreat, and subsequent acceleration. The precise mechanism(s) driving glacier instability, however, remain poorly understood, and while increasing evidence points to the importance of submarine melting, estimates of melt rates are uncertain. Here we estimate submarine melt rate by examining freeboard changes in the seasonal ice tongue of Kangiata Nunaata Sermia (KNS) at the head of Kangersuneq Fjord (KF), southwest Greenland. We calculate melt rates for March and May 2013 by differencing along-fjord surface elevation, derived from high-resolution TanDEM-X digital elevation models (DEMs), in combination with ice velocities derived from offset tracking applied to TerraSAR-X imagery. Estimated steady state melt rates reach up to 1.4 ± 0.5 m d−1 near the glacier grounding line, with mean values of up to 0.8 ± 0.3 and 0.7 ± 0.3 m d−1 for the eastern and western parts of the ice tongue, respectively. Melt rates decrease with distance from the ice front and vary across the fjord. This methodology reveals spatio-temporal variations in submarine melt rates (SMRs) at tidewater glaciers which develop floating termini, and can be used to improve our understanding of ice-ocean interactions and submarine melting in glacial fjords.

Satellite-derived submarine melt rates and mass balance (2011–2015) for Greenland's largest remaining ice tongues

Abstract. Ice-shelf-like floating extensions at the termini of Greenland glaciers are undergoing rapid changes with potential implications for the stability of upstream glaciers and the ice sheet as a whole. While submarine melting is recognized as a major contributor to mass loss, the spatial distribution of submarine melting and its contribution to the total mass balance of these floating extensions is incompletely known and understood. Here, we use high-resolution WorldView satellite imagery collected between 2011 and 2015 to infer the magnitude and spatial variability of melt rates under Greenland's largest remaining ice tongues – Nioghalvfjerdsbræ (79 North Glacier, 79N), Ryder Glacier (RG), and Petermann Glacier (PG). Submarine melt rates under the ice tongues vary considerably, exceeding 50 m a−1 near the grounding zone and decaying rapidly downstream. Channels, likely originating from upstream subglacial channels, give rise to large melt variations across the ice tongues. We compare the total melt rates to the influx of ice to the ice tongue to assess their contribution to the current mass balance. At Petermann Glacier and Ryder Glacier, we find that the combined submarine and aerial melt approximately balances the ice flux from the grounded ice sheet. At Nioghalvfjerdsbræ the total melt flux (14.2 ± 0.96 km3 a−1 w.e., water equivalent) exceeds the inflow of ice (10.2 ± 0.59 km3 a−1 w.e.), indicating present thinning of the ice tongue.

Nonlinear response of iceberg side melting to ocean currents

AbstractIcebergs calving into Greenlandic Fjords frequently experience strongly sheared flows over their draft, but the impact of this flow past the iceberg is not fully captured by existing parameterizations. We present a series of novel laboratory experiments to determine the dependence of submarine melting along iceberg sides on a background flow. We show, for the first time, that two distinct regimes of melting exist depending on the flow magnitude and consequent behavior of melt plumes (side‐attached or side‐detached), with correspondingly different meltwater spreading characteristics. When this velocity dependence is included in melt parameterizations, melt rates estimated for observed icebergs in the attached regime increase, consistent with observed iceberg submarine melt rates. We show that both attached and detached plume regimes are relevant to icebergs observed in a Greenland fjord. Further, depending on the regime, iceberg meltwater may either be confined to a surface layer or distributed over the iceberg draft.

A model for tidewater glacier undercutting by submarine melting

AbstractDynamic change at the marine‐terminating margins of the Greenland Ice Sheet may be initiated by the ocean, particularly where subglacial runoff drives vigorous ice‐marginal plumes and rapid submarine melting. Here we model submarine melt‐driven undercutting of tidewater glacier termini, simulating a process which is key to understanding ice‐ocean coupling. Where runoff emerges from broad subglacial channels we find that undercutting has only a weak impact on local submarine melt rate but increases total ablation by submarine melting due to the larger submerged ice surface area. Thus, the impact of melting is determined not only by the melt rate magnitude but also by the slope of the ice‐ocean interface. We suggest that the most severe undercutting occurs at the maximum height in the fjord reached by the plume, likely promoting calving of ice above. It remains unclear, however, whether undercutting proceeds sufficiently rapidly to influence calving at Greenland's fastest‐flowing glaciers.

Impact of a Localized Source of Subglacial Discharge on the Heat Flux and Submarine Melting of a Tidewater Glacier: A Laboratory Study

AbstractIdealized laboratory experiments have been conducted in a two-layer stratified fluid to investigate the leading-order dynamics that control submarine melting and meltwater export near a vertical ice–ocean interface as a function of subglacial discharge. In summer, the discharge of surface runoff at the base of a glacier (subglacial discharge) generates strong buoyant plumes that rise along the glacier front entraining ambient water along the way. The entrainment enhances the heat transport toward the glacier front and hence the submarine melt rate increases with the subglacial discharge rate. In the laboratory, the effect of subglacial discharge is simulated by introducing freshwater at freezing temperature from a point source at the base of an ice block representing the glacier. The circulation pattern observed both with and without subglacial discharge resembles those observed in previous observational and numerical studies. Buoyant plumes rise vertically until they find either their neutrally buoyant level or the free surface. Hence, the meltwater can deposit within the interior of the water column and not entirely at the free surface, as confirmed by field observations. The heat budget in the tank, calculated following a new framework, gives estimates of submarine melt rate that increase with the subglacial discharge and are in agreement with the directly measured submarine melting. This laboratory study provides the first direct measurements of submarine melt rates for different subglacial discharges, and the results are consistent with the predictions of previous theoretical and numerical studies.

Heat, Salt, and Freshwater Budgets for a Glacial Fjord in Greenland

AbstractIn Greenland’s glacial fjords, heat and freshwater are exchanged between glaciers and the ocean. Submarine melting of glaciers has been implicated as a potential trigger for recent glacier acceleration, and observations of ocean heat transport are increasingly being used to infer the submarine melt rates. The complete heat, salt, and mass budgets that underlie such methods, however, have been largely neglected. Here, a new framework for exploring glacial fjord budgets is developed. Building on estuarine studies of salt budgets, the heat, salt, and mass transports through the fjord are decomposed, and new equations for calculating freshwater fluxes from submarine meltwater and runoff are presented. This method is applied to moored records from Sermilik Fjord, near the terminus of Helheim Glacier, to evaluate the dominant balances in the fjord budgets and to estimate freshwater fluxes. Throughout the year, two different regimes are found. In the nonsummer months, advective transports are balanced by changes in heat/salt storage within their ability to measure; freshwater fluxes cannot be inferred as a residual. In the summer, a mean exchange flow emerges, consisting of inflowing Atlantic water and outflowing glacially modified water. This exchange transports heat toward the glacier and is primarily balanced by changes in storage and latent heat for melting ice. The total freshwater flux increases over the summer, reaching 1200 ± 700 m3 s−1 of runoff and 1500 ± 500 m3 s−1 of submarine meltwater from glaciers and icebergs in August. The methods and results highlight important components of fjord budgets, particularly the storage and barotropic terms, that have been not been appropriately considered in previous estimates of submarine melting.