Articles | Volume 16, issue 4
Research article
27 Apr 2022
Research article |  | 27 Apr 2022

Basal melt of the southern Filchner Ice Shelf, Antarctica

Ole Zeising, Daniel Steinhage, Keith W. Nicholls, Hugh F. J. Corr, Craig L. Stewart, and Angelika Humbert

Basal melt of ice shelves is a key factor governing discharge of ice from the Antarctic Ice Sheet as a result of its effects on buttressing. Here, we use radio echo sounding to determine the spatial variability of the basal melt rate of the southern Filchner Ice Shelf, Antarctica, along the inflow of Support Force Glacier. We find moderate melt rates with a maximum of 1.13 m a−1 about 50 km downstream of the grounding line. The variability of the melt rates over distances of a few kilometres is low (all but one < 0.15 m a−1 at < 2 km distance), indicating that measurements on coarse observational grids are able to yield a representative melt rate distribution. A comparison with remote-sensing-based melt rates revealed that, for the study area, large differences were due to inaccuracies in the estimation of vertical strain rates from remote sensing velocity fields. These inaccuracies can be overcome by using modern velocity fields.

1 Introduction

Filchner Ice Shelf (FIS), a West Antarctic ice shelf draining major East Antarctic ice streams (Bailey, Slessor, Recovery and Support Force glaciers), is thought to be vulnerable to a change in its basal mass balance within this century (Hellmer et al.2012) as a result of the possible penetration of relatively warm, off-shelf waters into the ocean cavity beneath the ice shelf. Subsequent thinning of the ice shelf would reduce its buttressing to inland glaciers, allowing them to speed up and thin and their grounding lines to retreat landward. If the stress perturbation is sufficiently large, then a positive ice-loss feedback may occur as the ice sheet's grounding line retreats across the deepening beds of the tributary ice streams (Schoof2012). The current discharge of ice across the grounding line at FIS is 106.3 ± 5.7  Gt a−1 (Rignot et al.2019), which is about 9.6 % of the discharge from East Antarctica, underlining the importance of understanding the current state of the ice shelf for assessing future change in basal melt. In addition, precise melt rates serve as validation for models projecting the future contribution of these ice streams to sea level change.

Basal melt rates can be derived from remote sensing data by solving the ice thickness evolution equation (Rignot et al.2013; Moholdt et al.2015; Berger et al.2017; Adusumilli et al.2020a). Although the Lagrangian approach adopted in recent years (Moholdt et al.2015) has led to improvements, major uncertainties from various factors remain. Hence, in situ observations of basal melt rates are required for assessing the reliability of remote sensing approaches. This is even more urgent, as remote-sensing-derived basal melt rates are used to construct parameterisations that diagnose basal melt rates from modelled sub-ice-shelf ocean conditions. These models are used to project the contribution of Antarctica to sea level change. Significant errors in observed distributions of basal melt rate therefore have a profound effect on the outcome of projections of future sea level rise, such as ISMIP6 (Seroussi et al.2020), as a result of their effect on the calibration of basal melt rate parameterisations (Jourdain et al.2020).

In recent years, the use of the phase-sensitive radio echo sounder (pRES) opened new possibilities for the precise determination of basal melt rates. Nicholls et al. (2015) and Stewart et al. (2019) presented basal melt rates from near Ross Island, Ross Ice Shelf, Antarctica, which were derived from 10 d of autonomous pRES (ApRES) measurements, and measurements from 78 stations, time-averaged between 2013 and 2014. Stewart et al. (2019) observed strong seasonal melt rate variability, with values up to 53 m a−1 within a 5 d period in January 2013 and an exponentially reducing mean annual basal melt rate with increasing distance from the calving front, with values up to 7.7 m a−1. Vaňková et al. (2020) presented a tidal melt and vertical strain analysis from 17 ApRES records across Filchner–Ronne Ice Shelf. They found the tidal vertical strain to be depth dependent only near the grounding line, with significant tidal melt measurable at some locations. The derived melt rates were used by Bull et al. (2021) to evaluate an ocean model. Marsh et al. (2016) investigated basal melt rates at 25 points at a melt channel near the grounding line of Ross Ice Shelf. They found basal melt rates decreasing from 22 m a−1 at the upstream end of the channel to 2.5 m a−1 40 km downstream. A strong seasonal variability in melt rate was recorded by Washam et al. (2019) on Petermann Glacier, Greenland, using an ApRES recording on the flank of a basal melt channel. In summer 2016, they found extreme melt rates equivalent to 80 m a−1 but most of the year the mean basal melt rate ranged from 0 to 10 m a−1.

Our survey is focused on the accessible southern part of FIS, which might be more susceptible to the potential inflow of warm waters (Hellmer et al.2012). Recent observations from hot-water-drilled boreholes through Filchner–Ronne Ice Shelf have revealed an interannual change in circulation mode starting in 2017, highlighting the variability in conditions within the sub-ice-shelf cavity (Hattermann et al.2021).

Here, we aim at understanding the magnitude and variation in basal melt over an area extending from the grounding line of Support Force Glacier, as far downstream as was feasible. In austral summer 2015/16, under the framework of the Filchner Ice Shelf Project (FISP), pRES measurements were carried out at a total of 94 locations and then repeated a year later. The stations were distributed along the central flow line of Support Force Glacier's extension on to FIS and along four cross sections, providing along-flow and across-flow melt rate distributions (Fig. 1). A further transect crossed the entire FIS south of Berkner Island. As far as safety allowed, we extended the profiles along the eastern margin towards the inland ice, to capture an area where gradients in the bathymetry were expected, steering the flow of water masses. With this observational design, we intended to measure the large-scale distribution of melt rates, but in addition we included more closely spaced stations to detect variations on short spatial scales. In the following, we first introduce the methodology and the data basis. We then present and discuss the derived basal melt rates and compare them with remote sensing data.

Figure 1(a) Map of the Ronne and Filchner ice shelves (BedMachine Antarctica; Morlighem2020; Morlighem et al.2020) with the marked study area (black box) near the Support Force Glacier (SFG). (b) Study area with derived basal melt rates (dots), grouped depending on their location on the central flow line (CFL) and five cross sections CSA, CSB, CSC, CSD and CSE. Nearby pRES measurements with distances < 2 km are shown by squares, and the location of an ApRES station (FSE1 in Vaňková et al.2020) is shown by a star. Multi-channel Coherent Radar Depth Sounder (MCoRDS) profiles (P1–P3) flown as part of NASA's Operation IceBridge (OIB) campaign in 2016 (Paden et al.2014) are shown by black lines (echograms are included in the Appendix, Fig. A1). Background colour shows the ice shelf thickness from BedMachine Antarctica (Morlighem2020; Morlighem et al.2020).

2 Materials and methods

Our estimation of basal melt rates is based on measurements using a pRES that is described in detail in Brennan et al. (2014) and Nicholls et al. (2015). The pRES transmits a frequency-modulated sweep (chirp) from 200 to 400 MHz over a period of 1 s via two skeleton slot antennas, separated by roughly 9 m. The exact locations were marked with bamboo for precise relocation a year later. After internal processing, only the difference in frequency between the transmitted and received signals, called the deramped frequency, is saved. Details of the internal processing are given by Brennan et al. (2014). By repeating the measurements after a time period, we are able to track changes in depth of internal reflectors within the ice, and of the basal echo, to a precision of millimetres. This allows the study of firn densification, vertical strain due to ice flow, and the (Lagrangian) change in ice shelf thickness. Being a Lagrangian measurement, no steady state assumption is required, and the basal melt rate can be separated from the overall change in ice thickness.

Our 94 measurement stations are grouped depending on their locations on the central flow line (CFL) and five cross sections (CS) A–E (for location, see Fig. 1). The time period between repeated measurements varied from 323 (18 January 2016–6 December 2016) to 356 d (31 December 2015–21 December 2016). To improve the signal-to-noise ratio, we recorded 100 chirps at each site. Correlations were calculated between each chirp and the 99 remaining chirps, and those chirps with a low average correlation were discarded. Those remaining chirps were averaged and then Fourier transformed to yield a complex (amplitude and phase) profile as a function of two-way-travel time. To convert the profile into a function of range, we calculated the velocity profile of the electromagnetic wave for each location by estimating the density–depth profile based on Herron and Langway (1980) with accumulation and mean annual temperature from RACMO 2.3/ANT (van Wessem et al.2014). Nevertheless, the uncertainty of the propagation velocity is about 1 % (Fujita et al.2000).

Using a procedure similar to that described by Corr et al. (2002) and Jenkins et al. (2006), we aligned the two radar profiles using a 6 m window below the firn–ice boundary by cross-correlating the amplitude profiles. This provided a datum within the ice column, removing the effects of instrument temperature change, firn densification, and snow accumulation and ablation.

The thickness change (DHi/Dt) in the solid ice below the aligned reflector is caused only by the dynamic ice thickness change due to vertical strain (Hiε˙zz) and by the basal melt rate ab:

(1) D H i D t = H i ε ˙ zz - a b ,

with Hi the solid-ice thickness below the aligned reflector and ε˙zz the vertical strain rate. In order to determine the vertical strain, the displacement between visits was calculated with a cross-correlation of the amplitude and phase information for each layer deeper than the aligned reflector. Under the plane-strain assumption the vertical strain is constant with depth; a least-squares method was used to calculate a linear fit of the shift of those layers that exhibited a high correlation value. The gradient of the linear fit is the vertical strain. The change in ice thickness below the aligned reflector is derived from the shift of the basal reflector, which was calculated in the same way as the shift of the internal layers. The largest error in the calculation comes from the alignment of the data because it is based only on the amplitude correlation. The uncertainty in the calculation of the phase shift is closely related to the signal-to-noise ratio of the reflectors. An additional uncertainty arises from the assumption of a linear displacement–depth relation, although this is generally thought reliable for plug flow (Greve and Blatter2009). These uncertainties add up to 0.03 m.

3 Results and discussion

3.1 Large-scale spatial variability

Of the 94 measurements, 79 were suitable for retrieval of basal melt rates. The main reasons for excluding the other 15 stations are (1) low correlation values in the depth of the firn–ice transition, which made it impossible to align the measurements; (2) changes in the shape of the basal reflector that prevented the reflections from being unequivocally matched; and (3) too few high correlation values for a linear fit to be used to calculate the vertical strain rate. An explanation for the low correlation values could be errors in operating the pRES, such as inaccuracies in the alignment of the antennas or incorrectly seated cables. In addition, changes in settings such as the attenuation affected the signal-to-noise ratio, thereby reducing the number of high correlation values.

Figure 2Distribution of pRES-derived results of (a) the basal melt rate ab, (b) the change in solid-ice thickness DHi/Dt and (c) the dynamic change in solid-ice thickness Hiε˙zz. The numbers in the upper right corner state the mean value and the standard deviation (SD).


For the remaining stations, we found a mean basal melt rate of 0.38 ± 0.26 m a−1 (mean + standard deviation; Fig. 2a) with a maximum of 1.13 ± 0.03 m a−1 at a location about 50 km downstream of the grounding line and freezing in the northernmost part of the central flow line (CFL). Of similar size but with different sign, the mean value of DHi/Dt is 0.38 ± 0.32 m a−1 (Fig. 2b), representing a thinning of the solid ice, whereas the mean value of Hiε˙zz (0.01 ± 0.20 m a−1) is close to zero (Fig. 2c).

This study focuses on the spatial variability of the melt rates, rather than the overall, annual average values, since we did not measure the interannual variability. However, the different sites were occupied for different periods, and thus, seasonality in basal melt would affect spatial variability if it exists. The time periods range from 365 minus 42 d to 365 minus 9 d with an increasing data acquisition interval southwards. Seasonality may affect the derived annual melt rate differently at the different sites. However, in the northern part of our study area an autonomous pRES (ApRES) station recorded for more than a year, including the period of the once-repeated pRES measurements (Fig. 1b, Vaňková et al.2020). For the time period between 18 January 2016 and 6 December 2016 (same period as the pRES measurements with the shortest time interval), we derived a melt rate of 0.02 ± 0.03 m a−1. In the period of 365 d (18 January 2016–17 January 2017), a slightly lower melt rate of 0.01 ± 0.03 m a−1 was derived. This indicates that no enhanced melting occurred at the location of the ApRES in summer 2016/2017. However, we cannot assess if melt rates are enhanced at other locations.

We present the distribution of ab in Fig. 1b, as well as an along-flow profile (CFL) and five cross sections (CSA-CSE) in Fig. 3. A total of 70 % of the estimated basal melt rates range between 0 and 0.50 m a−1. Higher melt rates were found for nine stations within 100 km of the grounding line at the CFL, CSA and CSB. All three stations with ab>1ma-1 are located in this part of the study area. The variation in ab along ice flow is weak and shows no clear trend of increasing melt towards the grounding line, despite the increasing ice draft (Fig. 3a). In the direction across ice flow (Fig. 3b–f) the largest variations in ab appear in the two southernmost cross sections (CSA, CSB). Apart from the southern part, higher basal melt rates of up to 0.82 ± 0.03 m a−1 occur only at CSD. The northernmost cross section, ranging from Berkner Island towards the inland ice (CSE), has a generally low ab. Three stations, all at the northernmost part of CFL, indicate freezing.

A key assumption made during pRES processing is that the phase shift on reflection at the ice–ocean interface remains constant. Although this is valid for a fresh ice–seawater interface, it is a poor assumption for the interface between fresh ice and possibly slushy marine ice, itself underlain by seawater. This means that if either of a pair of measurements is made during a period of freeze-on it is not possible to distinguish the change in the phase of the basal reflection that results from a change in its range, from the phase shift that results from the change in the nature of the basal interface. Thus, we can not determine the amount of accretion at the three sites.

In an ice shelf cavity where the water speeds are relatively high as a result either of tides, as in this case (Mueller et al.2018), or of strong buoyancy-driven flows, as in the case of a steeply inclined ice base over relatively warm water (Lazeroms et al.2019), basal melt rates are mainly controlled by three factors: the basal drag coefficient, the thermal driving and the water speed in the boundary layer. The thermal driving is the difference in the temperature of the water near the ice base and the freezing point of that water at the pressure of the ice base. The water speed and the basal drag generate the shear-driven turbulence that efficiently diffuses heat and salt towards the ice base (Holland and Jenkins1999). A fourth factor, discussed below, is the basal vertical temperature gradient in the ice.

In our study area the basal slopes are generally low (Morlighem2020; Morlighem et al.2020), as is the thermal driving. We therefore expect tidal speeds to dominate buoyancy-driven flows. Ice draft plays the role of modifying the thermal driving: lower basal pressures reduce the freezing point, thereby reducing the local freezing temperature. Mueller et al. (2018) find a strong increase in tidal speeds from the grounding line to the cross section CSE; this parallels a reducing basal draft, which will act to reduce the thermal driving. We expect these two tendencies to work together to modulate the large-scale spatial variation in basal melt rates.

The large-scale spatial variation in ab can also be influenced by changes in vertical gradients of the ice temperature. An ice shelf fed by a fast glacier typically contains a cold core as a result of ice advection, leading to larger vertical temperature gradients some distance from the grounding line. However, with melting over centuries, the ice temperature is more likely to approach a parabolic profile, with only moderate temperature gradients (Humbert2010).

Figure 3Variation in the basal melt rate (a) along the central flow line (CFL) of Support Force Glacier's extension on the FIS and (b–f) the cross sections CSA–CSE (left y axis). Locations are shown in Fig. 1. The pRES-derived values are shown in dark blue. The dark grey line represents remote-sensing-derived melt rates, and the light grey bounds display the uncertainty, both published by Adusumilli et al. (2020a). Uncertainties of the pRES-derived melt rates are 0.03 cm and therefore too small to visualise. The variation in the ice draft (in metres above sea level) is shown on the right y axis. For CFL, the distance refers to the grounding line (GL) of Support Force Glacier and for all cross sections to the CFL with positive distances on the eastern side.


3.2 Small-scale spatial variability

In order to assess the small-scale spatial variability of the basal melt rates and hence the representativeness of our point measurements over larger distances, we carried out 17 pRES measurements that were located within 2 km of another measurement (Fig. 1).

Overall, the local differences in melt rate (Δab) are small (Fig. 4). Except for one outlier, Δab ranges from 0 to 0.13 m a−1 with a median of only 0.02 m a−1. In the case of the outlier station pair (pRES060 and pRES061), the difference was 0.68 ± 0.06 m a−1 within a distance of 977 m (Fig. 3 and Fig. B1). Here, the higher melt rate was measured at site pRES060, which has a 10.5 m higher ice shelf draft, as derived from the surface elevation given by BedMachine (Morlighem2020; Morlighem et al.2020) and the pRES-derived ice thickness.

Variability at small spatial scales will not result from variations in tidal speed: in the absence of strong sea floor or ice base topography, strong horizontal gradients in tidal speed are not expected. As previously discussed, tides will dominate buoyancy-driven currents, which are therefore also unlikely to play a significant role in controlling local variations in melt rate. However, through its effect on thermal driving, even a quite modest local variation in basal depth is a candidate for causing variation in basal melting. A change in ice draft of, say, 10 m will change the thermal driving by about 0.007 C (e.g. Holland and Jenkins1999). Using the algorithm proposed by Jenkins et al. (2010) for the nearby Ronne Ice Shelf, and for an appropriate mean tidal speed of 0.1 m s−1 (Mueller et al.2018), a 0.007 C change in thermal driving would result in a melt rate difference of 0.17 m a−1. Given thermal stratification in the underlying water column, it is also possible for a locally increased draft to result in the ice base having contact with warmer waters, leading to a local increase in melt rates. Since such an increased draft of 10.5 m was found at pRES060 compared to pRES061, this could be a possible explanation for the large difference in melt rate of 0.68 ± 0.06 m a−1. Another candidate driver of local variability in melt rate is spatial variability in basal roughness. Differences in the drag coefficient at the ice base will directly affect melt rates (e.g. Holland and Jenkins1999). However, neither the ice roughness nor its spatial variation is known well enough to determine its importance in driving local spatial variation in melt rates.

Overall, this analysis gives evidence that individual measurements are representative of a large area on the scale of one to three ice thicknesses. Only minor variation is to be expected due to the specific choice of measurement location. Airborne radar echograms (Fig. A1), recorded within NASA's Operation IceBridge (OIB) with a Center for Remote Sensing of Ice Sheets (CReSIS) Multi-channel Coherent Radar Depth Sounder (MCoRDS) in 2016 (Paden et al.2014), show that our study area is mostly characterised by a smooth ice shelf base without terrace structures, supporting our interpretation that small-scale variability in basal melt rate is relatively low. One exception is a basal channel in the west with a height of approx. 50 m (Fig. A1e and f). At three locations around this channel, pRES measurements have been performed (western part of CSC in Fig. 3d). These show low variability, with a basal melt rate  0.2 m a−1 lower in the centre of the channel.

Figure 4Small-scale variability of basal melt rates. Difference in basal melt rate Δab as function of distance between two nearby stations (average distance: 1126 ± 296 m). The colour of each dot represents ab of the station with the larger basal melt rate.


4 Comparison with remote sensing

The analysis of remote-sensing-derived basal melt rates is based on precisely measured elevation changes of the ice shelf surface and on the correction of the surface mass balance, firn densification and dynamic change in ice thickness (e.g. Moholdt et al.2015; Adusumilli et al.2020a). The dynamic change in ice thickness and thus the vertical strain rate is often derived from the divergence of a satellite sensor-derived surface velocity field.

We used the pRES-derived vertical strain rates to assess the reliability of strain rates derived from different remote sensing velocity fields. Satellite-derived melt rates at FIS from Rignot et al. (2013), Moholdt et al. (2015) and Adusumilli et al. (2020a) were all based on the strain rates derived from the same early MEaSUREs velocities (Rignot et al.2011; Scheuchl et al.2012). However, this velocity field contained some significant data gaps in our study area that were not present in modern velocity fields such as the Landsat Ice Speed of Antarctica (LISA) product from which vertical strain rates were derived by Alley et al. (2018) or the newest MEaSUREs data set (Mouginot et al.2019a, b). Instead of comparing the vertical strain rate itself, we compared the dynamic ice thickness change (Hiε˙zz) that was derived from the vertical strain rate and the solid-ice thickness. The results reveal a significant improvement over the last decade in the accuracy of the determination of vertical strain rates from remote sensing.

While the average deviation between the pRES-derived product and that from Moholdt et al. (2015) was 0.40 ± 0.44 m a−1 (mean ± standard deviation; Fig. 5a and b and Fig. C1), there were much smaller deviations (0.01 ± 0.35 m a−1) from the product of Alley et al. (2018) (Fig. 5c and d and Fig. C1). The comparison with the dynamic ice thickness that we calculated using the latest MEaSUREs data set (Mouginot et al.2019a, b) also showed only minor deviations of 0.04 ± 0.17 m a−1 (Fig. 5e and f and Fig. C1). Here, similar to Moholdt et al. (2015), we applied a Gaussian filter with a 27 km× 27 km window to smooth the velocity data and calculated the divergence to obtain the vertical strain rate. The comparison highlights the recent improvement in the estimation of velocity fields for more accurate calculation of dynamic ice thickness changes and demonstrates good agreement between remote-sensing-derived strain rates and those from in situ measurements.

Remote-sensing-derived melt rates published by Rignot et al. (2013), Moholdt et al. (2015) and Adusumilli et al. (2020a) suggested a similar pattern of melt rates: southeast of Berkner Island, a freezing regime in the west switches to a melting regime eastwards, with melting persisting towards the south to the Support Force Glacier. However, a data gap in the velocity field meant that no melt rates could be determined by Rignot et al. (2013) for a large part of our study area.

The comparison with the results from Adusumilli et al. (2020a) reveals a broader distribution of the remote-sensing-derived melt rate (1.1–1.6 m a−1) at the pRES locations with an average deviation from the pRES-derived values of 0.35 ± 0.57 m a−1 (Figs. 3 and 5g and h), which is of a size similar to the deviation of the dynamic change in ice thickness. Another reason for the discrepancies can be the different measurement periods over which the basal melt rates were estimated: Adusumilli et al. (2020a) show that basal melt rates can vary at interannual timescales. In order to investigate whether different measurement periods contributed to the discrepancies between the results from the different methods, we compared the change in ice thickness DHi/Dt (Eq. 1) after the correction for the surface mass balance and firn densification (Fig. C2). Some of the differences occur because Adusumilli et al. (2020a) define Hi as the ice shelf thickness in units of metres of ice equivalent, which is slightly higher than the solid-ice thickness that we use for the pRES-based estimates. However, the comparison of DHi/Dt shows a good agreement, with an average difference of only 0.04 ± 0.24 m a−1 (Fig. 5i and j). Since variations in basal melt rate contribute to DHi/Dt and this only shows slight differences, a temporal variation in basal melting can be excluded as the reason for the significant discrepancies that we find. Furthermore, this indicates that the techniques derive consistent changes in ice thickness from their initial measurements after applying the corrections for the surface mass balance and the firn densification and that the large differences in basal melt rates result principally from differences in the strain rate, which can be improved by the use of modern surface velocity products.

Figure 5Comparison of remote sensing (grey) and pRES-derived (blue) results. The left column shows the distributions (a) of the dynamic change in ice thickness Hiε˙zz for the results published by Moholdt et al. (2015), (c) Alley et al. (2018) and (e) derived from the MEaSUREs product (Mouginot et al.2019a, b); (g) of the basal melt rate ab; and (i) of the change in ice thickness DHi/Dt, both in comparison with those from Adusumilli et al. (2020a). The right column (b, d, f, h and j) shows the distribution of the deviation between remote sensing and pRES-derived values according to (a, c, e, g and i). The numbers in the upper right corner state the mean value and the standard deviation (SD). Positive values refer to larger numbers derived from the remote-sensing-based method.


5 Conclusions

We have presented the first spatial distribution of basal melt rates in the southern Filchner Ice Shelf derived from repeated phase-sensitive radar measurements. In general the melt rates are moderate with maximum values in the centre of less than 1.13 m a−1. We tested the representativeness of individual measurements by assessing the variability over short distances. Spatial variability in ab is low, with occasional outliers possibly linked to large basal gradients. This gives us confidence that a small number of widely spaced measurements accurately represent the large-scale melt pattern. Temporal variability, however, is not captured.

Our in situ measurements reveal that inaccuracies in the estimation of dynamic ice thickness change negatively affected recent remote-sensing-derived melt rates at our study area at the Filchner Ice Shelf. A comparison with strain rates published by Alley et al. (2018) and with those derived from the newest MEaSUREs velocity field indicates that these inaccuracies can be overcome by using state-of-the-art velocity fields, in which data gaps could be closed. Our study demonstrates that satellite-derived basal melt rates hold great promise, but care needs to be taken, as modelling of the future contribution of Antarctica to sea level rise is currently calibrated using such products (Jourdain et al.2020). This highlights the need to obtain more data sets such as the one presented here, from across different ice shelves, and to conduct repeated field surveys to assess temporal variability.

Appendix A: Airborne radar echograms

Figure A1Airborne radar echograms (a) P1, (c) P2 and (e) P3 (location in Fig. 1), recorded with a Multi-channel Coherent Radar Depth Sounder (MCoRDS) as part of NASA's Operation IceBridge (OIB) campaign in 2016 (Paden et al.2014; Arnold et al.2020). (b, d and f) Insets showing enlarged basal section visualised by black box in (a), (c)and (e). The white dots mark the depth of the ice base derived from a nearby pRES measurement.


Appendix B: pRES echograms

Figure B1(a–d) Amplitude profiles of first (grey line) and repeated measurement at locations pRES060 (a, b; blue) and pRES061 (c, d; red). Insets in (b) and (d) showing enlarged basal section, visualised by black boxes in (a) and (c). (b, d) Vertical dashed lines mark the ice thickness and DHi the change in ice thickness between both visits. The correlation coefficients of the basal segments are 0.95 (pRES060) and 0.96 (pRES061). (e) Vertical displacements of internal segment (small dots) and of the basal segment (large dot) for pRES060 (blue) and pRES061 (red). The values given are the basal melt rates at both locations.


Appendix C: Comparison with remote sensing

Figure C1Variation in the dynamic ice thickness change Hiε˙zz (a) along the central flow line (CFL) of Support Force Glacier's extension on the FIS and (b–f) the cross sections CSA–CSE. Locations are shown in Fig. 1. pRES-derived values are shown in blue. Remote-sensing-derived values are represented by the solid grey line for results published by Moholdt et al. (2015), by a dashed line for results published by Alley et al. (2018) and by a dotted line for estimations derived from the MEaSUREs product (Mouginot et al.2019a, b). The bounds of the results from Moholdt et al. (2015) display the uncertainties. Derived errors of the pRES measurements are too small to visualise. For CFL, the distance refers to the grounding line (GL) of Support Force Glacier and for all cross sections to the CFL with positive distances on the eastern side.


Figure C2Variation in the ice thickness change DHi/Dt (a) along the central flow line (CFL) of Support Force Glacier's extension on the FIS and (b–f) the cross sections CSA–CSE. Locations are shown in Fig. 1. pRES-derived values are shown in blue. Remote-sensing-derived values are represented by the dark grey line for results published by Adusumilli et al. (2020a). The light grey bounds display the uncertainties. Derived errors of the pRES measurements are too small to visualise. For CFL, the distance refers to the grounding line (GL) of Support Force Glacier and for all cross sections to the CFL with positive distances on the eastern side.


Code and data availability

MATLAB routines for pRES processing are available from the corresponding author on request. Raw data of the pRES measurements and derived melt rates (, Zeising et al.2021) are available at the World Data Center PANGAEA. Echograms recorded with a Center for Remote Sensing of Ice Sheets (CReSIS) Multi-channel Coherent Radar Depth Sounder (MCoRDS) within NASA's Operation IceBridge (OIB) campaign in 2016 can be accessed at (Paden et al.2014). Basal melt rate data published by Adusumilli et al. (2020a) can be accessed at (Adusumilli et al.2020b). Ice shelf divergence and thickness data published by Moholdt et al. (2015) can be accessed at (Moholdt et al.2016). Strain rate data published by Alley et al. (2018) can be accessed through open ftp by contacting the National Snow and Ice Data Center (NSIDC) (last access: 25 June 2021). The MEaSUREs velocity product can be accessed at (Mouginot et al.2019a). The BedMachine Antarctica product can be accessed at (Morlighem2020) (last access: 12 April 2021).

Author contributions

OZ processed the data. OZ and AH analysed the basal melt rates and wrote the manuscript together with KWN. AH designed the study and conducted the field study together with DS. KWN supported the field study and contributed to melt rate analysis and its discussion together with HFJC and AH. CLS wrote the main parts of the MATLAB routines and supported data analysis. All authors reviewed various versions of the paper.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


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


This work was funded by the AWI strategy fund project FISP. Support for this work came from the UK Natural Environment Research Council large grant “Ice shelves in a warming world: Filchner Ice Shelf System” (NE/L013770/1). We acknowledge provision of data products by Susheel Adusumilli (Scripps Institution of Oceanography, University of California San Diego, La Jolla, CA, USA), that have been published in Adusumilli et al. (2020a). We want to thank Graham Niven and Bradley Morrell for their support in the field. We are grateful for discussions with Veit Helm and Niklas Neckel (AWI) on remote-sensing-derived vertical strain rates and our results. We thank the two anonymous reviewers for careful reading and very helpful suggestions.

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 Nicolas Jourdain and reviewed by two anonymous referees.


Adusumilli, S., Fricker, H. A., Medley, B., Padman, L., and Siegfried, M. R.: Interannual variations in meltwater input to the Southern Ocean from Antarctic ice shelves, Nat. Geosci., 13, 616–620,, 2020a. a, b, c, d, e, f, g, h, i, j, k, l

Adusumilli, S., Fricker, H. A., Medley, B. C., Padman, L., Siegfried, M. R.: Data from: Interannual variations in meltwater input to the Southern Ocean from Antarctic ice shelves, UC San Diego Library Digital Collections [data set],, 2020b. a

Alley, K. E., Scambos, T. A., Anderson, R. S., Rajaram, H., Pope, A., and Haran, T. M.: Continent-wide estimates of Antarctic strain rates from Landsat 8-derived velocity grids, J. Glaciol., 64, 321–332,, 2018. a, b, c, d, e, f

Arnold, E., Leuschen, C., Rodriguez-Morales, F., Li, J., Paden, J., Hale, R., and Keshmiri, S.: CReSIS airborne radars and platforms for ice and snow sounding, Ann. Glaciol., 61, 58–67,, 2020. a

Berger, S., Drews, R., Helm, V., Sun, S., and Pattyn, F.: Detecting high spatial variability of ice shelf basal mass balance, Roi Baudouin Ice Shelf, Antarctica, The Cryosphere, 11, 2675–2690,, 2017. a

Brennan, P. V., Nicholls, K., Lok, L. B., 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

Bull, C. Y. S., Jenkins, A., Jourdain, N. C., Vaňková, I., Holland, P. R., Mathiot, P., Hausmann, U., and Sallée, J.-B.: Remote Control of Filchner–Ronne Ice Shelf Melt Rates by the Antarctic Slope Current, J. Geophys. Res.-Oceans, 126, e2020JC016550,, 2021. a

Corr, H. F. J., Jenkins, A., Nicholls, K. W., and Doake, C. S. M.: Precise measurement of changes in ice-shelf thickness by phase-sensitive radar to determine basal melt rates, Geophys. Res. Lett., 29, 73–1–74–4,, 2002. 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, edited by: Hondoh, T., pp. 185–212, Hokkaido University Press, Sapporo, 2000. a

Greve, R. and Blatter, H.: Dynamics of Ice Sheets and Glaciers, Springer Science & Business Media, Berlin Heidelberg,, 2009. a

Hattermann, T., Nicholls, K. W., Hellmer, H. H., Davis, P. E. D., Janout, M. A., Østerhus, S., Schlosser, E., Rohardt, G., and Kanzow, T.: Observed interannual changes beneath Filchner–Ronne Ice Shelf linked to large-scale atmospheric circulation, Nat. Commun., 12, 2961,, 2021. a

Hellmer, H. H., Kauker, F., Timmermann, R., Determann, J., and Rae, J.: Twenty-first-century warming of a large Antarctic ice-shelf cavity by a redirected coastal current, Nature, 485, 225–228,, 2012. a, b

Herron, M. M. and Langway, C. C.: Firn densification: An empirical model, J. Glaciol., 25, 373–385, 1980. a

Holland, D. M. and Jenkins, A.: Modeling Thermodynamic Ice–Ocean Interactions at the Base of an Ice Shelf, J. Phys. Oceanogr., 29, 1787–1800,<1787:MTIOIA>2.0.CO;2, 1999. a, b, c

Humbert, A.: The temperature regime of Fimbulisen, Antarctica, Ann. Glaciol., 51, 56–64,, 2010. a

Jenkins, A., Corr, H. F. J., Nicholls, K. W., Stewart, C. L., and Doake, C. S. M.: Interactions between ice and ocean observed with phase-sensitive radar near an Antarctic ice-shelf grounding line, J. Glaciol., 52, 325–346,, 2006. a

Jenkins, A., Nicholls, K. W., and Corr, H. F. J.: Observation and Parameterization of Ablation at the Base of Ronne Ice Shelf, Antarctica, J. Phys. Oceanogr., 40, 2298–2312,, 2010. a

Jourdain, N. C., Asay-Davis, X., Hattermann, T., Straneo, F., Seroussi, H., Little, C. M., and Nowicki, S.: A protocol for calculating basal melt rates in the ISMIP6 Antarctic ice sheet projections, The Cryosphere, 14, 3111–3134,, 2020. a, b

Lazeroms, W. M. J., Jenkins, A., Rienstra, S. W., and van de Wal, R. S. W.: An Analytical Derivation of Ice-Shelf Basal Melt Based on the Dynamics of Meltwater Plumes, J. Phys. Oceanogr., 49, 917–939,, 2019. a

Marsh, O. J., Fricker, H. A., Siegfried, M. R., Christianson, K., Nicholls, K. W., Corr, H. F., and Catania, G.: High basal melting forming a channel at the grounding line of Ross Ice Shelf, Antarctica, Geophys. Res. Lett., 43, 250–255,, 2016. a

Moholdt, G., Padman, L., and Fricker, H. A.: Basal mass budget of Ross and Filchner–Ronne ice shelves, Antarctica, derived from Lagrangian analysis of ICESat altimetry, J. Geophys. Res.-Earth, 119, 2361–2380,, 2015. a, b, c, d, e, f, g, h, i, j, k

Moholdt, G., Padman, L., and Fricker, H. A.: Elevation change and mass budget of Ross and Filchner-Ronne ice shelves, Antarctica, Norwegian Polar Data Centre [data set],, 2016. a

Morlighem, M.: MEaSUREs BedMachine Antarctica, Version 2, NASA National Snow and Ice Data Center Distributed Active Archive Center, Boulder, Colorado USA [data set], (last access: 12 April 2021), 2020. a, b, c, d, e

Morlighem, M., Rignot, E., Binder, T., Blankenship, D., Drews, R., Eagles, G., Eisen, O., Ferraccioli, F., Forsberg, R., Fretwell, P., Goel, V., Greenbaum, J. S., Gudmundsson, H., Guo, J., Helm, V., Hofstede, C., Howat, I., Humbert, A., Jokat, W., Karlsson, N. B., Lee, W. S., Matsuoka, K., Millan, R., Mouginot, J., Paden, J., Pattyn, F., Roberts, J., Rosier, S., Ruppel, A., Seroussi, H., Smith, E. C., Steinhage, D., Sun, B., van den Broeke, M. R., van Ommen, T. D., van Wessem, M., and Young, D. A.: Deep glacial troughs and stabilizing ridges unveiled beneath the margins of the Antarctic ice sheet, Nat. Geosci., 13, 132–137,, 2020. a, b, c, d

Mouginot, J., Rignot, E., and Scheuchl, B.: MEaSUREs Phase-Based Antarctica Ice Velocity Map, Version 1, NASA National Snow and Ice Data Center Distributed Active Archive Center, Boulder, Colorado USA [data set], (last access: 13 April 2021), 2019a. a, b, c, d, e

Mouginot, J., Rignot, E., and Scheuchl, B.: Continent-Wide, Interferometric SAR Phase, Mapping of Antarctic Ice Velocity, Geophys. Res. Lett., 46, 9710–9718,, 2019b. a, b, c, d

Mueller, R. D., Hattermann, T., Howard, S. L., and Padman, L.: Tidal influences on a future evolution of the Filchner–Ronne Ice Shelf cavity in the Weddell Sea, Antarctica, The Cryosphere, 12, 453–476,, 2018. a, b, c

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

Paden, J., Li, J., Leuschen, C., Rodriguez-Morales, F., and Hale, R.: IceBridge MCoRDS L1B Geolocated Radar Echo Strength Profiles, Version 2. IRMCR1B, NASA National Snow and Ice Data Center Distributed Active Archive Center, Boulder, Colorado, USA [data set], (last access: 25 April 2021), 2014, updated 2019. a, b, c, d

Rignot, E., Mouginot, J., and Scheuchl, B.: Ice Flow of the Antarctic Ice Sheet, Science, 333, 1427–1430,, 2011. a

Rignot, E., Jacobs, S., Mouginot, J., and Scheuchl, B.: Ice-Shelf Melting Around Antarctica, Science, 341, 266–270,, 2013. a, b, c, d

Rignot, E., Mouginot, J., Scheuchl, B., van den Broeke, M., van Wessem, M. J., and Morlighem, M.: Four decades of Antarctic Ice Sheet mass balance from 1979–2017, P. Natl. Acad. Sci. USA, 116, 1095–1103,, 2019. a

Scheuchl, B., Mouginot, J., and Rignot, E.: Ice velocity changes in the Ross and Ronne sectors observed using satellite radar data from 1997 and 2009, The Cryosphere, 6, 1019–1030,, 2012. a

Schoof, C.: Marine ice sheet stability, J. Fluid Mech., 698, 62–72,, 2012. a

Seroussi, H., Nowicki, S., Payne, A. J., Goelzer, H., Lipscomb, W. H., Abe-Ouchi, A., Agosta, C., Albrecht, T., Asay-Davis, X., Barthel, A., Calov, R., Cullather, R., Dumas, C., Galton-Fenzi, B. K., Gladstone, R., Golledge, N. R., Gregory, J. M., Greve, R., Hattermann, T., Hoffman, M. J., Humbert, A., Huybrechts, P., Jourdain, N. C., Kleiner, T., Larour, E., Leguy, G. R., Lowry, D. P., Little, C. M., Morlighem, M., Pattyn, F., Pelle, T., Price, S. F., Quiquet, A., Reese, R., Schlegel, N.-J., Shepherd, A., Simon, E., Smith, R. S., Straneo, F., Sun, S., Trusel, L. D., Van Breedam, J., van de Wal, R. S. W., Winkelmann, R., Zhao, C., Zhang, T., and Zwinger, T.: ISMIP6 Antarctica: a multi-model ensemble of the Antarctic ice sheet evolution over the 21st century, The Cryosphere, 14, 3033–3070,, 2020. 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,, 2019. a, b

van Wessem, J. M., Reijmer, C. H., Morlighem, M., Mouginot, J., Rignot, E., Medley, B., Joughin, I., Wouters, B., Depoorter, M. A., Bamber, J. L., Lenaerts, J. T. M., van de Berg, W. J., van den Broeke, M. R., and van Meijgaard, E.: Improved representation of East Antarctic surface mass balance in a regional atmospheric climate model, J. Glaciol., 60, 761–770,, 2014. a

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

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., Steinhage, D., and Humbert, A.: Repeated pRES measurements and derived Lagrangian basal melt rates at Filchner Ice Shelf between 2015/16 and 2016/17, PANGAEA [data set],, 2021. a

Short summary
Remote-sensing-derived basal melt rates of ice shelves are of great importance due to their capability to cover larger areas. We performed in situ measurements with a phase-sensitive radar on the southern Filchner Ice Shelf, showing moderate melt rates and low small-scale spatial variability. The comparison with remote-sensing-based melt rates revealed large differences caused by the estimation of vertical strain rates from remote sensing velocity fields that modern fields can overcome.