Interannual Variation of Landfast Ice Using Ascending and Descending Sentinel-1 Images from 2019 to 2021: A Case Study of Cambridge Bay

: Landfast ice has undergone a dramatic decline in recent decades, imposing potential effects on ice travel for coastal populations, habitats for marine biota, and ice use for industries. The mapping of landfast ice deformation and the investigation of corresponding causes of changes are urgent tasks that can provide substantial data to support the maintenance of the stability of the Arctic ecosystem and the development of human activities on ice. This work aims to investigate the time-series deformation characteristics of landfast ice at multi-year scales and the corresponding inﬂuence factors. For the landfast ice deformation monitoring technique, we ﬁrst combined the small baseline subset approach with ascending and descending Sentinel-1 images to obtain the line-of-sight deformations for two ﬂight directions, and then we derived the 2D deformation ﬁelds comprising the vertical and horizontal directions for the corresponding periods by introducing a transform model. The vertical deformation results were mostly within the interval [ − 65, 23] cm, while the horizontal displacement was largely within the range of [ − 26, 78] cm. Moreover, the magnitude of deformation observed in 2019 was evidently greater than those in 2020 and 2021. In accordance with the available data, we speculate that the westerly wind and eastward-ﬂowing ocean currents are the dominant reasons for the variation in the horizontal direction in Cambridge Bay, while the factors causing spatial differences in the vertical direction are the sea-level tilt and ice growth. For the interannual variation, the leading cause is the difference in sea-level tilt. These results can assist in predicting the future deformation of landfast ice and provide a reference for on-ice activities.


Introduction
As an integral and sensitive component of the global climate system, sea ice can regulate the exchange of fluxes and radiative transfer between air and sea; it can also alter sea-surface temperature and salinity, contributing significantly to climate and environmental change in polar regions [1,2]. Landfast ice is a type of sea ice that is attached to a shore, to the front of an ice shelf, or to the perimeter of a grounded iceberg, and it is a significant component of the local physical and ecological systems [3]. Compared with floating sea ice, landfast ice plays a more crucial role in the assessment of a climate system given its characteristics of a long ice age, greater thickness, and notable seasonal trends [3,4]. Arctic landfast ice forms a rigid interface between the sea and air in nearshore areas, effectively buffering coastal erosion due to waves and, thus, influencing the nearshore sedimentation and shoreline development [5]. Landfast ice also serves as a key habitat for marine biota and as a platform for coastal populations to live on [6]. Over the last century, the expansion of transportation and resource exploitation has triggered increased human activities and further diversification of ice use in the Arctic [7]. As global temperatures have continued to rise since the 1970s, climatic conditions in high-altitude areas of the Arctic have dramatically changed and sea ice has been rapidly declining [8], exposing ice users to unknown risks, as well as potential hazards [9]. The continued warming of the Arctic and the accelerated melting of sea ice-particularly the remarkable reduction in the extent of landfast ice in coastal areas-have led to increased ship traffic and resource exploration [10]. Therefore, determining landfast ice deformation is particularly important for local populations and organisms to anticipate potential hazard events and provide effective early warning in real time. Awareness is growing that sea ice conditions for Arctic marine operations will be a challenge that will require considerable monitoring and advanced methods in the future. This improvement must be implemented on a large scale and a regional scale to assess environmental hazards and develop effective and timely preventive measures.
The accurate determination of the large-scale and long-term deformation of landfast ice while ensuring the immediate validity and feasibility of a method has been a challenge [11]. Observing sea ice conditions in situ is difficult due to the specificity and complexity of the terrain. Satellite remote sensing is an essential tool for detecting sea ice parameters. Optical remote sensing is always effective in investigating sea ice properties, such as the extent [12] and type classification [13], with the characteristics of large spatial coverage, high spatio-temporal resolution, and a huge accessible database. However, optical methods are ineffective in making timely measurements due to weather limitations, such as clouds and fog. By the early 1990s, synthetic aperture radar (SAR) was introduced to create a superior dataset for sea ice studies, given its advantages of performing in all weather, working on all days, higher resolution, and wider monitoring coverage [5]. Intensity information from SAR has been widely applied for sea ice classification [14] and the analysis of motion [15], thickness [16], concentration [17], roughness [18] and snow-cover extent [19][20][21]. However, SAR backscatter does not typically provide deformation information about landfast ice or transiently stable floating ice, as the internal motion of landfast ice is too small to be identified via change detection. The advent of the SAR interferometry (InSAR) technology has opened up the possibility of using SAR data to monitor landfast ice deformation by extracting the phase difference between two SAR images with similar imaging geometries. In this manner, topography can be retrieved if acquisitions are separated in space (nonzero perpendicular baseline), or surface movement along the line-of-sight (LOS) direction can be measured if it is separated in time (nonzero temporal baseline) [22,23]. For example, Li et al. first introduced InSAR into the detection of sea ice deformation [24]. Wang et al. combined ascending-and descending-pass images to obtain the vertical and horizontal deformations of landfast ice for a single period in the Baltic Sea [25]. The preceding studies only retrieved the LOS or a 2D deformation field for a single period. This information is evidently insufficient for characterizing the deformation of landfast ice. With the launches of the Sentinel-1A satellite in 2014 and the Sentinel-1B satellite in 2016, ascending and descending datasets that cover the same area can be utilized to derive a 2D deformation field of landfast ice over a long period. Although ascending and descending SAR data cannot be combined for interferometry, they can provide InSAR observations with different imaging geometries, and 3D surface deformation can be estimated by building a multi-source heterogeneous InSAR observation model [26,27].
The coherence of ice as a monitored object is highly susceptible to external factors, such as temperature, precipitation, and wind. The small baseline subset (SBAS) InSAR technique ensures the high coherence and reliability of the inversion results by selecting interferometric combinations with small spatio-temporal baselines that are as far as possible [26,28]. SBAS-InSAR is a specific type of multi-temporal InSAR that utilizes a small subset of images to create highly accurate results [29]. This method is ideal for areas where the availability of satellite data is limited or for monitoring large areas that would otherwise be impractical with traditional InSAR. One of the key benefits of SBAS-InSAR is its ability to detect even small changes in the Earth's surface over time. This makes it a valuable tool for monitoring the effects of natural disasters, such as earthquakes [30] and volcanic eruptions [31], as well as human-made changes, such as urban development and subsi-dence [32]. Additionally, SBAS-InSAR can be used to monitor ground deformation and subsidence in areas where oil and gas extraction is taking place [33], which helps to ensure the safety of nearby communities and infrastructure. Another advantage of SBAS-InSAR is that it is a non-invasive technique that does not require any physical access to the area being monitored. This makes it ideal for monitoring remote or difficult-to-reach areas, such as mountainous regions. Based on this, it has a great potential to be applied in monitoring natural disasters of the cryosphere, such as ice avalanches [34] and sea ice motion [35].
The current work uses ascending and descending Sentinel-1 images and the SBAS approach to derive LOS deformation results for both modes from January to April in 2019-2021. Then, the vertical and horizontal time-series deformations can be decomposed on the basis of a transform model. Finally, the temporal and spatial deformation characteristics of the test area and the factors that affect the deformation changes are analyzed in relation to the experimental results.

SBAS-InSAR Data Processing
SBAS-InSAR is a powerful remote sensing technique that is used to monitor changes in the Earth's surface over time [36]. This technology combines the use of SAR images with the principle of interferometry to produce highly accurate and detailed maps of the Earth's surface [37]. The procedure of SBAS-InSAR involves several steps, as shown in the pink frame of Figure 1.
(1) Data acquisition: The first step in SBAS-InSAR processing is to acquire the necessary SAR images. This typically involves using a SAR satellite to acquire multiple images of the same area at different times. The images should be acquired from similar viewing angles to ensure the highest accuracy.
(2) Image preprocessing: Once the SAR images have been acquired, they must be preprocessed to remove any noise or errors in the data. This may include correction for atmospheric effects, radiometric calibration, and speckle filtering [38].
(3) SBAS network: The next step is to select a small subset of the SAR images to use for the SBAS-InSAR analysis. This subset should be chosen based on the spatial and temporal distribution of the images and the desired accuracy of the results.
(4) Interferogram stacks: After the small baseline subset has been selected, an interferogram can be generated by applying interferometric techniques to the subset of images.
The interferogram provides information about the relative height differences between the surfaces in the two images and is used to generate the final SBAS-InSAR results.
(5) Phase unwrapping: The next step is to unwrap the phase information in the interferogram to remove any phase ambiguities. This is typically done by using algorithms such as the least-squares method or the multichannel algorithm [39].
(6) Time-series deformation inversion: A singular value decomposition (SVD) is utilized here to estimate the minimal norm least-squares solution of the average deformation rate of each period [40,41]; meanwhile, the times-series deformation can be calculated as well [42].

Calculation of Horizontal and Vertical Deformations
The time-series deformation results of ascending and descending orbits were derived via the SBAS method. In this work, the time-series deformation results were used in combination with the classical function model to monitor the 3D deformation of landfast ice in the study area.
The classical function model was established on the basis of the geometric relationship between the InSAR observations and 3D deformations. Figure 2 shows a schematic diagram of 3D deformations from different perspectives. In accordance with Figure 2a,b, the deformation vector along the LOS direction is located on a 2D plane composed of a ground direction (projection of the LOS direction on the horizontal plane) and a vertical direction.
The ground range and vertical direction are perpendicular to each other; thus, the following expression can be obtained according to the geometric relationship in Figure 2a [43]: where d u and d h define, respectively, the vertical and ground displacement vectors; θ A inc represents the pixel-based incidence angle of the radar pulse signal. The ground displacement vector is composed of the eastern and northern components, as illustrated in Figure 2b [30,44].
where d e and d n represent, respectively, the east-west and north-south displacement vectors; α A AZI is the orbit azimuth angle (positive clockwise from the north), and α A AZI − 3π/2 actually defines the angle between the ground range direction and the north. By substituting Equation (2) into Equation (1), the relationship between the LOS deformation and the 3D surface deformation vector can be established in the case of the ascending mode.
Similarly, the relationship between the LOS deformation and 3D surface deformation can be constructed in the descending mode as follows [45]: Sentinel-1 is a near-polar orbit satellite with a low sensitivity in the north-south direction. Therefore, performing high-precision detection of movement in this direction is difficult. We assume that the north-south deformation of landfast ice in the study area is zero. Then, 2D deformations comprising the east-west and vertical deformations can be decomposed from the combination of ascending and descending results of the same epoch. This issue can be resolved by combining Equations (3) and (4).

Cambridge Bay
The study site was located in the southern part of Cambridge Bay, as shown in Figure 3. Cambridge Bay is a hamlet located on Victoria Island in the Kitikmeot Region in Nunavut, Canada. Given its unique geographical situation, i.e., it is located on the western end of Queen Maud Gulf, where it narrows into Dease Strait, Cambridge Bay is a transportation and administrative hub for the Kitikmeot Region and the largest stop for passenger and polar research vessels traversing the Arctic Ocean's Northwest Passage [46]. Cambridge Bay has a polar climate, that is, no month has an average temperature of 10 • C or higher. In accordance with the regional sea ice charts available from the Canadian Ice Service Data Archive (CISDA), landfast ice dominates the ice type, with various thickness values at different dates during the winter [47].

Sentinel-1 Images and Reference Data
The Sentinel-1 mission comprises a constellation of two polar-orbiting satellites (Sentinel-1A and Sentinel-1B, which share the same orbital plane) that operate day and night to perform C-band SAR imaging, thus enabling them to offer reliable and repeated wide-area monitoring regardless of the weather [48]. To investigate the interannual deformation characteristics of landfast ice in the area of interest, the three-year (2019-2021) Sentinel-1 single-look complex (SLC) products over Cambridge Bay were used to retrieve the 2D deformation results (the specific characteristics are listed in Table 1). Considering the coherence level of the observed objects, SAR images from January to April in each year were selected for the experiments to avoid the incoherence caused by ice melting. The interferometric processing was accomplished with the GAMMA software, which included interferometry and phase unwrapping. In the process of SBAS-InSAR processing, the master image of each subset was connected to the two slave images whose acquisition dates were the closest to that of the master ( Figure A1) to minimize the probability of temporal decoherence and obtain ideal interferometric and unwrapping results. The coherence threshold of 0.3 was adopted to ensure relatively reliable high-quality point selection. The multi-look ratio was set to 10:2 in this work. In addition, a Goldstein filter was used to suppress the noise phase, and the commonly used minimum cost flow algorithm was applied for phase unwrapping. The reference point R of phase unwrapping was located in the residential area of Cambridge Bay, which was regarded as stable ( Figure 3).
The temperature, tide, and wind of two weather stations (Austin Bay and Cambridge Bay, the locations of which are shown in Figure 3) were collected from the official website of the Government of Canada to understand the ice conditions and interpret the deformation characteristics in Cambridge Bay. Moreover, the ice draft data available from the Ocean Networks Canada Data Archive were downloaded and are displayed in Figure 4.

Average Coherence and Interferograms
Coherence in InSAR is defined as the cross-correlation coefficient in the registration of two SLCs [49]. It is a widely used byproduct of InSAR [50,51]. It is also a significant indicator that is generally used to evaluate the quality of the interferometric phase. A higher coherence can effectively alleviate the problem of phase discontinuity, which introduces errors in the follow-up steps, yielding a more reliable displacement result. In the current work, the average coherence from C-band Sentinel-1 InSAR pairs was counted and is displayed in Figure A2 for the period of January-April 2019-2021. The majority of average coherence values were within the range of [0.3, 0.4]. The coherence level of landfast ice in the present work was evidently lower than those of traditional research objects, such as roads, bridges, mining areas, and earthquakes [52,53]. The factors for decorrelation could mostly be attributed to the physical properties of the ice, which was susceptible to temperature, wind, and precipitation. In addition, the ice growth exhibited the potential to decrease coherence by increasing freeboard, and the internal random motion (e.g., the dynamics of brine inclusions) of ice and snow could also lower it [22].
From Equation (5), we acquired two interferogram subsets (ascending and descending) for each year during the period of 2019-2021. Given the form of the fringes, they were unlikely to be caused solely by the displacement on the surface snow layer or ice growth. Instead, the inducer was likely to be a mechanical metamorphosis. We chose a representative selection of filtered interferograms from each subset for presentation ( Figure A3). Each color-cycle change indicates that the deformation change along the LOS direction is half of a wavelength in length. Dense interferometric fringes correspond to a large phase gradient. Rather, a sparse fringe distribution indicates that the region remained relatively stable. As a general rule, horizontal movement dominated the surface deformation when the fringes exhibit an evident parallel trend, while circular fringes indicate that vertical deformation played a major role [54]. From the spatial distribution shown in Figure A3, sparse suborbicular fringes were distributed near Cambridge Bay, while the quantity of fringes gradually increased with increasing distance from land, and the fringes were arranged in parallel.
This finding suggests that landfast ice closer to the shore was more stable under the effect of a buttressing force. In terms of time-series variation, high-density interferometric fringes were spread across the study area in January and February. With the passing of time, the quantity and shape of fringes changed markedly from March to April, with mostly sparse circular fringes, indicating that the deformation in this area was relatively significant from January to February, and the landfast ice gradually became stable from March to April, with enhanced interference immunity. Some fringe discontinuities could be observed in these interferograms; these were largely due to the influence of ice ridges on the ice roads caused by anthropogenic activities during winter. In addition, four feature points ( P 1 , P 2 , P 3 , and P 4 ) were marked for a better exploration of the deformation evolution pattern in this region. The concrete positions of these points are shown in Figure 3. P 1 is located on land, with the primary use of checking the reliability of the experimental results. P 2 is close to Cambridge Bay, where interferometric fringes are sparse and the magnitude of movement is relatively small. P 3 is located near the ice ridge. It can be used to investigate the deformation characteristics of landfast ice with support from external forces. P 4 is located on smooth landfast ice without the interference of visible factors. Numerous dense fringes are distributed near P 4 , indicating the existence of a large deformation.

Ascending and Descending LOS Deformations
As described in Section 2.1, the SBAS method was used to derive the time-series deformation from 2019 to 2011 on the basis of ascending and descending Sentinel-1 images. All InSAR measurements were calibrated as changes relative to a reference point, which is labeled with R in Figure 3. The corresponding LOS deformation results from the ascending and descending datasets are displayed in Figures A4-A9. A positive or negative sign represents movement closer to or away from the SAR sensor. In terms of the overall spatial distribution, the results for the ascending and descending modes reflected nearly opposite deformation trends, with a maximum accumulated movement of −52 cm for ascending images and 55 cm for descending images. The different geometries and polarizations (HH was used in the ascending mode, while VV was used in the descending mode) adopted in the deformation monitoring contributed to the discrepancies in the results of the ascending and descending images. The same deformation may exhibit the opposite direction to that of a deformation monitored under the opposite geometry. The time-series deformation results of the four feature points located in the different areas mentioned in Section 4.1 were extracted to further reveal the generated deformation characteristics, as shown in Figure 5. P 1 and P 2 remained relatively stable in terms of temporal variation, which was closely related to their geographical locations. P 1 was located on land, and point P 2 was close to Cambridge Bay. By contrast, the remaining two feature points had significantly larger deformation magnitudes and exhibited opposite deformation directions for the different datasets. P 3 experienced a downward displacement along the LOS direction at a total of −4.7 cm, −7 cm, and −9.6 cm from 2019 to 2021 with ascending images, and it presented the opposite trend at 21 cm, 15 cm, and 20 cm from 2019 to 2021 under the descending geometry. The deformation trend of P 4 was similar to that of P 3 , and they only differed in magnitude.

Horizontal and Vertical Deformations
The temporal coverage of the Sentinel-1 ascending images was not completely consistent with that of the descending images. The 2D time-series deformation results could only cover overlapping periods. The horizontal time-series displacement maps from 2019 to 2021 are displayed in Figures 6-8, and the vertical deformation results are presented in Figures 9-11. A positive or negative sign indicates that the movement was uplifted (or moved to the east) or sunk (or moved to the west). A vertical deformation within the range of [−68, 23] cm and a horizontal deformation within the range of [−26, 78] cm were derived in 2019. Meanwhile, the magnitude of vertical deformation in 2020 and 2021 was apparently lower, i.e., within the range of [− 45,12] cm. In summary, the movement of the coastal zone tended towards zero in both directions, and it was evidently smaller than that of the channel. The results showed a predominant subsidence on the western side, with an overall tendency for landfast ice to migrate eastwards. In terms of the time sequence, the period with a large deformation rate occurred mostly from early January to early March. Thereafter, the deformation rate in the study area significantly decreased. As for the horizontal deformation, the landfast ice retained similar variation characteristics from 2019 to 2021, within the range of [−26, 78] cm. Most areas showed a trend of eastward movement, especially in the western flank.
The time-series deformation results of the feature points mentioned above were extracted and are displayed in the first and second columns of Figure 5 and Table 2. In 2019, P 1 , P 2 , P 3 , and P 4 experienced settlement at 0.1 cm, 0.6 cm, 12.8 cm, and 34.9 cm, respectively. P 1 , P 2 , P 3 , and P 4 experienced eastward movement at 0.1 cm, 0.5 cm, 15.8 cm, and 35.3 cm, respectively. Similarly, the magnitude of vertical deformation in 2019 was greater than those in the other two years. In 2020, the cumulative maximum deformations of P 1 , P 2 , P 3 , and P 4 were 0.1 cm, −0.3 cm, −4.1 cm, and −14.8 cm, respectively, in the vertical direction and 0.1 cm, 0.4 cm, 17.4 cm, and 39.1 cm, respectively, in the horizontal direction. In 2021, the cumulative maximum deformations of P 1 , P 2 , P 3 , and P 4 were 0 cm, −0.3 cm, −7.6 cm, and −11.9 cm, respectively, in the vertical direction and 0.03 cm, 0.4 cm, 24.7 cm, and 37.5 cm, respectively, in the horizontal direction. The factors that affected the deformation characteristics will be discussed in Section 5. In addition, SAR images obtained at different times were slightly different under different SAR conditions, such as the atmospheric conditions and polarization [55]. This phenomenon led to changes in intensity, coherence, and measurement in the InSAR analysis, and these are not evaluated in the current work.

Validation
To compensate for the unavailability of external validation data over Cambridge Bay, the residual phase of each interferogram was calculated and selected as the accuracy index. In addition, we also compared the results of this work with those of similar published studies to verify the reliability of the experimental results.
The residual phase can be estimated by subtracting the low-pass deformation phase from the unwrapped phase, and this has the ability to evaluate the accuracy of deformation modeling in InSAR analysis [56]. A smaller residual phase indicates a higher accuracy of deformation inversion. The average residual phase for each subset of SBAS processing is provided in Table 3. Its value was less than 1 rad in any flight direction and acquisition year, demonstrating that the conventional linear model used in the SBAS method was sufficiently accurate for deformation estimation in the current work.  [57]. Only a single period of LOS deformation was obtained here due to the paucity of data, thus failing to recover the true 3D deformation field. Wang et al. combined ascending and descending data to retrieve the vertical and horizontal deformations of landfast ice in the Baltic Sea [25]. They extended the dimensionality of InSAR technology to landfast ice deformation monitoring. For long-term monitoring, Chen et al. utilized the multi-dimensional SBAS method to process time-series observations for retrieving 2D deformations that comprised the vertical and horizontal directions in Cambridge Bay [58]. At present, only Chen et al. has applied multi-temporal InSAR to the monitoring of landfast ice deformation in Cambridge Bay. However, they only investigated deformation characteristics for 1 year and did not systematically address the potential causes of deformation. In the current work, ascending and descending datasets were used to detect the time-series 2D deformation for 3 years (2019-2021), and the possible causes of deformation are summarized and discussed in a later section. Given the lack of external in situ measurements, we compared the results retrieved for the same geographical location with those of previous studies to verify the accuracy and improvements of the time-series 2D deformation in this work. The deformation results for the same acquisition period (2019) are listed in Table 4. A small deformation difference indicates that the two results maintained a high level of agreement, thus proving the reliability and authenticity of the experimental results in the current work.

Suggested Reasons for Deformation
In this work, 2D deformation results that included the east-west and vertical displacements of landfast ice in Cambridge Bay from January to April in 2019-2021 were obtained by combining ascending and descending images. The suggested reasons for the deformations are described from two perspectives. The factors that affect horizontal deformations include ocean currents, winds, and thermal expansion [25]. Meanwhile, the factors that affect vertical deformation include temperature, ice growth, and sea-level tilt [58], which are further described below.
The horizontal deformation of the landfast ice presented an eastward migration trend, except for landfast ice around the coastal area, which remained relatively stable. The cold ocean current from Beaufort Sea flows through Amundsen Gulf, Coronation Gulf, and Cambridge Bay from west to east [59], indicating the west-to-east current direction in this area. Thus, the cold ocean current possesses the ability to promote the eastward movement of landfast ice near Cambridge Bay in the horizontal direction. For the winds, their direction and magnitude exerted a certain degree of influence on the horizontal deformation. Figure 12 shows the wind information of the corresponding research period (2019-2021), including the speed and direction. Overall, the study area was dominated by westerly winds. The continuous winds from the west caused ice floes to drift eastward and compress the western boundary of the landfast ice, resulting in an eastward movement. Simultaneous dilatation was observed along the western boundary, which typically occurs under such forcing conditions [60]. In particular, the largest horizontal displacement rate for P 3 occurred between 4 January 2019 and 5 March 2019, with an accumulative deformation of 14 cm. As shown in Figure 12, the wind direction remained unchanged, while the speed rose from 17 km/h to 27 km/h, indicating an enhanced wind force. The combined actions of the wind speed and direction caused the most severe horizontal movement in this period.
Ice draft data collected from an ice profiler are shown in Figure 4. Ice growth began in late October every year and consistently maintained an upward trend, exhibiting some fluctuations before reaching the maximum value of approximately 1.8 m in the middle of or in late May. In general, a 9:1 proportional relationship existed between the ice draft and freeboard [61], that is, the thickness of the freeboard would increase by approximately 20 cm when the increment in the ice draft reaches 1.8 m. This assumption was inconsistent with the deformation magnitude observed in the current work. Therefore, except for the growth of ice, variations in sea level should be taken into account with regard to the vertical movement of landfast ice. Absolute sea-level changes cannot be reflected in an interferogram, as landfast ice is afloat, but sea-level tilt affects the distribution of vertical movement. Fringes can be affected by sea-level variation in accordance with ∆R = ϕ/2k , where ϕ defines the interferometric phase, and k = 2π/λ represents the wave number [62]. Horizontal movements could be estimated through the placement of oceanographic stations. Two oceanographic stations, Austin Bay and Cambridge Bay, could provide the required tidal data with a sampling frequency of 1 h. The corresponding spatial distribution of vertical deformation could be derived by measuring the differences in sea-level changes between the two meteorological stations during the same period. For example, the sea level in Austin fell from 18 cm to 8 cm from 31 December 2018 to 22 April 2019, whereas in Cambridge Bay, a reduction from 46 cm to 45 cm was observed. Therefore, the maximum relative change in the water level was 9 cm. According to the statistical results presented in Figure 13, the relative sea-level variation amplitude in the Austin area was higher than that in Cambridge Bay during most periods, leading to a phenomenon in which the settlement magnitude in the west was greater than that in the east. This phenomenon was consistent with the spatial distribution of the vertical deformation.

Possible Causes of Interannual Variation
In accordance with the interannual experimental results, the variability of landfast ice in the study area tended to be stable. The magnitude of deformation in 2019 was significantly higher than those of the other two years. By collating and analyzing the external data, we determined that the three-year average temperatures during the period covered by the SAR images were −27.4°C, −28.4°C, and −27.6°C, with slight variations. Meanwhile, the maximum temperature was considerably below 0°C. The possibility that melting caused the severe deformation of the ice is inexistent. According to the previous analysis, ice growth and sea-level tilt are the most likely factors that affect vertical movement in this area. The growth rate and maximum value of the ice draft in these 3 years were nearly equal, as shown in Figure 4. The sea-level tilt for each acquisition date in 2019 was evidently larger than those for the other two years. For example, the east-west tidal tilt from January to March 2019 was 0.5 cm, while those for the same periods in 2020 and 2021 were only 0.1 cm and 0.2 cm, respectively. We speculate that this phenomenon may have been one of the most important reasons for the vertical deformation of landfast ice in 2019. During the three studied years, the wind was dominated by westerlies, with minimal variation in wind speed. No extreme weather or ocean current events were reported in the study area in 2019-2021. Consequently, the deformation characteristics in the horizontal direction remained relatively immobile. In addition, the random internal motion of landfast ice would also have an influence on the deformation results. This influence was not considered in this work.

Conclusions
In the current work, the time-series 2D deformation that comprised the horizontal and vertical directions over the landfast ice area in Cambridge Bay was investigated by using ascending and descending datasets from Sentinel-1A and the multi-dimensional SBAS approach. This work is the first to utilize the multi-dimensional SBAS approach to explore the interannual variations (2019-2021) in landfast ice. We first obtained the LOS deformation of landfast ice from January to April in 2019-2021, and we discovered that the experimental results obtained from the ascending and descending datasets tended to exhibit opposite deformation directions due to differences in their imaging geometries. Then, the time-series 2D deformation could be decomposed from the ascending and descending LOS deformation by introducing a transform model. The derived results revealed that the vertical deformation in the study area was mostly within the interval of [−65, 23] cm, while the horizontal deformation was largely within the range of [−26, 78] cm. The maximum annual deformation rate of landfast ice occurred from early January to the end of March, while the movement remains relatively flat from April to May. Furthermore, the comparison of the results for the three years indicated that the magnitude of the timeseries deformation results obtained in 2019 was significantly larger than those in the other two years, particularly in the vertical direction.
With the collected external data, including tides, temperature, ice draft, and wind, the factors that affected deformation were classified into two categories in accordance with the deformation directions. The first category included the influencing factors associated with vertical deformation, such as ice growth, temperature, and sea-level tilt. The second category comprised the inducers acting on the horizontal deformation, such as ocean currents and winds. We concluded that wind and sea-level tilt dominated the variability in the deformation in Cambridge Bay, with the largest contributor to the magnitude anomalies of deformation in 2019 being the differences in sea-level tilt changes.
Overall, this work utilized the SBAS-InSAR method to determine interannual landfast ice deformation, which is particularly important for local populations and organisms in anticipating potential hazard events and providing effective early warnings in real time. In addition, this proved that the multi-temporal InSAR technology possesses the capacity of being introduced into surface deformation monitoring in the cryosphere. This work still has some weaknesses. Firstly, deriving the time-series deformation for a period of more than 5 or 10 years is difficult due to the limited data that are currently available. This issue may be addressed in the future with the launches of more SAR satellites, such as Sentinel-1C, ALOS-3, and Tandem-L. Secondly, the internal variability, such as the dynamics of brine inclusions and ice rheology, was not considered when analyzing the influencing factors. Therefore, we will be involved in using longer time series and multi-source data to conduct research and analyze the deformation patterns from a multidisciplinary perspective in the future. Acknowledgments: The authors would like to acknowledge the official website of the Government of Canada (https://www.canada.ca/en.html (accessed on 21 January 2023)) for providing the temperature, tide, and wind data. The ice draft data are available from the Ocean Networks Canada Data Archive (https://data.oceannetworks.ca/ (accessed on 21 January 2023).

Conflicts of Interest:
The authors declare no conflicts of interest.