1. INTRODUCTION
Almost three quarters of the Antarctic grounded ice is discharged into floating ice shelves (Bindschadler and others, Reference Bindschadler2011). Thinning of ice shelves, as observed by Paolo and others (Reference Paolo, Fricker and Padman2015), imprints the current mass loss by reducing the restraint (i.e. buttressing), which the ice shelves exert on the inland flow, subsequently causing acceleration of tributary glaciers (Scambos and others, Reference Scambos, Bohlander, Shuman and Skvarca2004; Dupont and Alley, Reference Dupont and Alley2005; Joughin and others, Reference Joughin, Smith and Medley2014; Mouginot and others, Reference Mouginot, Rignot and Scheuchl2014; De Rydt and others, Reference De Rydt, Gudmundsson, Rott and Bamber2015). Ice-shelf buttressing originates from lateral friction of embayments and/or from rigid obstacles emerging from the bathymetry and attaching to the ice-shelf from below.
Localized grounded areas within ice shelves form ice rumples and ice rises (Smith, Reference Smith1986; Matsuoka and others, Reference Matsuoka2015). Ice rumples, on the one hand, decelerate ice flow, but are still overrun by the ice-shelf (e.g. Schmeltz and others, Reference Schmeltz, Rignot and MacAyeal2001; Humbert and Steinhage, Reference Humbert and Steinhage2011). Ice rises, on the other hand, develop a local flow regime isolated from the surrounding ice shelf (e.g. Drews and others, Reference Drews2015; references therein). Distinguishing an ice rise from an ice rumple based on the (non-)locality of the flow pattern is not always straightforward, especially for smaller features that extend over only a few square kilometres; in those cases we will use the generic term pinning point.
Because ice rises (Goldberg and others, Reference Goldberg, Holland and Schoof2009; Favier and others, Reference Favier2014; Favier and Pattyn, Reference Favier and Pattyn2015) and ice rumples (Favier and others, Reference Favier, Gagliardini, Durand and Zwinger2012) affect the dynamics of the grounding line and that of the upstream catchment area, it is important to better understand what role these features play in defining the buttressing strength of ice shelves.
Although small pinning points (e.g. inset; Fig. 1) are common features of ice shelves in Dronning Maud Land (e.g. Matsuoka and others, Reference Matsuoka2015), they are often not covered by airborne or ground-based radars and hence not resolved in Antarctic-wide ice thickness datasets such as Bedmap2 (Fretwell and others, Reference Fretwell2013; Fürst and others, Reference Fürst2015). Moreover, their surface elevation is often erroneous since satellite-altimeters tracks are too widely spaced. Antarctic-wide ice velocity maps, on the other hand, typically show the deceleration of ice flow near pinning points, but are not necessarily trustworthy at the pinning points themselves (for instance, regrounding necessitates a specific calibration for interferometric synthetic aperture radar).
Overview of the Roi Baudouin Ice Shelf (Dronning Maud Land, East Antarctica). The map locates ground-truth velocity data and the satellite scenes. WP, DIR and WRG indicate the western promontory, Derwael Ice Rise and the West Ragnhild Glacier, respectively. The blue rectangle locates the inset that displays the pinning point and the dark blue curve shows the grounding line from Bindschadler and others (Reference Bindschadler2011), except in the inset where the curve designates the break in slope at the pinning point. The background images are from Radarsat mosaic (Jezek and RAMP-Product-Team, Reference Jezek2002) and Landsat 8 of September 2014 (inset).
Ice thickness, surface elevation and ice velocity are the primary input variables for ice flow model initialisation. A typical approach for initialising ice flow models is to invert observed surface velocities for poorly known parameters such as basal friction and ice rigidity (e.g MacAyeal, Reference MacAyeal1993; Fürst and others, Reference Fürst2015). As a result, if pinning points are insufficiently resolved in the observational datasets, their impact on ice-shelf dynamics may erroneously propagate to the ice-shelf rigidity during the inversion, leading to erroneous predictions in prognostic ice sheet modelling. To address this point, we analyse a small pinning point in the Roi Baudouin Ice Shelf (RBIS), Dronning Maud Land, East Antarctica, by combining radar remote sensing and inverse modelling.
The western part of the Roi Baudouin Ice Shelf (26°E and 71°S), referred to hereafter as the RBIS, extends over 8500 km2 and presents several ice-shelf channels (i.e. along-flow features of diminished thickness, expressed both at the base and at the top surface of the ice shelf). The RBIS is fed by the West Ragnhild Glacier and constrained by an ice promontory on its western side, Derwael Ice Rise towards the East, and an 8.7 km2 pinning point at the ice-shelf front (Fig. 1). This pinning point is absent in the bathymetry, ice thickness and surface elevation of Bedmap2 (Fretwell and others, Reference Fretwell2013). The West Ragnhild Glacier drains ~10% of the ice in Dronning Maud Land (13–14 Gt a−1; Callens and others, Reference Callens2014) and is therefore among the major contributors to the ice discharge of Dronning Maud Land. Compared with ice shelves of other major outlet glaciers (e.g. Jutulstraulmen and Shirase Glacier), the RBIS flows more slowly (maximum speed ≈345 m a−1) and ice is discharged through a comparatively wide gate at the grounding line (Callens and others, Reference Callens2014).
Using interferometric synthetic aperture radar (InSAR) and speckle tracking, we first derive an ice-flow map, which is calibrated and validated with Global Navigation Satellite System (GNSS) ground-truth data (Section 2.1). We then modify the bed topography and ice thickness from Bedmap2 to account for the missing pinning point. This modified dataset is further referred to as mBedmap2. To assess the impact of the pinning point on ice flow models, we perform three inverse experiments based on three different scenarios (Section 2.2): (1) the Standard scenario using Antarctic-wide velocities (Rignot and others, Reference Rignot, Mouginot and Scheuchl2011b) and ice thickness from Bedmap2 (Fretwell and others, Reference Fretwell2013); (2) the Intermediate scenario combining Rignot and others (Reference Rignot, Mouginot and Scheuchl2011b)'s velocities with mBedmap2; (3) the High-resolution scenario using the high-resolution flow field derived here and in mBedmap2. We validate and discuss our results (Sections 3 and 4) and close by drawing conclusions with respect to the role of small pinning points on ice-shelf buttressing (Section 5).
2. DATA AND METHODS
2.1. Derivation of surface velocities
We derived surface velocities performing InSAR and speckle tracking on data from the European Remote Sensing satellites (ERS-1 and 2 from 1996) and from the Advanced Land Observation Satellite – Phased Array type L-band synthetic aperture radar (ALOS-PALSAR from 2010), respectively. Both techniques compare two co-registered images of the same area, acquired at different times. InSAR evaluates phase differences to derive ice flow, whereas speckle tracking tracks random, but deterministic, amplitude variations. Assuming steady state, velocities were mosaicked and calibrated with ground-control points collected in 2012/13 (Fig. 1).
2.1.1. Interferometric SAR using ERS-1/2
The individual frames of the ERS-1/2 scenes (Fig. 1; Table 1) were pre-processed to single-look complex images with Gamma's Modular SAR processor (Werner and others, Reference Werner, Wegmüller, Strozzi and Wiesmann2000). Pairs with 1-day time intervals were co-registered and differenced in phase forming interferograms (Goldstein and others, Reference Goldstein, Engelhardt, Kamb and Frolich1993) with lines of constant phase-differences (a. k. a. fringes). The fringe pattern depends on the satellite orbits, surface topography and ice flow. Orbital effects were removed using the precise orbits provided by Technical University Delft (Scharroo and Visser, Reference Scharroo and Visser1998) and the topographic contribution was cancelled using a DEM (Bamber and others, Reference Bamber, Gomez-Dans and Griggs2009). Phase differences were unwrapped using a minimum cost-flow algorithm (Costantini, Reference Costantini1998), transformed into a line-of-sight velocity map, and calibrated with the ground-control points marked in red in Figure 1. For the latter, the floating parts were discriminated using the dense fringe pattern caused by the tidal uplift in the grounding zone (e.g. Goldstein and others, Reference Goldstein, Engelhardt, Kamb and Frolich1993; Gray and others, Reference Gray2002; Fricker and others, Reference Fricker2009; Rignot and others, Reference Rignot, Mouginot and Scheuchl2011a). Calibration on the floating areas presupposes only vertical movement by tides (i.e. without tilting).
Characteristics of the satellite data; ΔT, λ and B ⊥ are the temporal baseline, the wavelength of the sensor and the perpendicular spatial baseline between the master and slave images, respectively. The satellite frames are shown in Figure 1
2.1.2. Speckle tracking using ALOS-PALSAR
Because the available InSAR data do not entirely cover the RBIS, we also used speckle tracking (Strozzi and others, Reference Strozzi, Luckman, Murray, Wegmuller and Werner2002; Luckman and others, Reference Luckman, Murray, Jiskoot, Pritchard and Strozzi2003; Werner and others, Reference Werner, Wegmuller, Strozzi and Wiesmann2005) on ALOS-PALSAR images of 2010 (Fig. 1; Table 1). Each pair of PALSAR images was treated separately to generate independent flow fields. Co-registration was done solely on Derwael Ice Rise, the western promontory and the pinning point because ice flow in those areas is <15 m a−1 (Drews, Reference Drews2015; Drews and others, Reference Drews2015), leading to a negligible displacement of the ice-sheet surface in the 46-day interval of the PALSAR images.
Offsets caused by ice flow were tracked every 12 and 36 pixels in range and azimuth, using a normalised cross-correlation of image patches (64 × 192 pixels), with a signal-to-noise-ratio threshold of 7 and an over-sampling factor of 2 (cf. Rankl and others, Reference Rankl, Kienholz and Braun2014). Range and azimuth offsets were georeferenced and converted to horizontal velocity fields with 125 m spacing. Mismatched pixels were removed (following Mouginot and others, Reference Mouginot, Scheuchl and Rignot2012) and the flow fields were finally smoothed by averaging within 9 × 9 diamond-shaped windows. Because speckle tracking is insensitive to tides, only one of the InSAR ground-control points was used for the entire area (filled red triangle; Fig. 1).
2.1.3. Mosaicking
The four PALSAR flow fields (Table 1) were mosaicked and then blended with the ERS velocities. To reduce cutting edges, velocity maps were feathered by applying a linear taper ranging from 0 to 1 in the overlapping areas. Data gaps on the western promontory and Derwael Ice Rise were filled with Rignot and others (Reference Rignot, Mouginot and Scheuchl2011b)'s flow field and feathered over 4.5 km.
2.1.4. On-site ice-flow measurements
To measure ice flow, 3 m long markers were placed on the ice-shelf and revisited the next year (Fig. 1). The 42 stakes were positioned using geodetic, dual-phase GNSS receivers with an antenna mounted on top of the stakes, measuring for at least 30 min. The data were post-processed using Precise Point Positioning (PPP) from the Canadian Spatial Reference System. Measurements were collected over a 3 a period. The wider-spaced markers (average distance 5 km) were measured in 2012/13, and a denser network crossing an ice-shelf channel (defined in Section 1) was measured in 2013/14 (Drews, Reference Drews2015). In 2012/13 some markers were tilted, leading to a maximum uncertainty of ~2 m in the horizontal position. This error is larger than the uncertainty acquainted with the PPP, which is typically within cm (Kouba and Héroux, Reference Kouba and Héroux2001). Because the average ice velocity in this area is ~265 m a−1, this error is <1%. We used the GNSS-derived velocities to calibrate and validate the satellite-based surface velocities (Sections 2.1 and 3.1).
Another set of ice-flow measurements stems from a triangulation network of 74 stakes measured multiple times with theodolites between 1965 and 1967 (Derwael, Reference Derwael2014; Fig. 1). We use this set of measurements to investigate the steady state assumption (Section 4.1).
2.2. Numerical modelling
Ice-sheet response to environmental conditions depends on the ice rheology and the conditions at the bed. These variables, however, are difficult to estimate experimentally and often parametrised in ice sheet and ice-shelf models. Inverse methods enable us to infer these model parameters from observed surface velocities and ice sheet thickness.
2.2.1. Ice-flow model
We used the adaptive mesh finite-volume ice sheet model BISICLES (http://BISICLES.lbl.gov) (Cornford and others, Reference Cornford2013), which solves the Schoof-Hindmarsh approximation (L1L2) of the full-Stokes equations on an adaptive horizontal two-dimensional (2-D) grid rendered by the Chombo adaptive mesh refinement toolkit. The L1L2 solution is based on the Shallow Shelf Approximation (SSA); vertical shearing is computed from the Shallow Ice Approximation (SIA) and represented within the effective strain-rate in Glen's flow law. BISICLES compares well with the Elmer/Ice ice sheet model, which solves the full Stokes set of equations. Both models yield similar results for ideal (Pattyn and others, Reference Pattyn and L2013) and real-case Antarctic outlet glaciers (Favier and others, Reference Favier2014), both in terms of ice dynamics and inversion.
In the BISICLES model, the ice rheology is given by Glen's flow law:
$${\bi S} = 2\phi \eta {\dot {\bi \varepsilon}} $$
where
S
is the deviatoric stress tensor;
${\dot {\bi \varepsilon}} $
is the strain-rate tensor; η is the effective viscosity depending on the temperature field of the ice mass (updated from Pattyn, Reference Pattyn2010) and on the effective strain rate (Wright and others (Reference Wright2014) and Cornford and others (Reference Cornford2013) provide a detailed description of η); ϕ is a stiffening factor that can be tuned to account for effects such as ice damage and anisotropy, along with uncertainties for ice temperature. Equation (1) is the constitutive relation used by the forward model to calculate ice flow velocities from the stress-balance equations.
For the grounded part of the ice sheet, a linear friction law governs resistance to basal sliding:
$${\tau _{\rm b}} = - C{u_{\rm b}}$$
where τ b is the basal traction, u b is the velocity of ice at the bottom interface and C is an estimate of the friction coefficient. The magnitude of C depends on different factors such as subglacial hydrology, the presence of subglacial till, etc. (MacAyeal and others, Reference MacAyeal, Bindschadler and Scambos1995).
2.2.2. Inverse method
To solve inverse problems, BISICLES uses the control method described by Cornford and others (Reference Cornford2015), which is comparable with those of MacAyeal (Reference MacAyeal1993) and Morlighem and others (Reference Morlighem2010). A nonlinear conjugate gradient method is used to seek a minimum of the cost function J, by expressing its gradient in terms of the adjoint equations constructed from the stress-balance equations, and defined as:
$$J = {J_{\rm m}} + {J_{\rm p}}$$
where J m is the misfit between observed and modelled velocities. The Tikhonov penalty function J p is:
$${J_{\rm p}} = {\lambda _C}J_C^{reg} + {\lambda _\phi} J_\phi ^{reg} $$
where λ
C
and λ
ϕ
are the Tikhonov parameters and
$J_C^{reg} $
and
$J_\phi ^{reg} $
, respectively, represent the spatial gradients of C and ϕ integrated over the domain (Cornford and others (Reference Cornford2015) provide the detailed formalism). A L-curve analysis (as in Fürst and others, Reference Fürst2015) was performed on the High-resolution scenario to calibrate λ
C
and λ
ϕ
, subsequently used in all the scenarios. Figure 2 presents the results of the L-curve analysis and highlights the chosen values of λ
C
= 5 × 102 Pa−2 m6 a−4 and λ
ϕ
= 5 × 109 m4 a−2.
L-curve analysis to select the Tikhonov parameters λ
ϕ
and λ
C
: (a) 3-D scatter plot of the model-data misfit J
m as a function of the regularisation terms
$J_C^{reg} $
and
$J_\phi ^{reg} $
. (b) 2-D cross section for variable λ
ϕ
and λ
C
fixed to 5 × 102 Pa−2 m6 a−4. (c) Reverse case where λ
ϕ
is fixed to 5 × 109 m4 a−2 and λ
C
varies. The units of J
m and
$J_C^{reg} $
are m4a−2 and Pa2 m−2 a2, respectively.
$J_\phi ^{reg} $
has no unit.
2.2.3. Inverse experiments
We inverted simultaneously the stiffening factor (ϕ) and the basal friction coefficient (C) for the three scenarios as described in Section 1. For the mBedmap2 dataset, used in the Intermediate and High-resolution scenarios, we modified the surface elevation and bathymetry of Bedmap2 by incorporating a pinning point with a Gaussian shape, with dimensions based on radar and GNSS profiles crossing the observed pinning point on the RBIS (Drews, Reference Drews2015). Figure 4a shows the modifications for an along-flow cross section. We imposed a stress-free upper surface, water pressure at the calving front and the observed surface velocities at the other lateral boundaries. On the considered domain, all the inversions were performed at a 1 km spatial gridding.
2.2.4. Sensitivity of the inversion to the initial guesses
The inverse method is sensitive to the initial guesses of ϕ and C. To dampen this sensitivity, we chose physically justified initial fields for C and ϕ, called hereafter basis fields. The basis initial friction coefficient was computed from the SIA assuming balance between the driving stress and the basal friction. The basis initial stiffening factor was 1.
The former assumption is valid for most of the Antarctic ice sheet (Morlighem and others, Reference Morlighem, Rignot, Mouginot, Seroussi and Larour2014), but may not be applicable for ice streams where basal friction is low. We addressed this uncertainty by multiplying the SIA-based friction by 0.6, 0.8 and 1.2. As the inversion was robust under those different fields (<1% change of the RMS difference for the observed and modelled velocities), the model was initialised with the SIA-based friction fields.
The basis initial stiffening factor of 1 corresponds to an ice-creep parameter from Pattyn's (Reference Pattyn2010) updated temperatures only. In the ice sheet, no strong simplifications were applied (apart from neglecting horizontal diffusion) to calculate those temperatures, so we assumed they have sufficient physical meaning to be kept as initial conditions. However, in the ice-shelf, the temperature calculation neglected horizontal advection and internal strain heating, according to the analytical solution due to Holland and Jenkins (Reference Holland and Jenkins1999), so we investigated the sensitivity to an initial stiffening factor on the ice-shelf, ranging from 0.6 to 1.4. Increasing the initial value of ϕ decreased the velocity misfits (~20% difference between 0.6 and 1.4), and stiffened the ice shelf. Ice shelves are generally stiffer than an isotropic reference, because of their strain-induced anisotropy (Ma and others, Reference Ma2010). However, since the corresponding results are quantitatively almost similar and did not give qualitative differences with the highest values of the initial ϕ, we chose to initialise the stiffening factor to 1 for the ice shelf.
3. RESULTS
3.1. Satellite-derived velocities on the RBIS and strain rates
Velocities in Figure 3a show how Derwael Ice Rise and the pinning point form a flux gate through which most of the ice of the RBIS is directed. Ice-shelf velocities average 205 m a−1 and reach 345 m a−1 at the front. Velocities decrease down to 0.5 m a−1 on the pinning point, which is comparable with the flow velocities near the ice divides of Derwael Ice Rise and the Western Promontory.
(a) Horizontal surface velocities and (b) associated shear-strain rates. The grounding line (Bindschadler and others, Reference Bindschadler2011) is overlaid in black and the pinning point is marked in white. The labels a–e locate the areas of interest discussed in the text. The background image is from Radarsat mosaic (Jezek and RAMP-Product-Team, Reference Jezek2002).
The flow field deviates by −5.2 ± 4.5 m a−1 and 7.5 ± 16.5 m a−1 from the GNSS-derived velocities and Derwael's (Reference Derwael2014) triangulation network of 1965–67, respectively. The deviations for both ground-truth datasets do not show any spatial trend. Differences between the ERS and PALSAR mosaics decrease more or less linearly from positive values (~25 m a−1) in the northern part of the overlapping area to negative values (~−25 m a−1) in the southernmost parts. Mean ± standard deviation of the difference field is −5.6 ± 12.9 m a−1. This systematic error likely reflects uncertainties in the baseline estimation of the ERS satellites and/or a systematic phase residual from SAR focusing (Drews and others, Reference Drews, Rack, Wesche and Helm2009). Strain rates were computed for a local coordinate system aligned with the main flow directions (y along- and x across-flow). Horizontal shear-strain rates reveal a number of features (Fig. 3b), which are discussed in Section 4.2.
3.2. Inversion results
The three inversion schemes all adequately reproduce the observed velocities: the RMSD between the modelled and observed velocities are 11.51, 11.25 and 12.33 m a−1 for the Standard, Intermediate and High-resolution scenarios on the whole domain and fall to 8.48, 8.22 and 6.38 m a−1, respectively, on the ice-shelf.
Excluding the vicinity of the pinning point, the inverted basal friction displays the same general pattern for the three scenarios (Fig. 5b, d and f). The slowly flowing western promontory and Derwael Ice Rise are sticky at the bed (C peaks at around 6200 and 5900 Pa m−1a, respectively), whereas the fastest trunk of the West Ragnhild Glacier, upstream of the grounding line, is slippery with a basal friction coefficient ten times lower. The Standard, Intermediate, and High-resolution scenarios differ at the pinning point (insets in Fig. 5b, d and f). Because the pinning point is absent from the Standard scenario (Fig. 4a), the corresponding area is frictionless. In the Intermediate experiment, the pinning point is more slippery (mean C ≈ 485 Pa m−1a, for the pinning point) than the bed beneath the main stream of West Ragnhild Glacier, while in the High-resolution scenario it is on average 3.5 times stickier (mean C ≈ 1700 Pa m−1a, for the pinning point) than in the Intermediate scenario.
Profiles along the straight line shown in Figure 5b, d and f. (a) Bed and surface elevations relative to the EIGEN-GL04 geoid. (b) Comparison of the high-resolution flow field and velocities from Rignot and others (Reference Rignot, Mouginot and Scheuchl2011b). The red cross shows the ground-truth velocity measured on the pinning point. (c) Inverted friction coefficient.
All the experiments display four similar features in the stiffening factor ϕ (labelled 1–4 in Fig. 5a, c and e). Two wide patches of stiffened ice appear in the centre of the ice-shelf: one, a few km downstream of the grounding line (label 1 in Fig. 5; max (ϕ) ≃ 2.50) and the other, dozens of kilometres upstream of the calving front (label 4 – max (ϕ) ≃ 2.50). The most softened ice appears directly upstream of Derwael Ice Rise (label 2 – min (ϕ) ≃ 0.14) and, to a lesser extent, upstream of the pinning point (label 3 – min (ϕ) ≃ 0.40), in its shear band (Fig. 3b; label c). The Standard experiment stands out from the Intermediate and High-resolution cases and exhibits another patch with the most stiffened ice (label 5 – max (ϕ) ≈ 2.75). Compared with the Standard scenario, the ice-shelf is, on average, 4.0 and 4.6% softer in the Intermediate and High-resolution simulations, respectively.
(left: a, c and e) Stiffening factor ϕ and (right: b, d and f) friction coefficient C inverted for the Standard (top: a and b), Intermediate (middle: c and d) and High-resolution (bottom: e and f) scenarios. The grounding line (Bindschadler and others, Reference Bindschadler2011) is overlaid in black. A profile upstream/on top of the pinning point is shown in dotted, dashed and plain lines for the Standard, Intermediate and High-resolution scenarios, respectively. Theses profiles are plotted in Figure 4. The labels 1–5 locate areas of interest. Background images are from Radarsat Mosaic (Jezek and RAMP-Product-Team, Reference Jezek2002) and Landsat 8 of September 2014 (inset).
4. DISCUSSION
4.1. Is the Roi Baudouin Ice Shelf in steady state?
We mosaicked two different flow fields (1996 and 2010), both of which were calibrated on the ice-shelf with the same ground-control point collected in 2012/13 (Section 2.1). This approach holds if the RBIS is in steady state. To validate this hypothesis, we compared the mosaicked flow field with 74 ground-truth measurements collected in 1965–67 (Derwael, Reference Derwael2014; Section 2.1; Fig. 1). The deviations are not larger than the uncertainty of our high-resolution flow field and we conclude that the RBIS has not undergone prominent changes in average ice flow over the last five decades. This justifies our steady-state assumption in Section 2.1.
4.2. Velocities and strain rates
Our flow field (Fig. 3a) shows that ice virtually stagnates on the pinning point with a flow regime resembling more an ice rise than that of an ice rumple. This implies that the pinning point is stickier than what is suggested by the Antarctic-wide velocities (~150 m a−1 vs. 5 m a−1 for the ground truth; Fig. 4b) and that size alone is an insufficient criteria to classify ice rumples. Over a distance of only 5 km, ice from upstream is decelerated from 180 m a−1 to 0.5 m a−1. The horizontal shear band upstream of the pinning point (label c; Fig. 3b) separates the stagnant part on the west of the ice-shelf from the fast flowing centre and attests to the importance of buttressing on the Western part of the RBIS. Although it is difficult to discriminate the buttressing of the western promontory from the contribution of the pinning point, the alignment between the horizontal shear band and the location of the pinning point suggests a connection. Preliminary transient simulations indicate that unpinning repositions this shear band westwards for about two-thirds of the ice shelf upstream of the pinning point. In addition to horizontal shearing at the lateral boundaries of the RBIS (labels a and b; Fig. 3b), other remarkable features of the High-resolution flow field are the enhanced velocity gradients in the vicinity of ice-shelf channels (labels d and e; Fig. 3b). At the channels locations, the shearing reproduces the same pattern as that from the synthetic melting case of Drews (Reference Drews2015).
4.3. Inversion for the stiffening factor and basal friction
The High-resolution scenario has the highest averaged misfit between the observed and model velocities. However, the magnitudes of the misfits are primarily determined by the slow-flowing areas of Derwael Ice Rise and the Western Promontory, where the observations are correspondingly poor. Zooming in on the ice-shelf, the misfit of the High-resolution scenario is more than 20% lower compared with the two other scenarios.
The basal friction coefficient at the outlet of the West Ragnhild Glacier agrees with previous results (Callens and others, Reference Callens2014), based on a different inversion method. Callens and others (Reference Callens2014) showed that the downstream part of the West Ragnhild Glacier bed is wide, flat and covered by wet sediments, which induces basal sliding, thus a low coefficient C. At that place, the High-resolution scenario indicates a more slippery bed because of slightly higher (+50 m a−1) surface velocities between our flow field and that of Rignot and others (Reference Rignot, Mouginot and Scheuchl2011b). The sticky bed on Derwael Ice Rise and the western promontory concurs with low velocities and the negligible basal sliding in those areas.
The stiffening factor is used for fitting the modelled to the observed velocities. This means that it accounts for both physical approximations in the model and imperfections in the input data (e.g. missing factors such as damage, anisotropy and uncertainties in ice temperature); it is, however, not straightforward to distinguish between those different mechanisms. We will discuss hereafter the five patches of stiffened and softened ice (labels 1–5; Fig. 5a, c and e).
Patch 1 shows stiffened ice downstream from the grounding line, which we attribute to thermal effects. Our model uses the temperature field from Pattyn (Reference Pattyn2010), which neither accounts for horizontal advection nor horizontal shear-strain heating in the ice shelf (Holland and Jenkins, Reference Holland and Jenkins1999). This means that the local stiffening at least partially accounts for the cold ice that is advected from the tributary glacier to the ice shelf. The softened ice in patches 2 and 3 can be explained along the same lines, since in these areas (close to Derwael Ice Rise and upstream of the pinning point, respectively) horizontal strain heating probably occurs; moreover these areas are crevassed. The inverse model does not reveal weak margins of softer ice on the side of the ice stream, in contrast with what has been shown for faster ice shelves (e.g. Favier and others, Reference Favier2014; Larour and others, Reference Larour2014; Fürst and others, Reference Fürst2015). We have no conclusive evidence for patch 4, and a full thermomechanical coupling will be required to draw more solid conclusions.
Patches 1–4 occur in all model runs, whereas patch 5 only appears in the Standard scenario. In this scenario, the ice slows down a few kilometres upstream of the pinning point but the latter is ignored in the bathymetry. To match the observed velocities, the inverse method must stiffen the ice in patch 5, while in the two other scenarios the inverse method can also increase the friction at the pinning point. The value of the surface velocity at the pinning point is of secondary importance : the pinning point is stickier and the ice is less stiffened in the High-resolution scenario than in the Intermediate scenario. Changing the boundary condition of the ice leads to qualitative differences, while using different velocity values only quantitatively changes the results.
Predictive simulations of the Antarctic ice sheet generally determine initial state and parameters by inverting observational data. Therefore, unresolved pinning points in the latter can result in spurious ice dynamics for transient simulations of unpinning. In the case of the Standard scenario, for example, the unpinned RBIS is erroneously buttressed by the stiff ice around the pinning point. This is not the case for the other two scenarios, i.e. softer ice near the pinning point presumably leads to higher flow velocities in unpinning experiments. Our results show that it is more important to incorporate the pinning point in the bathymetry than to resolve the velocities at the pinning point (and correspondingly the inverted basal friction coefficient), even though both may influence transient simulations of the RBIS.
Pinning points such as that presented here are features small enough to be neglected/undersampled in observational datasets but important enough to impact the ice dynamics of their surroundings. Few kilometres-wide features can have larger-scale implications, although the gridding of Antarctic-wide datasets is too coarse to capture them in a modelling framework. For instance, the pinning point of the RBIS only extends over ~8 pixels when the spacing is 1 km (vs. 458 with 125 m) and can thus easily be missed out. The results derived here accord with the findings of Durand and others (Reference Durand, Gagliardini, Favier, Zwinger and le Meur2011), demonstrating that it is important to better resolve pinning points in the bathymetry around the Antarctic coastline.
5. CONCLUSION
We derived a flow field for the Roi Baudouin Ice Shelf, combining satellite-based InSAR and speckle tracking techniques. The velocities deviate by −5.2 ± 4.5 m a−1 compared with a set of 37 independent ground-control points. Comparison of our flow field with ground-truth measurements collected in 1965–67 provides no evidence for ice flow changes on decadal time scales. With its 125 m gridding, the derived flow field currently has the highest resolution available for the RBIS and shows that ice on top of a 8.7 km2 pinning point is virtually stagnant, contrary to what is indicated in other datasets. Despite its small size, this feature corresponds to a shear zone across the entire ice shelf.
Using different input scenarios of ice geometries/velocities we inverted for the stiffening factor (determining the ice rheology) and for the basal friction coefficient (determining the basal drag) using the BISICLES model. Surface velocities were adequately reproduced in all scenarios even in the one where the pinning point was omitted. However, including the pinning point in the bathymetry/ice thickness led to the best results on the ice shelf. The inversion for both ice rheology and basal conditions is the most sensitive to presence/absence of the pinning point in ice thickness and bathymetry. Neglecting basal drag at the pinning point is compensated by stiffened ice in its surrounding and overestimating the surface velocities at the pinning point results in too low basal friction coefficients. Both effects can cause errors in transient simulations. Our results emphasize that ice-flow models require high-resolution for observational datasets (especially bed topography) in coastal areas to fully capture the role of pinning points and their impact on the buttressing strength of ice shelves.
ACKNOWLEDGEMENTS
This paper forms a contribution to the Belgian Research Programme on the Antarctic (Belgian Federal Science Policy Office), Project SD/SA/06A Constraining ice mass changes in Antarctica (IceCon) as well as the FNRS-FRFC (Fonds de la Recherche Scientifique) project IDyRA. We thank the InBev Baillet Latour Antarctica Fellowship for partially financing fieldwork and the International Polar Foundation for providing all required logistics in the field. ERS1/2 and ALOS data were provided by the European Space Agency (Project C1P.10754) and further processed with the help of Wolfgang Rack at the Gateway Antarctica, Canterbury University, New Zealand. Ms. S. Berger is supported by a FRS-FNRS (Fonds de la Recherche Scientifique) PhD fellowship. Support and development of BISICLES from S. Cornford and D. Martin is greatly appreciated. We finally wish to thank A. Khazendar and an anonymous reviewer for constructive comments on this manuscript.
Open access
