Next Article in Journal
Comparing Water Indices for Landsat Data for Automated Surface Water Body Extraction under Complex Ground Background: A Case Study in Jilin Province
Previous Article in Journal
A Machine Learning Approach to Derive Aerosol Properties from All-Sky Camera Imagery
Previous Article in Special Issue
Multitemporal Glacier Mass Balance and Area Changes in the Puruogangri Ice Field during 1975–2021 Based on Multisource Satellite Observations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mass Balance of the Antarctic Ice Sheet in the Early 21st Century

School of Geospatial Engineering and Science, Sun Yat-sen University, Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519082, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(6), 1677; https://doi.org/10.3390/rs15061677
Submission received: 9 February 2023 / Revised: 15 March 2023 / Accepted: 18 March 2023 / Published: 20 March 2023
(This article belongs to the Special Issue Remote Sensing of Ice Loss Tracking at the Poles)

Abstract

:
Mass loss from the Antarctic Ice Sheet (AIS) is an important contributor to global sea level rise. To examine the recent ice loss, we estimated the mass budget of the AIS from 2000 to 2020 using multiple ice velocity datasets, state-of-the-art ice thickness datasets, and extended surface mass balance (SMB) records. The AIS lost mass at an average rate of −89 ± 99 Gt/yr over the study period. The East Antarctic Ice Sheet (EAIS) showed a slightly positive mass balance, while the West Antarctic Ice Sheet (WAIS) experienced a significant acceleration in mass loss. The ice discharge from the AIS increased from 1792 ± 47 Gt/yr in 2000 to 1940 ± 37 Gt/yr in 2017–2020, with the increase in the discharge from the WAIS being three to four times higher than that from the EAIS. Moreover, the average mass balance for 2017–2020 was −99 ± 93 Gt/yr, slightly more negative than the average for the early 21st Century. During this recent period, the ice discharge decreased in the East Indian Ocean sector, in contrast to its rapid increase from 2000 to 2013–2017. However, the discharge in the Amundsen Sea sector still greatly increased from 2013–2017 to 2017–2020. Overall, our results are in agreement with recent mass balance estimates for the AIS based on gravimetry and altimetry. Our assessments of the recent AIS mass balance with the mass budget method (input-output method) will contribute to the understanding of ice dynamic processes and provide insights into the stability of the AIS.

Graphical Abstract

1. Introduction

Approximately 70% of the world’s freshwater is stored as ice in the Antarctic Ice Sheet (AIS), and it holds enormous potential to raise global sea level by 58 m [1,2,3]. The AIS contributed 7.6 ± 3.9 mm to global mean sea level rise from 1992 to 2017 [4]. These freshwater ice changes can not only lead to sea level rise, but can also dramatically impact the ocean surface temperature, salinity, thermohaline circulation, and global carbon cycle [5,6,7]. Continuous monitoring of the mass balance of the AIS is essential for predicting sea level change. A more accurate understanding of the spatial and temporal distribution of the mass balance of the AIS, as well as how it responds to the patterns of climate change, is also essential for understanding the stability of the ice sheet.
Recent studies have shown that the AIS is losing mass mainly due to changes in ice dynamics caused by the intrusion of warmer, saltier sea water around Antarctica [8,9,10]. At present, various remote sensing-based methods have been used to estimate the mass balance of the AIS, including: (1) the mass budget method (or input-output method, IOM), which estimates the mass change by comparing the surface mass balance (input) of the AIS with the ice discharge (output) at the grounding line [11]; (2) the altimetry method, where the change in surface elevation of the entire AIS is measured by satellite altimeters, and the change in ice volume is translated into the mass change of the ice sheet [12]; and (3) the gravity-based method, which can obtain the time-varying signal of the Earth’s gravity field directly from gravity satellites, and then invert it into the relative change in ice sheet mass [13]. Each of these methods utilizes different data sources and inherent uncertainties, resulting in large differences in the results [4]. The IOM is more a parameterized and computationally complex method compared to the altimetry and gravity methods. However, it tends to utilize higher-resolution data (100 m to 1 km) [8]. Furthermore, the IOM is better suited for revealing the physical mass loss process, which is crucial for understanding the variations in ice mass changes across different basins of the AIS.
The IOM has long been used to assess the mass balance of the AIS [11,14,15]. Recently, Rignot et al. [8] estimated the mass balance for the entirety of Antarctica over nearly four decades from 1979–2017. An improvement to reduce the uncertainty in the discharge estimation was proposed by Gardner et al. [16]. They calculated the mass balance of the AIS for 2008–2015 using Antarctic-wide ice velocity products derived from Landsat images and analyzed the trends in velocity and discharge. Ice velocity products with a higher spatial resolution were used to estimate the mass balance of the AIS in 2008, 2014, and 2015 [17]. Other studies have used the IOM to estimate regional mass balance values at the drainage basin scale. Estimates of the mass balance of the Totten Glacier during 1989–2015 revealed a link between ice dynamic changes and ocean temperature, with a significant increase in ice velocity occurring in response to increasing ocean temperatures [18]. Assessments of the Amery Ice Shelf region show a long-term slight positive mass balance for this ice shelf basin [19,20,21]. A study of the Getz region in West Antarctica revealed that, from 1994 to 2008, two-thirds of the mass loss in this region was due to dynamic imbalance, which was likely related to ocean forcing [22]. From the results of the different ways to assess the mass balance published by the Ice sheet Mass Balance Inter-Comparison Exercise (IMBIE) in 2018, there are still relatively large differences among the three estimation methods. For example, for the mass balance results of the AIS from 2003–2010, the altimetry method yields −43 ± 21 Gt/yr, the gravimetry method yields −76 ± 20 Gt/yr, and the IOM yields −201 ± 82 Gt/yr [4]. Combining the results of the three methods while exploring the intrinsic relationships between them can help us better understand the causes of fluctuations in the mass balance of the AIS. The mass balance of the AIS has been estimated up to 2019 using the altimetry method [7] and up to 2020 using the gravimetry method [23]. However, the comprehensive evaluation of recent Antarctic-wide mass balance is currently lacking due to limitations in data and methodology. Estimates of the mass balance of the AIS using the IOM stopped in 2017 [8], preventing us from comparing the differences in the assessments of recent mass balance with the three methods.
In this study, we use the IOM to estimate the mass balance of the AIS up to the year 2020. We assess the mass balance variability of the AIS for individual basins from 2000 to 2020 using multiple surface velocity data products combined with surface mass balance (SMB) estimates. Additionally, in order to provide a comprehensive understanding of the mass balance of the AIS, we compare our results obtained via the IOM with the Gravity Recovery and Climate Experiment (GRACE) mass change results and altimetry-based estimates. We also investigate potential explanations for any discrepancies observed between the results obtained through different data and methods.

2. Data and Methods

2.1. Ice Velocity

The MEaSUREs Annual Antarctic Ice Velocity Maps [24] provide annual maps of Antarctic ice surface velocity, assembled using optical and synthetic aperture radar (SAR) satellite data. Our study used these annual ice surface velocity data from 2000–2020. The spatial resolution of these velocity data is 1 km, and the uncertainty is ~1–17 m/yr [24]. The Inter-Mission Time Series of Land Ice Velocity and Elevation (ITS_LIVE) ice velocity data derived from Landsat 4, 5, 7, and 8 imagery using the autonomous Repeat Image Feature Tracking (auto-RIFT) method were also used [16,25]. The spatial resolution of this dataset is 240 m, and the uncertainty is approximately 10 m/yr [16]. Other gaps in the coverage of the annual velocity product, especially for the regions south of 82.7°S, were filled using the MEaSUREs InSAR-Based Ice Velocity Maps of Central Antarctica [26] and MEaSUREs Phase-Based Antarctica Ice Velocity Map [27].
The MEaSUREs data products have a broad temporal coverage and use multiple satellite platforms, and these products form the majority of the ice velocity data used in this study. Although the MEaSUREs products have been able to represent the regional variation of the ice velocity well [28], there are still some missing values near the flux gate that need to be filled for the ice discharge calculation.
The Central Antarctica product was used to fill in missing velocity values located within the gap in the coverages of the satellite platforms at the pole. We assumed that the trend in velocity in these regions has been stable for a long time. The long-term velocity change rate was estimated from calculating the difference in Central Antarctica velocity between 2009 and 1997. Then, based on Central Antarctica velocity in 2009 and the long-term change rate, the missing values in each annual velocity map were calculated by a linear scaling method.
The ITS_LIVE annual products were further used to fill in missing values near the flux gate in the MEaSUREs annual velocity products with the nearest neighbor method prior to 2013. Note that after this step the annual ice velocity was already relatively complete. In addition, we used MEaSUREs Phase-Based Antarctica Ice Velocity Map to perform the final filling step to obtain a complete annual ice velocity mosaic. We divided the whole period into four epochs and listed ice velocity products used to calculate discharge in each of the epochs in Table 1. The ice velocity of each epoch was then estimated by averaging the annual ice velocity mosaics.

2.2. Ice Thickness, Flux Gate, Grounding Line, and Basin Boundary

Ice thickness has been considered the most significant source of error in estimating mass balance by the IOM [3]. The most accurate ice thickness data for the AIS come from airborne radio echo sounding (RES) [1]. However, RES data are very sparse over some of the AIS, and many regions are still not well-covered; in particular, more than 80% of the grounding line (GL) zones have not been sounded [16]. To improve the accuracy of the discharge estimation, we adapted the flux gate (FG) method proposed by Gardner et al. [16] (Figure 1). The FG was positioned close to the RES flight lines to include more RES data by sacrificing proximity to the GL. We derived the ice thickness along the FG from the BedMachine v2 dataset [3] (Figure 2a). This dataset includes most airborne RES data and introduces flow velocity and SMB data to refine the thickness estimates by correcting temporal changes in surface mass balance and ice thickness using the mass conservation method. Therefore, we expect it to have a low uncertainty in the fast-flowing regions and provide the most reliable estimates for our calculation.
We defined the area of each basin by combining an updated version of the GL data [16,30] with the Antarctic basin boundary from Zwally et al. [29] (Figure 1). The continent was divided into twenty-seven basins (B1–B27), including sixteen basins in the East Antarctica Ice Sheet (EAIS: B2–B17), seven basins in the West Antarctica Ice Sheet (WAIS: B1, B18–B23), and four basins on the Antarctic Peninsula (AP: B24–B27).

2.3. Ice Discharge

For the estimation of ice discharge, the conventional approach uses the position of GL as the FG, which means the ice flux (IF) is equal to the discharge, to directly calculate the discharge [2,8]. Because of the choice of FG, the IF and discharge are different in this study. We first discretize the FG line and then sort the FG nodes in clockwise ascending order to help determine the inflow or outflow of each FG unit (two consecutive FG nodes form an FG unit). Then, we use the ice velocity and thickness data at the location of the FG unit to estimate the IF. The calculation of the IF is expressed as follows:
I F = i = 1 n H i V x i W x i + V y i W y i ρ i c e
where n is the total number of FG units and Hi is the ice equivalent thickness. Vxi and Vyi denote the velocity components of the FG units in the x and y directions, respectively. Wxi is the FG projected in y coordinates and is perpendicular to Vxi. Wyi is the FG projected in x coordinates and is perpendicular to Vyi. ρice is the density of ice. Here, we use 917 kg/m3 with an uncertainty of 5 kg/m3 in the estimation [31].
We used the following equation based on the principle of mass conservation [32] to further estimate the discharge (D) from the IF:
D = I F + S M B G L F G d V d y n / d t
where SMBGL–FG is the unmeasured flux due to the (typically) positive surface mass balance in the region between the FG and the GL. dVdyn/dt is the dynamic volume change, which is the unmeasured flux due to ice dynamics between the FG and the GL. We assume that the ice elevation only changes due to dynamic changes, and dVdyn/dt can be regarded as a constant throughout the study period. We used the surface elevation change data provided by the CPOM (Centre for Polar Observation and Modelling) data portal, measured with ENVISAT and CryoSat-2 data, to calculate the dynamic volume change in each basin from 2002–2012 and 2013–2020 (Figure 2b). The CPOM elevation change data with 5 km spatial resolution are available at http://www.cpom.ucl.ac.uk/csopr/icesheets2/?icesheet_zone=antarctica accessed on 26 May 2021). The ice discharge of each epoch was estimated by averaging annual ice discharge results.

2.4. Surface Mass Balance

There are currently two conventional ways to assess the SMB of the AIS at the continental scale: either using a regional climate model (RCM) or using estimates from an atmospheric reanalysis product [33,34]. The Regional Atmospheric Climate Model (RACMO) is usually used to calculate the SMB in Antarctica [8,16,17]. Unfortunately, the latest RACMO Version 2.3p2 (RACMO2.3p2) is only available until 2017 [35]. Atmospheric reanalysis products usually include only two main SMB components: precipitation and evaporation (sublimation), with no components representing the liquid water balance [36]. Therefore, in this study, we use RACMO to assess the SMB of the AIS, and a combination of reanalysis data and RACMO to extend the SMB data to include recent years.
The spatial resolution of the outputs of the regional climate models RACMO2.3p1 and RACMO2.3p2 [35,37] is 27 km. RACMO2.3p2 is an updated version of RACMO2.3p1 that features changes to the topography, cloud scheme parameters, and snow properties [35]. We used RACMO2.3p1 to estimate the SMB values in the EAIS, as several studies have shown that RACMO2.3p1 performs better than RACMO2.3p2 in SMB estimates for the EAIS [8,38,39]. For the WAIS and the AP, we calculated the SMB using RACMO2.3p2. The SMB data product was first resampled to a grid size of 1 km. Then, the polygon vectors of individual basins (or the whole ice sheet) were masked to the grid.
The ERA5-Land product is a state-of-the-art reanalysis product, covering the period from 1950 to the present day [40]. The improvements in this recently released product include refined atmospheric physics and a revised snow scheme. We used ERA5-Land reanalysis data from 1979–2021 to extend the RACMO SMB data to 2020. We calculated the SMB from ERA5-Land by subtracting evaporation from precipitation.
For basins in the EAIS (B2–B17) in Figure 1, we used RACMO2.3p1, which only runs to 2015, and ERA5-Land data from 1979–2015 to estimate a scale factor. We calculated the total SMB from 1979–2015 for each individual basin, and the scale factor was estimated as the ratio of the SMB derived from those two datasets. Then, we scaled the SMB calculated from ERA5-Land reanalysis data to obtain the corrected SMB for the period 2015–2020. For the basins in the WAIS (B1, B18–B23) and AP (B24–B27), we performed the same calculation to extend the SMB record, but used the RACMO2.3p2 and ERA5-Land datasets from 1979–2017 to estimate the scale factor for each basin.

2.5. Mass Balance

The equation for calculating mass balance (MB) using the IOM is
M B = S M B D
We use the FG defined in Section 2.3 to reduce the uncertainty in calculating the discharge (D) of the AIS. In addition, we extend the SMB of the AIS with ERA5-Land reanalysis data. The MB of the AIS for 2000–2020 is then calculated.
The discharge and SMB in the northern AP have higher uncertainty than in other regions [8]. The region has a complex basal topography and steep slopes, and the ice velocity products also have poor quality. Due to the extreme surface topography, the annual variability in precipitation and melt is also large, and the measurement of the SMB is especially challenging [16]. Therefore, we use the net mass budget results estimated from the altimetry method during the 2003–2019 period [7] to replace our mass balance estimation for B25 and B26.

2.6. Uncertainty Assessment

Uncertainty assessment of the 27 basins of the AIS needs to consider both discharge and SMB components. The different parameters are assumed to be independent in the calculation. The uncertainty in the ice flux for each basin is as follows:
σ I F = σ H 2 + σ V 2 + σ ρ i c e 2
where the σIF is the uncertainty in ice flux. σH, σV, σρice are the uncertainty in ice flux that results from the uncertainty in ice thickness (H), ice velocity (V), and ice density (ρice), respectively.
The uncertainties in ice thickness, ice velocity, and ice density are calculated as follows:
σ H = i = 1 m ( u H i W i V i ρ i c e ) 2
σ V = i = 1 m ( u V i H i W i ρ i c e ) 2
σ ρ i c e = i = 1 m ( u ρ i c e V i H i W i ) 2
where m is the number of flux units of each basin. Hi and uHi represent the ice thickness and uncertainty in the ice thickness for each flux unit, respectively. Vi and uVi represent the ice surface velocity and its uncertainty, respectively. ρice and ice represent the ice density and its uncertainty, respectively. Wi is the length of each flux unit perpendicular to the velocity direction.
The uncertainty in discharge in each basin can be expressed as follows:
σ D = σ I F 2 + σ Δ S M B 2 + σ d V d y n / d t 2
Then, the uncertainty in Antarctic ice discharge is
σ D A I S = i = 1 n σ D i 2
where n is the number of basins. σΔSMB is the uncertainty in the SMB of the region between the GL and the FG, and σdVdyn/dt is the uncertainty in the dynamic volume change. In this study, 20% was used to assess the uncertainty in the SMB in each basin [16]. Uncertainty in the dynamic volume change cannot be rigorously quantified. Therefore, we conservatively assumed an uncertainty of 40% of the estimated dynamic volume change magnitude.
The uncertainty in the resulting AIS mass balance estimate, σMB, is calculated as follows:
σ M B = σ D A I S 2 + i = 1 n σ S M B i 2
where σSMBi is the uncertainty in each basin’s SMB.

3. Results

3.1. Changes in Ice Velocity and Ice Discharge from 2000 to 2020

We obtained ice velocity change (Figure 3) and ice discharge results change (Figure 4) for 2000, 2007–2011, 2013–2017, and 2017–2020. The AIS was further divided into six sea sectors to better demonstrate the spatial and temporal characteristics of ice velocity and ice discharge changes over the past 20 years.
Overall, our results show that the ice discharge has been increasing over time for the whole Antarctic in the four periods, with values of 1792 ± 47 Gt/yr, 1876 ± 51 Gt/yr, 1908 ± 51 Gt/yr, and 1940 ± 37 Gt/yr for the periods 2000, 2007–2011, 2013–2017, and 2017–2020, respectively (Table 2). The Amundsen Sea sector contributes ~70% of the discharge growth between the 2000 and 2017–2020 periods and is the region with the most rapidly increasing ice discharge in Antarctica (Figure 3). In addition, both the WAIS and the EAIS show increasing trends in ice discharge during the study period. From 2000 to 2017–2020, the discharge from the WAIS grew from 652 ± 32 Gt/yr to 753 ± 29 Gt/yr, an increase of up to 15%, and the discharge from the EAIS grew from 968 ± 32 Gt/yr to 1000 ± 22 Gt/yr. The growth rates of ice discharge vary from region to region over time. Below, we divided the AIS into six sectors to demonstrate the changes in detail.

3.1.1. Amundsen Sea Sector (B20–B22)

Our results show that the Amundsen Sea sector is the region with the most dramatic increase in Antarctic ice velocity and ice discharge, as previous studies have revealed [8,17]. The ice discharge in this sector was approximately 114 Gt/yr higher in 2017–2020 than in 2000, which is approximately ten times the change in other sea sectors. The enormous ice velocity change in B20 occurred at the Getz Ice Shelf (Getz IS). The ice velocity at the FG in this region increased up to 300 m/yr between 2000 and 2017–2020. The ice discharge in B20 increased by 22 Gt/yr during this period, with most of the discharge increases occurring between 2007–2011 and 2013–2017 due to significant ice velocity increases (R2 in Figure 4e).
The more representative region of B21 is the Thwaites Glacier, Smith Glacier, and Pope Glacier, which showed an average increase in ice velocity from 2000 to 2017–2020 of 83 m/yr. The ice velocity along the Thwaites Glacier FG increased by an average of 213 m/yr to a maximum of 1137 m/yr. The ice discharge in B21 increased ~20% from 2017–2020 compared with 2000. The ice velocity and discharge in this region showed a consistent increase in each period.
The main outlet glacier in B22 is Pine Island Glacier, one of the fastest-flowing glaciers in Antarctica. Between 2000 and 2017–2020, the average ice velocity increased by 365 m/yr at the FG of this glacier. The ice discharge in B22 increased by 41.1 Gt/yr, i.e., ~40%, the largest increase among all the basins in Antarctica. However, the ice discharge remained stable between 2007–2011 and 2013–2017 (Figure 3). This is similar to the conclusion of Rignot et al. [8], that the mass loss of the Pine Island Glacier stabilized after 2012. However, this glacier’s surface ice velocity increased again between 2013–2017 and 2017–2020 (R6 in Figure 4e). The ice velocity reached a maximum of just over 4 km/year, significantly increasing the mass loss from Pine Island Glacier.

3.1.2. Bellingshausen Sea Sector (B23–B25)

The overall trend of ice discharge growth in the Bellingshausen Sea sector was 13.7 Gt/yr from 2000 to 2017–2020. The largest contribution was 11.1 Gt/yr from B25, with a significant discharge increase found between 2000 and 2007–2011. The ice discharge in B24 maintained a slow growth trend during the study period, while the discharge in B23 exhibited almost no change and was relatively stable. The ice discharge for all the basins in this sector remained stable in the period from 2013–2017 to 2017–2020.

3.1.3. Weddell Sea Sector (B26, B27, and B1–B5)

There was a slight increase in ice discharge in the Weddell Sea sector during the study period, with an increase of approximately 10.6 Gt/yr from 2000 to 2017–2020. For B4 and B5 in this sector, the discharge changes mainly occurred before 2011. It is worth noting the significant increase in ice discharge in B1 between 2013–2017 and 2017–2020. The ice velocity in this basin also increased during this period (R5 in Figure 4e). The maximum velocity increases reached 97 m/yr in some regions. In addition, the Brunt Ice Shelf (Brunt IS) experienced a dramatic change in ice velocity (R4 in Figure 4e), as other studies have also revealed [42].

3.1.4. West Indian Ocean Sector (B6–B11)

During the study period, there was a weak increase in ice discharge in this sector, with a value of only approximately 9.9 Gt/yr. This increase was mainly due to the increase in ice discharge from B7, while the discharge of other basins remained stable during the study period. We also found that the ice velocity increased by a up to ~200 m/yr in some regions of B7.

3.1.5. East Indian Ocean Sector (B12–B14)

The increase in ice discharge in this sector between 2000 and 2017–2020 was approximately 18.9 Gt/yr. This increase mainly occurred between 2000 and 2013–2017, while the total ice discharge stopped increasing from 2013–2017 to 2017–2020.
The ice velocity along the FG in B13 showed an obvious increase between 2000 and 2007–2011, with the most significant increase occurring in the Vincennes Bay region with an average of 55 m/yr. As with the overall discharge changes in this sector, the ice discharge of B13 showed a slight decrease between 2013–2017 and 2017–2020, mainly resulting from the decrease in ice velocity of the two main glaciers in this basin. The ice velocity of Totten Glacier decreased by an average of 11 m/yr during this period, and that of Moscow University Glacier also showed a slight decline (R7 in Figure 4e).
In B12, the glaciers feeding the West Ice Shelf (West IS) and the Shackleton Ice Shelf (Shackleton IS) showed an average increase in ice velocity of approximately 17 m/yr, reaching up to 450 m/yr in some regions from 2000 to 2007–2011. Although glaciers feeding the West IS showed a deceleration from 2007–2011 to 2013–2017, the ice velocity of glaciers feeding the Shackleton IS continued to increase during this period (R1 in Figure 4e). There was a slight decrease in the ice velocity of glaciers feeding the West IS and the Shackleton IS from 2013–2017 to 2017–2020, similar to the trend of B13. The ice velocity of the Ninnis Glacier and the Mertz Glacier in B14 rapidly increased from 2007–2011 to 2013–2017 (R3 in Figure 4e), resulting in an ice discharge increase of ~6 Gt/yr during this period. However, again, the ice velocity remained almost unchanged between 2013–2017 and 2017–2020.

3.1.6. Ross Sea Sector (B15–B19)

The Ross Sea sector is the only region where ice discharge decreased during the study period, with a reduction of 17.0 Gt/yr from 2000 to 2017–2020. This is most likely due to the ice velocity decrease in B17 and B18. For the recent period from 2013–2017 to 2017–2020, the average annual ice discharge slightly increased by 7.0 Gt/yr.

3.2. SMB Calculation and Changes

We first performed a correlation test to demonstrate the similarity of the SMB results estimated from ERA5-Land and RACMO2.3 at the annual scale (Figure 5). In the EAIS, we used the SMB results of RACMO2.3p1 and ERA5-Land from 1979–2015 for the test, and they showed a strong correlation. The R2 is as high as 0.755 (Figure 5a), and the slope of the linear fit trend line is 0.924. In the WAIS and the AP, we used the results of RACMO2.3p2 and ERA5-Land from 1979–2017. The data for the WAIS and the AP also showed strong correlations over these 39 years, with R2 values of 0.877 (Figure 5b) and 0.929 (Figure 5c), respectively. In summary, ERA5-Land reanalysis results strongly correlated with the RACMO2.3 SMB outputs and could be used to extend the RACMO2.3 SMB results.
Using the method described in Section 2.4, we extended the RACMO2.3 SMB records to 2020 (Figure 6). The annual SMB results for all of Antarctica during the study period 2000–2020 range between 1694–2104 Gt, with a maximum interannual difference of up to 410 Gt. The interannual changes in the SMB for the EAIS, WAIS, and AP range from 994–1280 Gt, 486–687 Gt, and 159–292 Gt, respectively. The large interannual fluctuations in the SMB may directly contribute to the interannual change in the mass balance. However, the Antarctic SMB does not show a significant trend during our study period [43].

3.3. Changes in the Mass Balance

The mean value of the AIS mass balance for 2000–2020 was −89 ± 99 Gt/yr (Table 3). During the study period, almost all the years show a state of mass loss. This generally agrees with recent estimates of Antarctic mass change based on laser altimetry of −118 ± 24 Gt/yr for 2003–2019.
The mass balance of the EAIS during the study period was positive, which agrees well with most estimates for the EAIS [4,10,17]. The average mass accumulation of the EAIS from 2000 to 2020 was 53 ± 74 Gt/yr, and the average mass accumulation values in 2000, 2007–2011, 2013–2017, and 2017–2020 were 56 ± 62 Gt/yr, 13 ± 55 Gt/yr, 70 ± 58 Gt/yr, and 70 ± 55 Gt/yr, respectively (Table 3). The anomalously low SMB value in approximately 2007 led to a smaller mass balance in 2007–2011 than in other years. The basins with significant mass growth were B6–B8 (Queen Mary Land) and B17, which contributed 63 ± 34 Gt/yr to the mass accumulation of the EAIS. The mass of the Amery Ice Shelf region (B9-B11) was nearly balanced during the entire study period, which is consistent with most studies of the region within the stated errors (Table 3) [21,39]. In contrast, the East Indian Ocean sector (B12–B14) showed a significant acceleration in mass loss between 2000 and 2017–2020, which deserves more attention. The mass balance in this region changed from a weak mass equilibrium (1 ± 59 Gt/yr) in 2000 to a state of dramatic mass loss (−47 ± 55 Gt/yr) in 2017–2020. Considering that the ice discharge in this sector was almost unchanged between 2013–2017 and 2017–2020, the recent acceleration in mass loss was likely mainly associated with the decline in the SMB in the region in 2017–2020.
The WAIS is the region with the most severe mass loss, with a mean mass balance of −127 ± 56 Gt/yr (Table 3) over the period 2000–2020. In particular, the Amundsen Sea sector (B20–B22), which has undergone multiple ice shelf breakups and continued ice shelf thinning [44,45], had a mean mass loss of −155 ± 24 Gt/yr during the study period. The mass loss showed an acceleration between 2000 and 2013–2017, which is most likely attributable to the rapid increase in ice discharge in the Getz, Pine Island, Thwaites, and Crosson areas in this region (Figure 3). The WAIS maintained a high state of mass loss in 2017–2020 with no indication of slowing down, suggesting that this region will need continued monitoring.

4. Discussion

4.1. Comparison with Other Estimates

In order to evaluate the differences of mass balance assessments utilizing different methods, we conducted a comparative analysis of our estimates against those results from the gravity-based method, the altimetry-based method, and other IOM-based estimates (Table 4).

4.1.1. Comparison with Gravimetry-Derived Mass Change

We compared the results from our assessment with the mass balance calculated from gravimetry measurements spanning the years 2013–2020 [10]. We constructed a time series of mass balance data for the AIS from both our results and those obtained through the gravity calculations. During the 2013–2020 period, our estimated AIS mass loss was 595 Gt, while gravity results indicated a mass loss of 567 Gt for the period January 2013–September 2019 (Figure 7a). The results obtained through these two methods were found to be in close agreement. In the EAIS, both our results and the gravity inversions indicate that the EAIS was slowly gaining mass. However, our calculated mass increase (489 Gt) appears to be larger than the gravity results (419 Gt), which we speculate is due to basal melting of the ice sheet, as reported in [46]. Previous study has shown that the total basal mass balance over AIS decreased at an average rate of 21 ± 22 Gt/yr from 2003 to 2009, about twice as much mass change as the discrepancy observed in the EAIS. Although the time periods differ, the magnitude of annual total basal mass loss suggests that ice sheet basal melting is non-negligible, particularly when calculating the mass budget using IOM. In addition, recent study has indicated that current models significantly underestimate the basal melting rate in Antarctica [47], providing strong support for our speculation. The WAIS is the region that shows the best agreement between our results and the gravity results (Figure 7c). The AP also exhibits similar trends in our data and the gravity results (Figure 7d). Although the monthly estimates from gravity inversion contain more detail than our annual estimates, the results of the two methods agree well.

4.1.2. Comparison with Altimetry-Derived Mass Change

We conducted a comparison of our results with recent altimetry-based mass balance estimates of the AIS [7]. The altimetry-based results indicate that the mass change rate for the AIS, EAIS, WAIS, and AP between 2003 and 2019 was −118 ± 24 Gt/yr, 90 ± 21 Gt/yr, −169 ± 10 Gt/yr, and −39 ± 5 Gt/yr, respectively. Our estimates were found to be consistent with those of Smith et al. [7] within stated errors. Both results reveal mass gains in the EAIS and significant ice loss in the WAIS. However, the altimetry-based mass balance estimates were found to be more positive for the EAIS and more negative for the WAIS. We speculate that these differences are mainly due to the variations in the area of investigation. Our estimates were solely focused on the grounded ice sheet region, while the altimetry-based estimates also included peripheral islands. To enable a comparison of the two results on the same continental extent, we used our ice sheet boundary to mask the dm/dt results provided by Smith et al. [7]. Subsequently, the altimetry-based mass balance of the AIS, EAIS, WAIS, and AP were found to be −89 ± 24 Gt/yr, 89 ± 21 Gt/yr, −146 ± 10 Gt/yr, and −32 ± 5 Gt/yr, respectively. Considering the slight differences in the study period, our results agree quite well with the altimetry-based estimates.

4.1.3. Comparison with Other IOM-Based Results

The results derived from different IOM schemes have been found to exhibit substantial variation, primarily due to the different designs of the experiment, as reflected in the calculation of ice discharge. Our estimate of ice discharge of 1908 ± 51 Gt/yr during the 2013–2017 period differs from other estimates. The ice discharge values of the AIS estimated by Gardner et al. [16], Shen et al. [17], and Rignot et al. [8] were 1929 ± 40 Gt/yr (2015), 2136 ± 42 Gt/yr (2014–2015), and 2311 ± 142 Gt/yr (2013–2017), respectively. Our results were in good agreement with Gardner et al. [16], which may be due to the fact that we both estimated the AIS mass balance using a similar flux gate. On the other hand, differences in ice thickness data may be responsible for the discrepancy between our results and those reported by Shen et al. [17] and Rignot et al. [8]. Shen et al. [17] and Rignot et al. [8] primarily utilized data from Bedmap-2 [1] and Operation Ice Bridge, while we estimated the ice thickness from BedMachine. We speculate that the differences in the method and data used may have led to the discrepancy between different IOM studies.
Furthermore, it should be noted that the ice velocity products used in this study differ from previous studies. Thus, we carry out an analysis of the regional ice velocity differences from previous studies to determine if any variation observed are real or mere artifacts. In the EAIS, we note a continuous slight decrease in ice velocity from 2007–2011 to 2017–2020 on the Totten Glacier, which is distinct from the velocity increase recorded during the period of 2000–2007, as reported in a previous study [18]. This reduction in ice velocity has also been observed by Miles et al. [48], and is thought to result from intermittent contact between the ice shelf and bed obstacles. Note that an unusual increase in velocity at the north edge of Law Dome was observed in our results. We speculate this is just artifacts and not a real velocity variation. Furthermore, these artefacts have no impact on our ice discharge estimate as they are located far from the flux gate. On Ninnis Glacier and Mertz Glacier, we found acceleration mainly occurred at the floating ice shelves from 2007–2011 to 2017–2020, which was consistent with the results from Rignot et al. [28]. Additionally, a study of the calving cycle of the Ninnis Glacier found that its acceleration was closely related to the calving event and that the velocity at the grounding line was almost constant from 2013–2019 [49], which aligned with the results obtained in this study. Overall, the results on ice velocity changes in this study were consistent with most previous studies at the regional scale, although some velocity artifacts were present in the velocity products with little impact on the calculation of ice discharge.

4.2. Most Recent Changes

The WAIS exhibits pronounced changes in ice discharge, with an increase from 652 Gt/yr of the earlier 2000 period to 753 Gt/yr of the recent 2017–2020 period. About 70% of the discharge growth is attributed to the glacier speedup in the Amundsen Sea sector. As the region exhibiting the most dramatic changes, the Amundsen Sea sector merits further examination. The current ice velocity in the region has been found to be dramatically increasing (Figure 3), the ice thickness is decreasing (Figure 2), the rate of ice discharge is significantly increasing over time (Table 2), and the mass loss is significant (Table 3). The dramatic changes observed in this sector can be attributed to several factors. The primary cause is believed to be the wind-driven incursion of warm, salty, Circumpolar Deep Water (CDW) beneath the ice shelf [50], which leads to increased ice shelf basal melt rates and resultant changes in ice shelf buttressing [51,52]. In addition, the intrusion of maritime air masses over the ice sheet may result in extensive melt events [53].
B12–B14 in the East Indian Ocean sector have been identified as notable areas where numerous studies have reported mass losses [18,38]. Our study also found that this region exhibited a negative mass balance. In the period 2017–2020, we observed a slight decrease in ice velocity and a cessation of the continuous increase in ice discharge, as previously reported in another study [48]. However, the mass loss in this region remained higher than it was in the period 2013–2017 (Table 3), which suggests that a decrease in discharge alone is unlikely to have changed the state of the mass balance. The relatively low mass loss observed in 2013–2017 may have been related to the extreme El Niño that occurred during 2015–2016, which altered moisture transport patterns and resulted in increased precipitation in the coastal area [54]. The reasons for the decrease in the ice velocity in this region warrant further investigation.

5. Conclusions

In this study, we estimated the mass balance of the AIS for the period 2000 to 2020 using the IOM. We used interannual-scale Antarctic-wide ice velocity and state-of-the-art AIS thickness datasets and extended RACMO2.3 SMB records. We found that the overall discharge in the AIS showed an increasing trend, with the discharges in 2000, 2007–2011, 2013–2017, and 2017–2020 being 1792 ± 47 Gt/yr, 1876 ± 51 Gt/yr, 1908 ± 51 Gt/yr, and 1940 ± 37 Gt/yr, respectively. The discharge of the WAIS increased three to four times more than that of the EAIS in the same period. The Amundsen Sea sector was the most significant contributor to the increase in ice discharge. Additionally, for the East Indian Ocean sector (B12–B14), we found a slight decrease in ice discharge during 2017–2020 instead of a continuation of the upward trend of the previous period. The mass balance of the AIS was in a state of mass loss in the early 21st Century, with a mean value of –89 ± 99 Gt/yr over the study period. The EAIS showed a slightly positive mass balance, with some basins in the East Indian Ocean (B12–B14) showing mass losses. The mass loss in the WAIS was significant, accounting for approximately 70% of the total Antarctic mass loss and exhibiting a dramatic increasing trend. Furthermore, our mass balance estimates agreed with gravity inversion and altimetry-based results. In general, the discharges of the AIS have continued to grow in recent years (Figure 3). The discharge of B22 in the Amundsen Sea sector seems to have stopped growing between 2007–2011 and 2013–2017, but then increased again during 2017–2020. The East Indian Ocean sector (B12–B14) showed a decreased discharge during 2017–2020 (Table 2). The mass balance in 2017–2020 was –99 ± 93 Gt/yr, slightly above the average value of the study period.

Author Contributions

Conceptualization, Q.L.; methodology, Q.L. and T.Y.; software, T.Y.; validation, T.Y.; formal analysis, Q.L. and T.Y.; investigation, Q.L. and T.Y.; resources, X.C.; data curation, T.Y.; writing—original draft preparation, T.Y.; writing—review and editing, Z.C., F.H., X.C., T.L., L.Z. and Q.L.; visualization, T.Y.; supervision, Q.L.; project administration, Q.L.; funding acquisition, X.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (grant no. 42006192), Guangdong Basic and Applied Basic Research Foundation (2023A1515010925), and the Innovation Group Project of Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) (grant no. 311021008).

Data Availability Statement

Multiple ice velocity datasets, ice thickness, grounding lines, flux gates and basin boundary are available from the National Snow and Ice Data Center (NSIDC). The ERA5-Land are available from the European Centre for Medium-Range Weather Forecasts (ECMWF).

Acknowledgments

The authors thank the NSIDC for providing datasets. The authors thank the ECMWF for providing reanalysis product. The authors also thank J. M. van Wessem for providing the RACMO 2.3p1/p2 dataset.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fretwell, P.; Pritchard, H.D.; Vaughan, D.G.; Bamber, J.L.; Barrand, N.E.; Bell, R.; Bianchi, C.; Bingham, R.G.; Blankenship, D.D.; Casassa, G.; et al. Bedmap2: Improved Ice Bed, Surface and Thickness Datasets for Antarctica. Cryosphere 2013, 7, 375–393. [Google Scholar] [CrossRef] [Green Version]
  2. Bamber, J.L.; Westaway, R.M.; Marzeion, B.; Wouters, B. The Land Ice Contribution to Sea Level during the Satellite Era. Environ. Res. Lett. 2018, 13, 063008. [Google Scholar] [CrossRef]
  3. Morlighem, M.; Rignot, E.; Binder, T.; Blankenship, D.; Drews, R.; Eagles, G.; Eisen, O.; Ferraccioli, F.; Forsberg, R.; Fretwell, P.; et al. Deep Glacial Troughs and Stabilizing Ridges Unveiled beneath the Margins of the Antarctic Ice Sheet. Nat. Geosci. 2020, 13, 132–137. [Google Scholar] [CrossRef] [Green Version]
  4. The IMBIE team Mass Balance of the Antarctic Ice Sheet from 1992 to 2017. Nature 2018, 558, 219–222. [CrossRef] [Green Version]
  5. Thomas, E.R.; van Wessem, J.M.; Roberts, J.; Isaksson, E.; Schlosser, E.; Fudge, T.J.; Vallelonga, P.; Medley, B.; Lenaerts, J.; Bertler, N.; et al. Regional Antarctic Snow Accumulation over the Past 1000 Years. Clim. Past 2017, 13, 1491–1513. [Google Scholar] [CrossRef] [Green Version]
  6. Wadham, J.L.; Hawkings, J.R.; Tarasov, L.; Gregoire, L.J.; Spencer, R.G.M.; Gutjahr, M.; Ridgwell, A.; Kohfeld, K.E. Ice Sheets Matter for the Global Carbon Cycle. Nat. Commun. 2019, 10, 3567. [Google Scholar] [CrossRef] [Green Version]
  7. Smith, B.; Fricker, H.A.; Gardner, A.S.; Medley, B.; Nilsson, J.; Paolo, F.S.; Holschuh, N.; Adusumilli, S.; Brunt, K.; Csatho, B.; et al. Pervasive Ice Sheet Mass Loss Reflects Competing Ocean and Atmosphere Processes. Science 2020, 368, 1239–1242. [Google Scholar] [CrossRef] [PubMed]
  8. Rignot, E.; Mouginot, J.; Scheuchl, B.; van den Broeke, M.; van Wessem, M.J.; Morlighem, M. Four Decades of Antarctic Ice Sheet Mass Balance from 1979–2017. Proc. Natl. Acad. Sci. USA 2019, 116, 1095–1103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Schröder, L.; Horwath, M.; Dietrich, R.; Helm, V.; van den Broeke, M.R.; Ligtenberg, S.R.M. Four Decades of Antarctic Surface Elevation Changes from Multi-Mission Satellite Altimetry. Cryosphere 2019, 13, 427–449. [Google Scholar] [CrossRef] [Green Version]
  10. Velicogna, I.; Mohajerani, Y.; Landerer, F.; Mouginot, J.; Noel, B.; Rignot, E.; Sutterley, T.; Broeke, M.; Wessem, M.; Wiese, D. Continuity of Ice Sheet Mass Loss in Greenland and Antarctica From the GRACE and GRACE Follow-On Missions. Geophys. Res. Lett. 2020, 47, e2020GL087291. [Google Scholar] [CrossRef] [Green Version]
  11. Rignot, E.; Thomas, R.H. Mass Balance of Polar Ice Sheets. Science 2002, 297, 1502–1506. [Google Scholar] [CrossRef] [Green Version]
  12. Zwally, H.J.; Giovinetto, M.B.; Li, J.; Cornejo, H.G.; Beckley, M.A.; Brenner, A.C.; Saba, J.L.; Yi, D. Mass Changes of the Greenland and Antarctic Ice Sheets and Shelves and Contributions to Sea-Level Rise: 1992–2002. J. Glaciol. 2005, 51, 509–527. [Google Scholar] [CrossRef] [Green Version]
  13. Tapley, B.D.; Bettadpur, S.; Ries, J.C.; Thompson, P.F.; Watkins, M.M. GRACE Measurements of Mass Variability in the Earth System. Science 2004, 305, 503–505. [Google Scholar] [CrossRef] [Green Version]
  14. Rignot, E.; Bamber, J.L.; van den Broeke, M.R.; Davis, C.; Li, Y.; van de Berg, W.J.; van Meijgaard, E. Recent Antarctic Ice Mass Loss from Radar Interferometry and Regional Climate Modelling. Nat. Geosci. 2008, 1, 106–110. [Google Scholar] [CrossRef] [Green Version]
  15. Rignot, E.; Velicogna, I.; van den Broeke, M.R.; Monaghan, A.; Lenaerts, J.T.M. Acceleration of the Contribution of the Greenland and Antarctic Ice Sheets to Sea Level Rise: Acceleration of Ice Sheet Loss. Geophys. Res. Lett. 2011, 38, L05503. [Google Scholar] [CrossRef] [Green Version]
  16. Gardner, A.S.; Moholdt, G.; Scambos, T.; Fahnstock, M.; Ligtenberg, S.; van den Broeke, M.; Nilsson, J. Increased West Antarctic and Unchanged East Antarctic Ice Discharge over the Last 7 Years. Cryosphere 2018, 12, 521–547. [Google Scholar] [CrossRef] [Green Version]
  17. Shen, Q.; Wang, H.; Shum, C.K.; Jiang, L.; Hsu, H.T.; Dong, J. Recent High-Resolution Antarctic Ice Velocity Maps Reveal Increased Mass Loss in Wilkes Land, East Antarctica. Sci. Rep. 2018, 8, 4477. [Google Scholar] [CrossRef] [Green Version]
  18. 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. 2016, 43, 6366–6373. [Google Scholar] [CrossRef]
  19. Zhou, C.; Liang, Q.; Chen, Y.; Lei, H.; Fu, Z.; Zheng, L.; Liu, R. Mass Balance Assessment of the Amery Ice Shelf Basin, East Antarctica. Earth Space Sci. 2019, 6, 1987–1999. [Google Scholar] [CrossRef]
  20. Cui, X.; Du, W.; Xie, H.; Sun, B. The Ice Flux to the Lambert Glacier and Amery Ice Shelf along the Chinese Inland Traverse and Implications for Mass Balance of the Drainage Basins, East Antarctica. Polar Res. 2020, 39, 3582. [Google Scholar] [CrossRef]
  21. Xu, D.; Tang, X.; Yang, S.; Zhang, Y.; Wang, L.; Li, L.; Sun, B. Revisiting Ice Flux and Mass Balance of the Lambert Glacier–Amery Ice Shelf System Using Multi-Remote-Sensing Datasets, East Antarctica. Remote Sens. 2022, 14, 391. [Google Scholar] [CrossRef]
  22. Selley, H.L.; Hogg, A.E.; Cornford, S.; Dutrieux, P.; Shepherd, A.; Wuite, J.; Floricioiu, D.; Kusk, A.; Nagler, T.; Gilbert, L.; et al. Widespread Increase in Dynamic Imbalance in the Getz Region of Antarctica from 1994 to 2018. Nat. Commun. 2021, 12, 1133. [Google Scholar] [CrossRef]
  23. Diener, T.; Sasgen, I.; Agosta, C.; Fürst, J.J.; Braun, M.H.; Konrad, H.; Fettweis, X. Acceleration of Dynamic Ice Loss in Antarctica From Satellite Gravimetry. Front. Earth Sci. 2021, 9, 741789. [Google Scholar] [CrossRef]
  24. Mouginot, J.; Rignot, E.; Scheuchl, B.; Millan, R. Comprehensive Annual Ice Sheet Velocity Mapping Using Landsat-8, Sentinel-1, and RADARSAT-2 Data. Remote Sens. 2017, 9, 364. [Google Scholar] [CrossRef] [Green Version]
  25. Gardner, A.S.; Fahnestock, M.A.; Scambos, T.A. MEaSUREs ITS_LIVE Landsat Image-Pair Glacier and Ice Sheet Surface Velocities: Version 1; National Snow and Ice Data Center: Boulder, CO, USA, 2019. [Google Scholar] [CrossRef]
  26. Scheuchl, B.; Mouginot, J.; Rignot, E. MEaSUREs InSAR-Based Ice Velocity Maps of Central Antarctica: 1997 and 2009, Version 1; NASA National Snow and Ice Data Center Distributed Active Archive Center: Boulder, CO, USA, 2012. [Google Scholar] [CrossRef]
  27. Rignot, E.; Mouginot, J.; Scheuchl, B. MEaSUREs InSAR-Based Antarctica Ice Velocity Map, Version 2. [Indicate Subset Used]; NASA National Snow and Ice Data Center Distributed Active Archive Center: Boulder, CO, USA, 2017. [Google Scholar] [CrossRef]
  28. Rignot, E.; Mouginot, J.; Scheuchl, B.; Jeong, S. Changes in Antarctic Ice Sheet Motion Derived From Satellite Radar Interferometry Between 1995 and 2022. Geophys. Res. Lett. 2022, 49, e2022GL100141. [Google Scholar] [CrossRef]
  29. Zwally, H.J.; Giovinetto, M.B.; Beckley, M.A.; Saba, J.L. Antarctic and Greenland Drainage Systems. GSFC Cryospheric Sciences Laboratory. 2012. Available online: http://icesat4.gsfc.nasa.gov/cryo_data/ant_grn_drainage_systems.php (accessed on 30 May 2022).
  30. Depoorter, M.A.; Bamber, J.L.; Griggs, J.A.; Lenaerts, J.T.M.; Ligtenberg, S.R.M.; van den Broeke, M.R.; Moholdt, G. Calving Fluxes and Basal Melt Rates of Antarctic Ice Shelves. Nature 2013, 502, 89–92. [Google Scholar] [CrossRef]
  31. Griggs, J.A.; Bamber, J.L. Antarctic Ice-Shelf Thickness from Satellite Radar Altimetry. J. Glaciol. 2011, 57, 485–498. [Google Scholar] [CrossRef] [Green Version]
  32. Morlighem, M.; Rignot, E.; Seroussi, H.; Larour, E.; Ben Dhia, H.; Aubry, D. A Mass Conservation Approach for Mapping Glacier Ice Thickness: BALANCE THICKNESS. Geophys. Res. Lett. 2011, 38. [Google Scholar] [CrossRef] [Green Version]
  33. Agosta, C.; Amory, C.; Kittel, C.; Orsi, A.; Favier, V.; Gallée, H.; van den Broeke, M.R.; Lenaerts, J.T.M.; van Wessem, J.M.; van de Berg, W.J.; et al. Estimation of the Antarctic Surface Mass Balance Using the Regional Climate Model MAR (1979–2015) and Identification of Dominant Processes. Cryosphere 2019, 13, 281–296. [Google Scholar] [CrossRef] [Green Version]
  34. Gossart, A.; Helsen, S.; Lenaerts, J.T.M.; Broucke, S.V.; van Lipzig, N.P.M.; Souverijns, N. An Evaluation of Surface Climatology in State-of-the-Art Reanalyses over the Antarctic Ice Sheet. J. Clim. 2019, 32, 6899–6915. [Google Scholar] [CrossRef]
  35. Wessem, J.M.V.; van de Berg, W.J.; Noël, B.P.Y.; van Meijgaard, E.; Amory, C.; Birnbaum, G.; Jakobs, C.L.; Krüger, K.; Lenaerts, J.T.M.; Lhermitte, S.; et al. Modelling the Climate and Surface Mass Balance of Polar Ice Sheets Using RACMO2—Part 2: Antarctica (1979–2016). Cryosphere 2018, 12, 1479–1498. [Google Scholar] [CrossRef] [Green Version]
  36. Lenaerts, J.T.M.; Medley, B.; van den Broeke, M.R.; Wouters, B. Observing and Modeling Ice Sheet Surface Mass Balance. Rev. Geophys. 2019, 57, 376–420. [Google Scholar] [CrossRef] [Green Version]
  37. Wessem, J.M.V.; Reijmer, C.H.; Morlighem, M.; Mouginot, J.; Rignot, E.; Medley, B.; Joughin, I.; Wouters, B.; Depoorter, M.A.; Bamber, J.L.; et al. Improved Representation of East Antarctic Surface Mass Balance in a Regional Atmospheric Climate Model. J. Glaciol. 2014, 60, 761–770. [Google Scholar] [CrossRef] [Green Version]
  38. Mohajerani, Y.; Velicogna, I.; Rignot, E. Mass Loss of Totten and Moscow University Glaciers, East Antarctica, Using Regionally Optimized GRACE Mascons. Geophys. Res. Lett. 2018, 45, 7010–7018. [Google Scholar] [CrossRef] [Green Version]
  39. Mohajerani, Y.; Velicogna, I.; Rignot, E. Evaluation of Regional Climate Models Using Regionally Optimized GRACE Mascons in the Amery and Getz Ice Shelves Basins, Antarctica. Geophys. Res. Lett. 2019, 46, 13883–13891. [Google Scholar] [CrossRef] [Green Version]
  40. Muñoz-Sabater, J.; Dutra, E.; Agustí-Panareda, A.; Albergel, C.; Arduini, G.; Balsamo, G.; Boussetta, S.; Choulga, M.; Harrigan, S.; Hersbach, H.; et al. ERA5-Land: A State-of-the-Art Global Reanalysis Dataset for Land Applications. Earth Syst. Sci. Data 2021, 13, 4349–4383. [Google Scholar] [CrossRef]
  41. Haran, T.; Klinger, M.; Bohlander, J.; Fahnestock, M.; Painter, T.; Scambos, T. MEaSUREs MODIS Mosaic of Antarctica 2013–2014 (MOA2014) Image Map, Version 1. [Surface Morphology]; NASA National Snow and Ice Data Center Distributed Active Archive Center: Boulder, CO, USA, 2018. [Google Scholar] [CrossRef]
  42. Cheng, Y.; Xia, M.; Qiao, G.; Lv, D.; Li, Y.; Hai, G. Imminent Calving Accelerated by Increased Instability of the Brunt Ice Shelf, in Response to Climate Warming. Earth Planet. Sci. Lett. 2021, 572, 117132. [Google Scholar] [CrossRef]
  43. Mottram, R.; Hansen, N.; Kittel, C.; van Wessem, J.M.; Agosta, C.; Amory, C.; Boberg, F.; van de Berg, W.J.; Fettweis, X.; Gossart, A.; et al. What Is the Surface Mass Balance of Antarctica? An Intercomparison of Regional Climate Model Estimates. Cryosphere 2021, 15, 3751–3784. [Google Scholar] [CrossRef]
  44. Conway, H.; Rasmussen, L.A. Recent Thinning and Migration of the Western Divide, Central West Antarctica. Geophys. Res. Lett. 2009, 36, L12502. [Google Scholar] [CrossRef] [Green Version]
  45. Deb, P.; Orr, A.; Bromwich, D.H.; Nicolas, J.P.; Turner, J.; Hosking, J.S. Summer Drivers of Atmospheric Variability Affecting Ice Shelf Thinning in the Amundsen Sea Embayment, West Antarctica. Geophys. Res. Lett. 2018, 45, 4124–4133. [Google Scholar] [CrossRef]
  46. Kang, J.; Lu, Y.; Li, Y.; Zhang, Z.; Shi, H. Antarctic Basal Water Storage Variation Inferred from Multi-Source Satellite Observation and Relevant Models. Remote Sens. 2022, 14, 2337. [Google Scholar] [CrossRef]
  47. Artemieva, I.M. Antarctica Ice Sheet Basal Melting Enhanced by High Mantle Heat. Earth-Sci. Rev. 2022, 226, 103954. [Google Scholar] [CrossRef]
  48. Miles, B.W.J.; Stokes, C.R.; Jamieson, S.S.R.; Jordan, J.R.; Gudmundsson, G.H.; Jenkins, A. High Spatial and Temporal Variability in Antarctic Ice Discharge Linked to Ice Shelf Buttressing and Bed Geometry. Sci. Rep. 2022, 12, 10968. [Google Scholar] [CrossRef] [PubMed]
  49. Cheng, Y.; Xia, M.; Qiao, G.; Li, Y.; Hai, G.; Lv, D. Calving Cycle of Ninnis Glacier over the Last 60 Years. Int. J. Appl. Earth Obs. Geoinf. 2021, 105, 102612. [Google Scholar] [CrossRef]
  50. Thoma, M.; Jenkins, A.; Holland, D.; Jacobs, S. Modelling Circumpolar Deep Water Intrusions on the Amundsen Sea Continental Shelf, Antarctica. Geophys. Res. Lett. 2008, 35, L18602. [Google Scholar] [CrossRef] [Green Version]
  51. Jacobs, S.S.; Jenkins, A.; Giulivi, C.F.; Dutrieux, P. Stronger Ocean Circulation and Increased Melting under Pine Island Glacier Ice Shelf. Nat. Geosci. 2011, 4, 519–523. [Google Scholar] [CrossRef]
  52. MacGregor, J.A.; Catania, G.A.; Markowski, M.S.; Andrews, A.G. Widespread Rifting and Retreat of Ice-Shelf Margins in the Eastern Amundsen Sea Embayment between 1972 and 2011. J. Glaciol. 2012, 58, 458–466. [Google Scholar] [CrossRef] [Green Version]
  53. Scott, R.C.; Nicolas, J.P.; Bromwich, D.H.; Norris, J.R.; Lubin, D. Meteorological Drivers and Large-Scale Climate Forcing of West Antarctic Surface Melt. J. Clim. 2019, 32, 665–684. [Google Scholar] [CrossRef]
  54. Bodart, J.A.; Bingham, R.J. The Impact of the Extreme 2015–2016 El Niño on the Mass Balance of the Antarctic Ice Sheet. Geophys. Res. Lett. 2019, 46, 13862–13871. [Google Scholar] [CrossRef]
Figure 1. Ice velocity, flux gate (FG), and basin boundaries (B1–B27). The ice velocity is from the MEaSUREs annual ice velocity for 2019–2020 [24]. The red dashed line represents the FG [16], and the solid black line represents the basin boundaries and the grounding line (GL). The numbers 1–27 represent different basins [29].
Figure 1. Ice velocity, flux gate (FG), and basin boundaries (B1–B27). The ice velocity is from the MEaSUREs annual ice velocity for 2019–2020 [24]. The red dashed line represents the FG [16], and the solid black line represents the basin boundaries and the grounding line (GL). The numbers 1–27 represent different basins [29].
Remotesensing 15 01677 g001
Figure 2. (a) The BedMachine ice thickness of the AIS [3]. (b) Surface elevation changes of the AIS during 2013–2020.
Figure 2. (a) The BedMachine ice thickness of the AIS [3]. (b) Surface elevation changes of the AIS during 2013–2020.
Remotesensing 15 01677 g002
Figure 3. Ice discharge changes in different sea sectors and different basins. The bar charts show the anomaly from the mean value (left side of the bar chart) for different periods (shown in different colors). The circles show the discharge changes for different sea sectors between 2017–2020 and 2000. The bar length and the circle size are proportional to the absolute magnitude of the change in discharge. The circle colors represent the negative (blue) or positive (red) state. The base map is the MODIS mosaic of Antarctica [41], and the black lines are the basin boundaries [29].
Figure 3. Ice discharge changes in different sea sectors and different basins. The bar charts show the anomaly from the mean value (left side of the bar chart) for different periods (shown in different colors). The circles show the discharge changes for different sea sectors between 2017–2020 and 2000. The bar length and the circle size are proportional to the absolute magnitude of the change in discharge. The circle colors represent the negative (blue) or positive (red) state. The base map is the MODIS mosaic of Antarctica [41], and the black lines are the basin boundaries [29].
Remotesensing 15 01677 g003
Figure 4. Antarctic-wide ice velocity changes in four different periods: (a) from 2000 to 2007–2011; (b) from 2007–2011 to 2013–2017; (c) from 2013–2017 to 2017–2020; (d) from 2000 to 2017–2020. The black rectangles show the locations of R1–R7 in (e). The black line is the basin boundary produced with the flux gate. (e) Seven regions (R1–R7) with significant changes.
Figure 4. Antarctic-wide ice velocity changes in four different periods: (a) from 2000 to 2007–2011; (b) from 2007–2011 to 2013–2017; (c) from 2013–2017 to 2017–2020; (d) from 2000 to 2017–2020. The black rectangles show the locations of R1–R7 in (e). The black line is the basin boundary produced with the flux gate. (e) Seven regions (R1–R7) with significant changes.
Remotesensing 15 01677 g004
Figure 5. Comparison of annual SMB estimated from ERA5-Land and RACMO2.3 output. (a) The SMB of ERA5-Land against the SMB from RACMO2.3p1 in the EAIS from 1979 to 2015. (b) The SMB of ERA5-Land against the SMB from RACMO2.3p2 in the WAIS from 1979 to 2017. (c) The SMB of ERA5-Land against the SMB from RACMO2.3p2 in the AP from 1979 to 2017.
Figure 5. Comparison of annual SMB estimated from ERA5-Land and RACMO2.3 output. (a) The SMB of ERA5-Land against the SMB from RACMO2.3p1 in the EAIS from 1979 to 2015. (b) The SMB of ERA5-Land against the SMB from RACMO2.3p2 in the WAIS from 1979 to 2017. (c) The SMB of ERA5-Land against the SMB from RACMO2.3p2 in the AP from 1979 to 2017.
Remotesensing 15 01677 g005
Figure 6. Comparison of annual SMB estimated from ERA5-Land and RACMO2.3 output. (a) The SMB of ERA5-Land against the SMB from RACMO2.3 in the AIS. (b) The SMB of ERA5-Land against the SMB from RACMO2.3p1 in the EAIS. (c) The SMB of ERA5-Land against the SMB from RACMO2.3p2 in the WAIS. (d) The SMB of ERA5-Land against the SMB from RACMO2.3p2 in the AP.
Figure 6. Comparison of annual SMB estimated from ERA5-Land and RACMO2.3 output. (a) The SMB of ERA5-Land against the SMB from RACMO2.3 in the AIS. (b) The SMB of ERA5-Land against the SMB from RACMO2.3p1 in the EAIS. (c) The SMB of ERA5-Land against the SMB from RACMO2.3p2 in the WAIS. (d) The SMB of ERA5-Land against the SMB from RACMO2.3p2 in the AP.
Remotesensing 15 01677 g006
Figure 7. Comparison of our mass balance estimates with Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow-On (FO) mass change estimates for 2013–2020 in (a) the AIS, (b) the EAIS, (c) the WAIS, and (d) the AP [10]. The red shaded part in the figure represents the uncertainty range of our calculation, and the red dashed line and light blue dotted line represent the trends of our study’s results and the gravity-based results from 2013 to 2020, respectively. Gray vertical bar indicates the gap between the end of GRACE and the start of GRACE-FO.
Figure 7. Comparison of our mass balance estimates with Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow-On (FO) mass change estimates for 2013–2020 in (a) the AIS, (b) the EAIS, (c) the WAIS, and (d) the AP [10]. The red shaded part in the figure represents the uncertainty range of our calculation, and the red dashed line and light blue dotted line represent the trends of our study’s results and the gravity-based results from 2013 to 2020, respectively. Gray vertical bar indicates the gap between the end of GRACE and the start of GRACE-FO.
Remotesensing 15 01677 g007
Table 1. List of the main ice velocity products and filling ice velocity products used in each of the epochs.
Table 1. List of the main ice velocity products and filling ice velocity products used in each of the epochs.
EpochMEaSUREs Annual Velocity Products [24]Fill Missing Values Products
Central [26]ITS_LIVE [16]Phase [27]
20002000–20012000
2007–20112007–2008,
2008–2009, 2010–2011
2007, 2008, 2010
2013–20172013–2014, 2014–2015
2015–2016, 2016–2017
//
2017–20202017–2018,
2018–2019, 2019–2020
//
Table 2. The area of the 27 basins, ice discharge for the periods of 2000, 2007–2011, 2013–2017, and 2017–2020, and the mean SMB for the study years.
Table 2. The area of the 27 basins, ice discharge for the periods of 2000, 2007–2011, 2013–2017, and 2017–2020, and the mean SMB for the study years.
BasinArea (km2)Discharge (Gt/yr)Mean of SMB (Gt/yr)
20002007–20112013–20172017–2020
1479,840114 ± 25115 ± 20115 ± 20120 ± 7133 ± 27
2780,48850 ± 349 ± 849 ± 450 ± 455 ± 11
31,570,43265 ± 866 ± 567 ± 466 ± 275 ± 15
4244,69640 ± 843 ± 943 ± 1043 ± 548 ± 10
5186,88431 ± 435 ± 534 ± 735 ± 335 ± 7
6606,04059 ± 858 ± 1257 ± 1358 ± 670 ± 14
7493,85757 ± 963 ± 1263 ± 1167 ± 888 ± 18
8159,82921 ± 523 ± 822 ± 522 ± 335 ± 7
9145,67416 ± 418 ± 617 ± 317 ± 218 ± 4
10919,36243 ± 243 ± 242 ± 242 ± 140 ± 8
11257,47715 ± 314 ± 614 ± 314 ± 115 ± 3
12719,472114 ± 14119 ± 18124 ± 15124 ± 11121 ± 24
131,109,722208 ± 14212 ± 18218 ± 19216 ± 16191 ± 38
14711,661127 ± 11128 ± 16134 ± 15133 ± 14129 ± 26
15124,08131 ± 726 ± 730 ± 630 ± 628 ± 6
16265,50116 ± 315 ± 314 ± 215 ± 112 ± 2
171,843,82478 ± 1171 ± 1463 ± 1368 ± 982 ± 16
18266,47111 ± 111 ± 29 ± 210 ± 129 ± 6
19372,69648 ± 447 ± 647 ± 447 ± 344 ± 9
20182,833148 ± 15148 ± 10168 ± 17170 ± 15124 ± 25
21210,994155 ± 10175 ± 10183 ± 11192 ± 11106 ± 21
22213,74498 ± 4131 ± 3131 ± 4139 ± 387 ± 17
2375,15978 ± 880 ± 677 ± 1075 ± 873 ± 15
24100,86583 ± 885 ± 588 ± 788 ± 690 ± 18
2534,21248 ± 761 ± 860 ± 959 ± 871 ± 14
2641,51928 ± 430 ± 327 ± 428 ± 344 ± 9
2752,13613 ± 312 ± 112 ± 212 ± 121 ± 4
EAIS10,139,000968 ± 32980 ± 26991 ± 311000 ± 221042 ± 64
WAIS1,801,737652 ± 32707 ± 42730 ± 39753 ± 29596 ± 49
AP228,732172 ± 12188 ± 10187 ± 12187 ± 10226 ± 25
Total12,169,4691792 ± 471876 ± 511908 ± 511940 ± 371864 ± 85
Table 3. Mass balance in the four periods 2000, 2007–2011, 2013–2017, and 2017–2020 for the 27 basins.
Table 3. Mass balance in the four periods 2000, 2007–2011, 2013–2017, and 2017–2020 for the 27 basins.
BasinMass Balance (Gt/yr)
20002007–20112013–20172017–2020Mean
113 ± 3615 ± 3321 ± 3414 ± 2817 ± 32
23 ± 114 ± 137 ± 128 ± 126 ± 12
311 ± 176 ± 1510 ± 1611 ± 169 ± 16
41 ± 115 ± 134 ± 149 ± 125 ± 13
5−1 ± 7−3 ± 84 ± 103 ± 81 ± 9
613 ± 173 ± 1718 ± 2013 ± 1512 ± 18
720 ± 1820 ± 2123 ± 2131 ± 2124 ± 21
89 ± 811 ± 1013 ± 916 ± 913 ± 9
9−3 ± 40 ± 70 ± 52 ± 41 ± 6
10−7 ± 70 ± 9−5 ± 8−2 ± 8−2 ± 8
11−1 ± 41 ± 61 ± 51 ± 31 ± 5
123 ± 2814 ± 32−6 ± 29−12 ± 25−1 ± 28
1312 ± 46−27 ± 41−18 ± 44−40 ± 39−24 ± 42
14−14 ± 25−20 ± 2710 ± 325 ± 31−2 ± 30
150 ± 9−2 ± 80 ± 9−1 ± 8−1 ± 8
16−2 ± 4−4 ± 4−4 ± 30 ± 3−3 ± 3
1712 ± 215 ± 2113 ± 2025 ± 2114 ± 20
1817 ± 615 ± 621 ± 721 ± 619 ± 6
19−5 ± 9−7 ± 10−1 ± 10−1 ± 10−3 ± 10
20−13 ± 31−37 ± 26−44 ± 30−38 ± 30−37 ± 29
21−28 ± 27−74 ± 23−82 ± 23−81 ± 25−75 ± 24
224 ± 21−45 ± 17−46 ± 18−53 ± 18−43 ± 17
233 ± 18−2 ± 17−7 ± 17−7 ± 16−5 ± 17
243 ± 194 ± 199 ± 21−7 ± 173 ± 19
25 *−10 ± 13−10 ± 13−10 ± 13−10 ± 13−10 ± 13
26 *−16 ± 10−16 ± 10−16 ± 10−16 ± 10−16 ± 10
278 ± 57 ± 410 ± 59 ± 49 ± 4
EAIS56 ± 6213 ± 5570 ± 5870 ± 5553 ± 56
WAIS−9 ± 73−135 ± 75−137 ± 77−145 ± 71−127 ± 74
AP−15 ± 25−15 ± 28−7 ± 29−24 ± 26−15 ± 32
AIS33 ± 99−137 ± 97−75 ± 101−99 ± 93−89 ± 99
* The final mass balance for the northern AP (B25 and B26) is determined by the mass balance estimate of the altimetric method [7].
Table 4. Results of mass balance from our IOM-based estimate, the gravity-based method, and the altimetry-based method.
Table 4. Results of mass balance from our IOM-based estimate, the gravity-based method, and the altimetry-based method.
MethodRegionSpan (Years)Mean of Mass Balance (Gt/yr)Uncertanity
(Gt/yr)
IOMEAIS2013–20207066
WAIS2013–2020−14161
AP2013–2020−1424
AIS2013–2020−8593
GravimetryEAIS2013–202050/
WAIS2013–2020−118/
AP2013–2020−13/
AIS2013–2020−81/
Altimetry *EAIS2003–20199021
WAIS2003–2019−16910
AP2003–2019−395
AIS2003–2019−11824
* Note that the values are calculated combining our ice cover boundary as a mask with the dm/dt results provided by Smith et al. [7].
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, T.; Liang, Q.; Zheng, L.; Li, T.; Chen, Z.; Hui, F.; Cheng, X. Mass Balance of the Antarctic Ice Sheet in the Early 21st Century. Remote Sens. 2023, 15, 1677. https://doi.org/10.3390/rs15061677

AMA Style

Yang T, Liang Q, Zheng L, Li T, Chen Z, Hui F, Cheng X. Mass Balance of the Antarctic Ice Sheet in the Early 21st Century. Remote Sensing. 2023; 15(6):1677. https://doi.org/10.3390/rs15061677

Chicago/Turabian Style

Yang, Tian, Qi Liang, Lei Zheng, Teng Li, Zhuoqi Chen, Fengming Hui, and Xiao Cheng. 2023. "Mass Balance of the Antarctic Ice Sheet in the Early 21st Century" Remote Sensing 15, no. 6: 1677. https://doi.org/10.3390/rs15061677

APA Style

Yang, T., Liang, Q., Zheng, L., Li, T., Chen, Z., Hui, F., & Cheng, X. (2023). Mass Balance of the Antarctic Ice Sheet in the Early 21st Century. Remote Sensing, 15(6), 1677. https://doi.org/10.3390/rs15061677

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop