Front calving is a primary mechanism through which Antarctic ice shelves discharge ice mass into the Southern Ocean. It is an important process that influences ice shelf stability and thus, impacts the Antarctic Ice sheet's contribution to global sea level rise. Mélange dynamics inside rifts is recognized to potentially influence the rift propagation and subsequent iceberg calving. However, large-scale, high-resolution three-dimensional (3D) observations are scarce, which leads to their inability to capture small scale rift dynamics. Ultimately, the lack of knowledge in 3D rift structural changes and mélange dynamics hinders our understanding of the role of mélange in ice shelf retreat and mechanisms underlying the weakening of ice shelf stability. We propose an innovative multi-temporal digital elevation model (DEM) adjustment model (MDAM) to build a multi-satellite DEM time series from meter-level resolution small DEMs. It removes biases across large Antarctic ice shelves, as large as ∼ 6 m in elevation, caused by tides, ice flow dynamics, and observation errors. Using 30 Reference Elevation Model of Antarctica (REMA) and ZY-3 sub-DEMs, we establish a cross-shelf DEM time series from 2014 to 2021 for the Filchner-Ronne Ice Shelf, the second largest ice shelf in Antarctica. This unified and integrated DEM series, with an unprecedented submeter elevation accuracy, reveals quantitative 3D structural and mélange features of two ∼ 50 km long rifts, including rift lips, flank surface, pre-mélange cavities, and mélange elevations. For the first time, we have observed that while the mélange elevation in the rifts decreased by 2.7 ± 0.6 m from 2014 to 2021, the mélange within the rifts experienced a rapid expansion of (10.31 ± 0.03) × 109 m3, or 139 %. This expansion is attributed to newly calved shelf ice from rift walls, associated rift widening, and other factors related to rift-mélange interactions. The developed MDAM system and the 3D mélange dynamics analysis methods can be applied for quantifying ice shelf instability and the future contribution of the Antarctic Ice Sheet to global sea level rise.
National Natural Science Foundation of China42394131Ministry of Science and Technology of the People's Republic of China2024YFF0808302Introduction
The Antarctic Ice Sheet (AIS) has shown a persistent pattern of mass loss and has contributed to global sea level rise since the beginning of the satellite earth observation era in the 1960s (Jacobs et al., 1992; Allison et al., 2009; Rignot et al., 2019; Smith et al., 2020; Li et al., 2022, 2023a; Otosaka et al., 2023). The lost ice mass enters the Southern Ocean from ice shelves mainly through two processes, namely shelf front calving and basal melting, each accounting for ∼ 50 % (Depoorter et al., 2013; Liu et al., 2015; Smith et al., 2020; Greene et al., 2022). In contrast to basal melting that occurs on the bottom of ice shelves, calving events are directly observed by using satellite remote sensing in the past decades. Examples of large iceberg calving include a major calving event in the Filchner Ice Shelf in 1986 with an area loss of ∼ 11 500 km2 (Ferrigno and Gould, 1987), release of an iceberg of ∼ 5800 km2 from Larsen C Ice Shelf in 2017 (Larour et al., 2021; Wang et al., 2022), and break out of a part of the “loose tooth” of ∼ 1680 km2 from the Amery Ice Shelf in 2019 (Francis et al., 2021; Walker et al., 2021). Although the mechanism of ice shelf caving is not fully understood, calving processes are clearly related to the evolution of rifts as they advect to the shelf front and cut icebergs (Joughin and MacAyeal, 2005; Hulbe et al., 2010; Walker and Gardner, 2019; Greene et al., 2022). Mélanges inside rifts, which consist of shelf ice, snow, sea ice, and water, have been investigated in relation to ice shelf fracturing in Antarctica and glacier calving in Greenland (Rignot and MacAyeal, 1998; Larour et al., 2004; Cassotto et al., 2021). Specifically, reductions in mélange thickness have been observed during rift widening on the Amery, Ronne, and Larsen C ice shelves (Fricker et al., 2005; Walker et al., 2021; Larour et al., 2021). Furthermore, modelling results indicate that a mélange may have the capability of transmitting stresses across rift flanks and thus, influences rift propagation and shelf front calving (Larour et al., 2004, 2021). However, high-resolution observations covering large regions are scarce and small-scale mélange textures and 3D rift dynamics cannot be captured with a necessary accuracy. This, in turn, hinders our understanding of the role of mélange in ice shelf retreating and the mechanism of ice shelf stability weakening.
(a) High-resolution small-coverage DEMs (e.g., REMA strips of ∼ 16–18 km ×∼ 110–120 km) are unified and integrated for accurate 3D rift structural and mélange dynamic monitoring in a large Antarctic ice shelf environment (e.g., ∼ 165 km wide Filchner Ice Shelf). The LIMA mosaic (Bindschadler et al., 2008) is draped on the RAMP DEM (Liu et al., 2001) with a vertical exaggeration factor of 40 times. (b) Concept of the MDAM system to connect sub-DEMs across the ice shelf through “tie points” (TP, represented by circles) to create a DEM time series. The “ground control points” (GCP, represented by triangles) are used to fix the combined DEM model to grounded regions at ice shelf margins. The closely located red and blue circles (triangles) are from adjacent sub-DEMs and form corresponding TP (GCP) pairs.
In recent decades, altimetric and stereo mapping satellites have been used to acquire 3D data of the Earth surface features, which are further used to derive digital elevation models (DEMs) where 3D rift structure and mélange information is inherently available. DEMs from satellite altimetry data at the Antarctic continental scale have relatively lower resolutions, for example, the 500 m ICESat (DiMarzio et al., 2007) and ICESat-2 (Shen et al., 2022) DEMs, and the 1000 m CryoSat-2 DEM (Slater et al., 2018). The DEMs from optical and SAR satellite stereo mapping data are of higher resolutions and have the potential for geometric modeling and analysis of rifts and mélange features, including the 90 m TanDEM-X DEM (Wessel et al., 2021), the 30 m ASTER GDEM (Tachikawa et al., 2011), and especially, the 2 m REMA DEM (Howat et al., 2019). However, these high-resolution DEMs are derived from relatively narrow swath satellite data and have smaller extents. Thus, the smaller DEMs of different timespans are mosaiced together to cover large rifts and large ice shelves, in which significant discrepancies between individual DEMs are induced by influences of tides, ice flow dynamics, and observation errors (Fig. 1a). Existing methods for co-registration of DEMs are mostly developed based on global fitting and used in a land environment, e.g., a seven-parameter DEM transformation for mountain glacier mass balance estimation (Li et al., 2023b) and a global optimization process for computing three parameters between DEMs using all grid points (Nuth and Kääb, 2011). Furthermore, REMA DEMs are registered to a control data set, such as airborne altimetric data or satellite altimetric data of ICESat and CryoSat-2 (Howat et al., 2019; Shean et al., 2019; Zinck et al., 2023). While existing methods can be used to register DEMs, a mathematical model can help formalize the registration process and provide a consistent framework for quantifying and propagating uncertainties. Such a model-based approach may facilitate systematic comparison between different regions or time periods, and can offer clearer insight into the relationship between registration parameters and mélange dynamics.
We propose an innovative multi-temporal DEM adjustment model (MDAM) with the purpose to remove biases between adjacent sub-DEMs and to establish a unified and integrated multi-satellite DEM time series with a submeter elevation accuracy. We introduce tie points (TPs, see detailed description in Sect. 2.3) to connect the sub-DEMs and to formulate a DEM time series across a large ice shelf, which is then constrained to grounded regions using ground control points (GCPs). Since the GCPs are selected where ICESat-2 altimetric data are available, the connected DEM time series are then geometrically controlled at a centimeter elevation accuracy (Markus et al., 2017). We develop a data-driven approach to remove large offsets between adjacent sub-DEMs, namely the tide-induced elevation displacements and ice flow-induced horizontal discrepancies, by using a tide model and velocity maps, respectively. The residuals are then minimized in the MDAM system in a least-squares process to obtain the 3D correction parameters of each sub-DEM. The validated MDAM is applied to the Filchner Ice Shelf to establish a cross-shelf DEM time series from 2014–2021. The Filchner Ice Shelf (Fig. 1a) is one of the dual ice shelves of Filchner-Ronne Ice Shelf that is the second largest ice shelf in Antarctica. Berkner Island and Coats Land are located at the western and eastern margins of the ice shelf, respectively. The largest calving event recorded for this ice shelf, with an area loss of 11 500 km2, occurred in 1986 due to a rapid propagation of a prominent rift known as “Grand Chasm” (Ferrigno and Gould, 1987), producing three giant icebergs that greatly impacted the circulation and hydrography in the nearby ice shelf-ocean system (Grosfeld et al., 2001) (Fig. A1). Two large rifts, T1 and T2 (Fig. 3), were detected in a Landsat satellite image in 1996 and an ARGON satellite image in 1963, respectively (Li et al., 2017a; Walker and Gardner, 2019; Lv et al., 2022). These rifts have a combined length of ∼ 100 km and exhibit a similar propagation pattern as Grand Chasm. We demonstrate the capability of the MDAM by quantitatively characterizing 3D structural and mélange features of these two ∼ 50 km long rifts, including rift lips, flank surface, pre-mélange cavities, and mélange elevations. Finally, we demonstrate the observed mélange changes and estimated volumetric changes in relation to rift widening as the rifts advect toward the ice shelf front during the study period.
Data and MethodData
The sub-DEMs used in this study are REMA strips and a Ziyuan-3 (ZY-3) DEM which are generated from stereo satellite images of WorldView (Anderson and Marchisio, 2012) and ZY-3 (Wang et al., 2014), respectively (Table A1). Both are formed by the along-track stereo mechanism (Li, 1998). REMA DEMs are available for the entire Antarctica and are produced from WorldView images using an open-source automatic photogrammetric processing software system, Surface Extraction from TIN-based Searchspace Minimization (SETSM) (Noh and Howat, 2015). The WorldView data set includes WorldView-1, -2 and, -3, and GeoEye-1 images with a small swath of 13.1–17.6 km and a high resolution of 0.3–0.5 m (Howat et al., 2019). Thus, each derived REMA DEM provides a high level of ice surface details, but with a relatively small extent, resulting in difficulties in covering large rifts or the entire ice shelf. More specifically, we use REMA DEM strips of 2016–2017 (Fig. 3a) and 2020–2021 (Fig. 3b) with a grid spacing of 2 m. On the other hand, ZY-3, launched in 2012, is the first Chinese high-resolution natural resource mapping satellite (Wang et al., 2014). It is equipped with a three-line mapping system that forms quasi-real-time along-track stereo image pairs using fore (3.5 m), nadir (2.1 m), and aft (3.5 m) views. The fore and aft combination has a large intersection angle of 47°, facilitating a high elevation accuracy for DEM generation. Its wide swath of ∼ 52 km ensures a complete coverage of large rifts, such as T1 and T2 (Fig. 3a). The ZY-3 DEM used in this study has a time tag of 28 February 2014 and a grid spacing of 5 m. It has a relatively narrow strip because of cloud coverage, but it covers both T1 and T2. The DEM is registered to the ICESat-2 altimetry data and has an elevation accuracy of 0.3 m. The ZY-3 DEM of 2014 has a full coverage of T1 and T2 and can provide a reference for comparison with the adjusted sub-DEMs of REMA. The combined ZY-3 and REMA sub-DEMs form a seven-year long time series from 2014 to 2021 with a time interval of 3–4 years, which can be used to analyze mélange elevation and volume changes.
We use elevation points of the ICESat-2 ATL06 data from 2019 to 2021 (Markus et al., 2017) to register the ZY-3 DEM. The GCPs are also selected on or close to ICESat-2 ATL06 elevation points to establish a unified reference system in MDAM. To ensure that the GCPs are “stable” during the timespan between the ZY-3 DEM (REMA sub-DEMs) and ICESat-2, these GCPs are further required to be on grounded features with a low velocity (<10–20 m yr−1). In a special case we also select GCPs on the ice shelf where both the DEM images and ICESat-2 data were acquired within one day. The ICESat-2 data are also used for validation of the bias-corrected DEM time series. Along each of the six beam tracks, ICESat-2 ATL06 data provide aggregated elevation points with a 20 m spacing and a revisit cycle of 91 d. Field validation studies have demonstrated that the ICESat-2 data reached an elevation accuracy within 2–4 cm in Antarctica (Brunt et al., 2019; Magruder et al., 2020; Li et al., 2021).
Data processing principle of the multi-temporal DEM adjustment model (MDAM).
Multi-temporal DEM adjustment model (MDAM)
As illustrated in Fig. 1b, a set of sub-DEMs (DEM1, DEM2, …, DEMk, …, DEMN-1, DEMN) are connected through TPs in the overlapping areas between adjacent sub-DEMs. The two sub-DEMs at the two shelf margins (DEM1 and DEMN) cover grounded regions, while those in between (from DEM2 to DEMN-1) cover the floating ice shelf. GCPs are used to fix the entire DEM time series to the grounded regions. They are input into the multi-temporal DEM adjustment model (MDAM) for removing biases between adjacent sub-DEMs (Fig. 2). Initially, relatively large vertical offsets, up to ∼ 1.2 m, caused by tides are corrected for all sub-DEMs by using a tide model, CATS2008, with an accuracy of ∼ 8 cm (Padman et al., 2008; King et al., 2011). Then we manually select TPs in the overlapping areas between sub-DEMs on the floating ice to cover the entire ice shelf (see Sect. 2.3 for details). The relatively large horizontal displacements between corresponding TP pairs (Fig. 1b), up to hundreds of meters, induced by ice flow dynamics are compensated by using ITS_LIVE velocity maps (Gardner et al., 2019). At the same time, we tie the GCP locations to satellite altimetry data of ICESat-2 (see Sect. 2.3 for details). In the model inconsistencies at TP and GCP pairs are first reduced by the tide and velocity corrections, leaving the uncertainty originated from the sub-DEM production as residuals, including photogrammetric measurement errors and ephemeris data errors. Thereafter, we use the least-squares adjustment method (McGlone et al., 2004) to estimate the unknown bias corrections of all sub-DEMs by minimizing the residuals (Li, 1998; McGlone et al., 2004; Shean et al., 2019). Finally, the corrected sub-DEMs are applied for 3D analysis and modeling of large rifts and mélange dynamics.
For each pair of TPs Pik,k+1 and Pik+1,k (i=1, …, ik,k+1) (Fig. 1b), we measure 3D coordinates (XPik, YPik, ZPik) from DEMk and (XPik+1, YPik+1, ZPik+1) from DEMk+1. We define 3D translational parameters of (dXk, dYk, dZk) as unknowns for bias removal. Thus, we establish the 3D positional relationships at tie point Pik,k+1 between DEMk and DEMk+1 as observation equations for each pair of TPs:
XPik+dXk+VelXik×tk+1-tk=XPik+1+dXk+1+εxikYPik+dYk+VelYik×tk+1-tk=YPik+1+dYk+1+εyikZPik+dZk=ZPik+1+dZk+1+εzik,
where k=1, …, N-1; and i=1, …, ik,k+1. tk and tk+1 are the known data acquisition times for DEMk and DEMk+1. VelXik and VelYik are the ice flow speed values in X- and Y-directions at the TPs given in a velocity map. Here, ITS_LIVE image-pair velocity maps are selected to spatially match the extent of each sub-DEM and to temporally align as closely as possible with the acquisition times of two adjacent sub-DEMs. Similarly, (dXk+1, dYk+1, dZk+1) are unknown bias correction parameters for DEMk+1. (εxik, εyik, εzik) represent residuals at TPs in the observation equations, including photogrammetric measurement errors and ephemeris data errors.
In a similar way, observation equations for ground control point PGCPik are given as:
XGCPik+dXk=XPiGCP,k+εxiGCP,kYGCPik+dYk=YPiGCP,k+εyiGCP,kZGCPik+dZk=ZPiGCP,k+εziGCP,k
where k=1 and N; and i=1, …, iGCP,k. (XGCPik, YGCPik, ZGCPik) are 3D coordinates measured at the GCPs from DEMk. (XPiGCP,k, YPiGCP,k, ZPiGCP,k) are 3D coordinates at the same point from the control data set (e.g., ICESat-2). (εxiGCP,k, εyiGCP,k, εziGCP,k) represent residuals at GCPs in the observation equations, describing inconsistencies between the integrated DEM and the outside control data (e.g., ICESat-2). Furthermore, it is noticed that in the study area, elevation disagreements between adjacent REMA DEMs present a linear trend in the Y-direction (mainly ice flow direction). Thus, we use a linear term of aok+a1kYOLik to replace dZk in Eqs. (1) and (2). Here, aok and a1k are coefficients as unknowns, and YOLik is the Y coordinate defined in the local central coordinate system in the overlapping region, with the origin defined at the center of the overlapping area.
Equations (1) and (2) are combined and populated with all TPs and GCPs of the sub-DEMs. They are further formulated as an observation equation system in a matrix form. Finally, the unknowns, namely bias correction parameters (dXk, dYk, aok, and a1k) (k=1, …, N), are solved by an optimization procedure using the least-squares principle (McGlone et al., 2004). Uncertainties of the estimated bias correction parameters (σxik, σyik, σzik, σxiGCP,k, σyiGCP,k, σziGCP,k) are computed through an error propagation within the optimization procedure.
The Filchner Ice Shelf used as a validation and application site for the MDAM system. Colored boxes are REMA DEM coverages. Green lines at Berkner Island, Coats Land and on the ice shelf are ICESat-2 tracks. Blue circles are TPs and red triangles are GCPs. (a) REMA DEMs of 2016–2017, ZY-3 DEM of 2014 (light gray lines), and rifts T1 and T2. (b) REMA DEMs of 2020–2021. GCPs are also selected on the ice shelf where both the DEM images and ICESat-2 data were acquired within one day. Background elevation information is from the REMA DEM (Howat et al., 2019).
Implementation considerationsGround Control Point Selection
Selection of high quality GCPs is of great importance to the accuracy of the entire DEM time series. We choose GCPs on the grounded regions at both margins of the ice shelf (Figs. 1b and 3) with considerations of stability, low ice velocity, and local peaks (Feng et al., 2022). Meanwhile, the GCPs are also chosen to be on or close to the ICESat-2 ATL06 points from which their 3D coordinates can be measured. Despite the absence of outcrops and blue ice in the grounded regions of this study, we are able to identify relatively stable features for GCP selection, aided by multisource data, e.g., REMA DEMs, ICESat-2 data, ice flow velocity maps (Rignot et al., 2011), and bed topography (Morlighem et al., 2020).
Local 3D view for GCP candidate (red triangle) selection: shaded relief map (a) and velocity map (b) draped on the REMA DEM. 3D close-up view for GCP finalization (red cross): local maximum selected with the help of the ICESat-2 elevation profile (blue line) (c) and velocity map (d).
Figure 4 illustrates an example of GCP selection on the ice surface of Coats Land. Within a local area we drape a shaded relief map and a velocity map (Rignot et al., 2011) on the REMA DEM to make sure that a GCP candidate is selected not on a steep slope (<5°) and with a low ice flow speed (<10 m yr−1) (Fig. 4a and b). The candidate is approximately localized at a local peak (red triangle). In a close-up view (Fig. 4c and d), an ICESat-2 elevation profile (blue line) is further used to finalize the GCP position by adjusting the candidate position to the local maximum of the terrain (red cross). Finally, 3D coordinates of the GCP are measured on the ICESat-2 profile as knowns (XPiGCP,N, YPiGCP,N, ZPiGCP,N) and from the sub-DEM as observations (XGCPiN, YGCPiN, ZGCPiN) that are used in Eq. (2). The above GCP selection method is implemented as a manual procedure.
TP selection on rift lip: 3D views from two shaded relief maps of adjacent REMA DEMs (a) and (b), and verification using a 3D Landsat image view (c). TP selection on snow bridge: 3D views from two shaded relief maps of adjacent REMA DEMs (d) and (e); verification using a 3D Landsat image view (f).
Tie Point Selection
The key consideration for TP selection is to choose features that are identifiable on both adjacent sub-DEMs and unchanged during the time interval between the acquisitions of the sub-DEMs (Figs. 1b and 3). Further specific rules are followed: (1) if applicable, TPs be selected at distinct features, such as crevasses, rifts, and other relatively “stable” ice surface features; (2) elevated locations are preferred to avoid influences of shadows casted by different lighting conditions and obstacles in 3D visualization; and (3) TPs be approximately evenly distributed in the overlapping area. For example, to select a TP candidate on a rift lip, two shaded relief maps are produced from the adjacent sub-DEMs and used for 3D visualization (Fig. 5a and b) to check whether the feature (orange and blue dots in Fig. 5) is changed during the time interval between the two sub-DEMs. The TP is finalized (red cross) by a verification using a 3D Landsat image view with a different illumination geometry (Fig. 5c). Figure 5d–f show another example of TP selection on a snow bridge that is formed by an ice surface process when snow is blown into a crevasse and a rift or accumulated about tens of meters above the shelf surface (Williams et al., 2014; Thompson et al., 2020). Similarly, Fig. 5d and e show the 3D views from shaded relief maps with lower elevation angles produced from the adjacent sub-DEMs. After a verification using a Landsat image (Fig. 5f), the peak of the snow bridge is defined as the TP. In both cases, 3D coordinates (XPik, YPik, ZPik) and (XPik+1, YPik+1, ZPik+1) are measured and used in Eq. (1). The TP selection method is also implemented as a manual procedure.
Mélange Volume Estimation
To compute the mélange volumetric change, we apply the hydrostatic equilibrium principle to obtain the unknown thickness H of the mélange in each transect (Fricker et al., 2005):
H=ρwaterρwater-ρmélangeh-hsea
where h is the average mélange elevation along the transect and hsea is the average sea level (-15.2 m) estimated from the ICESat-2 altimetry data within a range of ∼ 10 km from the ice shelf front from March 2019 to December 2021 during which the ICESat-2 data become stable after the initial stage of the mission (Magruder et al., 2021). ρwater and ρmélange are densities of seawater and mélange (1028 and 865 kg m−3), respectively (King, 1994; Fricker et al., 2005; Adusumilli et al., 2020). The estimated change in mélange thickness H is relevant to ice sheet modeling. Consequently, the mélange volume of the transect is the product of the thickness H, rift width W, and transect interval, which is integrated over all transects to compute the mélange volume of the entire rift: Σ (H×W× transect interval).
Bias correction
The proposed MDAM system estimates the bias corrections (dXk, dYk, dZk) for all sub-DEM DEMk (k=1,2, …, N) through the least-squares procedure where their uncertainties are provided in the covariance matrix (McGlone et al., 2004). The REMA DEMs of 2016–2017 consists of 14 sub-DEMs that cover the entire front part of Filchner Ice Shelf (Fig. 3a). We use 12 GCPs that are measured by using ICESat-2 data to tie the connected DEMs to the grounded regions of Berkner Island and Coats Land where the average ice flow speed at the GCPs is ∼ 8 m yr−1. In total 51 TPs are used to connect the sub-DEMs on the floating part of the ice shelf, with 3–4 evenly distributed TPs in each overlapping area. Horizontal displacements at the TPs caused by the ice flow are on average ∼ 90 m and are corrected by using the velocity map (Gardner et al., 2019). We establish the observation equations in Eqs. (1) and (2) using different weights. The weights for TPs are computed as inverse distances from the TPs to the nearest GCPs on grounded regions. Hence, the weights for GCPs are set to 1. Those TPs that are farther away from the grounded regions have smaller weights.
We process the second set of REMA DEMs of 2020–2021 (Fig. 3b) in the same way in the least-squares process. In addition, the ZY-3 DEM of 2014 is reconstructed as a cross-ice-shelf DEM and does not need to go through this bias-correction process. The DEM is co-registered to the ICESat-2 ATL06 data of 2019 through a bundle adjustment procedure (McGlone et al., 2004; Li, 1998; Li et al., 2017b) using GCPs that are selected from “stable” features in the same way for those used in REMA sub-DEM co-registration.
Illustration of uncertainties of the estimated bias corrections (dXk, dYk, dZk) of the REMA DEMs of 2016–2017 in Table A2.
ResultsAdjusted DEM time series
The bias corrections (dXk, dYk, dZk) for all sub-DEM DEMk (k=1, 2, …, 14) in the REMA DEMs of 2016–2017 are listed in Table A2. As illustrated in Fig. 6, uncertainties of the corrections in both horizontal and elevation directions are effectively controlled for sub-DEMs close to grounded regions in Berkner Island and Coats Land through GCPs, leaving those in the middle section of the ice shelf relatively higher. The average uncertainty of the elevation corrections is 0.2 m, with the maximum of 0.4 m. The bias corrections (dXk, dYk, dZk) for the REMA DEMs of 2020–2021 are also listed in Table A2. The bias-corrected DEMs have an average uncertainty of 0.2 m for the elevation bias corrections. After a co-registration of the ZY-3 DEM of 2014 to the ICESat-2 ATL06 data of 2019, an elevation accuracy of 0.3 m is achieved.
Elevation offsets between adjacent sub-DEMs at TPs and GCPs in the REMA DEMs of 2016–2017 before bias-correction (black arrows, tide and ice velocity compensated) and thereafter (red arrows).
These corrections are then applied to each sub-DEM to produce the bias-corrected REMA DEMs of 2016–2017. The effectiveness of the MDAM system is clearly demonstrated by a dramatic reduction of the elevation offset between adjacent sub-DEMs at the 51 TPs and 12 GCPs (Fig. 7), in average, from 12.7 m before bias-correction (black arrows) to 0.2 m thereafter (red arrows, also called residuals). Residuals vary from 0.1 to 0.5 m.
Verification of bias-corrected REMA DEMs of 2020–2021: (a) data acquisition times and footprints of REMA DEMs in different colors, and times and two ground tracks of ICESat-2 data in green color; (b) and (c) elevation profiles of REMA DEMs of 2020–2021 prior to (gray) and post (orange) bias-correction, and ICESat-2 ATL06 data of 2021 (blue), with their corresponding differences at bottom.
To validate the bias-correction results using an external dataset, we directly compare the bias-corrected REMA DEMs of 2020–2021 with two tracks of ICESat-2 ATL06 data acquired from January to February 2021 (green lines in Fig. 8a). The ICESat-2 data set is chosen because of its high elevation accuracy and along-track point density. Furthermore, we select two tracks of the ICESat-2 data that were acquired on the same times (within one day) as the REMA sub-DEMs to avoid displacements caused by the shelf ice advection and internal mélange changes. The timespans of other two sets of DEMs, REMA DEMs 0f 2016–2017 and ZY-3 DEM, fall outside the ICESat-2 temporal coverage. Along Orbit 0375 Ground track GT3L, REMA DEM points prior to bias correction (gray) are relatively close to the ICESat-2 ATL06 elevation points (blue points), with a smaller average difference of 1.1 m (Fig. 8b). However, the average difference of 11.2 m along Orbit 0634 Ground track GT3L is significantly larger (Fig. 8c). After bias-corrections, the differences are effectively controlled within 0.1 m on these two tracks and 0.2 m on mélange top.
Multi-satellite DEM time series of two large rifts (T1 and T2) on Filchner Ice Shelf from 2014 to 2021 in 3D: (a) and (d) are reconstructed ZY-3 DEM of 2014, (b) and (e) are bias-corrected REMA DEMs of 2016–2017, and (c) and (f) are bias-corrected REMA DEMs of 2020–2021. The boxes in (b) and (e) indicate two rift sections where details of 3D structure and mélange are presented in Fig. 10.
3D mélange dynamics from 2014–2021
Based on the registered ZY-3 DEM of 2014 and bias-corrected sub-DEMs of REMA DEMs of 2016–2017 and 2020–2021, we establish a multi-satellite DEM time series of T1 and T2, from 2014 to 2021 (Fig. 9, Table A1). The 2D maps of the DEM time series are presented in Fig. A2. Each of these periodical unified and integrated DEMs presents the unprecedented 3D structural information of these ∼ 50 km long rifts. The unique DEM time series allows us to exploit 3D structure characteristics, including ice surface flank, rift lips, and mélange top and thickness (Fig. 10c and g). For the first time, we can accurately compute mélange elevation and volume changes within the entire rifts as they propagate and advect towards the shelf front during the period of 2014–2021 (Fig. 10d and h). To quantify 3D dynamic mélange changes within the rifts, we define a transect system for each rift based on the ZY-3 DEM of 2014 with transects spaced at an interval of 1 km along the rift centerline from the west tip to the east tip, which is also applied to the bias-corrected REMA DEMs of 2016–2017 and 2020–2021.
3D sectional structure and mélange: (a) and (e) are sectional shaded relief maps of REMA DEMs of 2016–2017 (4 January 2017 for T1 and 25 October 2016 for T2) indicated by the boxes in Fig. 9; (b) and (f) are WorldView images (0.5 m resolution) of 25 February 2016 for T1 and 2 March 2016 for T2; (c) and (g) are elevation profiles along the pink lines in (a) and (e), and 3D rift and mélange structure parameters; and (d) and (h) are rift and mélange changes along the profiles from 2014 to 2021 (2014 in blue, 2016–2017 in green, and 2020–2021 in red). Elevation displayed in (a) and (e) are exaggerated by 10 times.
The lengths of Rifts T1 and T2 are 47.1 and 48.3 km (maximum widths of 1.5 and 1.5 km), respectively, measured on the ZY-3 DEM of 2014. Based on the bias-corrected DEM time series from 2014 to 2021 (Fig. 9), we found that T1 and T2 propagated consistently during the period, resulting in an accelerated widening by 646 m (227 m yr−1) and 412 m (145 m yr−1), respectively, during 2014–2017; and 1109 m (268 m yr−1) and 834 m (203 m yr−1), respectively, during 2016–2021. Correspondingly, T1 and T2 have lengthened by 2393 m (840 m yr-1) and 1701 m (597 m yr−1), respectively, during 2014–2017; and 1082 m (262 m yr−1) and 1667 m (405 m yr−1), respectively, during 2016–2021. Based on the REMA DEMs of 2016–2017, the rift lips on the seaward sides are mostly higher than that on the landward side in the middle sections of the rifts (Table A3), with an average difference of 0.15 m. This phenomenon of the greater heights of seaward rift lips based on the precise measurements from the bias-corrected sub-DEMs is consistent with the results for some rifts on other Antarctic ice shelves presented in Walker and Gardner (2019).
Elevation (a for T1 and c for T2) and volume (b for T1 and d for T2) of ice mélange inside the rifts from the multi-satellite DEM series from 2014 to 2021. Average elevation and volume of each transect that are separated every 1 km along the rift centerline are illustrated from the west tip to the east tip.
Figure 11a and c illustrate the average mélange top elevations inside T1 and T2 from the ZY-3 DEM of 2014, REMA DEMs of 2016–2017 and REMA DEMs of 2020–2021. Mélange top elevation profiles in the central sections of the rifts are relatively flat where distinct linear features, called Seracs in Neuburg et al. (1959) and Alley et al. (2023), appear floating with a dominant orientation parallel to the rift centerline in both shaded relief map of DEM (Fig. 10a) and WorldView satellite image (Fig. 10b). We believe that the central sections of the rifts cut through the ice shelf and call the mélange “mature mélange”. In contrast, the mélange in the sections close to both rift tips has relatively rugged surface. More specifically, we found cavities that are freshly opened during the fracturing process in the unmatured mélange section close to the west tip of T1 (Fig. 11a). These pre-mélange cavities are verified in ZY-3 stereo images and reconstructed ZY-3 DEM (Fig. A3). They are initially as deep as 24.3 m in the ZY-3 DEM of 2014 (Fig. 11a). As the rift propagates in both length and width, new shelf ice in this section breaks apart from rift walls and fills in the cavities, resulting in a 5.7 m reduction in the cavity depth in the bias-corrected REMA DEM of 2016–2017 (Fig. 11a), and finally, disappearing during 2020–2021. During this process, we believe that the rift in this section cuts the ice shelf deeper and deeper, and the premature mélange becomes gradually mature.
The established DEM time series unveils an overall decrease of 1.8 ± 0.6 m, at a rate of 0.4 ± 0.1 m yr−1, in mélange elevation inside the two rifts from 2014 to 2021 (3.1 ± 0.4 and 0.4 ± 0.4 m inside T1 (excluding cavity) and T2, respectively) (between blue and red lines in Fig. 11a and c). Correspondingly, the mélange thickness decreased by 10.0 ± 3.3 m (1.6 ± 0.4 m yr−1) during the period (19.7 ± 2.5 and 2.8 ± 2.2 m for T1 and T2, respectively). This mélange thickness reduction observed during rift widening on Filchner Ice Shelf coincides with the earlier findings on Amery Ice shelf, Ronne Ice Shelf, and Larsen C Ice Shelf (Fricker et al., 2005; Walker et al., 2021; Larour et al., 2021). It may also indicate that the stress propagation between flanks may be compromised (Larour et al., 2004). More importantly, our results indicate that despite the significant decrease in mélange elevation from 2014 to 2021 (Fig. 11a and c), the inferred periodical mélange volumes reveal a rapid expansion of the mélange by (10.15 ± 0.03) × 109 m3 (151 %) during the period, (7.93 ± 0.03) × 109 m3 for T1 and (2.22 ± 0.01) × 109 m3 for T2, respectively (Fig. 11), mainly attributed to newly vacated space due to partial collapse of rift flanks, associated rift widening, and other rift-mélange interaction factors.
Discussion
The developed multi-temporal DEM adjustment model MDAM has shown its effectiveness in removing biases between adjacent multi-satellite sub-DEMs and establishing a unified and integrated DEM time series. We demonstrate the full 3D mapping capability for characterizing rift structure (rift lips, pre-mélange cavities, mélange seracs, etc.) and estimating dynamic mélange volume changes from the ICESat-2 controlled DEM time series, extending from rift topography and mélange thickness estimated along ground tracks from ICESat and ICESat-2 measurements in previous studies (Fricker et al., 2005; Walker et al., 2021). Such high-resolution mélange dynamic observations allow us to understand the mélange movement inside a closed environment of large Antarctic rifts, and to further study its role in rapid rift propagation and iceberg calving. To make this MDAM model working in a more dynamic open ocean environment, such as Pine Island and Thwaites ice shelves where ice velocity is significantly higher, modifications need to be carried out to address the rapid calving process with incoherent mélange changes in ice shelf front. For example, tie points between adjacent sub DEMs may be selected with decreased time intervals to minimize the uncertainty caused by the compensation of horizontal displacement due to ice dynamics using velocity maps.
A thick mélange layer can effectively “freeze” a rift, enabling mechanical stress transmission between its flanks and ultimately suppressing rift propagation (Rignot and MacAyeal, 1998; Larour et al., 2021). In this study, we show that the mélange thickness in T1 and T2 decreased as the rifts widened. Thus, the mélange layer may have thinned for a relatively long period, passing a possible stage of rift freeze. We further propose that the mélange production itself is related to the rift widening process. The seracs inside the rifts are formed from partial collapsing of rift flanks. We suggest that the bottom part of the rift wall is first excavated through interactions by tides (Padman et al., 2002, 2008), melting caused by intrusion of warmer sea water into the rift (Poinelli et al., 2023a), and other factors. This, in turn, causes the collapse of the upper part of the rift flank due to the removal of the bottom support. This process, like coastal bluff erosion, repeats itself and becomes one of the mechanisms that widen the rift and increase the mélange volume. We find that the increased mélange volume in a rift promotes its widening rate.
Although the icebergs calved from the shelf front in 1986 are mostly located inside the passive shelf ice area (Fig. A1), meaning no significant buttress reduction from the ice shelf (Doake et al., 1998; Fürst et al., 2016), the combined rifts of T1 and T2 propagated rapidly recently and have already covered ∼ 58 % of the ice shelf laterally from Coats Land to Berkner Island. They have the potential to cause a calving before reaching the passive shelf ice boundary (green line in Fig. A1). Furthermore, warm water observed near Berkner Island (Davis et al., 2022) may accelerate the propagation processes of T1 and T2 and destabilize the ice shelf through a mechanism like that proposed for the Larsen C Ice Shelf (Poinelli et al., 2023a, b).
Conclusions
We have developed an innovative multi-temporal DEM adjustment model (MDAM). This system successfully removed biases, as large as ∼ 6 m in elevation, between 50 sub-DEMs across the Filchner Ice Shelf in East Antarctica. The bias-corrected multi-satellite DEM time series achieved an accuracy of 0.1 m for the 2020–2021 REMA DEM assessed by ICESat-2 validation, and 0.2 m for the 2016–2017 REMA DEM and 0.3 m for the 2014 ZY-3 DEM estimated through the error propagation. Using this unified and integrated high-resolution DEM time series from 2014 to 2021, for the first time, we unveiled the 3D structural and mélange parameters of the ∼ 50 km long rifts of T1 and T2, including the rift lip height, pre-mélange cavity, and mélange elevation and thickness. Our findings indicate that while the mélange elevation decreased by 2.7 ± 0.6 m from 2014 to 2021, the mélange within the rifts experienced a rapid expansion by (10.31 ± 0.03) × 109 m3, or 139 %, during the period. This expansion is attributed to newly calved shelf ice from rift walls, associated rift widening, and other rift-mélange interaction factors. These new insights into 3D mélange dynamics are important for studying the mechanisms and the role of mélange in rift propagation.
Rifts of T1 and T2 (pink) in 2021, Grand Chasm (blue) in 1986; passive shelf ice (PSI) boundary (green) (Fürst et al., 2016); ice bergs calved in 1986; background is Landsat image of 11 February 2021.
Multi-satellite DEM time series of the two large rifts (T1 and T2) on Filchner Ice Shelf from 2014 to 2021 in 2D: (a) and (d) are reconstructed ZY-3 DEM of 2014, (b) and (e) are bias-corrected REMA DEMs of 2016–2017, and (c) and (f) are bias-corrected REMA DEMs of 2020–2021. The boxes in (b) and (e) indicates two rift sections of the rift where details of 3D structure and mélange are presented in Fig. 10.
Anaglyph stereo pair of 2014 ZY-3 satellite images (a) and reconstructed ZY-3 DEM (b) showing cavities close to the west tip of rift T1.
Information of images used for ZY-3 DEM reconstruction and REMA sub-DEMs (Howat et al., 2022).
Bias corrections (dXk, dYk, dZk) estimated for sub-DEMs of the REMA DEMs of 2016–2017 and the REMA DEMs of 2020–2021, ordered from Berkner Island to Coats Land.
Heights of rift lips measured on the land-ward (L) side and sea-ward (S) side, and their difference (S–L) at each transect along the centerlines of Rift T1 and T2 (REMA DEMs of 2016–2017).
Rift T1: Transect ID from west tip toLand-ward lipSea-ward lipDifferenceeast tip (1 km separation)height (L) (m)height (S) (m)(S–L) (m)0-1.70.52.213.64.51.022.74.72.030.81.70.944.41.5-2.950.71.71.061.82.00.272.32.60.381.31.40.290.11.71.710-1.61.83.411-0.71.11.7121.63.41.8132.83.00.2141.62.61.0151.62.30.8161.82.10.3171.51.90.5182.52.2-0.3191.61.90.2201.91.8-0.1212.31.6-0.7222.01.4-0.6231.70.9-0.7241.71.0-0.7250.91.00.1261.90.3-1.6271.00.1-0.9280.71.50.8292.43.41.0301.42.10.7313.33.70.5323.33.90.6333.72.1-1.6343.74.81.1352.01.90.0360.41.20.837-0.3-2.3-2.0382.83.50.7392.51.3-1.2402.41.7-0.6412.01.1-0.9420.81.50.7432.62.60.0444.64.2-0.4453.12.7-0.546-0.80.11.0471.81.90.148-1.4-1.8-0.4
Average lip height difference (s–l) within the section of mature mélange of T1: 0.1 m (max 1.8 m at 12 km from west tip).
Continued.
Rift T2: Transect ID from west tip toLand-ward lipSea-ward lipDifferenceeast tip (1 km separation)height (L) (m)height (S) (m)(S–L) (m)00.43.02.510.60.70.12-0.11.61.732.13.00.944.54.90.452.32.90.762.00.8-1.272.02.70.78-2.2-1.80.49-2.6-0.91.7100.10.50.511-0.70.10.8120.91.20.3130.3-0.9-1.314-0.3-3.2-2.915-0.3-2.4-2.2160.3-0.9-1.3170.21.20.9180.13.23.019-1.64.66.120-11.51.4–21-7.8-11.5-3.722-2.83.05.9230.62.72.2241.23.62.5250.7-16.7–26-18.6-5.24–27-0.3-1.9-1.6282.50.0-2.5291.81.6-0.1301.11.10.1311.21.1-0.1320.6-0.1-0.7331.91.5-0.4343.32.1-1.2352.92.80.0363.63.5-0.1374.74.4-0.3383.93.3-0.6394.33.5-0.8403.22.7-0.6413.02.0-1.0423.41.4-2.0433.61.7-1.9443.52.0-1.5453.51.7-1.9464.31.6-2.7475.31.7-3.6486.14.7-1.4495.06.21.2
Average lip height difference (s–l) within the section of mature mélange of T2: 0.2 m (max 5.9 m at 22 km from west tip).
Data availability
REMA DEMs used in this study are available at 10.7910/DVN/X7NDNY (Howat et al., 2019). ITS_LIVE ice flow velocity maps are at 10.5067/IMR9D3PEI28U (Gardner et al., 2019). The MDAM code (including TPs and GCPs) is available at https://github.com/menglianxia/MDAM (last access: 3 December 2025; Xia et al. 2025a). The produced DEMs (including adjusted REMA DEMs of 2016-2017 and 2020-2021, and ZY-3 DEM of 2014) are available at 10.5281/zenodo.15260323 (Xia et al. 2025b).
Author contributions
RL led the study and designed the MDAM system. MX programmed and implemented the model. ZL performed accuracy verification. RL, MX, and MS were involved in data analysis and presentation. RL, MS, LA, and GQ edited the manuscript.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Acknowledgements
We thank the editor, Professor Bert Wouters, for his guidance and suggestions during the review process. We sincerely thank Dr. Ann-Sofie Priergaard Zinck and Dr. Mattia Poinelli for their constructive comments and suggestions, which greatly improved the manuscript. We thank the National Snow and Ice Data Center (NSIDC) for the REMA DEMs, the United States Geological Survey (USGS) for the Landsat images, and the Chinese Ministry of Natural Resources for the ZY-3 images. The WorldView images are purchased from DigitalGlobe.
Financial support
This research has been supported by the National Natural Science Foundation of China (grant no. 42394131) and the Ministry of Science and Technology of the People's Republic of China, National Key Research and Development Program of China (grant no. 2024YFF0808302).
Review statement
This paper was edited by Bert Wouters and reviewed by Ann-Sofie Priergaard Zinck and Mattia Poinelli.
References Adusumilli, S., Fricker, H. A., Medley, B., Padman, L., and Siegfried, M. R.: Interannual variations in meltwater input to the Southern Ocean from Antarctic ice shelves, Nat. Geosci., 13, 616–620, 10.1038/s41561-020-0616-z, 2020.Alley, R. B., Cuffey, K. M., Bassis, J. N., Alley, K. E., Wang, S., Parizek, B. R., Anandakrishnan, S., Christianson, K., and DeConto, R. M.: Iceberg Calving: Regimes and Transitions, Annual Review of Earth and Planetary Sciences, 51, 189–215, 10.1146/annurev-earth-032320-110916, 2023.Allison, I., Alley, R. B., Fricker, H. A., Thomas, R. H., and Warner, R. C.: Ice sheet mass balance and sea level, Antarctic Sci., 21, 413–426, 10.1017/S0954102009990137, 2009.Anderson, N. T. and Marchisio, G. B.: WorldView-2 and the evolution of the DigitalGlobe remote sensing satellite constellation: introductory paper for the special session on WorldView-2, in: Algorithms and Technologies for Multispectral, Hyperspectral, and Ultraspectral Imagery XVIII, 2012/05/14, 166–180, 10.1117/12.919756, 2012.Bindschadler, R., Vornberger, P., Fleming, A., Fox, A., Mullins, J., Binnie, D., Paulsen, S. J., Granneman, B., and Gorodetzky, D.: The Landsat image mosaic of Antarctica, Remote Sens. Environ., 112, 4214–4226, 10.1016/j.rse.2008.07.006, 2008.Brunt, K. M., Neumann, T. A., and Smith, B. E.: Assessment of ICESat-2 Ice Sheet Surface Heights, Based on Comparisons Over the Interior of the Antarctic Ice Sheet, Geophys. Res. Lett., 46, 13072–13078, 10.1029/2019GL084886, 2019. Cassotto, R. K., Burton, J. C., Amundson, J. M., Fahnestock, M. A., and Truffer, M.: Granular decoherence precedes ice mélange failure and glacier calving at Jakobshavn Isbræ, Nat. Geosci., 14, 417–422, 2021.Davis, P. E. D., Jenkins, A., Nicholls, K. W., Dutrieux, P., Schröder, M., Janout, M. A., Hellmer, H. H., Templeton, R., and McPhail, S.: Observations of Modified Warm Deep Water Beneath Ronne Ice Shelf, Antarctica, From an Autonomous Underwater Vehicle, J. Geophys. Res. Oceans, 127, e2022JC019103, 10.1029/2022JC019103, 2022.Depoorter, M. A., Bamber, J. L., Griggs, J. A., Lenaerts, J. T. M., Ligtenberg, S. R. M., van den Broeke, M. R., and Moholdt, G.: Calving fluxes and basal melt rates of Antarctic ice shelves, Nature, 502, 89–92, 10.1038/nature12567, 2013.DiMarzio, J., Brenner, A., Schutz, R., Shuman, C. A., and Zwally, H. J.: GLAS/ICESat 500 m laser altimetry digital elevation model of Antarctica, National Snow and Ice Data Center, 10.5067/K2IMI0L24BRJ, 2007.Doake, C. S. M., Corr, H. F. J., Rott, H., Skvarca, P., and Young, N. W.: Breakup and conditions for stability of the northern Larsen Ice Shelf, Antarctica, Nature, 391, 778–780, 10.1038/35832, 1998.Feng, T., Li, Y., Wang, K., Qiao, G., Cheng, Y., Yuan, X., Luo, S., and Li, R.: A hierarchical network densification approach for reconstruction of historical ice velocity fields in East Antarctica, J. Glaciol., 1–20, 10.1017/jog.2022.58, 2022.Ferrigno, J. G. and Gould, W. G.: Substantial changes in the coastline of Antarctica revealed by satellite imagery, Polar Record, 23, 577–583, 10.1017/S003224740000807X, 1987.Francis, D., Mattingly, K. S., Lhermitte, S., Temimi, M., and Heil, P.: Atmospheric extremes caused high oceanward sea surface slope triggering the biggest calving event in more than 50 years at the Amery Ice Shelf, The Cryosphere, 15, 2147–2165, 10.5194/tc-15-2147-2021, 2021.Fricker, H. A., Young, N. W., Coleman, R., Bassis, J. N., and Minster, J.-B.: Multi-year monitoring of rift propagation on the Amery Ice Shelf, East Antarctica, Geophys. Res. Lett., 32, 10.1029/2004GL021036, 2005.Fürst, J. J., Durand, G., Gillet-Chaulet, F., Tavard, L., Rankl, M., Braun, M., and Gagliardini, O.: The safety band of Antarctic ice shelves, Nature Climate Change, 6, 479–482, 10.1038/nclimate2912, 2016.Gardner, A. S., Fahnestock, M. A., and Scambos, T. A.: MEaSUREs ITS_LIVE Landsat Image-Pair Glacier and Ice Sheet Surface Velocities, Version 1, Boulder, Colorado USA, NASA National Snow and Ice Data Center Distributed Active Archive Center, 10.5067/IMR9D3PEI28U, 2019.Greene, C. A., Gardner, A. S., Schlegel, N.-J., and Fraser, A. D.: Antarctic calving loss rivals ice-shelf thinning, Nature, 609, 948–953, 10.1038/s41586-022-05037-w, 2022.Grosfeld, K., Schröder, M., Fahrbach, E., Gerdes, R., and Mackensen, A.: How iceberg calving and grounding change the circulation and hydrography in the Filchner Ice Shelf-Ocean System, J. Geophys. Res. Oceans, 106, 9039–9055, 10.1029/2000JC000601, 2001.Howat, I. M., Porter, C., Smith, B. E., Noh, M.-J., and Morin, P.: The Reference Elevation Model of Antarctica, The Cryosphere, 13, 665–674, 10.5194/tc-13-665-2019, 2019.Howat, I., Porter, C., Noh, M.-J., Husby, E., Khuvis, S., Danish, E., Tomko, K., Gardiner, J., Negrete, A., Yadav, B., Klassen, J., Kelleher, C., Cloutier, M., Bakker, J., Enos, J., Arnold, G., Bauer, G., and Morin, P.: The Reference Elevation Model of Antarctica – Strips, Version 4.1 (V1), Harvard Dataverse, 10.7910/DVN/X7NDNY, 2022.Hulbe, C. L., LeDoux, C., and Cruikshank, K.: Propagation of long fractures in the Ronne Ice Shelf, Antarctica, investigated using a numerical model of fracture propagation, J. Glaciol., 56, 459–472, 10.3189/002214310792447743, 2010.Jacobs, S. S., Helmer, H. H., Doake, C. S. M., Jenkins, A., and Frolich, R. M.: Melting of ice shelves and the mass balance of Antarctica, J. Glaciol., 38, 375–387, 10.3189/S0022143000002252, 1992.Joughin, I. and MacAyeal, D. R.: Calving of large tabular icebergs from ice shelf rift systems, Geophys. Res. Lett., 32, 10.1029/2004GL020978, 2005.King, E. C.: Observations of a rift in the Ronne Ice Shelf, Antarctica, J. Glaciol., 40, 187–189, 10.3189/S0022143000003968, 1994.King, M. A., Padman, L., Nicholls, K., Clarke, P. J., Gudmundsson, G. H., Kulessa, B., and Shepherd, A.: Ocean tides in the Weddell Sea: New observations on the Filchner-Ronne and Larsen C ice shelves and model validation, J. Geophys. Res. Oceans, 116, 10.1029/2011JC006949, 2011.Larour, E., Rignot, E., and Aubry, D.: Modelling of rift propagation on Ronne Ice Shelf, Antarctica, and sensitivity to climate change, Geophys. Res. Lett., 31, 10.1029/2004GL020077, 2004.Larour, E., Rignot, E., Poinelli, M., and Scheuchl, B.: Physical processes controlling the rifting of Larsen C Ice Shelf, Antarctica, prior to the calving of iceberg A68, Proceedings of the National Academy of Sciences, 118, e2105080118, 10.1073/pnas.2105080118, 2021. Li, R.: Potential of high-resolution satellite imagery for national mapping products, Photogrammetric Engineering and Remote Sensing, 64, 1165–1170, 1998.Li, R., Xiao, H., Liu, S., and Tong, X.: A Systematic Study of the Fracturing of Ronne - Filchner Ice Shelf, Antarctica, Using Multisource Satellite Data from 2001 to 2016, The Cryosphere Discuss. [preprint], 10.5194/tc-2017-178, 2017a.Li, R., Ye, W., Qiao, G., Tong, X., Kong, F., Liu, S., and Ma X.: A new analytical method for estimating Antarctic ice flow in 1960s from historical optical satellite imagery, IEEE Transactions on Geoscience and Remote Sensing, 55, 2771–2785, 10.1109/TGRS.2017.2654484, 2017b.Li, R., Li, H., Hao, T., Qiao, G., Cui, H., He, Y., Hai, G., Xie, H., Cheng, Y., and Li, B.: Assessment of ICESat-2 ice surface elevations over the Chinese Antarctic Research Expedition (CHINARE) route, East Antarctica, based on coordinated multi-sensor observations, The Cryosphere, 15, 3083–3099, 10.5194/tc-15-3083-2021, 2021.Li, R., Cheng, Y., Cui, H., Xia, M., Yuan, X., Li, Z., Luo, S., and Qiao, G.: Overestimation and adjustment of Antarctic ice flow velocity fields reconstructed from historical satellite imagery, The Cryosphere, 16, 737–760, 10.5194/tc-16-737-2022, 2022.Li, R., Cheng, Y., Chang, T., Gwyther, D. E., Forbes, M., An, L., Xia, M., Yuan, X., Qiao, G., Tong, X., and Ye, W.: Satellite record reveals 1960s acceleration of Totten Ice Shelf in East Antarctica, Nat. Commun., 14, 4061, 10.1038/s41467-023-39588-x, 2023a.Li, T., Hu, Y., Liu, B., Jiang, L., Wang, H., and Shen, X.: Co-registration and residual correction of digital elevation models: a comparative study, The Cryosphere, 17, 5299–5316, 10.5194/tc-17-5299-2023, 2023b.Liu, H., Jezek, K., Li, B., and Zhao, Z.: Radarsat Antarctic Mapping Project digital elevation model version 2, Boulder, Colorado USA, NASA National Snow and Ice Data Center Distributed Active Archive Center, 10.5067/8JKNEW6BFRVD, 2001.Liu, Y., Moore, J. C., Cheng, X., Gladstone, R. M., Bassis, J. N., Liu, H., Wen, J., and Hui, F.: Ocean-driven thinning enhances iceberg calving and retreat of Antarctic ice shelves, Proceedings of the National Academy of Sciences, 112, 3263–3268, 10.1073/pnas.1415137112, 2015.Lv, D., Cheng, Y., Xiao, H., and Hai, G.: A Framework for Fracture Extraction Under Glaciological Property-Based Constraints: Scientific Application on the Filchner–Ronne Ice Shelf of Antarctica, Earth and Space Science, 9, e2022EA002293, 10.1029/2022EA002293, 2022.Magruder, L. A., Brunt, K. M., and Alonzo, M.: Early ICESat-2 on-orbit Geolocation Validation Using Ground-Based Corner Cube Retro-Reflectors, Remote Sensing, 12, 3653, 10.3390/rs12213653, 2020.Magruder, L., Neumann, T., and Kurtz, N.: ICESat-2 Early Mission Synopsis and Observatory Performance, Earth and Space Science, 8, e2020EA001555, 10.1029/2020EA001555, 2021.Markus, T., Neumann, T., Martino, A., Abdalati, W., Brunt, K., Csatho, B., Farrell, S., Fricker, H., Gardner, A., Harding, D., Jasinski, M., Kwok, R., Magruder, L., Lubin, D., Luthcke, S., Morison, J., Nelson, R., Neuenschwander, A., Palm, S., Popescu, S., Shum, C., Schutz, B. E., Smith, B., Yang, Y., and Zwally, J.: The Ice, Cloud, and land Elevation Satellite-2 (ICESat-2): Science requirements, concept, and implementation, Remote Sens. Environ., 190, 260–273, 10.1016/j.rse.2016.12.029, 2017. McGlone, J., Mikhail, E., Bethel, J., and Mullen, R.: Manual of photogrammetry, 5th ed., American Society for Photogrammetry and Remote Sensing, ISBN 978-1570830716, 2004.Morlighem, M., Rignot, E., Binder, T., Blankenship, D., Drews, R., Eagles, G., Eisen, O., Ferraccioli, F., Forsberg, R., and Fretwell, P.: Deep glacial troughs and stabilizing ridges unveiled beneath the margins of the Antarctic ice sheet, Nat. Geosci., 13, 132–137, 10.1038/s41561-019-0510-8, 2020.Neuburg, H. A. C., Thiel, E., Walker, P. T., Behrendt, J. C., and Aughenbaugh, N. B.: The Filchner Ice Shelf, Annals of the Association of American Geographers, 49, 110–119, 10.1111/j.1467-8306.1959.tb01603.x, 1959.Noh, M.-J. and Howat, I. M.: Automated stereo-photogrammetric DEM generation at high latitudes: Surface Extraction with TIN-based Search-space Minimization (SETSM) validation and demonstration over glaciated regions, GIScience & Remote Sensing, 52, 198–217, 10.1080/15481603.2015.1008621, 2015.Nuth, C. and Kääb, A.: Co-registration and bias corrections of satellite elevation data sets for quantifying glacier thickness change, The Cryosphere, 5, 271–290, 10.5194/tc-5-271-2011, 2011.Otosaka, I. N., Shepherd, A., Ivins, E. R., Schlegel, N.-J., Amory, C., van den Broeke, M. R., Horwath, M., Joughin, I., King, M. D., Krinner, G., Nowicki, S., Payne, A. J., Rignot, E., Scambos, T., Simon, K. M., Smith, B. E., Sørensen, L. S., Velicogna, I., Whitehouse, P. L., A, G., Agosta, C., Ahlstrøm, A. P., Blazquez, A., Colgan, W., Engdahl, M. E., Fettweis, X., Forsberg, R., Gallée, H., Gardner, A., Gilbert, L., Gourmelen, N., Groh, A., Gunter, B. C., Harig, C., Helm, V., Khan, S. A., Kittel, C., Konrad, H., Langen, P. L., Lecavalier, B. S., Liang, C.-C., Loomis, B. D., McMillan, M., Melini, D., Mernild, S. H., Mottram, R., Mouginot, J., Nilsson, J., Noël, B., Pattle, M. E., Peltier, W. R., Pie, N., Roca, M., Sasgen, I., Save, H. V., Seo, K.-W., Scheuchl, B., Schrama, E. J. O., Schröder, L., Simonsen, S. B., Slater, T., Spada, G., Sutterley, T. C., Vishwakarma, B. D., van Wessem, J. M., Wiese, D., van der Wal, W., and Wouters, B.: Mass balance of the Greenland and Antarctic ice sheets from 1992 to 2020, Earth Syst. Sci. Data, 15, 1597–1616, 10.5194/essd-15-1597-2023, 2023.Padman, L., Fricker, H. A., Coleman, R., Howard, S., and Erofeeva, L.: A new tide model for the Antarctic ice shelves and seas, Ann. Glaciol., 34, 247–254, 10.3189/172756402781817752, 2002.Padman, L., Erofeeva, S. Y., and Fricker, H. A.: Improving Antarctic tide models by assimilation of ICESat laser altimetry over ice shelves, Geophys. Res. Lett., 35, 10.1029/2008GL035592, 2008.Poinelli, M., Schodlok, M., Larour, E., Vizcaino, M., and Riva, R.: Can rifts alter ocean dynamics beneath ice shelves?, The Cryosphere, 17, 2261–2283, 10.5194/tc-17-2261-2023, 2023a.Poinelli, M., Nakayama, Y., Larour, E., Vizcaino, M., and Riva, R.: Ice-Front Retreat Controls on Ocean Dynamics Under Larsen C Ice Shelf, Antarctica, Geophys. Res. Lett., 50, e2023GL104588, 10.1029/2023GL104588, 2023b.Rignot, E. and MacAyeal, D. R.: Ice-shelf dynamics near the front of the Filchner–Ronne Ice Shelf, Antarctica, revealed by SAR interferometry, J. Glaciol., 44, 405–418, 10.3189/S0022143000002732, 1998.Rignot, E., Mouginot, J., and Scheuchl, B.: Ice flow of the Antarctic ice sheet, Science, 333, 1427–1430, 10.1126/science.1208336, 2011.Rignot, E., Mouginot, J., Scheuchl, B., Van Den Broeke, M., Van Wessem, M. J., and Morlighem, M.: Four decades of Antarctic Ice Sheet mass balance from 1979–2017, Proceedings of the National Academy of Sciences, 116, 1095–1103, 10.1073/pnas.1812883116, 2019.Shean, D. E., Joughin, I. R., Dutrieux, P., Smith, B. E., and Berthier, E.: Ice shelf basal melt rates from a high-resolution digital elevation model (DEM) record for Pine Island Glacier, Antarctica, The Cryosphere, 13, 2633–2656, 10.5194/tc-13-2633-2019, 2019.Shen, X., Ke, C.-Q., Fan, Y., and Drolma, L.: A new digital elevation model (DEM) dataset of the entire Antarctic continent derived from ICESat-2, Earth Syst. Sci. Data, 14, 3075–3089, 10.5194/essd-14-3075-2022, 2022.Slater, T., Shepherd, A., McMillan, M., Muir, A., Gilbert, L., Hogg, A. E., Konrad, H., and Parrinello, T.: A new digital elevation model of Antarctica derived from CryoSat-2 altimetry, The Cryosphere, 12, 1551–1562, 10.5194/tc-12-1551-2018, 2018.Smith, B., Fricker, H. A., Gardner, A. S., Medley, B., Nilsson, J., Paolo, F. S., Holschuh, N., Adusumilli, S., Brunt, K., Csatho, B., Harbeck, K., Markus, T., Neumann, T., Siegfried, M. R., and Zwally, H. J.: Pervasive ice sheet mass loss reflects competing ocean and atmosphere processes, Science, 368, 1239–1242, 10.1126/science.aaz5845, 2020.Tachikawa, T., Hato, M., Kaku, M., and Iwasaki, A.: Characteristics of ASTER GDEM version 2, 2011 IEEE international geoscience and remote sensing symposium, Vancouver, BC, Canada, 3657–3660, 10.1109/IGARSS.2011.6050017, 2011.Thompson, S. S., Cook, S., Kulessa, B., Winberry, J. P., Fraser, A. D., and Galton-Fenzi, B. K.: Comparing satellite and helicopter-based methods for observing crevasses, application in East Antarctica, Cold Regions Science and Technology, 178, 103128, 10.1016/j.coldregions.2020.103128, 2020.Walker, C. and Gardner, A.: Evolution of ice shelf rifts: Implications for formation mechanics and morphological controls, Earth Planet. Sci. Lett., 526, 115764, 10.1016/j.epsl.2019.115764, 2019.Walker, C. C., Becker, M. K., and Fricker, H. A.: A High Resolution, Three-Dimensional View of the D-28 Calving Event From Amery Ice Shelf With ICESat-2 and Satellite Imagery, Geophys. Res. Lett., 48, e2020GL091200, 10.1029/2020GL091200, 2021.Wang, S., Liu, H., Jezek, K., Alley, R. B., Wang, L., Alexander, P., and Huang, Y.: Controls on Larsen C Ice Shelf Retreat From a 60-Year Satellite Data Record, J. Geophys. Res.-Earth Surf., 127, e2021JF006346, 10.1029/2021JF006346, 2022.Wang, T., Zhang, G., Li, D., Tang, X., Jiang, Y., Pan, H., Zhu, X., and Fang, C.: Geometric Accuracy Validation for ZY-3 Satellite Imagery, IEEE Geosci. Remote Sens. Lett., 11, 1168–1171, 10.1109/LGRS.2013.2288918, 2014.Wessel, B., Huber, M., Wohlfart, C., Bertram, A., Osterkamp, N., Marschalk, U., Gruber, A., Reuß, F., Abdullahi, S., Georg, I., and Roth, A.: TanDEM-X PolarDEM 90 m of Antarctica: generation and error characterization, The Cryosphere, 15, 5241–5260, 10.5194/tc-15-5241-2021, 2021.Williams, R. M., Ray, L. E., Lever, J. H., and Burzynski, A. M.: Crevasse Detection in Ice Sheets Using Ground Penetrating Radar and Machine Learning, IEEE J. Sel. Topics Appl. Earth Observ. Remote Sens., 7, 4836–4848, 10.1109/JSTARS.2014.2332872, 2014. Xia, M., Li, R., Scaioni, M., An, L., Li, Z. and Qiao, G.: MDAM, GitHub [code], https://github.com/menglianxia/MDAM (last access: 2 December 2025), 2025a.Xia, M., Li, R., Scaioni, M., An, L., Li, Z. and Qiao, G.: Dataset belonging to the article ”Building multi-satellite DEM time series for insight into mélange inside large rifts in Antarctica”, Zenodo [data set], 10.5281/zenodo.15260323, 2025b.Zinck, A.-S. P., Wouters, B., Lambert, E., and Lhermitte, S.: Unveiling spatial variability within the Dotson Melt Channel through high-resolution basal melt rates from the Reference Elevation Model of Antarctica, The Cryosphere, 17, 3785–3801, 10.5194/tc-17-3785-2023, 2023.