Articles | Volume 15, issue 7
Research article
06 Jul 2021
Research article |  | 06 Jul 2021

Indication of high basal melting at the EastGRIP drill site on the Northeast Greenland Ice Stream

Ole Zeising and Angelika Humbert

The accelerated ice flow of ice streams that reach far into the interior of the ice sheets is associated with lubrication of the ice sheet base by basal meltwater. However, the amount of basal melting under the large ice streams – such as the Northeast Greenland Ice Stream (NEGIS) – is largely unknown. In situ measurements of basal melt rates are important from various perspectives as they indicate the heat budget, the hydrological regime and the relative importance of sliding in glacier motion. The few previous estimates of basal melt rates in the NEGIS region were 0.1 m a−1 and more, based on radiostratigraphy methods. These findings raised the question of the heat source, since even an increased geothermal heat flux could not deliver the necessary amount of heat. Here, we present basal melt rates at the recent deep drill site EastGRIP, located in the centre of NEGIS. Within 2 subsequent years, we found basal melt rates of 0.19±0.04 m a−1 that are based on analysis of repeated phase-sensitive radar measurements. In order to quantify the contribution of processes that contribute to melting, we carried out an assessment of the energy balance at the interface and found the subglacial water system to play a key role in facilitating such high melt rates.

1 Introduction

Ice sheet models are used to quantify the contribution of the Greenland Ice Sheet (GrIS) to future sea-level rise under different climatic scenarios. In these simulations, the distinctive extent of Greenland's largest ice stream – the Northeast Greenland Ice Stream (NEGIS, Fig. 1) – can only be reproduced well if a higher-order approximation is considered for the momentum balance and initial states are based on inversion (Goelzer et al.2018) or involve subglacial hydrological models (Smith-Johnsen et al.2020a). Primarily, this is due to the model's inability to accurately represent lubrication and thus the subsequent sliding at the ice stream base that occurs.

The NEGIS is the only large ice stream in Greenland, extending from a distance of 100 km from the ice divide over a length of about 700 km towards the coast (Fahnestock et al.1993, 2001b; Joughin et al.2001). It drains about 12 % of Greenland's ice through three major outlet glaciers Nioghalvfjerdsbrae, Zachariæ Isstrøm and Storstrømmen Glacier (Rignot and Mouginot2012). Loss of the floating tongue of Zachariæ Isstrøm has already led to ice flow acceleration and increased mass loss (Mouginot et al.2015). Consequently, it is expected and projected that NEGIS will contribute significantly to sea-level rise in the future (Khan et al.2014), highlighting the importance of understanding the general ice-flow dynamics and its driving mechanisms.

One hypothesis for the genesis of NEGIS is locally increased basal melting at the onset area that enables and enhances basal sliding (Fahnestock et al.2001a; Christianson et al.2014; Franke et al.2021) and forms a subglacial hydrological system. The coupling with basal sliding is facilitated via the water pressure, so that the sliding velocity rises with increasing water pressure (e.g. Beyer et al.2018; Smith-Johnsen et al.2020a). However, little is known about the amount of subglacial water below the up to ∼3300 m thick ice sheet.

First estimates of basal melt rates by Fahnestock et al. (2001a) and later by Keisling et al. (2014) and MacGregor et al. (2016) are based on the interpretation of chronology in radiostratigraphy. All three studies found melt rates of 0.1 m a−1 and more – which is extremely large for inland ice. However, these estimates may be prone to limited validity given the assumptions about the flow regime and constant accumulation rate. The cause for such intensive melt was attributed to a high geothermal heat flux, which possibly originates from the passage of Greenland over the Icelandic hot spot (Fahnestock et al.2001a; Rogozhina et al.2016; Martos et al.2018; Alley et al.2019).

In order to directly observe, among other things, flow regimes and basal conditions of ice streams, an ice core is being drilled as part of the East Greenland Ice-Core Project (EastGRIP) near the onset of the NEGIS. Here, surface velocities reach about 57 m a−1 (Hvidberg et al.2020) and the NEGIS widens (Fig. 1). Smith-Johnsen et al. (2020a) forced an ice model with a locally increased heat flux below the EastGRIP drill site. They found that a heat flux of 0.97 W m−2 (corresponding to a basal melt rate of 0.1 m a−1; Fahnestock et al.2001a) is necessary to reasonably reproduce the velocities of NEGIS. By utilising a coupled subglacial hydrology and ice sheet model, Smith-Johnsen et al. (2020b) demonstrated the large impact of an uncertainty in geothermal heat flux on the flow of NEGIS arising from the subglacial hydrological system and hence from basal melting and water pressure, as well as from friction.

However, measurements with an adequate accuracy are still required to narrow down the basal melt rates further. Here, we present the first estimates of basal melt rates from repeated in situ phase-sensitive radar measurements from the EastGRIP drill site and consider the contribution of different heat sources at the ice base.

Figure 1Surface ice flow velocity map of the Greenland Ice Sheet (Joughin et al.2018). The box in the overview map (upper left corner) marks the boundaries of the main figure showing northeast Greenland and the Northeast Greenland Ice Stream (NEGIS), which drains into the three major outlet glaciers, namely Nioghalvfjerdsbrae (79 N Glacier, 79NG), Zachariæ  Isstrøm (ZI) and Storstrømmen Glacier (SG). The location of the EastGRIP drill site is denoted by the black triangle.

2 Data and methods

2.1 Instrument, data acquisition and processing

The autonomous phase-sensitive radio-echo sounder (ApRES; Brennan et al.2014; Nicholls et al.2015) is a low-power, ground-based radar that allows measurements to be carried out autonomously with a selected interval over long periods of time. By analysing the phase shifts of the return signals, the vertical displacements of internal reflections and of the base can be precisely determined in the millimetre range. Thus, the ApRES is often used to determine Lagrangian basal melt rates and their temporal variability of ice shelves (e.g. Stewart et al.2019; Washam et al.2019; Vaňková et al.2020). In order to derive an annual mean basal melt rate at the EastGRIP drill site, we deployed an ApRES within a near-surface trench (Fig. 1). The ApRES performed a measurement once a day during winter from August 2017–April 2018 and August 2018–May 2019.

In the following, we shortly describe the theory of operation of the ApRES and the signal processing, whereas a more detailed description is given by Brennan et al. (2014), Nicholls et al. (2015) and Stewart et al. (2019). Within a single measurement, the ApRES transmits a sequence of 100 chirps, each with a duration of 1 s in which the frequency of the transmitted electromagnetic wave is increased from 200 to 400 MHz. After reflection, the received signal is mixed with a replica of the transmitted signal and sampled with 40 kHz (Nicholls et al.2015). Since the frequencies of the resulting deramped signal are related to the two-way travel time, a spectral analysis needs to be done in the processing to obtain depth profiles of the amplitude and phase. For the conversion from travel time to depth, we used a vertical propagation velocity of 168 914 km s−1 according to the relative permittivity of εr=3.15 (Fujita et al.2000). Prior to the spectral analysis, the performed chirps were averaged to increase the signal-to-noise ratio. However, due to weak reflections in the lower part of the ice, the signal-to-noise ratio is reduced. As a consequence, no reliable analysis of the data is possible below the noise level depth limit of ∼1450 m, with the exception of the basal return. Next, we present how we derive basal melt rates from the radar data.

2.2 Ice thickness evolution

The method we use to derive a basal melt rate is based on the ice thickness evolution equation that is valid in both the Eulerian and Lagrangian reference systems:

(1) H t = - div Q + a s - a b ,

with the ice thickness H, the time t, the volume flux Q, the surface mass balance as and the basal melt rate ab (positive for melting) (e.g. Greve and Blatter2009). Equation (1) states that a temporal change in ice thickness is caused by a changing volume flux arising from deformation and accumulation or ablation at the ice surface and base. It is worth noting that a basal melt rate larger than the accumulation rate only leads to thinning of the glacier, if the volume flux cannot supply sufficient ice to balance this out. The volume flux Q is defined as the vertically integrated horizontal velocities vx,vy(x,y,z,t):

(2) Q = Q x Q y = 0 H v x d z 0 H v y d z ,

(Greve and Blatter2009), and it represents the ice thickness change due to deformation and sliding and thus stretching or compression in the horizontal direction. This may, for example, be due to changes in basal velocities or ice creeping across a bedrock undulation. Using the continuity equation for incompressible materials, div v=0, and Leibniz's integral rule, we can rewrite div Q as

(3) div Q = x 0 H v x d z + y 0 H v y d z = 0 H v x x + v y y d z = 0 H v z z d z = 0 H ε ˙ z z d z

with ε˙zz the vertical strain rate ε˙zz=vz/z.

The recorded ApRES time series allows for a precise estimation of changes in ice thickness ΔH from the vertical displacement of the basal reflector and of internal layers from consecutive measurements. However, applying the ice thickness evolution equation (Eq. 1) to the ApRES measurements requires some modifications. Since the ApRES is located within a trench below the surface, the “measured ice thickness” H is defined as the range between the ApRES and the ice base. The total ice thickness – the range from the surface to the ice base – is about 7 to 8 m thicker and includes the upper firn and snow layers. Thus, ΔH is independent of the surface mass balance, as=0 m a−1, but influenced by firn densification that significantly affects the vertical displacement in the upper ∼100 m. As this is not considered in Eq. (1), we add the term ΔHf/Δt to correct for the densification process below the ApRES.

Equation (3) states that the divergence of the volume flux in Eq. (1) can be expressed by the depth-integrated vertical strain rate. However, we derive vertical displacements uz from ApRES measurements instead of vertical velocities vz. Thus, we can calculate strain εzz=uz/z for a time period of Δt. Therefore, Eq. (3) needs to be reformulated as

(4) 0 H ε ˙ z z d z = 1 Δ t 0 H ε z z d z = Δ H ε Δ t ,

with the change in ice thickness due to vertical strain ΔHε. Finally, the modified ice thickness evolution equation can be written as

(5) Δ H Δ t = Δ H f Δ t + Δ H ε Δ t - a b .

All three quantities ΔH, ΔHf and ΔHε, which are needed to derive ab, are described by vertical displacements and hence by the radar measurement itself in a consistent manner.

2.3 Derivation of basal melt rates

In order to derive vertical displacements of internal layers and of the basal return from the ApRES time series, we slightly modified the processing of Vaňková et al. (2020) (details below). Both methods are based on phase differences estimated from cross-correlation of the repeated measurements.

Firstly, we divided the depth profile into 6 m segments with a 3 m overlap from a depth of 20 m below the antennas to 20 m above the ice base and a wider segment of 10 m (9 to +1 m) around the basal return, characterised by a strong increase in amplitude. In order to derive vertical displacements, each depth segment of the first measurement (t1) was cross-correlated with the same segment of each repeated measurement (ti). This is in contrast to Vaňková et al. (2020), who derived displacements from pairwise time-consecutive measurements (ti−1ti). The lag of the minimum mean phase difference obtained from the cross-correlation gives the cumulative displacement at the given depth. The range of expected lag was limited by the estimation to the previous measurement (t1ti−1). This results in a time series of displacements for each segment individually. The vertical displacement of the basal segment is the change in the measured ice thickness ΔH.

Next, we estimate the vertical strain εzzobs and quantify ΔHf as well as ΔHε based on a regression analysis of the vertical displacements. To avoid influences of firn densification on the determination of εzzobs, we excluded all segments above a depth of 250 m (∼9 % of all segments). In addition, segments below the noise-level depth limit of h≈1450 m (where noise prevents an unambiguous estimation) were excluded (∼45 % of all segments). Furthermore, outliers were filtered out (∼7 %). We found a linear fit uz(z) of

(6) u z ( z ) = ε z z obs z + Δ H f , 250 m z h

that best matches the cumulative vertical displacements of the remaining  400 segments within the ice. The gradient of this fit is εzzobs, and the shift between the intercept at the depth of the ApRES and ΔH is ΔHf. However, εzz for zh is unknown. Here, we used two scenarios to estimate ΔHε (Fig. 2, Appendix Fig. A1). First, we assumed that εzz is constant with depth:

(7) ε z z const ( z ) = ε z z obs , 0 z H .

As a second scenario, we used a vertical strain distribution (εzzsim) obtained from an ice sheet model based on inverse surface flow velocities (Rückamp et al.2020). Here, εzzsim increases with depth and reaches values of roughly twice εzzobs at the base.

In order to be less dependent on a single measurement, we compute for each of the last 65 d (records; roughly 25 % of the measurements) of a year an annual melt rate and compute from these 65 melt rate estimates a mean annual value by averaging. Finally, ΔHε was derived from Eq. (4) for the two vertical strain distributions (ΔHεconst, ΔHεsim), and the basal melt rate ab from Eq. (5). Given errors are based on the standard deviation of the estimates based on the considered 65 measurements and a 1 % uncertainty in the signal propagation speed in ice (Fujita et al.2000). For visualisation, we calculated the cumulative vertical displacement referenced to the ice base (Fig. 2).

3 Results

The analysis of the 2017/18 ApRES time series revealed a measured ice thickness (distance between radar and ice base) of roughly 2668 m at the EastGRIP drill site with an annual mean change of −0.471 m a−1 (Table 1, Fig. 2). The firn densification – the intercept (see Fig. 2) of the linear fit at z=0 m (the elevation of the ApRES) – occurring below the radar is 0.074 m a−1. We derived a vertical strain of εzz=-0.068×10-3 from reliable estimates of vertical displacements feasible to a depth of 1450 m. The dynamic thinning of the ice derived from the two scenarios ranges from −0.181 m a−1 (ΔHεconst) to −0.194 m a−1 (ΔHεsim). This results in a basal melt rate of 0.210 ± 0.015 m a−1. The numbers derived from the time series recorded in 2018/19 differ slightly (Table 1, Appendix Fig. A2). The annual mean change in measured ice thickness is 27 mm (or 6 %) lower and the firn densification 15 mm (or 20 %) larger compared to the values derived in 2017/18. The resulting basal melt rate of 0.167 ± 0.018 m a−1 is ∼20 % lower than the year before. Finally, we derive an averaged melt rate over both years of 0.19 ± 0.04 m a−1.

Table 1Results for measured ice thickness change (ΔH), firn densification (ΔHf), vertical strain (εzz), dynamic ice thickness change obtained from a constant vertical strain (ΔHεconst) and a simulation (ΔHεsim), and basal melting (ab) for both time series projected to 365 d. Negative values contribute to the thinning of the ice column, whereas a positive melt rate represents melting.

Download Print Version | Download XLSX

Figure 2Derived vertical displacements uz of the depth of the ApRES (ΔH; blue dot) and of selected internal layers referenced to the ice base from the 2017/18 ApRES time series. Derived displacements used for melt rate estimations are marked by red dots and of layers within the firn by grey dots. The estimated displacements between a depth of 250 m and h are used to calculate a linear fit (solid black line), the gradient of which is the vertical strain. Extrapolations to the bottom are shown by the dashed lines. The offset at the ice base is caused by basal melting, and the difference between the intercept of the linear fit at z=0 m and ΔH is the firn compaction.


4 Discussion

We used estimated vertical displacements from the upper half of the ice column to estimate the dynamic thinning, since noise prevents an unambiguous estimation of the vertical strain for the lower half. To cover a range of variations in the dynamic thinning, we used two different scenarios for vertical strain distribution. The resulting dynamic thinning of the simulated vertical strain and the constant strain differ only slightly. However, we cannot exclude the possibility that larger strain values are reached at the base, which would lead to an overestimation of the basal melt rates. In the case of a non-existing melt rate, the dynamic thinning of the lower half of the ice column would be, on average, more than 4 times as large as the one of the upper half. However, a strong increase is not found in higher-order ice sheet simulations (Rückamp et al.2020). A frequently used strain distribution (e.g. Fahnestock et al.2001a; Keisling et al.2014; MacGregor et al.2016) that takes into account deviating strain within a shear zone is the Dansgaard–Johnsen distribution model (Dansgaard and Johnsen1969). As this model assumes a linearly decreasing strain in the shear zone that reaches zero at the ice base, the resulting basal melt rate at EastGRIP would be even larger. However, the Dansgaard–Johnsen model represents a no-slip boundary condition at the ice base. As this is an unrealistic assumption in an ice stream, we did not consider the Dansgaard–Johnsen model further.

The derived vertical strain is based on more than 300 vertical displacements estimated between the firn–ice transition and about 1450 m. In contrast, the estimation of the displacement of the basal return is based on the phase shift of only one segment around the basal return, slightly above the noise level. This makes the determination more prone to errors. Instead of comparing the first measurement (t1) with all repeated measurements (ti), the pairwise comparison of time-consecutive measurements (ti−1 and ti), as shown by Vaňková et al. (2020), leads to a lower thinning rate of ΔH in 2017/18 than in 2018/19 (-0.441±0.004 m a−1 in 2017/18, -0.467±0.009 m a−1 in 2018/19). Thus, the variability found is not necessarily a variability of the ice sheet system but can rather be influenced by the methodology.

A variation in the selected depth limit of densification, to exclude segments affected by densification, causes slight changes in vertical strain and thus in basal melt rate on the order of millimetres. However, we observed an increased densification rate within the considered 65 records. The increased densification can possibly be a result of increased load from the camp at the surface.

Our derived basal melt rate of 0.19 ± 0.04 m a−1 is above previous estimates from airborne radar measurements. Fahnestock et al. (2001a) and MacGregor et al. (2016) found melt rates on the order of 0.1 m a−1 in the vicinity of the EastGRIP drill site but larger melt rates of >0.15 m a−1 further upstream in the onset region of NEGIS. Both studies used a constant vertical strain over depth where basal melting occurs. Smith-Johnsen et al. (2020a) found that basal melt rates of 0.1 m a−1, derived from a heat flux of 0.97 W m−2, are needed at the location of the EastGRIP drill site to reproduce the NEGIS in an ice sheet model.

4.1 Considerations of the energy balance at the ice base

In order to constrain the heat flux required to sustain the basal melt rates ab derived in this study, we consider the energy balance at the ice base. As for any surface across which a physical quantity may not be continuous, a jump condition is formulated. In a typical continuum mechanical formulation, the jump ([[ψ]]) of a quantity ψ is defined as [[ψ]]=ψ+-ψ-, meaning the difference in the quantity ψ across the interface (Greve and Blatter2009). The jump condition of the energy at the ice base reads as


with the heat flux q, the velocity v, the velocity of the singular surface w, the normal vector n pointing outwards from the ice body, the Cauchy stress t, the ice density ρi and the internal energy u (Greve and Blatter2009). The jump of the heat flux [[qn]] becomes (qgeo+qsw)n-κ(T)gradT, with qgeo the geothermal heat flux and qsw the heat flux from subglacial water with a temperature above pressure melting point, T temperature and κ thermal conductivity. For the jump in work of surface forces we find

(9) [ [ v t n ] ] = v sw t sw n - v b i t i n ,

with tsw the Cauchy stress of the subglacial water side of the singular surface, vbi the ice velocity and ti the stress field of the ice at the base.

Figure 3Sketch of the energy balance at the ice base. Thermal components are coloured in orange, and mechanical components are in blue-purple colour. For a detailed explanation see the main text.


We split the traction vector of the subglacial water into a normal and tangential component, with the water pressure psw and the stress in the normal direction. Following the same approach as at an ice shelf base (Greve and Blatter2009), we employ an empirical relation

(10) t sw n = - p sw n + C i / sw ρ sw | v sw | 2 e t

with et=vsw/|vsw| and etn. The drag coefficient at the underside of the ice is Ci/sw, similar to the Manning roughness taken into account in subglacial conduits. So the part of the subglacial water becomes


with vsw and vsw the normal and tangential velocity of the subglacial water, respectively. This formulation is quite similar to the treatment of the jump condition at an ice shelf base. For the traction vector at the ice base, we follow the same procedure and find

(12) t i n = - N n + τ b e t ,

with N the normal component and τb the component in the tangential plane. For vbitin we find

(13) v b i t i n = - N v b i n + τ b v b i e t .

With the jump of the internal energy [[u]]=L, we can reformulate Eq. (8) to


The tangential components Ci/swρsw|vsw|3 and τbvbiet are frictional heating and dominate the contribution of heat arising from work of surface forces. They need to be seen as two end members of the system: either the ice is only in contact with a thick subglacial hydrological system, in which case Ci/swρsw|vsw|3 is active, or the subglacial hydrological system is permanently in contact with a lubricated base, in which case the second term τbvbiet is governing. The components are visualised in Fig. 3.

Next, we aim at constraining the individual terms for which we use the following material parameters: ρi=910 kg m−3, the latent heat of fusion, L=335 kJ kg−1 and the thermal conductivity for ice at the pressure melting point of 270.81 K κ(270.81 K)=2.10 W m−1 K−1 (Greve and Blatter2009).

We consider three scenarios: (i) there is only temperate ice that is melting, (ii) heat is required to warm the ice to the pressure melting point and (iii) friction at the base contributes significantly to basal melting. (i) For temperate ice and no heat arising from work of surface forces, we find a melt rate of at least 0.19 m a−1 to correspond to a heat flux of 1.84 W m−2. (ii) Considering grad T to be less than 10−1 K m−1, this increases the required heat flux from scenario (i) by up to 0.21 W m−2, as this additional heat is required to warm the ice to the pressure melting point. (iii) Heat arising from work of the surface forces may, however, reduce the required heat flux into the ice to melt this amount of ice. To this end, we need to estimate the magnitude of the components of the stress tensors.

We assume that the normal stress component N is hydrostatic and bridging stresses to be negligible. With a mean density of ice of 910 kg m−3 we find pi=23.8 MPa. The normal velocity is of the order of the basal melt rate vb-0.2 m a−1 by assuming the velocity of the interface (w) to be zero. The normal component of the ice side is then on the order of 0.15 W m−2. For the tangential components of the ice side, we consider the shear stress at the base to be |τb|1 to 100 kPa. This compares to basal shear stress found by Rückamp et al. (2020) of 50 kPa. To constrain the sliding velocity, we assume it to be maximum the surface velocity of 57 m a−1 and minimum half of the surface velocity. This leads to a tangential component on the ice side to be up to 0.15 W m−2 (Fig. 4).

Next, we constrain the normal component of the subglacial water pswvsw. A water pressure of 10 to 23 MPa is consistent with subglacial hydrological modelling (Beyer et al.2018; Smith-Johnsen et al.2020a). Assuming the normal velocity to be at most as large as the basal melt rate, we find the range of this term to be between 0.05 and 0.12 W m−2 (Fig. 4). The tangential component Ci/swρw|vsw|3 needs an assumption on the roughness Ci/sw, for which we consider a range from the roughness of the ice shelf base of 10−3 to a maximum roughness 10 times as large.

Figure 4Magnitude of terms in the energy jump condition at the ice base. (a) Contribution of work of surface forces in the normal direction from water pressure (lines) and ice overburden pressure (dots). (b) Contribution of work of surface forces in the tangential direction from friction on the subglacial water. (c) Contribution of work of surface forces in the tangential direction from friction on the ice side.


The motivation for this is that ice shelf roughness is governed by convection cells at the interface, whereas in the inland ice, the interaction with the bedrock may lead to a larger roughness. As nothing is known about the shape of the subglacial conduit, the range of velocity cannot be constrained well. We consider a speed similar to the one of the ocean 0.1 m s−1, but as surface rivers easily reach 1.0 m s−1, we take this as an upper limit (Fig. 4). Thus, the contribution of friction to the energy available for basal melting may account for at least ∼0.20 W m−2, with the potential to be far larger based on the assumptions we made.

To summarise, the jump in the tangential component (friction) has the potential to govern the heat budget, depending on flow speeds in subglacial water and roughness of the ice base, as can be seen in Fig. 4. However, assuming the geothermal heat flux to be on the order of 𝒪(qgeo)≈0.25 W m−2 makes evident that the key player in facilitating such high melt rates is the subglacial water system, which may supply the ice base with an additional heat flux.

We have focused our consideration on the interface between a subglacial water layer and the ice, as this drives the basal melt rate. However, observations of Christianson et al. (2014) highlight the existence of a wet till layer beneath the ice stream. Depending on the thickness of the water layer, the velocity and pressure of the water, and the porosity of the till layer, complex interaction between the till and water may arise, too. Kutscher et al. (2019) present high-resolution simulations of a comparable system that highlight the importance of studying this interface as well. To date, it is unclear which vertical extent of the water layer is required to decouple the interaction of a water–till interface from the ice–water interface and thus the influence on the basal melt rate.

Large basal melting mainly affects basal sliding, as it increases the effective normal pressure. Considering sliding to be the dominant part of ice stream motion, large basal melt rates are also plausible in the respect that the subglacial hydrological system needs to be sustained over time, and hence creep closure of cavities or conduits needs to be balanced by melting.

Future measurements at EastGRIP after successful completion of the drilling to the ice base will shed more light on the sliding speed and may also provide more information on the characteristics of the subglacial hydrological system. This will enable the community to put our melt rate estimates into further context.

5 Conclusions

We estimated annual mean basal melt rates at the EastGRIP drill site from time series of high-precision phase-sensitive radar measurements. We derived the change in the measured ice thickness, thinning from firn densification occurring below the instrument and the vertical strain in the upper 1450 m of the roughly 2668 m thick ice. Two different scenarios for vertical strain distribution were used to quantify a plausible range of dynamic thinning. Thus, we derived an averaged melt rate of 0.19 ± 0.04 m a−1. We are aware that these melt rates require an extremely large amount of heat that we suggest to arise from the subglacial water system and the geothermal heat flux. However, these melt rates are based on measurements with a modern ice-penetrating radar whose penetration depth is limited due to transmitting power. Thus, no assumptions on past accumulation rates or other uncertainties in age reconstruction are involved. Our major uncertainty is the vertical strain in the lower part of the ice stream. This could be overcome if a more powerful radar with a similar vertical resolution could be operated autonomously over several months.

Appendix A: Additional figures of ApRES processing

Figure A1Observed vertical strain distribution (solid black line) and different extensions to the ice base (dashed lines) for two scenarios based on a constant (dashed black line), and a simulated (Rückamp et al.2020, dashed orange line) distribution.


Figure A2Derived vertical displacements uz of the depth of the ApRES (ΔH; blue dot) and of selected internal layers referenced to the ice base from 2018/19 ApRES time series. Derived displacements used for melt rate estimations are marked by red dots and for layers within the firn by grey dots. The estimated displacements between a depth of 250 m and h are used to calculate a linear fit (solid black line), the gradient of which is the vertical strain. Extrapolations to the bottom are shown by the dashed lines. The offset at the ice base is caused by basal melting, and the difference between the intercept of the linear fit at z=0 m and ΔH is the firn compaction.


Code availability

Our own developed MATLAB routines for ApRES processing are available from the corresponding author on request.

Data availability

Raw data of the ApRES measurements (; Zeising and Humbert2021) are available at the World Data Center PANGAEA.

Author contributions

OZ and AH conducted part of the field study. OZ processed the ApRES data and estimated the resulting basal melt rates. AH designed the study and conducted the analysis of the energy jump condition. OZ and AH discussed the time series. OZ and AH wrote the manuscript jointly.

Competing interests

The authors declare that they have no conflict of interest.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


Data were acquired at the EastGRIP camp that kindly hosted this activity as an associate project. EastGRIP is directed and organised by the Centre for Ice and Climate at the Niels Bohr Institute, University of Copenhagen. It is supported by funding agencies and institutions in Denmark (A. P. Møller Foundation, University of Copenhagen), the United States (US National Science Foundation, Office of Polar Programs), Germany (Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research), Japan (National Institute of Polar Research and Arctic Challenge for Sustainability), Norway (University of Bergen and Trond Mohn Foundation), Switzerland (Swiss National Science Foundation), France (French Polar Institute Paul-Emile Victor, Institute for Geosciences and Environmental research), Canada (University of Manitoba) and China (Chinese Academy of Sciences and Beijing Normal University). We would like to thank Johanna Kerch, Sepp Kipfstuhl and Daniel Steinhage for their support in the deployment and the dismount of the measuring device. We are grateful for discussions with Thomas Kleiner and Martin Rückamp (AWI) and Ralf Greve (ITLS, Japan) on vertical strain rates and our results.

Financial support

The article processing charges for this open-access publication were covered by the Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research (AWI).

Review statement

This paper was edited by Alexander Robinson and reviewed by Basile de Fleurian and two anonymous referees.


Alley, R. B., Pollard, D., Parizek, B. R., Anandakrishnan, S., Pourpoint, M., Stevens, N. T., MacGregor, J. A., Christianson, K., Muto, A., and Holschuh, N.: Possible Role for Tectonics in the Evolving Stability of the Greenland Ice Sheet, J. Geophys. Res.-Earth, 124, 97–115,, 2019. a

Beyer, S., Kleiner, T., Aizinger, V., Rückamp, M., and Humbert, A.: A confined–unconfined aquifer model for subglacial hydrology and its application to the Northeast Greenland Ice Stream, The Cryosphere, 12, 3931–3947,, 2018. a, b

Brennan, P. V., Lok, L. B., Nicholls, K., and Corr, H.: Phase-sensitive FMCW radar system for high-precision Antarctic ice shelf profile monitoring, IET Radar Sonar Nav., 8, 776–786,, 2014. a, b

Christianson, K., Peters, L. E., Alley, R. B., Anandakrishnan, S., Jacobel, R. W., Riverman, K. L., Muto, A., and Keisling, B. A.: Dilatant till facilitates ice-stream flow in northeast Greenland, Earth Planet. Sc. Lett., 401, 57–69,, 2014. a, b

Dansgaard, W. and Johnsen, S.: A flow model and a time scale for the ice core from Camp Century, Greenland, J. Glaciol., 8, 215–223,, 1969. a

Fahnestock, M., Bindschadler, R., Kwok, R., and Jezek, K.: Greenland Ice Sheet Surface Properties and Ice Dynamics from ERS-1 SAR Imagery, Science, 262, 1530–1534,, 1993. a

Fahnestock, M. A., Abdalati, W., Joughin, I., Brozena, J., and Gogineni, P.: High geothermal heat flow, basal melt, and the origin of rapid ice flow in central Greenland, Science, 294, 2338–2342,, 2001a. a, b, c, d, e, f

Fahnestock, M. A., Joughin, I., Scambos, T. A., Kwok, R., Krabill, W. B., and Gogineni, S.: Ice-stream-related patterns of ice flow in the interior of northeast Greenland, J. Geophys. Res.-Atmos., 106, 34035–34045,, 2001b. a

Franke, S., Jansen, D., Beyer, S., Neckel, N., Binder, T., Paden, J., and Eisen, O.: Complex Basal Conditions and Their Influence on Ice Flow at the Onset of the Northeast Greenland Ice Stream, J. Geophys. Res.-Earth, 126, e2020JF005689,, 2021. a

Fujita, S., Matsuoka, T., Ishida, T., Matsuoka, K., and Mae, S.: A summary of the complex dielectric permittivity of ice in the megahertz range and its applications for radar sounding of polar ice sheets, in: Physics of ice core records, Hokkaido University Press, 185–212, 2000. a, b

Goelzer, H., Nowicki, S., Edwards, T., Beckley, M., Abe-Ouchi, A., Aschwanden, A., Calov, R., Gagliardini, O., Gillet-Chaulet, F., Golledge, N. R., Gregory, J., Greve, R., Humbert, A., Huybrechts, P., Kennedy, J. H., Larour, E., Lipscomb, W. H., Le clec'h, S., Lee, V., Morlighem, M., Pattyn, F., Payne, A. J., Rodehacke, C., Rückamp, M., Saito, F., Schlegel, N., Seroussi, H., Shepherd, A., Sun, S., van de Wal, R., and Ziemen, F. A.: Design and results of the ice sheet model initialisation experiments initMIP-Greenland: an ISMIP6 intercomparison, The Cryosphere, 12, 1433–1460,, 2018. a

Greve, R. and Blatter, H.: Dynamics of Ice Sheets and Glaciers, Springer Science & Business Media,, 2009. a, b, c, d, e, f

Hvidberg, C. S., Grinsted, A., Dahl-Jensen, D., Khan, S. A., Kusk, A., Andersen, J. K., Neckel, N., Solgaard, A., Karlsson, N. B., Kjær, H. A., and Vallelonga, P.: Surface velocity of the Northeast Greenland Ice Stream (NEGIS): assessment of interior velocities derived from satellite data by GPS, The Cryosphere, 14, 3487–3502,, 2020. a

Joughin, I., Fahnestock, M., MacAyeal, D., Bamber, J. L., and Gogineni, P.: Observation and analysis of ice flow in the largest Greenland ice stream, J. Geophys. Res.-Atmos., 106, 34021–34034,, 2001. a

Joughin, I., Smith, B. E., and Howat, I. M.: A complete map of Greenland ice velocity derived from satellite data collected over 20 years, J. Glaciol., 64, 1–11,, 2018. a

Keisling, B. A., Christianson, K., Alley, R. B., Peters, L. E., Christian, J. E., Anandakrishnan, S., Riverman, K. L., Muto, A., and Jacobel, R. W.: Basal conditions and ice dynamics inferred from radar-derived internal stratigraphy of the northeast Greenland ice stream, Ann. Glaciol., 55, 127–137,, 2014. a, b

Khan, S. A., Kjær, K. H., Bevis, M., Bamber, J. L., Wahr, J., Kjeldsen, K. K., Bjørk, A. A., Korsgaard, N. J., Stearns, L. A., Van Den Broeke, M. R., et al.: Sustained mass loss of the northeast Greenland ice sheet triggered by regional warming, Nat. Clim. Change, 4, 292–299,, 2014. a

Kutscher, K., Geier, M., and Krafczyk, M.: Multiscale simulation of turbulent flow interacting with porous media based on a massively parallel implementation of the cumulant lattice Boltzmann method, Comput. Fluids, 193, 103733,, 2019. a

MacGregor, J. A., Fahnestock, M. A., Catania, G. A., Aschwanden, A., Clow, G. D., Colgan, W. T., Gogineni, S. P., Morlighem, M., Nowicki, S. M., Paden, J. D., Price, S. F., and Seroussi, H.: A synthesis of the basal thermal state of the Greenland Ice Sheet, J. Geophys. Res.-Earth, 121, 1328–1350,, 2016. a, b, c

Martos, Y. M., Jordan, T. A., Catalán, M., Jordan, T. M., Bamber, J. L., and Vaughan, D. G.: Geothermal heat flux reveals the Iceland hotspot track underneath Greenland, Geophys. Res. Lett., 45, 8214–8222,, 2018. a

Mouginot, J., Rignot, E., Scheuchl, B., Fenty, I., Khazendar, A., Morlighem, M., Buzzi, A., and Paden, J.: Fast retreat of Zachariæ Isstrøm, northeast Greenland, Science, 350, 1357–1361,, 2015.  a

Nicholls, K. W., Corr, H. F., Stewart, C. L., Lok, L. B., Brennan, P. V., and Vaughan, D. G.: A ground-based radar for measuring vertical strain rates and time-varying basal melt rates in ice sheets and shelves, J. Glaciol., 61, 1079–1087,, 2015. a, b, c

Rignot, E. and Mouginot, J.: Ice flow in Greenland for the international polar year 2008–2009, Geophys. Res. Let., 39, L11501,, 2012. a

Rogozhina, I., Petrunin, A. G., Vaughan, A. P., Steinberger, B., Johnson, J. V., Kaban, M. K., Calov, R., Rickers, F., Thomas, M., and Koulakov, I.: Melting at the base of the Greenland ice sheet explained by Iceland hotspot history, Nat. Geosci., 9, 366–369,, 2016. a

Rückamp, M., Goelzer, H., and Humbert, A.: Sensitivity of Greenland ice sheet projections to spatial resolution in higher-order simulations: the Alfred Wegener Institute (AWI) contribution to ISMIP6 Greenland using the Ice-sheet and Sea-level System Model (ISSM), The Cryosphere, 14, 3309–3327,, 2020. a, b, c, d

Smith-Johnsen, S., de Fleurian, B., Schlegel, N., Seroussi, H., and Nisancioglu, K.: Exceptionally high heat flux needed to sustain the Northeast Greenland Ice Stream, The Cryosphere, 14, 841–854,, 2020a. a, b, c, d, e

Smith-Johnsen, S., Schlegel, N.-J., de Fleurian, B., and Nisancioglu, K. H.: Sensitivity of the Northeast Greenland Ice Stream to Geothermal Heat, J. Geophys. Res.-Earth, 125, e2019JF005252,, 2020b. a

Stewart, C. L., Christoffersen, P., Nicholls, K. W., Williams, M. J., and Dowdeswell, J. A.: Basal melting of Ross Ice Shelf from solar heat absorption in an ice-front polynya, Nat. Geosci., 12, 435–440,, 2019. a, b

Vaňková, I., Nicholls, K. W., Corr, H. F., Makinson, K., and Brennan, P. V.: Observations of tidal melt and vertical strain at the Filchner-Ronne Ice Shelf, Antarctica, J. Geophys. Res.-Earth, 125, e2019JF005280,, 2020. a, b, c, d

Washam, P., Nicholls, K. W., Münchow, A., and Padman, L.: Summer surface melt thins Petermann Gletscher Ice Shelf by enhancing channelized basal melt, J. Glaciol., 65, 662–674,, 2019. a

Zeising, O. and Humbert, A.: Autonomous phase-sensitive radar (ApRES) measurements at EastGRIP drill site, Greenland, PANGAEA,, 2021. a

Short summary
Greenland’s largest ice stream – the Northeast Greenland Ice Stream (NEGIS) – extends far into the interior of the ice sheet. Basal meltwater acts as a lubricant for glaciers and sustains sliding. Hence, observations of basal melt rates are of high interest. We performed two time series of precise ground-based radar measurements in the upstream region of NEGIS and found high melt rates of 0.19 ± 0.04 m per year.