We present a 1000 km transect of phase-sensitive radar measurements of ice thickness, basal reflection strength, basal melting and ice-column deformation across the Ross Ice Shelf (RIS). Measurements were gathered at varying intervals in austral summer between 2015 and 2020, connecting the grounding line with the distant ice shelf front. We identified changing basal reflection strengths revealing a variety of basal conditions influenced by ice flow and by ice–ocean interaction at the ice base. Reflection strength is lower across the central RIS, while strong reflections in the near-front and near-grounding line regions correspond with higher basal melt rates, up to 0.47 ± 0.02 m a−1 in the north. Melting from atmospherically warmed surface water extends 150–170 km south of the RIS front. Melt rates up to 0.29 ± 0.03 m a−1 and 0.15 ± 0.03 m a−1 are observed near the grounding lines of the Whillans and Kamb Ice Stream, respectively. Although troublesome to compare directly, our surface-based observations generally agree with the basal melt pattern provided by satellite-based methods but provide a distinctly smoother pattern. Our work delivers a precise measurement of basal melt rates across the RIS, a rare insight that also provides an early 21st-century baseline.

Ice tongues are unconfined by land on their lateral margins and are sensitive to external forcing from the ocean. They are found sporadically around the Antarctic coast but are common in the western Ross Sea. Lateral flexure creates bending stresses within these ice tongues which is likely to contribute to their fragility and may restrict their spatial distribution. A three-year time series (2017–2020) of synthetic aperture radar interferometry (InSAR) and differential interferometry (DInSAR) is used to observe the lateral flexure of the 10 km long Erebus Ice Tongue as a result of ocean currents. The fast ice area around the ice tongue was mapped during the same period. It was found that when fast ice was absent (34.7% of the time), the average lateral movement of the ice tongue was twice as high (0.44 m) as when it was embedded in fast ice (0.19 m). A significant correlation (0.45) between flexure and tidal currents was found when fast ice was absent. An analytical model tuned to observations suggests that even without sea ice for stabilisation, the lateral bending stresses induced by the ocean are insufficient to cause calving without additional amplifying factors.

The only complete inventory of New Zealand glaciers was based on aerial photography starting in 1978. While there have been partial updates using 2002 and 2009 satellite data, most glaciers are still represented by the 1978 outlines in contemporary global glacier databases. The objective of this project is to establish an updated glacier inventory for New Zealand. We have used Landsat 8 OLI satellite imagery from February and March 2016 for delineating clean glaciers using a semi-automatic band ratio method and debris-covered glaciers using a maximum likelihood classification. The outlines have been checked against Sentinel-2 MSI data, which have a higher resolution. Manual post processing was necessary due to misclassifications (e.g. lakes, clouds), mapping in shadowed areas, and combining the clean and debris-covered parts into single glaciers. New Zealand glaciers cover an area of 794 ± 34 km2 in 2016 with a debris-covered area of 10%. Of the 2918 glaciers, seven glaciers are >10 km2 while 71% is <0.1 km2. The debris cover on those largest glaciers is >40%. Only 15 glaciers are located on the North Island. For a selection of glaciers, we were able to calculate the area reduction between the 1978 and 2016 inventories.

The Darwin–Hatherton Glacial system (DHGS) connects the East Antarctic Ice Sheet (EAIS) with the Ross Ice Shelf and is a key area for understanding past variations in ice thickness of surrounding ice masses. Here we present the first detailed measurements of ice thickness and grounding zone characteristics of the DHGS as well as new measurements of ice velocity. The results illustrate the changes that occur in glacier geometry and ice flux as ice flows from the polar plateau and into the Ross Ice Shelf. The ice discharge and the mean basal ice shelf melt for the first 8.5 km downstream of the grounding line amount to 0.24 ± 0.05 km3 a−1 and 0.3 ± 0.1 m a−1, respectively. As the ice begins to float, ice thickness decreases rapidly and basal terraces develop. Constructed maps of glacier geometry suggest that ice drainage from the EAIS into the Darwin Glacier occurs primarily through a deep subglacial canyon. By contrast, ice thins to <200 m at the head of the much slower flowing Hatherton Glacier. The glaciological field study establishes an improved basis for the interpretation of glacial drift sheets at the link between the EAIS and the Ross Ice Sheet.

A variety of data are used to investigate Larsen B, which is at present the northernmost section of the Larsen Ice Shelf, Antarctic Peninsula. Recently declassified U.S. Argon satellite photographs from 1963 Kosmos photographs from 1975, Landsat images from 1986,1988 and 1990, ERS-1/2 synthetic aperture radar images from 1992-97, Radarsat from 1998 and field surveys are used to analyze the areal extent, surface characteristics and dynamic behaviour of this ice-shelf section over more than three decades. Visible and radar imagery together with field observations are used synergistically to describe the ice-shelf morphology, including meltwater features and rifts. In contrast to the retreat of the ice-shelf sections in the north, Larsen B advanced steadily from 1963 to early 1995, when the area decreased significantly due to a major calving event. Analysis of different satellite images indicates that melting is proceeding further south in accordance with the regional warming trend. In addition? fracturing processes and rapid development of new rifts are observed, associated with recent acceleration of ice motion close to the front. All observations predict major calving events for this ice-shelf section in the near future.

The motion field of the northern Larsen Ice Shelf, Antarctic Peninsula, was analyzed using radar interferometry in combination with field measurements and finite-element model calculations. The ice shelf between Jason Peninsula and Seal Nunataks has been in steady retreat since January 1995. Model calculations suggest that the ice shelf is in a stage of irreversible retreat since the last calving events in summer 1998/99. The interferometric analysis is based on synthetic-aperture radar data of the tandem mission of European remote-sensing satellites ERS-1 and –2 in austral spring 1995. The phase contributions due to tidal motion were estimated from the vertical displacement at those parts of the grounding zone where the horizontal motion is close to zero" in order to separate the vertical and horizontal motion components over the ice shelf. Satellite-derived velocities compare well with the long-term field measurements along a transverse and a longitudinal profile. The real interferograms and synthetic interferograms, calculated from model velocities, show reasonable agreement over the main parts of the ice shelf, but differ in the boundary zones where the details are not resolved by the model.

The retreat of the Larsen B ice shelf, Antarctic Peninsula, and the collapse of its northern section are analyzed using satellite images acquired between January 1995 and May 2003. Over 1 week during March 2002, after a period of steady retreat since 1995, 2300 km2 of the ice shelf broke up into many small icebergs. This rapid collapse occurred at the end of an exceptionally warm summer, and after a multi-year period of decreasing surface net mass balance, ice thinning, flow acceleration and widening of rifts. The ice-shelf area decreased from 11 512 km2 in January 1995 to 3463 km2 in March 2002, and 2667 km2 in April 2003. ERS synthetic aperture radar (SAR) images were used to identify ice-shelf zones with different surface morphology, which generated icebergs of different sizes and shapes. The pattern of retreat and break-up, similar to that of Larsen A in 1995, suggests that fracturing enhanced by abundant surface melt played a key role. In addition, the recent changes of grounded and residual floating ice north of Larsen B are analyzed by means of Envisat advanced synthetic aperture radar (ASAR) images up to summer 2003, showing significant loss of grounded ice upstream of those ice-shelf sections which disintegrated in 1995 and 2002.

Changes of Larsen Ice Shelf, Antarctica, and the surrounding glaciers after its collapse in 1995 were investigated using satellite radar imagery, with emphasis on changes in the glaciers which previously nourished the ice shelf north of Seal Nunataks and now calve directly into the sea. The large glaciers retreated several kilometres inland of the previous grounding line. The velocity field of Drygalski Glacier, the largest glacier in this area, was mapped by means of interferograms derived from pairs of European Remote-sensing Satellite synthetic aperture radar images from 1995 and 1999. The main part of the glacier showed a significant acceleration of flow over these 4 years, with an increase of velocity up to three-fold at the terminus. Similar accelerations were observed by means of interferometry on several other grounded glaciers, suggesting that the removal of ice shelves could lead to an effect on eustatic sea level. For Larsen B, the northernmost surviving part of Larsen Ice Shelf, the retreat of the ice front to October 2000 is documented.

We characterize the basal mass balance of the Ekström Ice Shelf, Dronning Maud Land, Antarctica, using interferometrically derived surface velocities and ice thickness measurements from radio-echo sounding (RES). The surface velocities are based on data from European Remote-sensing Satellites-1 and -2 (ERS-1/2) during 1994–97. The ice thickness grid consists of 136 RES profiles acquired between 1996 and 2006. Mass fluxes are calculated along selected RES profiles where possible, to reduce uncertainties from ice thickness interpolation. Elsewhere large-scale mass fluxes are calculated using interpolated ice thickness data. Themass flux into the Ekström Ice Shelf from the main grounded drainage basins is estimated to be 3.19±0.4Gt a–1. The mass flux near the ice shelf front is 2.67±0.3Gt a–1. Assuming steady state, and based on the equation of continuity, we interpret the residual mass flux as a combined effect of snow accumulation and subglacial melting/refreezing. Using net snow accumulation rates from previous studies, we link the mass flux divergence in irregular-shaped polygons to processes beneath the ice shelf. The highest subglacial melt rates of ~1.1ma–1 are found near the grounding zone of two main inflow glaciers, and around the German station Neumayer III. The detection of unlikely refreezing in a small area ~15 km west of Neumayer III is attributed to both dataset inaccuracies and a (possibly past) violation of the steady-state assumption. In general, the method and input data allow mapping of the spatial distribution of basal melting and the results are in good agreement with several previous studies.

Based on ground-penetrating radar profiles, we analyze area-wide spatial and temporal characteristics of accumulation rate in the vicinity of the EPICA (European Project for Ice Coring in Antarctica) deep-drilling site in Dronning Maud Land, Antarctica (EDML). An area of 1600km2 is covered by 500km of radar profiles, organized in a star-like pattern with eight 20–25km legs and a 10km grid with 1–3km spacing, each pattern being centred on the EDML drilling location. Distributions of density, cumulative mass, age and the electromagnetic wave speed with depth are available from physical ice-core records. Nine internal reflection horizons are continuously tracked within the upper 110m of ice over the whole area, yielding a spatial picture of accumulation rate history over >1000 years. The mean accumulation rate over the last 153 years varies between 50 and 75 kgm–2 a–1 over 50 km perpendicular to the ice divide; the spatial average is ~61 kgm–2 a–1. This general pattern is overlain by small-scale variations of accumulation rate on the order of 10% of the mean. Maximum local gradients in accumulation rate are ~2–3 kgm–2 a–1 km–1, about five times the regional accumulation rate gradient. Comparison of topography and accumulation rate along a 20 km profile in the direction of the mean winds indicates that variations in accumulation rate over short distances are linked to surface undulations. Our results show that advected spatial variations of accumulation rate are on the same order and even exceed temporal changes over the investigated periods. Ice flow and upstream effects therefore might influence accumulation rates reconstructed from the EDML ice core.

The motion field of the northern Larsen Ice Shelf, Antarctic Peninsula, was analyzed, based on Landsat data from 1986 to 1989, Earth Resources Satellite (ERS) synthetic aperture radar (SAR) data from 1992 to 1997, and comparative field measurements along three transects. During this period the northern sections of the ice shelf showed steady retreat, which culminated in the disintegration of the two ice-shelf sections north of Seal Nunataks in January 1995. Velocities of these two sections were derived by cross-correlation, using SAR images of 1 year time intervals and Landsat images of 1-3 year intervals. A slight increase of velocity was observed as crevasses and rifts opened before the final disintegration. In addition, an interferometric motion analysis was carried out for the ice shelf around and south of Seal Nunataks, based on an image pair from the ERS-1/2 Tandem Mission in 1995. This analysis reveals a complex pattern of tidal flexure in the grounding zones, as well as rifting and shear zones on the ice shelf. In addition, the motion of the main input glaciers was derived.

We investigate snowpack properties at a site in west-central Greenland with ground-penetrating radar (GPR), supplemented by stratigraphic records from snow pits and shallow firn cores. GPR data were collected at a validation test site for CryoSat (T05 on the Expéditions Glaciologiques Internationales au Groenland (EGIG) line) over a 100 m × 100 m grid and along 1 km sections at frequencies of 500 and 800 MHz. Several internal reflection horizons (IRHs) down to a depth of 10 m were tracked. IRHs are usually related to ice-layer clusters in vertically bounded sequences that obtain their initial characteristics near the surface during the melt season. Warm conditions in the following melt season can change these characteristics by percolating meltwater. In cold conditions, smaller melt volumes at the surface can lead to faint IRHs. The absence of simple mechanisms for internal layer origin emphasizes the need for independent dating to reliably interpret remotely sensed radar data. Our GPR-derived depth of the 2003 summer surface of 1.48 m (measured in 2004) is confirmed by snow-pit observations. The distribution of IRH depths on a 1 km scale reveals a gradient of increasing accumulation to the northeast of about 5 cm w.e. km−1. We find that point measurements of accumulation in this area are representative only over several hundred metres, with uncertainties of about 15% of the spatial mean.

Antarctic tabular icebergs are important active components of the ice–ocean system. To investigate the relevance of inherent ice dynamics to iceberg evolution, we developed a numerical model based on the fundamental equations of ice-shelf flow and heat transfer, forced by environmental parameters of the ice–ocean–atmosphere system. Model experiments with idealized icebergs of constant density show that the strain thinning rate for a typical iceberg with a thickness of 250 m and a temperature of −15°C is about 1 m a−1. Sensitivity studies for different scenarios of environmental conditions confirmed the reliability of our model. A 5 year simulation of the evolution of iceberg A-38B yielded a mean decrease in thickness from 220 m to 106.3 m, 95% of which was caused by basal melting, 1% by surface melting and 4% by strain thinning. We found iceberg spreading decelerating by about 75%, and ice temperatures being strongly affected by progressive erosion of the relatively warm basal layers and warming in the uppermost part. According to the model results, basal melting is the primary cause of change of iceberg geometry during drift, whereas strain thinning is only relevant in cold areas where basal melting is low.

Tidal flexure in ice shelf grounding zones has been used extensively in the past to determine grounding line position and ice properties. Although the rheology of ice is viscoelastic at tidal loading frequencies, most modelling studies have assumed some form of linear elastic beam approximation to match observed flexure profiles. Here we use density, radar and DInSAR measurements in combination with full-Stokes viscoelastic modelling to investigate a range of additional controls on the flexure of the Southern McMurdo Ice Shelf. We find that inclusion of observed basal crevasses and density dependent ice stiffness can greatly alter the flexure profile and yet fitting a simple elastic beam model to that profile will still produce an excellent fit. Estimates of the effective Young's modulus derived by fitting flexure profiles are shown to vary by over 200% depending on whether these factors are included, even when the local thickness is well constrained. Conversely, estimates of the grounding line position are relatively insensitive to these considerations for the case of a steep bed slope in our study region. By fitting tidal amplitudes only, and ignoring phase information, elastic beam theory can provide a good fit to observations in a wide variety of situations. This should, however, not be taken as an indication that the underlying rheological assumptions are correct.

Ice-thickness measurements in Antarctic ice-shelf grounding zones are necessary for calculating the mass balance of individual catchments, but remain poorly constrained for most of the continent. We describe a new inverse modelling optimization approach to estimate ice thickness in the grounding zone of Antarctic outlet glaciers and ice shelves using spatial patterns of tide-induced flexure derived from differential interferometric synthetic aperture radar (InSAR). We demonstrate that the illposedness of the inverse formulation of the elastic-plate equations for bending can be controlled by regularization. In one dimension, the model recreates smooth, synthesized profiles of ice thickness from flexure information to within 1–2%. We test the method in two dimensions and validate it in the grounding zone of Beardmore Glacier, a major outlet glacier in the Transantarctic Mountains, using interferograms created from TerraSAR-X satellite imagery acquired in 2012. We compare our results with historic and modern ice-thickness data (radio-echo sounding from 1967 and ground-penetrating radar from 2010). We match both longitudinal and transverse thickness transects to within 50 m root-mean-square error using an effective Young’s modulus of 1.4 GPa. The highest accuracy is achieved close to the grounded ice boundary, where current estimates of thickness based on surface elevation measurements contain a systematic bias towards thicker ice.

We present a diagnostic glacier flowline model parameterized and constrained by new velocity data from ice-surface GPS installations and speckle tracking of TerraSAR-X satellite images, newly acquired airborne-radar data, and continental gridded datasets of topography and geothermal heat flux, in order to better understand two outlet glaciers of the East Antarctic ice sheet. Our observational data are employed as primary inputs to a modelling procedure that first calculates the basal thermal regime of each glacier, then iterates the basal sliding coefficient and deformation rate parameter until the fit of simulated to observed surface velocities is optimized. We find that the two glaciers have both frozen and thawed areas at their beds, facilitating partial sliding. Glacier flow arises from a balance between sliding and deformation that fluctuates along the length of each glacier, with the amount of sliding typically varying by up to two orders of magnitude but with deformation rates far more constant. Beardmore Glacier is warmer and faster-flowing than Skelton Glacier, but an up-glacier deepening bed at the grounding line, coupled with ice thicknesses close to flotation, lead us to infer a greater vulnerability of Skelton Glacier to grounding-line recession if affected by ocean-forced thinning and concomitant acceleration.

Grounding zones are vital to ice-sheet mass balance and its coupling to the global ocean circulation. Processes here determine the mass discharge from the grounded ice sheet, to the floating ice shelves. The response of this transition zone to tidal forcing has been described by both elastic and viscoelastic models. Here we examine the validity of these models for grounding zone flexure over tidal timescales using field data from the Southern McMurdo Ice Shelf (78° 15′S, 167° 7′E). Observations of tidal movement were carried out by simultaneous tiltmeter and GPS measurements along a profile across the grounding zone. Finite-element simulations covering a 64 d period reveal that the viscoelastic model fits best the observations using a Young's modulus of 1.6 GPa and a viscosity of 1013.7 Pa s (≈ 50.1 TPa s). We conclude that the elastic model is only well-constrained for tidal displacements >35% of the spring-tidal amplitude using a Young's modulus of 1.62 ± 0.69 GPa, but that a viscoelastic model is necessary to adequately capture tidal bending at amplitudes below this threshold. In grounding zones where bending stresses are greater than at the Southern McMurdo Ice Shelf or ice viscosity is lower, the threshold would be even higher.

The areal changes of the northern Larsen Ite Shelf (LIS), Antarctic Peninsula, between March 1986 and March 1997 have been analyzed, based on synthetic aperture radar images of the European remote-sensing satellites ERS-1 and ERS-2 and on Landsat images. This analysis is complemented by data on ice motion and surface mass balance which have been obtained during several field campaigns since the early 1980s. After a period of retreat, coinciding with atmospheric warming and with decreasing net accumulation at the surface due to melt losses, the two northernmost sections of LIS disintegrated almost completely within a few days in January 1995. Recent observations of the ice-shelf section north of Jason Peninsula, which is presently the northernmost section of LIS, show increased summer melt and intensification of the rifting processes, probably causing accelerated retreat of this section in the near future. The retreat and the disintegration event of LIS indicate high sensitivity of ice shelves to prolonged perturbations of the mass balance.