Coupled modelling of subglacial hydrology and calving-front melting at Store Glacier, West Greenland

Abstract. We investigate the subglacial hydrology of Store Glacier in West Greenland, using the open-source, full-Stokes model Elmer/Ice in a novel 3D application that includes a distributed water sheet, as well as discrete channelised drainage, and a 1D model to simulate submarine plumes at the calving front. At first, we produce a baseline winter scenario with no surface meltwater. We then investigate the hydrological system during summer, focussing specifically on 2012 and 2017, which provide examples of high and low surface-meltwater inputs, respectively. We show that the common assumption of zero winter freshwater flux is invalid, and we find channels over 1 m2 in area occurring up to 5 km inland in winter. We also find that the production of water from friction and geothermal heat is sufficiently high to drive year-round plume activity, with ice-front melting averaging 0.15 m d−1. When the model is forced with seasonally averaged surface melt from summer, we show a hydrological system with significant distributed sheet activity extending 65 and 45 km inland in 2012 and 2017, respectively; while channels with a cross-sectional area higher than 1 m2 form as far as 55 and 30 km inland. Using daily values for the surface melt as forcing, we find only a weak relationship between the input of surface meltwater and the intensity of plume melting at the calving front, whereas there is a strong correlation between surface-meltwater peaks and basal water pressures. The former shows that storage of water on multiple timescales within the subglacial drainage system plays an important role in modulating subglacial discharge. The latter shows that high melt inputs can drive high basal water pressures even when the channelised network grows larger. This has implications for the future velocity and mass loss of Store Glacier, and the consequent sea-level rise, in a warming world. </jats:p

A Full-Stokes 3-D Calving Model Applied to a Large Greenlandic Glacier

Iceberg calving accounts for around half of all mass loss from both the Greenland and Antarctic ice sheets. The diverse nature of calving and its complex links to both internal dynamics and climate make it challenging to incorporate into models of glaciers and ice sheets. Here, we present results from a new open-source 3D full-Stokes calving model developed in Elmer/Ice. The calving model implements the crevasse depth criterion, which states that calving occurs when surface and basal crevasses penetrate the full thickness of the glacier. The model also implements a new 3D rediscretization approach and a time-evolution scheme which allow the calving front to evolve realistically through time. We test the model in an application to Store Glacier, one of the largest outlet glaciers in West Greenland, and find that it realistically simulates the seasonal advance and retreat when two principal environmental forcings are applied. These forcings are 1) submarine melting in distributed and concentrated forms, and 2) ice mélange buttressing. We find that ice mélange buttressing is primarily responsible for Store Glacier’s seasonal advance and retreat. Distributed submarine melting prevents the glacier from forming a permanent floating tongue, while concentrated plume melting has a disproportionately large and potentially destabilizing effect on the calving front position. Our results also highlight the importance of basal topography, which exerts significant control on calving, explaining why Store Glacier has remained stable during a period when neighboring glaciers underwent prolonged interannual retreat

Effect of near-terminus subglacial hydrology on tidewater glacier submarine melt rates

Submarine melting of Greenlandic tidewater glacier termini is proposed as a possiblemechanism driving their recent thinning and retreat. We use a general circulation model, MITgcm, tosimulate water circulation driven by subglacial discharge at the terminus of an idealized tidewater glacier.We vary the spatial distribution of subglacial discharge emerging at the grounding line of the glacier andexamine the effect on submarine melt volume and distribution. We find that subglacial hydrology exerts animportant control on submarine melting; under certain conditions a distributed system can induce a factor5 more melt than a channelized system, with plumes from a single channel inducing melt over only alocalized area. Subglacial hydrology also controls the spatial distribution of melt, which has the potential tocontrol terminus morphology and calving style. Our results highlight the need to constrain near-terminussubglacial hydrology at tidewater glaciers if we are to represent ocean forcing accurately

Scalings for Submarine Melting at Tidewater Glaciers from Buoyant Plume Theory

Rapid dynamic changes at the margins of the Greenland Ice Sheet, synchronous with ocean warming, have raised concern that tidewater glaciers can respond sensitively to ocean forcing. Understanding of the processes encompassing ocean forcing nevertheless remains embryonic. The authors use buoyant plume theory to study the dynamics of proglacial discharge plumes arising from the emergence of subglacial discharge into a fjord at the grounding line of a tidewater glacier, deriving scalings for the induced submarine melting. Focusing on the parameter space relevant for high discharge tidewater glaciers, the authors suggest that in an unstratified fjord the often-quoted relationship between total submarine melt volume and subglacial discharge raised to the ⅓ power is appropriate regardless of plume geometry, provided discharge lies below a critical value. In these cases it is then possible to formulate a simple equation estimating total submarine melt volume as a function of discharge, fjord temperature, and calving front height. However, once linear stratification is introduced—as may be more relevant for fjords in Greenland—the total melt rate discharge exponent may be as large as ¾ (⅔) for a point (line) source plume and display more complexity. The scalings provide a guide for more advanced numerical models, inform understanding of the processes encompassing ocean forcing, and facilitate assessment of the variability in submarine melting both in recent decades and under projected atmospheric and oceanic warming

Iceberg melting substantially modifies oceanic heat flux towards a major Greenlandic tidewater glacier

Fjord dynamics influence oceanic heat flux to the Greenland ice sheet. Submarine iceberg melting releases large volumes of freshwater within Greenland’s fjords, yet its impact on fjord dynamics remains unclear. We modify an ocean model to simulate submarine iceberg melting in Sermilik Fjord, east Greenland. Here we find that submarine iceberg melting cools and freshens the fjord by up to ~5 °C and 0.7 psu in the upper 100-200 m. The release of freshwater from icebergs drives an overturning circulation, resulting in a ~10% increase in net up-fjord heat flux. In addition, we find that submarine iceberg melting accounts for over 95% of heat used for ice melt in Sermilik Fjord. Our results highlight the substantial impact that icebergs have on the dynamics of a major Greenlandic fjord, demonstrating the importance of including related processes in studies that seek to quantify interactions between the ice sheet and the ocean

Coupled modelling of subglacial hydrology and calving-front melting at Store Glacier, West Greenland

Abstract. We investigate the subglacial hydrology of Store Glacier in West Greenland, using the open-source, full-Stokes model Elmer/Ice in a novel 3D application that includes a distributed water sheet, as well as discrete channelised drainage, and a 1D model to simulate submarine plumes at the calving front. At first, we produce a baseline winter scenario with no surface meltwater. We then investigate the hydrological system during summer, focussing specifically on 2012 and 2017, which provide examples of high and low surface-meltwater inputs, respectively. We show that the common assumption of zero winter freshwater flux is invalid, and we find channels over 1 m2 in area occurring up to 5 km inland in winter. We also find that the production of water from friction and geothermal heat is sufficiently high to drive year-round plume activity, with ice-front melting averaging 0.15 m d−1. When the model is forced with seasonally averaged surface melt from summer, we show a hydrological system with significant distributed sheet activity extending 65 and 45 km inland in 2012 and 2017, respectively; while channels with a cross-sectional area higher than 1 m2 form as far as 55 and 30 km inland. Using daily values for the surface melt as forcing, we find only a weak relationship between the input of surface meltwater and the intensity of plume melting at the calving front, whereas there is a strong correlation between surface-meltwater peaks and basal water pressures. The former shows that storage of water on multiple timescales within the subglacial drainage system plays an important role in modulating subglacial discharge. The latter shows that high melt inputs can drive high basal water pressures even when the channelised network grows larger. This has implications for the future velocity and mass loss of Store Glacier, and the consequent sea-level rise, in a warming world. </jats:p

BedMachine v3: Complete Bed Topography and Ocean Bathymetry Mapping of Greenland From Multibeam Echo Sounding Combined With Mass Conservation

Greenland's bed topography is a primary control on ice flow, grounding line migration, calving dynamics, and subglacial drainage. Moreover, fjord bathymetry regulates the penetration of warm Atlantic water (AW) that rapidly melts and undercuts Greenland's marine-terminating glaciers. Here we present a new compilation of Greenland bed topography that assimilates seafloor bathymetry and ice thickness data through a mass conservation approach. A new 150 m horizontal resolution bed topography/bathymetric map of Greenland is constructed with seamless transitions at the ice/ocean interface, yielding major improvements over previous data sets, particularly in the marine-terminating sectors of northwest and southeast Greenland. Our map reveals that the total sea level potential of the Greenland ice sheet is 7.42 ± 0.05 m, which is 7 cm greater than previous estimates. Furthermore, it explains recent calving front response of numerous outlet glaciers and reveals new pathways by which AW can access glaciers with marine-based basins, thereby highlighting sectors of Greenland that are most vulnerable to future oceanic forcing