Abstract
Wilkes Land and Totten Glacier (TG) in East Antarctica (EA) have been losing ice mass significantly since 1989. There is a lack of knowledge of long-term mass balance in the region which hinders the estimation of its contribution to global sea level rise. Here we show that this acceleration trend in TG has occurred since the 1960s. We reconstruct ice flow velocity fields of 1963â1989 in TG from the first-generation satellite images of ARGON and Landsat-1&4, and build a five decade-long record of ice dynamics. We find a persistent long-term ice discharge rate of 68â±â1 Gt/y and an acceleration of 0.17â±â0.02 Gt/y2 from 1963 to 2018, making TG the greatest contributor to global sea level rise in EA. We attribute the long-term acceleration near grounding line from 1963 to 2018 to basal melting likely induced by warm modified Circumpolar Deep Water. The speed up in shelf front during 1973â1989 was caused by a large calving front retreat. As the current trend continues, intensified monitoring in the TG region is recommended in the next decades.
Similar content being viewed by others
Introduction
The Antarctica Ice Sheet (AIS) has been losing ice mass and contributing to global sea level rise (SLR) at an accelerated pace over the last few decades1. This has been determined by estimates from altimetric, gravimetric, and optical and SAR imaging satellite observations, with the main mass loss from West Antarctica (WA) and the Antarctic Peninsula (AP)1,2,3. A reconciled solution shows a mass change rate with a large uncertainty (5â±â46 Gt/y) in East Antarctica (EA)4. However, a recent finding indicates that Wilkes Land sets itself apart from this overall trend in EA by showing a high mass loss rate of 51â±â13 Gt/y (refs. 1, 5, 6). Furthermore, studies1,7 have noted that the Totten Glacier (TG), fed by the largest drainage basin in Wilkes Land8, has undergone consistent mass loss since 1989, losing a total of ~175 Gt (or 0.5âmm SLR) at a rate of 7â±â2 Gt/y from 1989 to 2015. Limited by the quality and availability of satellite data, there is generally a lack of mass balance estimates for the AIS and this region prior to 1989 (ref. 9). Consequently, the lack of long-term mass balance knowledge hinders our ability to understand earlier glaciological responses in EA and the entire AIS to climate change.
Totten Glacier is marine based, mostly grounded below sea level and the bedrock slopes downwards away from the ocean10. Thus, the glacier is subject to marine ice sheet instability (MISI) and vulnerable to potential rapid ice sheet retreat due to the intrusion of the relatively warm modified Circumpolar Deep Water (mCDW) from the Southern Ocean11,12,13,14,15. With a catchment of ~570,000âkm2 (ref. 12), the TG basin has the capacity to raise the global sea level by 3.9âm (ref. 16) if completely melted. Its ice shelf has been experiencing increased basal melting at a rate from ~9.1 m/y (1992â2007) (ref. 13) to ~18 m/y (2005â2011) (ref. 17). The grounding line retreated up to ~3 km during 1996â2013 (ref. 16). Between 1989 and 2015, the ice flow acceleration was found to be as large as 18% (refs. 7, 18). However, a longer satellite observation record is needed to put the recent variability in a better context. Although recent satellite data of InSAR, altimetric and gravimetric sensors are not available before 1989, early optical images collected on films by the ARGON intelligence satellites during the 1960s (ref. 19) and the Landsat satellites in the 1970s and 1980s can be used to recover historical glacier topography and ice flow velocity fields9,20,21 with comparable uncertainties of earlier and recent velocity products1,22,23. Furthermore, these ice velocity data along with other supporting data will allow us to estimate mass balance during historical periods at both glacial and continental scales through the inputâoutput (IO) method1,5.
To quantify the early state of ice dynamics and mass balance in the TG region, we developed and applied an innovative method of hierarchical network densification to systematically map Antarctic ice velocity from historical optical satellite images. Here, we used images in the TG region acquired by the first-generation satellites of ARGON and early Landsat missions in order to reconstruct the ice velocity fields of multiple periods from 1963 to 1989. We found ice flow acceleration and high-level ice discharge during the period that are linked to a large calving front retreat during 1973â1985 and increased basal melting of the ice shelf. Finally, using our historical estimates alongside recent results, we built a nearly 60-year record and proved that there has been a persistent, long-term ice discharge increase driven by continued ice shelf basal melting over 6 decades in this marine-based region, making TG the glacier with the highest mass loss in EA. We suggest that the current acceleration and mass loss process in the TG region had already started by the 1960s.
Results
Reconstructed velocity maps
A historical velocity map of the TG region from 1963 to 1989 is generated (Fig. 1a), reconstructing the velocity field using 71,442 velocity vectors derived from ARGON and Landsat MSS and TM satellite images. This overall map is produced by weight-averaging three periodic velocity maps with shorter timespans (1963â1973, 1973â1989, and 1989) which are used to analyze ice flow dynamics and estimate the mass balance of the drainage basin over the 26-year period (Fig. S1). Displacements of the velocity vectors were computed by tracking ice flow features in image pairs using the innovative hierarchical matching and network densification method developed for glacier mapping using historical Antarctic images21,23. Unique processing techniques for overcoming difficulties inherent in the historical images were developed and applied, including film deformation, large format lens distortion, and georeferencing errors exceeding 10âkm (refs. 21, 24,25,26). For instance, we adopted a semi-automatic algorithm for recognition and measurement of damaged fiducial marks on ARGON films; exterior orientation parameters were initially estimated from the ephemeris data and ground features, and then refined through a bundle adjustment procedure. Consequently, uncertainties of the produced velocity maps range from 4âm to 79âm (Table S1), which were estimated from the orthorectification error of the images, feature identification and matching errors, and timespan according to the error propagation law described in the âMethodsâ section.
a Ice flow velocity field from 1963 to 1989 with glacier and ice shelf centerlines marked AAâ² and BBâ², respectively, Box 1 in shelf front (~50âkm from shelf front), Box 2 near grounding line (~4âkm from grounding line) and shelf fronts showing the largest retreat during 1973â1989. The gray dashed line is the dividing line between the velocity maps generated in this study and the regional velocity map29, 30. b Velocity difference map (1973â1989 minus 1963â1973), illustrating an ice shelf-wide increase of 60â±â11âm/y, due to the shelf front retreat. Red denotes acceleration and blue denotes deceleration. c Velocity difference map (1989 minus 1973â1989), showing that the ice velocity in 1989 remained high, with a difference of only â15â±â55âm/y. d Ice velocity changes from 1963 to 1989 along centerline AAâ² of the main trunk of TG. e Ice velocity changes from 1963 to 1989 along centerline BBâ² of the eastern tributary glacier. The backgrounds of aâc are from the Landsat image mosaic of Antarctica (LIMA mosaic)69. The ice shelf fronts in a are digitized from the ARGON and Landsat orthoimages. The grounding line is from ref. 70. The Passive Shelf Ice (PSI) boundary (green line) in a is from ref. 33. Bathymetry data (100âm isobath) on the continental shelf in a are from ref. 59. Source data are provided in the Source Data file.
The long timespan velocity maps of 1963â1973 and 1973â1989 have the benefit of covering the grounded regions and most of the ice shelf at a low uncertainty of 4â17 m/y (Table S1), especially the low velocity areas (<10â80 m/y), such as the Law Dome (Figs. S1a, b and 1a). However, the fast-flow shelf front area (up to ~2300 m/y) cannot be mapped due to lost tracking features that were calved into the ocean over timespans of up to 16 years. This data gap was effectively filled by the 8-month velocity map of 1989 (Fig. S1c) at an uncertainty of 79 m/y. The overall velocity map of the TG from 1963 to 1989 with a grid spacing of 500âm (Fig. 1a) was generated by weight-averaging the velocity points from the three periodic maps within a natural neighbor27,28. The extended area (outside the gray dashed line in Fig. 1a) is covered by using a regional velocity map29,30.
Overestimation (OE) of velocity due to ice flow acceleration was found in long timespan velocity maps (e.g., over 3 years)9,31. If not corrected, the implication is that the overestimated historical velocity, when combined with the recent unbiased velocity, may cause an underestimation of the long-term velocity change used to study climate impact on a decadal scale. We developed an innovative Lagrangian velocity-based method for OE correction without the use of field observations or additional image data9. The OEs in this region are mainly distributed along the main trunk, in the grounding zone, and at the shelf front (Fig. S2). The average of the OE corrections is 50â±â39 m/y for the velocity maps of 1963â1973 and 1973â1989. The velocity map of 1989 did not need OE correction because of the short timespan of 8 months.
Ice flow kinematics
The velocity of grounded ice in the slow-flowing area (<80 m/y, Fig. 1a) of the TG changed little, on average 6â±â16 m/y (~0%), from 1963 to 1989 as indicated by the two velocity difference maps (Fig. 1b, c), indicating insignificant variation of discharge from the upper part of the grounded ice which would contribute to sea-level rise. However, acceleration was mainly found on the ice shelf, with an average velocity increase of 60â±â11 m/y (~7%) from 1963â1973 to 1973â1989 for the floating ice. An increase of up to 135â±â9 m/y (~10%) occurred near the shelf front during the same period. This shelf-wide acceleration was in the same time range as a large calving front retreat between 1973 and 1989 (Fig. 1a), with an area loss of ~645 km2 (12%). We further narrowed the retreat period to 1973â1985 using images of the Advanced Very High-Resolution Radiometer (AVHRR), during which the retreat may have resulted from a major calving event or a series of smaller calving activities initiated by numerous rifts32. This retreat is found to be the largest for the Totten Ice Shelf (TIS) since the region was photographed by the earliest ARGON satellite in 1963. In particular, a large portion of the lost area crossed the âpassive shelf-iceâ (PSI) boundary33 in the western shelf margin (Fig. 1a). Loss of this âactive shelf-iceâ area (104 km2 or ~2%) caused velocity responses throughout the ice shelf32,34. Ice velocity on the ice shelf remained high in 1989, at least 4 years after the retreat, without significant changes from the 1973â1989 map (Fig. 1c). Furthermore, the acceleration from 1963â1973 to 1973â1989 is clear on the ice shelf and in a 20 km zone upstream and downstream from the grounding line along centerline AAâ² (Fig. 1a, d). Similarly, velocity increased along centerline BBâ² of the eastern tributary glacier (Fig. 1a) by up to 207â±â55 m/y (~14%) near the shelf front (Fig. 1e). However, the velocity on grounded ice and near grounding line showed no significant changes, indicating no variations in discharge and direct contribution to sea-level rise from this part of the glacier. It is also suggested that the strong calving front retreat in the western margin may have not affected the ice dynamics in the tributary on the other side of the ice shelf.
To link the reconstructed ice dynamics during our study period with recent acceleration in the TG, we built a nearly 6-decade record of velocity in two locations, one near the shelf front and the other near grounding line (Boxes 1 and 2 in Fig. 1a), using our historical velocity maps and recent velocity products35,36,37. We calculate changes in percentage during our study period relative to the long-term average from 1963 to 2018. The velocity near the shelf front was initially high, at 1328â±â9 m/y during 1963â1973 (Fig. 2a), which is 28â±â19 m/y (~2%) above the long-term average velocity. Subsequent acceleration due to the large calving front retreat caused an increase of 113â±â9 m/y (~11%), making the velocity of 1973â1989 the highest over the 55 years. After 1989, the velocity fluctuated mostly within a range of ~50 m/y. Note that the acceleration was, therefore, not constant near the shelf front throughout the time period. The long-term ice velocity trend (Fig. 2a) appears to be consistent with that of the area loss and gain caused by calving front retreat and advance (Fig. 2b), with R2â=â0.8 estimated between the two time series using timespan as weights (Table S3). Therefore, we suggest that the long-term ice velocity trend near the shelf front of TIS has mainly been controlled by calving activities. The ice velocity during our study period was historically high due to the largest calving front retreat between 1973 and 1985.
a Ice velocity near the shelf front (Box 1 in Fig. 1a) during the study period of 1963â1989 shows acceleration caused by the largest calving front retreat and the highest velocity in nearly 6 decades; b shelf area change (loss or gain) due to calving front retreat or advance indicates correlation with the ice velocity trend near the shelf front in a (R2â=â0.8); c ice velocity near grounding line (Box 2 in Fig. 1a) from 1963 to 2018 shows a long-term increasing trend (black dashed line) consistent with basal melt rate in d, in addition to the acceleration induced by the calving front retreat during 1973â1989; d simulated area-averaged basal melt rate of TIS from 1960 to 2007 (ref. 38) with an increasing trend (black dashed line). Source data are provided in the Source Data file.
Conversely, the velocity near grounding line (Box 2 in Fig. 1a) was initially low during 1963â1973, 123â±â21 m/y (~15%) below the long-term average velocity (Fig. 2c). There is a long-term increasing trend of velocity over the period of 1963â2018 which is consistent with that of the modeled melt rate from 1960 to 2007 (ref. 38), simulated with COREv2 forcing-39 (Fig. 2d). The simulated melt rate shows lower basal melting of ~5 to 6 m/y. Thereafter melt rate estimated from multi-mission altimetry observations increased to 11.5â±â2 m/y from 1994 to 2018 (ref. 40) and 17.9â±â1.2 m/y from 2005 to 2011 (ref. 17). We reconstructed a velocity time series from that in Fig. 2c by subtracting its long-term linear trend (black dashed line) estimated by a regression, which then shows a timespan weighted correlation of R2â=â0.6 with the shelf front area change in Fig. 2b. Thus, the acceleration during 1973â1989 in Fig. 2c, with a velocity increase of 124â±â8 m/y (~0%), appears to be caused by the largest calving front retreat during the same period, whose area loss intruded into the âactiveâ shelf ice region inside the PSI boundary and influenced ice kinematics in the grounding zone. The ice velocity declines as elsewhere on the ice shelf in 1989, but then picks back up again. This recent acceleration after 1989, also reported in refs. 1, 5, 6, 16, was induced by ice shelf basal melting caused by intrusion of mCDW from the continental shelf into the ice shelf cavity through the trough on the seafloor7,8,41. The combined long-term melt rates from 1960 to 2018 based on modeling and satellite observations suggest that ice shelf basal melting may have existed in TIS as early as 1960. Here we demonstrate that the long-term acceleration trend in the grounding region induced mainly by ice shelf basal melting (Fig. 2c) started in 1963, 26 years earlier than reported.
Mass balance
Due to the ice flow acceleration in the grounding zone from 1963 to 1989 caused by both calving front retreat and basal melting, ice discharge in the TG was at a high level, close to the reference surface mass balance (SMB) that is the long-term average SMB from 1979 to 2016 (Fig. 3a). Consequently, a total of 1774â±â16âGt (66â±â1 Gt/y) ice mass was discharged across the grounding line during the period (Fig. 3b). However, SMB in this period was higher, on average ~3âGt/y above the reference SMB (Fig. 3a). The excessive cumulative SMB of 1904â±â22 Gt (71â±â1 Gt/y) in the TG basin overrode the discharge and resulted in a total net mass gain of 130â±â27 Gt (5â±â1 Gt/y) (Fig. 3b). However, as the SMB decreased and discharge increased in the TG basin, the net mass gain reduced to about zero at the end of the period (Table 1). Therefore, during 1963â1989, the TG basin experienced a high level of ice discharge, transiting from positive to negative mass balance.
a Ice discharge with uncertainty (blue line and shaded margin) and SMB from RACMO2.3 p2 (red line and shaded margin) in the TG basin during 1963â1989. Reference SMB (black dashed line) is averaged during 1979â2016. Average SMB (red dashed line) is calculated over the period of 1979â1989. Zoom inset âa1â shows details for 1972â1974. b Cumulative discharge (blue line), cumulative SMB (red line) and total mass change (net mass gain, black line) for the same period. Zoom inset âb1â shows net mass gain with uncertainty (shaded margin) of the last two years. Source data are provided in the Source Data file.
Discussion
Bindoff et al.42 and Williams et al.43 observed mCDW on the continental shelf adjacent to TIS in the austral summer of 1995â1996 and austral winter of 2007, respectively. Intrusion of the on-shelf warm mCDW water towards the TIS resulted in circulation of the cyclonic gyre and westward Antarctic Slope Current42,44, modulated by polynya activities near the calving front45. A major entrance, a ~5âkm wide trough, which allows the mCDW water access to the ice shelf cavity, was found in the western part of the TIS calving front using inverted bathymetry from airborne gravity and magnetic data collected between 2008 and 2012 (ref. 8). The trough entrance was more precisely described, by shipboard bathymetry data acquired in January 2015, as being ~10âkm wide, ~600âm deep, and at a depth of ~1097âm (maximum) from the sea surface41. Li et al.7 suggested that ice flow acceleration in the TG between 1989 and 2015 was linked to ocean temperature changes. Furthermore, Greene et al.14 and Gwyther et al.38 indicated that the troughs underneath the ice shelf provided effective pathways for the mCDW water to intrude into the cavity and grounding zone, causing increased ice shelf basal melting and thinning at a rate from ~9 m/y for 1992â2008 (refs. 13, 46) to ~18 m/y for 2005â2011 (ref. 17).
Despite the spatial and temporal sparsity of oceanographic observations in our study area, we are able to use six conductivity, temperature, and depth (CTD) profiles (IDs 1â6) in an extended region (Figure S3), mostly collected during austral summer of 1995â1996 (ref. 42), from the World Ocean Database (WOD) supported by the NOAA Climate and Global Change Program (https://www.ncei.noaa.gov/). Limited by fast ice and grounded icebergs, the CTD profiles were not deployed close to the shelf front. They reveal the presence of mCDW of relatively high temperature (>~0â°C) and salinity (>~34.5 PSU) at ~300â600 m depth, near the continental shelf break and in front of TIS (Fig. S3b, c). We further performed a regional COREv2 simulation from 1960 to 2007 (Fig. S3d) and show that the shelf-wide pattern (Fig. 2d) is primarily attributed to the high melt signals in the grounding line region (~11 m/y) and margin regions (~5 m/y; Fig. S3e) where ice velocity and ice-flux can be altered. The melt rate in the shelf front region is low at ~2 m/y. Given the same bathymetry of the continental shelf, and the entrances to and troughs underneath the ice shelf cavity, we suggest that during our study period warm mCDW water may have crossed the continental shelf and intruded into the cavity, resulting in increased basal melting and changing in driving stress at the grounding line or resistive stresses along shear margins47,48,49,50,51. Therefore, the mCDW-induced basal melting process is likely responsible for the acceleration during our study period from 1963 to 1989 and further for the long-term acceleration trend of nearly 6 decades in the grounding line region of the TG (Fig. 2c).
Our ice shelf modeling results (âMethodâ section and Fig. S6) reveal that the shelf front retreat during 1973â1985, which caused loss of ice shelf contact with a long near-front section of western margin, induced significant speed up (>300 m/y) of shelf flow in proximity to the western margin (Fig. S6F). This instantaneous response explains the observed acceleration detected in Box 1 using the velocity maps (Fig. 2a). After that, margin contact is regained; ice speed declines and then remains around the long-term average. Furthermore, the modified boundary conditions due to the shelf front retreat impacted the entire ice shelf and immediate grounded areas of grounding line with an average velocity increase of ~84 m/y (Fig. S6F). This explains the observed shelf-wide acceleration from 1963â1973 to 1973â1989 (Fig. 1b). The modeled velocity response in the grounding line region to the calving front retreat is <10m/y, lower than the observed increase in Box 2 (Fig. 2c).
On the decadal scale, the mCDW-induced basal melting indirectly drove the persistent ice discharge increase from 1963 to 2018 (Fig. S4a). The TG basin discharged ice mass under the long-term average SMB (reference SMB) until reaching an equilibrium point in 1989. After this transition, the glacier flowed at speeds above that required to maintain a state of mass balance. TG discharged a total of 3830â±â31 Gt over the entire period from 1963 to 2018. Using the transition point of 1989 as the starting time and reference SMB (67.7â±â4.1 Gt/y) as the basis for forward and backward cumulative mass change computation, we demonstrate that TG discharged 66â±â1 Gt/y on average during the first 27 years, compensating for the high SMB by an additional amount of 54â±â12 Gt (2â±â1 Gt/y) and reaching equilibrium in 1989 (Fig. S4b). Although we cannot show short timespan (e.g., annual) variations, the overall discharge increase trend during the period is clearly demonstrated. Another 92â±â28 Gt was discharged excessively in the last 29 years (3â±â1 Gt/y) to accelerate mass loss. The overall ice discharge acceleration during the 56 years is 0.17â±â0.02 Gt/y2. The cumulative mass balance of the TG basin was dominated by ice discharge, but it was modulated by SMB ranging from -26 Gt to 27 Gt and two significant decreases during 1980â1989 and 2005â2018. Over the entire period of more than 5 decades, the cumulative mass balance in the TG basin changed from a net mass gain of 130â±â27 Gt during 1963â1989 to a net mass loss of 136â±â36 Gt during 1990â2018.
TG is a marine-based glacier16 and has the eastern lobe of its grounding line along the main trunk, sitting on a retrograded bed, away from the ocean (Fig. S3d). Thus, it is subject to marine ice sheet instability (MISI) and has been retreating faster than the western lobe, located on a forward slope. Given the current grounding line position, the retreat rate of 0.15 km/y (refs. 16, 52) and the bed topography, we suggest that the eastern lobe was at the beginning of the current retro-slope during our study period (1963â1989). The mCDW-induced grounding line retreat16 may have begun as early as the 1960s. We further suggest that the TG is currently experiencing a process that occurred earlier at Pine Island Glacier, another marine-based and rapidly changing glacier on the AIS53,54,55 where basal melting induced grounding line retreat and speedup occurred for several decades prior to sustained large scale shelf front retreat was observed since 201756. Analogously, since 1990, the observed TIS calving front positions have been changing modestly, mostly between the shelf fronts of 1973 and 1989 (Fig. S5), and likely with influence from sea ice dynamics14. Meanwhile the mCDW-induced basal melting caused grounding line retreat and acceleration in the TIS for over three decades. Therefore, more rapid calving activities may be expected as basal melting persistently weakens the stability of the ice shelf, resulting in an imbalanced shelf front retreat.
Using historical images of the first-generation satellites of ARGON and Landsat-1 and -4, we reconstructed velocity fields in the TG region, EA from 1963 to 1989. The developed velocity mapping techniques effectively handled the lower quality satellite images and corrected the velocity overestimates caused by the combined effects of ice flow acceleration and long timespan between image pairs, which would have otherwise introduced underestimated long-term velocity changes and mass balance. We found an ice shelf-wide velocity increase trend during 1963â1989. The acceleration in ice shelf front for 1973â1989 is attributed to the large calving front retreat, while that in the grounding line region is mainly caused by the mCDW-induced basal melting, indirectly leading to an increase in ice discharge. The reconciled satellite observations and modeling results reveal that the continued basal melting in TIS drove the long-term acceleration near grounding line and associated ice discharge over the period from 1963 to 2018. With the results of the accelerated ice discharge and decreased SMB, the mass balance state of the TG changed from a mass gain of 5â±â1 Gt/y during 1963â1989 to a mass loss of â5â±â2 Gt/y during 1989â2018. We suggest that recently reported ice flow acceleration and the mass loss trend in the TG basin and the Wilkes Land sector of East Antarctica since the 1980s may have started in the 1960s. As this trend continues, intensified monitoring of ice-air-water interactions in the TIS region is recommended in the next decades.
Methods
Velocity map uncertainty. We estimate velocity map uncertainty \({{{{{{\rm{\sigma }}}}}}}_{{vel}}\) from the orthorectification error \({\sigma }_{{ortho}}\), feature identification error \({\sigma }_{{ident}}\), matching error \({{{{{{\rm{\sigma }}}}}}}_{{match}}\), and timespan \(\Delta t:\,{\sigma }_{{vel}}=\frac{1}{\Delta t}\sqrt{{{\sigma }_{{ortho}}}^{2}+{{\sigma }_{{ident}}}^{2}+{{\sigma }_{{match}}}^{2}}\). \({\sigma }_{{ortho}}\) is calculated using the displacements of check points between the known positions and those measured from the ARGON and Landsat orthoimages (Table S1). \({\sigma }_{{ident}}\) is the error for identifying a feature on the reference (first) image and is set as 0.5 pixels. \({{{{{{\rm{\sigma }}}}}}}_{{match}}\) is also set as 0.5 pixels. We used 2 pixels for ARGON \({\sigma }_{{ident}}\) and 1 pixel for MSS \({{{{{{\rm{\sigma }}}}}}}_{{match}}\) in the 1963â1973 maps because of their lower image quality. The computed uncertainties of the three velocity maps are listed in Table S1. Details of the uncertainty estimation method are presented in ref. 23.
Mass balance. Since RACMO2.3 p2 data are available after 1979, we used the average SMB of 1979â2016 (67.7â±â4.1 Gt/y) as the long-term reference SMB. We used the reference SMB to fill the data gap of 1963â1978 (Fig. 3a). SMB uncertainties are estimated according to the elevation SMB bias scheme7,57,58. To compute the ice flux, we used the ice thickness dataset BedMachine Antarctica59. The flux gate (thick black dashed line in Fig. 1a) is set based on FG2 in ref. 5, but improved by moving toward locations with ice thickness errors close to 30 m59. The average ice thickness uncertainty in this region is reduced from 61 m along FG2 to 32 m along the improved flux gate. Furthermore, an adjustment computed as the SMB between the flux gate and grounding line using RACMO2.3 p2 is applied to obtain the flux (discharge) across the grounding line. The MB of the TG basin is finally estimated for each velocity map period by subtracting the discharge from the SMB (Table 1)5,23. The accumulative discharge, SMB, and MB (Fig. 3b) are calculated, with regards to the reference SMB, as an integration of each variable from 1963 to 1989.
Basal melting 1960â2007. The TG is simulated with the Regional Ocean Modelling System60 framework, which has modifications to include ice-ocean interaction (following refs. 61, 62) using the three-equation parameterization63. The model domain was from 104.5°Eâ130°E with 1/15° resolution and 60°Sâ68°S with 1/30° resolution. Model bathymetry is based on RTopo64 but with the addition of the Totten ice shelf cavity including the western trough cavity entrance8. Lateral boundary forcing and wind forcing are from the COREv2 1949â2007 reanalysis dataset39, while surface heat and salt fluxes are composed of a combination of SSM/I sea ice formation observations65 with COREv2 evaporation minus precipitation fields. The combination allows for long-term estimates of polynya activity and sea ice formation. Further information is given in ref. 38.
Calving modeling. We used Ice Sheet System Model (ISSM)66, a finite element model, to simulate the velocity response of the Totten glacier and ice shelf to calving front retreat. This model is created using velocity and geometry prior to 1973 (Fig. S6). The model construction procedure includes inverting for an ice rigidity parameter, B, within the shelf and a friction coefficient, C, for grounded ice through a minimization of the difference between observed ice flow velocity of 1963â1973 and modeled velocity (Fig. S6B, C). The fidelity of the inferred basal friction structure, common to both domains, is not of interest so long as the ice flow is well described at larger (10s of ice-thicknesses). Consequently, the friction coefficient is obtained during different minimization iterations, without the use of a regularization term. Regularization is included and calibrated with an L-curve analysis67,68, for the inversion of ice rigidity within the shelf where fine scale variations may impact the analysis outcomes (Fig. S6A). Finally, the model was used to simulate the velocity of 1985 with the 1985 shelf front (Fig. S6E) to test the model response to the calving front retreat during 1973â1985.
Data availability
All of the data used in this study are available online: ARGON, Landsat and LIMA images from the United States Geological Survey (USGS) (https://earthexplorer.usgs.gov/); RACMO 2.3 p2 from Institute for Marine and Atmospheric Research, Utrecht University (https://www.projects.science.uu.nl/iceclimate/publications/data/2018/); MEaSUREs BedMachine Antarctica, Version 2 from NASA National Snow and Ice Data Center (NSIDC) (https://nsidc.org/data/NSIDC-0756/versions/2). Velocity maps generated in this study have been deposited in the Dryad database (https://datadryad.org/stash/share/qJNkqLJ5acd1pKeka5u6jFV5cxg5g4XL60W8GfREN6E). Source data are provided in the Source Data file of this paper. Source data are provided with this paper.
Code availability
Code for velocity overestimation calculation is available from corresponding author upon request.
References
Rignot, E. et al. Four decades of Antarctic Ice Sheet mass balance from 1979â2017. Proc. Natl Acad. Sci. USA 116, 1095â1103 (2019).
Smith, B. et al. Pervasive ice sheet mass loss reflects competing ocean and atmosphere processes. Science 368, 1239â1242 (2020).
Velicogna, I. et al. Continuity of ice sheet mass loss in Greenland and Antarctica from the GRACE and GRACE follow-on missions. Geophys. Res. Lett. 47, e2020GL087291 (2020).
The IMBIE team. Mass balance of the Antarctic Ice Sheet from 1992 to 2017. Nature 558, 219â222 (2018).
Gardner, A. S. et al. Increased West Antarctic and unchanged East Antarctic ice discharge over the last 7 years. Cryosphere 12, 521â547 (2018).
Shen, Q. et al. Recent high-resolution Antarctic ice velocity maps reveal increased mass loss in Wilkes Land, East Antarctica. Sci. Rep. 8, 4477 (2018).
Li, X., Rignot, E., Mouginot, J. & Scheuchl, B. Ice flow dynamics and mass loss of Totten Glacier, East Antarctica, from 1989 to 2015. Geophys. Res. Lett. 43, 6366â6373 (2016).
Greenbaum, J. S. et al. Ocean access to a cavity beneath Totten Glacier in East Antarctica. Nat. Geosci. 8, 294â298 (2015).
Li, R. et al. Overestimation and adjustment of Antarctic ice flow velocity fields reconstructed from historical satellite imagery. Cryosphere 16, 737â760 (2022).
Young, D. A. et al. A dynamic early East Antarctic Ice Sheet suggested by ice-covered fjord landscapes. Nature 474, 72â75 (2011).
Weertman, J. Stability of the junction of an ice sheet and an ice shelf. J. Glaciol. 13, 3â11 (1974).
Roberts, J. L. et al. Refined broad-scale sub-glacial morphology of Aurora Subglacial Basin, East Antarctica derived by an ice-dynamics-based interpolation scheme. Cryosphere 5, 551â560 (2011).
Gwyther, D. E., Galton-Fenzi, B. K., Hunter, J. R. & Roberts, J. L. Simulated melt rates for the Totten and Dalton ice shelves. Ocean Sci. 10, 267â279 (2014).
Greene, C. A., Blankenship, D. D., Gwyther, D. E., Silvano, A. & van Wijk, E. Wind causes Totten Ice Shelf melt and acceleration. Sci. Adv. 3, e1701681 (2017).
Bell, R. E. & Seroussi, H. History, mass loss, structure, and dynamic behavior of the Antarctic Ice Sheet. Science 367, 1321â1325 (2020).
Li, X., Rignot, E., Morlighem, M., Mouginot, J. & Scheuchl, B. Grounding line retreat of Totten Glacier, East Antarctica, 1996 to 2013. Geophys. Res. Lett. 42, 8049â8056 (2015).
Liu, Y. et al. Ocean-driven thinning enhances iceberg calving and retreat of Antarctic ice shelves. Proc. Natl Acad. Sci. USA 112, 3263â3268 (2015).
Roberts, J. et al. Ocean forced variability of Totten Glacier mass loss. Geol. Soc. Lond. Spec. Publ. 461, 175â186 (2018).
Bindschadler, R. & Seider, W. Declassified Intelligence Satellite Photography (DISP) Coverage of Antarctica (NASA Goddard Space Flight Center, 1998).
Wang, S., Liu, H., Yu, B., Zhou, G. & Cheng, X. Revealing the early ice flow patterns with historical Declassified Intelligence Satellite Photographs back to 1960s. Geophys. Res. Lett. 43, 5758â5767 (2016).
Li, R. et al. A new analytical method for estimating Antarctic ice flow in the 1960s from historical optical satellite imagery. IEEE Trans. Geosci. Remote Sens. 55, 2771â2785 (2017).
Frezzotti, M., Tabacco, I. E. & Zirizzotti, A. Ice discharge of eastern Dome C drainage area, Antarctica, determined from airborne radar survey and satellite image analysis. J. Glaciol. 46, 253â264 (2000).
Feng, T. et al. A hierarchical network densification approach for reconstruction of historical ice velocity fields in East Antarctica. J. Glaciol. https://doi.org/10.1017/jog.2022.58 (2022).
Zanter, K. Landsat Collection 1 Level 1 Product Definition (United States Geological Survey, 2017).
Li, R. Potential of high-resolution satellite imagery for national mapping products. Photogramm. Eng. Remote Sens. 64, 1165â1169 (1998).
Ye, W. et al. Improved geometric modeling of 1960s KH-5 ARGON satellite images for regional Antarctica applications. Photogramm. Eng. Remote Sens. 83, 477â491 (2017).
Childs, C. Interpolating Surfaces in ArcGIS Spatial (Environmental Systems Research Institute, Inc. (ESRI), Redlands, California, 2004).
Gold, C. M., Remmele, P. R. & Roos, T. Algorithmic Foundations of Geographic Information Systems (eds. van Kreveld, M., Nievergelt, J., Roos, T. & Widmayer, P.) vol. 1340. p. 21â35 (Springer Berlin Heidelberg, 2006).
Cheng, Y. et al. Ice flow velocity mapping of East Antarctica from 1963 to 1989. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. XLII-2/W13, 1735â1739 (2019).
Luo, S. et al. Ice flow velocity mapping in East Antarctica using historical images from 1960s to 1980s: recent progress. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. XLIII-B3-2021, 491â496 (2021).
Berthier, E., Raup, B. & Scambos, T. New velocity map and mass-balance estimate of Mertz Glacier, East Antarctica, derived from Landsat sequential imagery. J. Glaciol. 49, 503â511 (2003).
Cook, S. et al. Modelled fracture and calving on the Totten Ice Shelf. Cryosphere 12, 2401â2411 (2018).
Fürst, J. J. et al. The safety band of Antarctic ice shelves. Nat. Clim. Change 6, 479â482 (2016).
Hogg, A. E. & Gudmundsson, G. H. Impacts of the Larsen-C Ice Shelf calving event. Nat. Clim. Change 7, 540â542 (2017).
Mouginot, J., Scheuchl, B. & Rignot, E. MEaSUREs Annual Antarctic Ice Velocity Maps, Version 1 https://doi.org/10.5067/9T4EPQXTJYW9 (2017).
Gardner, A. S., Fahnestock, M. A. & Scambos, T. A. MEaSUREs ITS_LIVE Landsat Image-Pair Glacier and Ice Sheet Surface Velocities: Version 1. https://doi.org/10.5067/IMR9D3PEI28U (2019).
Rignot, E., Scheuchl, B. & Mouginot, J. MEaSUREs Multi-year Reference Velocity Maps of the Antarctic Ice Sheet, Version 1. https://doi.org/10.5067/FB851ZIZYX5O (2022).
Gwyther, D. E., OâKane, T. J., Galton-Fenzi, B. K., Monselesan, D. P. & Greenbaum, J. S. Intrinsic processes drive variability in basal melting of the Totten Glacier Ice Shelf. Nat. Commun. 9, 3141 (2018).
Large, W. G. & Yeager, S. G. The global climatology of an interannually varying airâsea flux data set. Clim. Dyn. 33, 341â364 (2009).
Adusumilli, S., Fricker, H. A., Medley, B., Padman, L. & Siegfried, M. R. Interannual variations in meltwater input to the Southern Ocean from Antarctic ice shelves. Nat. Geosci. 13, 616â620 (2020).
Rintoul, S. R. et al. Ocean heat drives rapid basal melt of the Totten Ice Shelf. Sci. Adv. 2, e1601610 (2016).
Bindoff, N. L., Rosenberg, M. A. & Warner, M. J. On the circulation and water masses over the Antarctic continental slope and rise between 80 and 150°E. Deep Sea Res. Part II Top. Stud. Oceanogr. 47, 2299â2326 (2000).
Williams, G. D. et al. Late winter oceanography off the Sabrina and BANZARE coast (117â128°E), East Antarctica. Deep Sea Res. Part II Top. Stud. Oceanogr. 58, 1194â1210 (2011).
Nakayama, Y. et al. Antarctic slope current modulates ocean heat intrusions towards Totten Glacier. Geophys. Res. Lett. 48, e2021GL09414 (2021).
Silvano, A. et al. Freshening by glacial meltwater enhances melting of ice shelves and reduces formation of Antarctic Bottom Water. Sci. Adv. 4, eaap9467 (2018).
Khazendar, A. et al. Observed thinning of Totten Glacier is linked to coastal polynya variability. Nat. Commun. 4, 2857 (2013).
Pritchard, H. D. et al. Antarctic ice-sheet loss driven by basal melting of ice shelves. Nature 484, 502â505 (2012).
Meyer, C. R. & Minchew, B. M. Temperate ice in the shear margins of the Antarctic Ice Sheet: Controlling processes and preliminary locations. Earth Planet. Sci. Lett. 498, 17â26 (2018).
De Rydt, J., Reese, R., Paolo, F. S. & Gudmundsson, G. H. Drivers of Pine Island Glacier speed-up between 1996 and 2016. Cryosphere 15, 113â132 (2021).
Feldmann, J., Reese, R., Winkelmann, R. & Levermann, A. Shear-margin melting causes stronger transient ice discharge than ice-stream melting in idealized simulations. Cryosphere 16, 1927â1940 (2022).
Summers, P. T., Elsworth, C. W., Dow, C. F. & Suckale, J. Migration of the shear margins at Thwaites Glacier: dependence on basal conditions and testability against field data. J. Geophys. Res. Earth Surf. 128, e2022JF006958 (2023).
Li, T., Dawson, G. J., Chuter, S. J. & Bamber, J. L. Grounding line retreat and tide-modulated ocean channels at Moscow University and Totten Glacier ice shelves, East Antarctica. Cryosphere Discuss https://doi.org/10.5194/tc-2022-129 (2022).
Rignot, E., Vaughan, D. G., Schmeltz, M., Dupont, T. & Macayeal, D. Acceleration of Pine Island and Thwaites Glaciers, West Antarctica. Ann. Glaciol. 34, 189â194 (2002).
Payne, A. J., Vieli, A., Shepherd, A. P., Wingham, D. J. & Rignot, E. Recent dramatic thinning of largest West Antarctic ice stream triggered by oceans. Geophys. Res. Lett. 31, L23401 (2004).
Christianson, K. et al. Sensitivity of Pine Island Glacier to observed ocean forcing: PIG response to ocean forcing. Geophys. Res. Lett. 43, 10,817â10,825 (2016).
Joughin, I., Shapero, D., Smith, B., Dutrieux, P. & Barham, M. Ice-shelf retreat drives recent Pine Island Glacier speedup. Sci. Adv. 7, eabg3080 (2021).
Van Wessem, J. M. et al. Improved representation of East Antarctic surface mass balance in a regional atmospheric climate model. J. Glaciol. 60, 761â770 (2014).
van Wessem, J. M. et al. Modelling the climate and surface mass balance of polar ice sheets using RACMO2 â Part 2: Antarctica (1979â2016). Cryosphere 12, 1479â1498 (2018).
Morlighem, M. et al. Deep glacial troughs and stabilizing ridges unveiled beneath the margins of the Antarctic ice sheet. Nat. Geosci. 13, 132â137 (2020).
Shchepetkin, A. F. & McWilliams, J. C. The regional oceanic modeling system (ROMS): a split-explicit, free-surface, topography-following-coordinate oceanic model. Ocean Model. 9, 347â404 (2005).
Dinniman, M. S., Klinck, J. M. & Smith, W. O. Influence of sea ice cover and icebergs on circulation and water mass formation in a numerical circulation model of the Ross Sea, Antarctica. J. Geophys. Res. 112, C11013 (2007).
Galton-Fenzi, B. K., Hunter, J. R., Coleman, R., Marsland, S. J. & Warner, R. C. Modeling the basal melting and marine ice accretion of the Amery Ice Shelf. J. Geophys. Res. Oceans 117, C09031 (2012).
Holland, D. M. & Jenkins, A. Modeling thermodynamic iceâocean interactions at the base of an ice shelf. J. Phys. Oceanogr. 29, 1787â1800 (1999).
Timmermann, R. et al. A consistent data set of Antarctic ice sheet topography, cavity geometry, and global bathymetry. Earth Syst. Sci. Data 2, 261â273 (2010).
Tamura, T., Ohshima, K. I. & Nihashi, S. Mapping of sea ice production for Antarctic coastal polynyas. Geophys. Res. Lett. 35, L07606 (2008).
Larour, E., Seroussi, H., Morlighem, M. & Rignot, E. Continental scale, high order, high spatial resolution, ice sheet modeling using the Ice Sheet System Model (ISSM). J. Geophys. Res. Earth Surf. 117, F01022 (2012).
Morlighem, M., Seroussi, H., Larour, E. & Rignot, E. Inversion of basal friction in Antarctica using exact and incomplete adjoints of a higher-order model. J. Geophys. Res. Earth Surf. 118, 1746â1753 (2013).
Wolovick, M., Humbert, A., Kleiner, T. & Rückamp, M. Regularization and L-curves in ice sheet inverse models, a case study in the Filchner-Ronne catchment. EGUsphere [Preprint]. https://doi.org/10.5194/egusphere-2023-741 (2023).
Bindschadler, R. et al. The Landsat image mosaic of Antarctica. Remote Sens. Environ. 112, 4214â4226 (2008).
Rignot, E., Mouginot, J. & Scheuchl, B. Antarctic grounding line mapping from differential satellite radar interferometry. Geophys. Res. Lett. 38, L10504 (2011).
Acknowledgements
We appreciate constructive comments and suggestions from referees and editors. This research has been supported by the National Key Research and Development Program of China (No. 2017YFA0603100, R.L., G.Q., and X.T.; 2021YFB3900105, L.A. and G.Q.), the National Science Foundation of China (41730102, R.L.), and the Fundamental Research Funds for the Central Universities (R.L., X.T., L.A., and G.Q.).
Author information
Authors and Affiliations
Contributions
R.L. designed and lead research; Y.C., T.C., D.E.G., M.F., L.A., M.X., Y.H., and W.Y. performed research; R.L., Y.C., D.E.G., M.F., L.A., G.Q., and X.T. analyzed data; and R.L. wrote the paper.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Li, R., Cheng, Y., Chang, T. et al. Satellite record reveals 1960s acceleration of Totten Ice Shelf in East Antarctica. Nat Commun 14, 4061 (2023). https://doi.org/10.1038/s41467-023-39588-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-023-39588-x
This article is cited by
-
Status and trends in the stability of the three largest ice shelves in Antarctica
Science China Earth Sciences (2024)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.