Mapping pan-Arctic landfast sea ice stability using Sentinel-1 interferometry

Arctic landfast sea ice has undergone substantial changes in recent decades, affecting ice stability and including potential impacts on ice travel by coastal populations and on industry ice roads. We present a novel approach for evaluating landfast sea ice stability on a pan-Arctic scale using Synthetic Aperture Radar Interferometry (InSAR). Using Sentinel-1 images from spring 2017, we discriminate between bottomfast, stabilized, and nonstabilized landfast ice over the main marginal seas of the Arctic Ocean (Beaufort, Chukchi, East Siberian, Laptev, and Kara seas). This approach draws on the evaluation of relative changes in interferometric fringe patterns. This first comprehensive assessment of Arctic bottomfast sea ice extent has revealed that most of the bottomfast sea ice is situated around river mouths and coastal shallows. The Laptev and East Siberian seas dominate the aerial extent, covering roughly 4100 and 5100 km2, respectively. These seas also contain the largest extent of stabilized and nonstabilized landfast ice, but are subject to the largest uncertainties surrounding the mapping scheme. Even so, we demonstrate the potential for using InSAR for assessing the stability of landfast ice in several key regions around the Arctic, providing a new understanding of how stability may vary between regions. InSAR-derived stability may serve for strategic planning and tactical decision support for different uses of coastal ice. In a case study of the Nares Strait, we demonstrate that interferograms may reveal early-warning signals for the breakup of stationary sea ice.

Abstract. Arctic landfast sea ice has undergone substantial changes in recent decades, affecting ice stability and including potential impacts on ice travel by coastal populations and on industry ice roads. We present a novel approach for evaluating landfast sea ice stability on a pan-Arctic scale using Synthetic Aperture Radar Interferometry (InSAR). Using Sentinel-1 images from spring 2017, we discriminate between bottomfast, stabilized, and nonstabilized landfast ice over the main marginal seas of the Arctic Ocean (Beaufort, Chukchi, East Siberian, Laptev, and Kara seas). This approach draws on the evaluation of relative changes in interferometric fringe patterns. This first comprehensive assessment of Arctic bottomfast sea ice extent has revealed that most of the bottomfast sea ice is situated around river mouths and coastal shallows. The Laptev and East Siberian seas dominate the aerial extent, covering roughly 4100 and 5100 km 2 , respectively. These seas also contain the largest extent of stabilized and nonstabilized landfast ice, but are subject to the largest uncertainties surrounding the mapping scheme. Even so, we demonstrate the potential for using In-SAR for assessing the stability of landfast ice in several key regions around the Arctic, providing a new understanding of how stability may vary between regions. InSAR-derived stability may serve for strategic planning and tactical decision support for different uses of coastal ice. In a case study of the Nares Strait, we demonstrate that interferograms may reveal early-warning signals for the breakup of stationary sea ice.

Landfast sea ice stability and stakeholder dependence
Sea ice is an important component of Arctic ecosystems and provides important functions as a climate regulator (Screen and Simmonds, 2010), habitat for marine biota (Thomas, 2017), and a platform for coastal populations . During the last century, an expansion of transportation and resource extraction have led to increased human presence in the Arctic and further diversification of ice use (Eicken et al., 2009). The recent retreat of sea ice observed over the past several decades (Stroeve et al., 2012;Comiso and Hall, 2014;Meier et al., 2014) has already resulted in widespread consequences for ice users (ACIA, 2004;Aporta and Higgs, 2005;Fienup-Riordan and Rearden, 2010;Orviku et al., 2011;Druckenmiller et al., 2013), as well as increasing hazards (Ford et al., 2008;Eicken and Mahoney, 2015). At the same time, the related increased accessibility to Arctic waters (Stephenson et al., 2011) is leading to more ship traffic and resource exploration (Lovecraft and Eicken, 2011;Eguíluz et al., 2016). It is further recognized that sea ice conditions for future Arctic marine operations will be challenging, and will require substantial monitoring and improved observations (Arctic Council, 2009). This improvement will require observations on local and regional scales in order to provide an assessment of environmental hazards and effective emergency responses . Most of the Arctic Ocean is dominated by drifting pack ice, whereas stationary landfast ice occupies much of the  , derived from sea ice charts based on optical instruments and SAR. The data for this figure were obtained from the National Snow and Ice Data Center (Yu et al., 2014).
Arctic coastlines roughly between November and June, depending on location ( Fig. 1; Yu et al., 2014). Sections of landfast ice, often from several kilometers to hundreds of kilometers wide, are held in place by grounded ridges, islands, or coastline morphology, such as embayments or fjords. Similar to drifting pack ice, landfast ice has declined significantly during the last few decades, particularly in terms of delayed freeze-up in the Beaufort (Mahoney et al., 2014) and Laptev (Selyuzhenok et al., 2015) seas, as well as to a significantly reduced extent in the Chukchi Sea (Mahoney et al., 2014). Later freeze-up critically impacts stakeholders through reduced stability of landfast ice in response to fewer grounded ridges capable of withstanding wind, ocean, or ice forcing (Dammann, 2017). Previous research suggests that landfast ice stability can be expressed in terms of combined frictional resistance provided by relevant grounding or attachment points (e.g., islands and grounded ridges; Mahoney et al., 2007;Druckenmiller, 2011). Although landfast ice is stationary, it deforms at the centimeter to meter scale, on timescales of days to months due to forcing from wind, currents, and drifting ice (Dammann et al., 2016). Its stability in part determines the rate at which the ice deforms, and ultimately, the severity of breakout events or magnitude of structural defects. We suggest that landfast sea ice can be further categorized into three regimes, defined through their respective stabilities: (1) bottomfast ice, (2) floating ice enclosed in lagoons or fjords, or sheltered by point features such as grounded ridges or islands, and (3) floating ice extensions (Table 1). A typical landfast ice regime is illustrated in Fig. 2, where the stability of the landfast ice area decreases from the coast toward the open ocean (Dammann et al., 2016).
Bottomfast sea ice can grow laterally to kilometer scale during winter, depending on local bathymetry Stevens et al., 2010). This bottomfast ice allows for heat loss from the seafloor and is therefore an integral part of aggregating and maintaining subsea permafrost Stevens et al., 2008Stevens et al., , 2010Stevens, 2011), as well as controlling coastal stability and morphology (Are and Reimnitz, 2000;Eicken et al., 2005). Bottomfast ice is also relevant for fish, as it reduces habitable shallow waters during winter (Hirose et al., 2008), and for on-ice operations, as it can support a much larger load than floating ice. High to moderately stable landfast ice is of relevance to industrial (Potter et al., 1981) and subsistence ice use (Druckenmiller et al., 2013), as well as for habitats (Tibbles et al., 2018). Ringed seals, for instance, are dependent on stable landfast ice for denning (Smith, 1980). Low-stability ice is potentially relevant for ocean-based operations, such as shipping through trans-Arctic passages close to the coast, where patches of landfast ice occasionally break off and drift into nearby shipping lanes, potentially causing hazards. Even areas hundreds of kilometers from landfast ice can be impacted through the failure of ice arches.
Ice arches may be considered an additional zone of "temporarily stabilized pack ice". Ice arches form when ice moving through a narrow passage experiences flow stoppage as a result of confining pressure and begins to behave like landfast ice, though potentially without cohesive strength between individual floes (Hibler et al., 2006). Ice arches typically form between November and March (Moore and McNeil, 2018), and can block the export of ice through straits as wide as 100 km (Melling, 2002). When formed, such arches represent a significant obstacle to marine traffic, due to the high confining pressures that make icebreaking impossible for all but the most powerful vessels. Arches can in some locations prevail into the following season (Melling, 2002), but typ- ically collapse in July-August (Kwok, 2005). Conversely, their breakup can lead to the advection of large amounts of thick multiyear ice into high-traffic shipping routes (Barber et al., 2018), creating a well-known hazard (Bailey, 1957;Wilson et al., 2004;Howell et al., 2013). Recent and ongoing sea ice decline is leading to an increasing presence of thinner ice in the Canadian Archipelago (Haas and Howell, 2015), and weaker ice due to warmer temperatures (Melling, 2002) may lead to earlier breaching of ice arches. This could result in a larger amount of advected ice with potentially longer travel paths, increasing the severity of such events (Melling, 2002;Barber et al., 2018). One location of particular interest is the Nares Strait, situated between Greenland and Ellesmere Island, featuring a seasonal ice arch (Kwok, 2005;Kwok et al., 2010) with important implications for the multiyear ice budget of the Arctic Ocean (Kwok et al., 2010). This stability is also relevant for destinational cargo shipping in the Arctic, as less stable, thinner ice is easier to break through, resulting in opportunities for docking in areas of landfast ice. For navigating through landfast ice, stabilization through ridging is also important to identify, as ridges can be problematic to navigate and are often associated with hazards (Hui et al., 2017).

Remote sensing of landfast ice stability
Satellite remote sensing is an important tool for measuring ice conditions in the Arctic, including the mapping of landfast ice (Muckenhuber and Sandven, 2017). Optical and thermal satellite data, such as from the Advanced Very High Resolution Radiometer (AVHRR), were used to produce operational ice charts until the early 1990s, when SAR was introduced into the charting production (Yu et al., 2014) as a superior data set, thanks to its independence from light and weather conditions and its higher (∼ 100 m) resolution -both advantageous to stakeholders . Different techniques exist to map the boundaries of land-fast sea ice, typically derived by evaluating unchanged sections of ice between consecutive SAR backscatter scenes (Johannessen et al., 2006;Giles et al., 2008;Mahoney et al., 2014). In addition to its use in the mapping of landfast ice, SAR backscatter can also discriminate between multiyear and first-year ice (Onstott, 1992) and identify different roughness regimes . SAR has also been used to estimate the advection of ice through straits in the Canadian Archipelago (Melling, 2002;Kwok, 2006;Howell et al., 2013). However, SAR backscatter typically does not give information pertaining to the stability of landfast ice or temporarily stabilized pack ice, since the internal movement of the landfast ice is too small (mm day −1 ) to be identified with change detection.
SAR interferometry (InSAR) is a signal-processing technique, which extracts the phase difference between two SAR images acquired with similar viewing geometries. This phase difference (typically referred to as interferometric phase) can either signify topography if acquisitions are separated in space (i.e., nonzero perpendicular baseline) or measure the line-of-sight motion if acquisitions are separated in time (nonzero temporal baseline; Bamler and Hartl, 1998;Ferretti et al., 2007). InSAR has been successfully used to map landfast ice edges by identifying areas of slow motion . InSAR has also provided information pertaining to landfast ice dynamics (Li et al., 1996;Morris et al., 1999;Vincent et al., 2004;Marbouti et al., 2017) and topography Dierking et al., 2017) by evaluating the phase change between acquisitions. InSAR has also been shown to reveal plausible rheologies for landfast ice (Dammert et al., 1998) and has been used to determine the origin of internal ice stresses (Berg et al., 2015). Combined with inverse modeling, InSAR enables us to determine ice deformation modes (Dammann et al., 2016), rates, and the associated stress and fracture potentials (Dammann et al., 2018b). These studies have demonstrated (1) the potential of In-SAR as a tool for assessing landfast ice dynamics and stability through local case studies and (2) its utility as a planning tool for on-ice operations (Dammann et al., 2018a, b). They have also laid the foundation for applying InSAR on a larger scale, potentially as a means for generating operational information products and evaluating long-term trends. The coverage and access to InSAR-compatible SAR scenes has been an obstacle in the past, but has improved significantly since the launch of Sentinel-1 (the suitability of Sentinel-1 for automatic SAR processing has been shown, e.g., in Meyer et al., 2015). In this study, we explore InSAR as a tool for providing pan-Arctic information on ice stability, which is relevant to subsea permafrost, biological habitats, and sea ice use. The goal of this work is to determine Sentinel-1 interferometric data availability along substantial parts of the circumpolar coastlines and to explore their applicability for consistently mapping landfast sea ice stability zones in different geographic regions. We further explore the limitations of the technology and possible applications.
2 Data and methods

InSAR principles
The interferometric phase may be related to lateral (e.g., thermal contraction or displacement due to compressional or shear forces) or vertical (e.g., through buckling or tidal displacement) sea ice motion that occurs in between the acquisition times for the two SAR images. A phase signature can sometimes also be attributed to factors not related to surface motion, such as atmospheric phase delay. Of the phase change attributed to motion, only displacement in line-ofsight direction ( r LOS ) results in a phase change disp , according to disp = 4π r LOS /λ. The observed phase is measured within the wrapped interval of [0; 2π [. The interferogram is a series of fringes representing the projection of the true three-dimensional ice motion onto the line-of-sight vector. The orientation of fringes can be used to interpret the direction of the three-dimensional motion field, and fringe spacing is an indicator of the deformation rate. The interpretation of observed fringe patterns is, however, not straightforward, and it typically requires the use of an inverse model (Dammann et al., 2016). The interferometric phase values will only be useful if scattering elements remain largely unchanged throughout the time interval bracketed by the image pairs used in processing. InSAR phase stability, referred to as InSAR coherence, depends largely on topography coupled with a perpendicular baseline, as well as the temporal stability of the scatterers on the ground surface. Coherence ranges between 0 (pure noise) and 1 (no noise), and serves as a measure of the quality of the interferogram. Coherence is generally high if scatterers remain unchanged and low if there is significant change in the scattering medium .

Sentinel-1 data
This study uses Sentinel-1, a constellation of two C-band SAR systems (Sentinel-1A and B) in operation since 2014 and 2016, with a repeat-pass interval of 6 to 12 days, depending on whether both or just one of the satellites acquire data. Thanks to a free-and-open data policy and large spatial coverage, we obtained Sentinel-1 acquisitions for five marginal seas of the Arctic Ocean, enabling mapping of landfast sea ice on a pan-Arctic scale. All images used were captured in interferometric wideswath (IW) mode, with a single-look resolution of roughly 3 m × 22 m in slant range and azimuth, and a ∼ 250 km swath width. Images were acquired almost exclusively between March and May 2017 (see Supplement for full list of images used). We acquired over 100 SAR images, covering nearly the entire continental coastlines of the Beaufort, Chukchi, East Siberian, Laptev, and Kara seas. To reduce computational costs, we omitted Greenland, some island groups, and the Canadian Archipelago, which are characterized by extensive coastline lengths. In this work, we focused on the Alaska and Russian marginal seas of the Arctic Ocean. These coastlines have high economic significance for the shipping and natural resource industries, and also feature dynamically diverse ice regimes. Large areas of bottomfast ice are expected in these regions. Except for one approximately 50 km long section of coast in the Kara Sea and the eastern Laptev Sea, multiple InSAR compatible pairs were available for the specified time frame. This allowed us to select interferograms centered around the end of April, when most Arctic landfast ice is at its maximum extent and thickness.
In addition to images obtained for the large-scale mapping of stability zones, a series of six consecutive image pairs were acquired covering the Nares Strait and the breakup of an ice arch during spring 2017. This image sequence featured a 6-day temporal baseline covering a time span of 36 days.

Data processing
All the complex Sentinel-1 data were processed to obtain backscatter in order to visually interpret features (e.g., landfast ice edge, fracturing, ice roughness and types). Data were further processed for interferometry. Depending on the perpendicular baseline, sea ice topography can have a modest impact on the phase difference. Due to tight baseline limits (∼ 50 m standard deviation) for the Sentinel-1 constellation, and since sea ice topography rarely exceeds 10 m, impacts for interferogram interpretation are minimal for the data shown here. In this work we predominately employed acquisitions with a temporal lag of either 12 or 24 days, depending on data availability. For this timespan, the coherence over landfast ice was found to be generally high, due to its stationary nature.
However, coherence loss was evident in some areas -in particular in the Chukchi Sea, such as in the Kotzebue Sound. This was likely predominately due to surface melt, as air temperature reached above freezing between SAR acquisitions. Other possible contributing reasons for coherence loss in this region could include ice motion, subsurface ice thinning from river runoff, and low signal-to-noise ratio. Significant decorrelation can also occur in late spring, as the onset of melt causes substantial changes in the scattering medium. In this work, we obtained images as close to late April as possible. This time frame was found to be ideal for our purpose, as ice thickness is near its maximum, leading to maximum stability and minimizing impacts from the onset of melt. To ensure a realistic representation of an operationally produced synoptic, contiguous pan-Arctic interferogram, we did not attempt to derive alternative interferograms (i.e., different time periods) in cases of low coherence.
All backscatter images and interferograms in this work were produced using a standard Sentinel-1 workflow in Gamma software (Werner et al., 2000). The IW images initially consist of independent bursts and swaths, which we combined to utilize the full extent of the acquisition. We further coregistered pairs of acquisitions to ensure images cover exactly the same area with subpixel accuracy. Images were then multilooked, averaging 10 pixels in range and 2 pixels in azimuth, resulting in reduced speckle and a final pixel spacing of roughly 23 × 28 m. Next, spectral filtering was performed to ensure both images comprise the same spectral range, reducing phase noise in the final interferogram. The interferometric phase was calculated for each pixel of the coregistered and filtered images. Furthermore, the expected phase ramp in cross-track direction from a stationary flat earth surface was removed. The phase noise of the final interferogram was reduced using an adaptive phase filter (Goldstein and Werner, 1998). The result was 52 interferograms, covering almost the entire coastline between the Canadian Archipelago and the Barents Sea.

Mapping landfast ice stability zones
In this work, we evaluated relative ice stability based on fringe spacings within individual interferograms. This allowed us to identify variations within an area imaged under largely the same conditions (e.g., close to the same wind and ocean forcing). Trends from higher to lower fringe density will likely correspond to increasing ice stability. Therefore, interferograms can provide information related to the spatial variations in stability. Meyer et al. (2011) demonstrated that interferometry can be used to map the landfast ice based on a coherent phase response. Their work also suggested that fringe patterns are significantly impacted by grounded ridges and reduced fringe density. Dammann et al. (2018c) further showed that bottomfast ice can be mapped based on a nearzero phase change where ice is frozen to the seafloor. We built on this work by suggesting that InSAR can be used to map three different zones of relative stability: bottomfast ice, stabilized ice, and nonstabilized ice (Table 1). These three zones are subjectively and manually mapped without the use of specific threshold values.
Bottomfast ice is identified with near-zero phase change in the interferogram. It can often be distinguished from adjacent floating ice commonly featuring a nonzero phase change or low coherence (Fig. 3a). The shoreward boundary of bottomfast ice is difficult to obtain from the interferogram alone, since the phase signatures over bottomfast ice and low-lying coastal areas are similar. The use of the landmask (Wessel and Smith, 1996) is not ideal, since subtle coastal features such as sediment bars are often not captured. We therefore delineated the coastline (i.e., shoreward boundary of bottomfast ice) using the backscatter signature in a composite image with backscatter and phase (Fig. 3b). Plotting bottomfast ice with the landmask can thus give the wrongful appearance that (1) areas of near-zero phase should have been mapped as bottomfast and (2) bottomfast ice appears in sporadic areas along the coast separated by floating ice.
We could often identify stabilized ice by a stark fringe discontinuity between different fringe densities (Fig. 3c). However, in some regions, changes in stability are more gradual between zones (Figs. 6c and 7c). Mapping of such regions are therefore more subjective and possibly less exact. In cases lacking stark fringe discontinuities, stabilized ice is also mapped in regions featuring a very slight phase response with no clear fringe patterns (Fig. 3d).
Nonstabilized ice is identified as landfast ice (i.e., areas featuring relatively high interferometric coherence) otherwise not marked as bottomfast or stabilized ice. Nonstabilized ice commonly features clear fringe patterns (Fig. 3e).

Validation areas and data
The Beaufort Sea coast of Alaska was used for validation, as landfast sea ice in this area includes all three stability zones (bottomfast, stabilized, and nonstabilized), and ample validation data are available from previous studies. Alaska's Beaufort Sea coast is of major interest in the context of local and indigenous ice use, as well as industry resource exploration and extraction. Some areas along this coastline feature similar landfast ice extent over timescales from months to years. It was found that these regions ("nodes") of consistent landfast ice extent are often tied to the location of the 20 m isobath, a water depth associated with the grounding of pressure ridges stabilizing the landfast ice (Mahoney et al., 2014). Indigenous knowledge and a field study also indicate persistent grounded ridges in the location of the node closest to Utqiaġvik, Alaska . We also evaluated our results near Stolbovoy Island in the Laptev Sea. This area features a shoal of < 10 m water depth, leading to earlier formation of fast ice than the surrounding areas (Selyuzhenok et al., 2015), likely due to the formation of grounded ridges on the shoal and resulting in increased stability.  (Mahoney et al., 2004) or can break-off from the stabilized ice. 3 Results

Evaluating landfast ice stability zones
We constructed a series of Sentinel-1 interferograms along the coastlines of five marginal seas in the Arctic Ocean during 2017: the Beaufort, Chukchi, East Siberian, Laptev, and Kara seas. As seen in the interferograms , landfast sea ice varies substantially between the seas in terms of the extent and interferometric fringe density.
The landfast sea ice extent in the Beaufort Sea ranges from almost zero up to 100 km (Fig. 4a). River outlets such as the Colville and Mackenzie deltas feature extensive regions of bottomfast ice several kilometers wide (Fig. 4b, c). Bottomfast ice is also prominent in many lagoons along the coast. Much of the floating ice along the coast from Prudhoe Bay to Point Barrow is stabilized. This stabilized ice can be identified by a stark fringe discontinuity separating regions of different fringe density and stability (Fig. 4a, b). The line of discontinuity features several seaward points, an expected pat- tern surrounding grounded ridges. This is because grounded ridges result in a shoreward increase in stability that does not extend to areas immediately to the side of the ridges (the alongshore direction). Examples of likely grounding points are indicated with white arrows in Fig. 4b, and similar patterns are also apparent near the Mackenzie Delta (Fig. 4c). The landfast ice in the eastern part of the Beaufort Sea also consists of large areas of stabilized ice. Here, the landfast ice is noticeably sheltered by land features, resulting in lowerdensity fringes (Fig. 4d).
Landfast sea ice in the Chukchi Sea is generally less extensive than in the Beaufort Sea, particularly along the Russian coast (Fig. 5a). Bottomfast ice in the Chukchi is constrained mostly to lagoons. Some of these lagoons, such as the Kasegaluk, consist almost exclusively of bottomfast ice (Fig. 5b). Only a few areas of landfast ice appear to be stabilized, including the northern coast of Alaska near Peard Bay (Fig. 5c) and the southern Chukchi Sea near Shishmaref (Fig. 5d). The Chukchi Sea consists predominantly of nonstabilized ice, with the most extensive region of landfast ice situated off the shore of the village of Shishmaref (Fig. 5d). The Chukchi Sea features coherence loss in several regions such as the Kotzebue Sound (Fig. 5a).
The landfast ice in the East Siberian Sea is more extensive than in the Chukchi and Beaufort seas and can extend over 100 km from the shore (Fig. 6a). Bottomfast ice is also more extensive than in the Beaufort and Chukchi seas. The bottomfast ice in the East Siberian Sea follows several sections of coastline even tens of kilometers away from major rivers, though most of the bottomfast ice is situated near the Kolyma and Indigirka deltas (Fig. 6b). In contrast to the Beaufort and Chukchi seas, stabilized ice extends several tens of kilometers offshore without being sheltered by coastline morphology or islands (Fig. 6c). These large areas also lack clear indications of the presence of grounded ridges (Fig. 6d).
Landfast ice in the Laptev Sea, similar to the East Siberian Sea, extends upwards of 100 km from the shore (Fig. 7a).
Here, most of the bottomfast ice is situated around river outlets and in particular near the Lena Delta, extending tens of kilometers from shore (Fig. 7b). This delta features a large amount of small, low-lying land areas (e.g., gravel islands) only partly covered by the landmask. This has made it problematic to delineate all areas of bottomfast ice and led to more approximate delineations than in the other deltas mapped. On the east side of the Lena Delta and south of the Great Lyakhovsky Island, there are extensive sections of stabilized ice (Fig. 7c). Some regions of the eastern Laptev Sea lack a clear discontinuity, but at the same time feature locally reduced fringe density, indicative of stabilized ice (Fig. 7c). We also considered these areas to be stabilized (Fig. 7c), though possibly as a result of different ice types or thicknesses, rather than through grounding or sheltering. However, one offshore area is clearly identified as stable by a lack of consistent fringe patterns and a clear discontinuity, likely due to grounded ridges (Fig. 7d). Landfast ice in the Kara Sea features a much smaller ice extent than the other Russian seas (Fig. 8a). Bottomfast ice is also much less prevalent and largely situated near the Pyasina River (Fig. 8b). The landfast ice extends tens of kilometers from shore, predominately in areas supported by offshore islands and archipelagos (Fig. 8c, d). In these archipelagos, the ice confined by islands is largely stabilized (Fig. 8c).
Interferograms have enabled the mapping of landfast ice stability zones based on subjective interpretations of interferometric fringes (Fig. 9). The resulting stability map allows for a large-scale comparison and analysis of bottomfast, stabilized, and nonstabilized landfast ice, within and between the different seas. For this comparison, we have calculated the area of each stability zone (Table 2). However, it is important to note these area calculations are not complete, as the analysis omitted some island groups and included some data gaps.
Most areas with extensive bottomfast ice reaching several kilometers from shore are located either in the vicinity of river deltas or within lagoons. The East Siberian Sea and its three large river systems (the Indigirka, Bogdashkina, and Kolyma rivers) contain the most bottomfast ice of the regions considered here. The Laptev Sea also contains a large area of bottomfast sea ice. Together, the Laptev and East Siberian seas contain over half (∼ 57 %) of the total areal extent of bottomfast ice calculated, while the Chukchi Sea features the lowest extent of bottomfast ice of the regions considered here. Bottomfast ice is predominately situated in lagoons.
Stabilized ice was found in all marginal seas (Fig. 9), though its relative contributions to overall landfast ice extent varied widely. The largest extent of stabilized landfast ice in our study region was found in the Laptev and East Siberian seas. These regions feature particularly large continuous areas of stabilized ice labeled A-F in Fig. 9. Even so, as we  The bottomfast ice zone is constrained between its outer extent (interpreted from the phase) and the coast (as interpreted from the backscatter scenes). The stabilized zone is constrained between its outer extent (as interpreted from the phase) and the bottomfast ice or the landmask (Wessel and Smith, 1996). The nonstabilized ice is constrained between the outer extent of nonzero coherence and the bottomfast ice, stabilized ice, or the landmask. delineated here, the Beaufort Sea is the only sea that features more stabilized ice than nonstabilized ice. This is likely attributed to the large grounded sections, as well as areas sheltered by coastal morphology. The Laptev Sea also features large areas confined by coastlines. However, in the Laptev sea, these regions also commonly feature nonstabilized ice. Meanwhile a large part of the landfast ice in the Kara Sea is mapped as stabilized, largely due to the fraction of landfast ice situated between islands and archipelagos. With a relatively narrow landfast ice extent compared to other seas and the absence of regions of sheltered ice, the Chukchi Sea contains the lowest total extent of stabilized ice.
In the Chukchi Sea, we identified the vast majority of landfast ice as nonstabilized (Fig. 9), resulting in our largest areal fraction (nonstabilized ice vs. stabilized ice). Though the largest total areas of nonstabilized ice can be found in the Laptev and East Siberian seas. Here, the distinction between stabilized and nonstabilized landfast ice is not as straightforward as in the Beaufort and Chukchi seas, due to a lack of clear boundaries between areas of different fringe densities. Even so, it is clear that landfast ice extent in the East Siberian and Laptev seas is dominated by vast areas of nonstabilized ice. However, unlike the Chukchi Sea, we also identified significant areas of stabilized landfast ice along these two seas. The Kara Sea features predominately nonstabilized ice along the coast and along the outer margins of archipelagos.

Evaluating stability of temporarily stabilized pack ice
Sentinel-1 SAR backscatter imagery captured the location and breakup of the ice arch in Nares Strait in 2017 (Fig. 10). This breakup event occurred relatively early compared with past events (Kwok, 2005), partly in response to thinner ice conditions and northerly winds (Moore and McNeil, 2018). The arch appeared stable on 6 May (Fig. 10b), before eventually failing sometime before 12 May (Fig. 10c) as seen in the SAR backscatter images. The interferograms revealed the ice deformation around the location of fracture up until the failure event. As seen in the interferograms, the ice arch features various levels of centimeter-to meter-scale deformation and fractures prior to breakup, resulting in fringe discontinu-ities ( Fig. 11a) most pronounced near the arch terminus to the south. Near the failure line, there was no sign of a fringe discontinuity up until 12 April (Fig. 11a), when the interferogram displays near-cross-track parallel fringes, indicating compression towards the terminus (Fig. 11b). There is a significant contrast in fringe density on either side of the line, which may be indicative of a fracture, where the ice to the west is being compressed more rapidly than the ice close to the coast. The interferogram between 18 and 24 April features widespread coherence loss, possibly due to continued compression (Fig. 11c). Deformation is less severe from 24 April, when the fringe density is significantly reduced. However, we did notice a fringe discontinuity to the east of the failure line, featuring perpendicular intermediate fringe patterns following late April (see arrows in Fig. 11d). These patterns develop further into circular patterns often associated with vertical lifts and depressions (Fig. 11e) whole arch appears to fail through shear motion along this same fault (Fig. 11f).

Validating stability zones with areas of known ice stability
The InSAR technique used to map bottomfast sea ice was thoroughly validated in several regions by Dammann et al. (2018c). The high stability of these regions can be inferred from the ice resting on the seafloor. However, other stability zones (i.e., stabilized and nonstabilized ice) are based on relative stability, in terms of whether the ice is anchored or sheltered. Determining absolute stability (i.e., whether an area is stable enough for a specific use, such as ice roads) would be problematic using fringe density alone. This is because there are many factors that affect fringe density in addition to stability, including changing wind and ocean currents, satellite viewing geometry, and the prevalent mode of ice deformation (Dammann et al., 2016). A measure of whether ice is practically stable would also depend on specific stakeholders and their dependence on stability. For example, on shorter timescales, industry ice roads would be able to accommodate less strain than community ice trails, due to different modes of transportation and user-specific needs. Further steps to identify such thresholds are outlined in Dammann et al. (2018a).
There is limited information that can be used to validate these stability classes -namely, the separation between stabilized-and nonstabilized ice. Even so, we compared our mapping approach here with one region in the Beaufort and one in the Laptev Sea with areas of known stabilization points. We examined a backscatter mosaic from three images (8, 15, and 17 April) along the Beaufort Sea coast. These images exhibit, in certain locations, a sharp discontinuity in backscatter, which can identify the location of the landfast ice edge (see white arrows in Fig. 12a). The landfast ice edge identified using backscatter is consistent with the three nodes (A-C) identified by Mahoney et al. (2007Mahoney et al. ( , 2014our Fig. 12b). These nodes signify a persistent landfast ice edge, believed to be a result of reoccurring grounded ice features (Mahoney et al., 2014). The ice shoreward of these three nodes is expected to be stabilized, because grounded ridges are known to stabilize landfast ice, leading to reduced strain shoreward of the grounding points (Mahoney et al., 2007;Druckenmiller, 2011). Interferograms exhibit a phase response, suggesting stabilized ice directly shoreward of nodes A and C (Fig. 12c). Here, node A is known to correspond to the location of large grounded ridges, offering stability to the ice cover . Nodes B and C are also expected to be regions of persistent grounded ridges since the nodes coincide roughly with the 20 m isobath (Mahoney et al., 2014). However, ice directly shoreward of node B appears nonstabilized, with stabilization only further in. This may be due to the reduced keel depth of ridges in 2017 or the possibly reduced grounding strength of ridges present in Node B. Certain sections of the border between stabilized and nonstabilized ice extend relatively far from the coast (see black arrows in Fig. 12d). At these points, the stability is higher than adjacent areas with the same distance from shore. This is consistent with increased stability behind grounded ridges.
Although the landfast ice edge can in some instances be mapped using a single backscatter image, stabilized ice cannot easily be discriminated from nonstabilized ice. This is apparent when comparing grounding locations as obtained with InSAR with backscatter images (see black arrows in Fig. 12a). It is also worth noting that relying on backscatter to discriminate landfast or drifting ice only works in some cases. There must be noticeable differences in backscatter between landfast and drifting ice, or a severely deformed landfast ice edge as a result of shear interaction with the pack ice (Druckenmiller et al., 2013).
Similar patterns indicating grounded ridges were found in the Laptev Sea, where an April interferogram shows a section of stabilized ice roughly 100 km offshore (see A in Fig. 13). The full area extent of the stabilized ice cannot be deter- mined due to limited data availability in the region, and one interferogram had to be acquired as early as February, before this region had stabilized. Stabilized ice is expected in this region, which features a large shoal, earlier ice formation than surrounding areas, and grounded ridges (Selyuzhenok et al., 2015). The location of this large shoal, along with smaller ones, are obtained from Jakobsson et al. (2012) and displayed in Fig. 13b. Here, it is apparent that even some of the smaller shoals are associated with stabilized ice (see B and C in Fig. 13b). It is also clear that the extensive stabilized ice that stretches out halfway between Great Lakhovsky Island and Stolbovoy Island is potentially anchored between the coast and the shallow areas (see D in Fig. 13b).

Methodological limitations for mapping stability zones
There are a number of sources of uncertainty that affect our map of landfast ice and its relative stability. Dammann et al. (2018c) have determined that, in some instances, bottomfast ice has to be approximated on the subkilometer scale due to ambiguities associated with low fringe density or fringes parallel to the bottomfast ice edge. We also acknowledge that small islands or sandbars not represented by our landmask may be erroneously identified as bottomfast ice. We have reduced such errors by not mapping areas that appear to be low-lying land in the SAR backscatter images. However, discriminating between ice and low-lying land can be difficult based on strictly SAR. Here, other remote-sensing systems such as optical systems could be applied to further reduce biases from coastline errors. In areas where the landmask does not appear to fit the coastline due to errors or coastline changes, mapping intricate coastal morphology can be a time-consuming task -hence mapping on a pan-Arctic scale will inevitably contain inaccuracies. It is also worth mentioning that the other stability zones are mapped against the landmask, also likely resulting in errors. However, as the extent of these zones are larger, the relative contribution of such errors will be much smaller.
In this work, we did not apply strict mapping thresholds to distinguish between stabilized and nonstabilized ice but rather made subjective determinations based on fringe patterns. This approach works well in the Chukchi and Beaufort seas, where regions of low fringe density lie adjacent to the coast or bottomfast ice and can be easily distinguished from regions of higher fringe density. However, in some regions, especially in the Russian Arctic, there is often a lack of distinct boundaries between regions of different fringe spacing, introducing ambiguities between stabilized and nonstabilized ice on scales from kilometers to even tens of kilometers (Figs. 6c and 7c). The difficulty of making distinctions between these two zones may result from reduced pack ice interaction along the Russian shelf, given the predominately divergent ice regimes (Reimnitz et al., 1994;Alexandrov et al., 2000;Jones et al., 2016). This generally results in reduced ice forcing and landfast ice strain in contrast to the western Arctic. These regions are expected to feature reduced dynamically induced strain (and therefore fewer interferometric fringes) in nonstabilized ice, making it appear more stable. This is visible in the different fringe densities of the nonstabilized ice in Figs. 4d and 6d. Additionally, the greater extent of landfast ice on the shoreward side of the grounding points provides a greater fetch, which may cause stabilized ice on the Russian Shelf to exhibit higher fringe densities than in the Chukchi or Beaufort seas. This suggests that there is likely a spectrum of landfast ice stability. Additional zones may be necessary to fully characterize landfast ice regimes in different regions for different ice uses or research aims. Expanding on the classes presented here would likely require a different set of evaluation criteria for fringes, depending on regions. Additional data such as bathymetry would also likely strengthen this analysis.
We have focused on some examples with possibly suboptimal classification. One potential candidate for reclassification is landfast ice in sheltered bays, such as the Khatanga Gulf in the western Laptev Sea, which exhibited predominantly high fringe densities (Fig. 7a). Hence, the Khatanga Gulf was largely identified as nonstabilized, despite being nearly landlocked (Fig. 9). Due to the shallow water in this region, it is likely that the high fringe density is caused in part by vertical motion associated with tides and coastal setup. Since vertical motion has less impact on stability in well-confined landfast ice, such examples suggest the potential need for an additional zone of stability, allowing higher fringe densities in coastally confined regions. Such additional classification would depend on other data sets such as a landmask or bathymetry to identify the level of restricted ice movement in response to likely forcing conditions. Another, larger-scale example is the eastern Laptev sea, which is an area of landfast ice sheltered by the New Siberian Islands and is typically considered stable (Eicken et al., 2005). However, based on relatively high fringe density, particularly offshore of the Lena Delta, we classified much of the landfast ice in this region as nonstabilized (Fig. 9). This suggests the criteria for stabilized ice used in this analysis is different than in Eicken et al. (2005) and can provide new information related to stability in the region. Based on the overall fringe counts and patterns, the majority of the phase response is due to lateral displacement and potentially only partially due to vertical displacement (circular fringe patterns with low density -see Dammann et al., 2016) due to tidal motion. It is possible that landfast ice in this region may be less stable than previously thought, and that a partially stabilized zone may be appropriate. This would be consistent with a recent SAR backscatter analysis of landfast ice in the Laptev Sea (Selyuzhenok et al., 2017), which showed that areas identified as landfast ice in operational ice charts may actually contain pockets of partly mobile ice. This was shown for the month after initial landfast ice formation, but could possibly result in more dynamic ice throughout spring due to reduced ice thickness.
Sensitivity to specific atmospheric and oceanographic conditions during the time period between SAR acquisitions may place a limitation on the number of stability zones that can be mapped. For example, in the absence of dynamic interaction with pack ice, there may be little difference in fringe spacing between landfast ice seaward and shoreward of stabilizing anchor points. Without evaluating the phase response for each area of interest in detail during different forcing scenarios, it may be difficult to understand under what conditions the ice remains stable. Classification of stability based on relative differences in fringe density is also complicated by the use of nonsimultaneous interferograms to provide complete coverage of a region. The interferograms used here were obtained as close to maximum ice extent and sta-bility as possible (roughly late April), but once had to be obtained as early as February. Fringe density tends to decrease over the winter as the ice thickens (Dammann et al., 2016). Hence, the use of interferograms based on different dates can aid interpretation by confirming consistent fringe patterns and discontinuities that identify temporal changes. Temporal changes result in phase discontinuities at the image stitchings that are not related to different stability zones, which further complicates the mapping process. Sentinel-1 IW imagery is predominantly acquired over land, so it is likely not possible to construct interferograms away from the coast, for extensive landfast ice approaching the 250 km IW swath, such as that in the East Siberian Sea. Data availability further restricts the temporal baseline between images to a minimum of 12 days, though this now represents a shorter period than past work identifying landfast ice (Mahoney et al., 2004;Meyer et al., 2011;Dammann et al., 2016). Further studies should investigate the effect of different temporal baselines on the stability product. A shorter baseline will result in higher temporal resolution. However, with a shorter baseline (e.g., Sentinel-1 6-day baseline), mapping of the seaward landfast ice edge may incorporate stationary pack ice. A longer baseline will result in lower interferometric coherence. With a 12-day baseline, some regions, such as the Kotzebue Sound region, already feature consistent coherence loss. Such regions can most often be identified through a spatially inconsistent progression, from high to a complete loss of coherence. In such cases, the mapping of landfast ice type boundaries is not possible. It is worth mentioning that this technique can only be used before the onset of melt, when widespread coherence loss occurs. Therefore, it is not possible to evaluate the retreat of bottomfast ice or the reduction of ice stability in response to melt.

Conclusion
In a time of rapidly changing sea ice conditions and continued interest in the Arctic from a range of stakeholders, we stress the need for new assessment strategies to enable safe and efficient use of sea ice. InSAR is gaining growing attention in the sea ice scientific community, and here we demonstrate its value for identifying zones of landfast ice stability. We are also highlighting the application of InSAR for the development of operational sea ice information products, for both long-term strategic planning and short-term tactical decisions. Using interferograms generated by a standardized workflow, we show that three stability zones of landfast ice can be identified based on fringe density and continuity, indicative of differential ice motion occurring between SAR acquisitions. Along the Beaufort Sea coast of Alaska, we find that the landfast ice regime can be well described with three stability zones: bottomfast ice, where the sea ice is frozen to or resting on the seabed; stabilized ice, which is floating but sheltered by coastlines or anchored by islands or grounded ridges; and nonstabilized ice, which represent floating extensions seaward of any anchoring points. These findings are supported by comparison with the location of stable nodes, identified through analysis of hundreds of landfast ice edge positions over the period 1996-2008. Not only does this provide some validation of our results, but it demonstrates the ability of InSAR to capture useful information in just two snapshots, compared to previously requiring analysis over many years.
Based on our findings, it is likely that InSAR-derived maps could provide substantial value as a standalone product for some regions such as the Beaufort Sea. With that said, the stability zones in the Beaufort Sea and the Russian Arctic appear to be qualitatively different. This makes it challenging to directly adapt the proposed scheme to the East Siberian and Laptev seas, which are associated with substantial uncertainties. Even so, we demonstrate the data availability and application of this InSAR-based approach, which can provide added value to ice charts and other products. In ice charting, multiple information products are evaluated with local knowledge to create final products. Similarly, the value of In-SAR may be greatly enhanced by linking it with other products (e.g., InSAR time series analysis, SAR-based and optical remote-sensing products, local knowledge, coastal morphology and bathymetry, and atmospheric and ocean forcing data).
The use of a standardized workflow facilitates large-scale application of this approach, which we demonstrate on a near-pan-Arctic scale using 52 Sentinel-1 acquisition pairs during spring 2017. This has allowed us to map the same zones of landfast ice in the Beaufort, Chukchi, East Siberian, Laptev, and Kara seas. To our knowledge, our results represent the first mapping of bottomfast ice extent at this scale and the first attempt to map the extents of different landfast ice stability zones on any scale. It also enabled us to estimate and compare the total area covered by each stability zone in each marginal sea. However, we note that these comparisons are based on the assumptions that the landfast ice regimes in all these seas can be well described by the same three stability zones. Although the delineation of different zones can be subjective, in particular in the Russian Arctic, our results clearly show that not all landfast ice is equally stable. Here, InSAR is potentially able to detect small-scale motion up to hundreds of kilometers from the shore that have previously been overlooked. In addition, there are uncertainties associated with the exact mapping of stability zones -in particular in terms of the exact delineation between stabilized and nonstabilized ice in the East Siberian and Laptev seas. Here, the boundaries between stabilized and nonstabilized ice are more difficult to discriminate, likely due to fewer pinning points where the ice is grounded or supported. Therefore, what we present here is not an operational ice chart, but the ability and application to discriminate stability classes on a pan-Arctic scale using InSAR.
The method presented in this work has a broad set of potential applications for monitoring, including subsea permafrost, biological habitats both beneath and above the ice surface, and ice use by a range of stakeholders. Bottomfast ice is important because it helps the aggregation of subsea permafrost, which serves to constrain the location of permafrost-rich shorelines. Utilizing InSAR, it is likely possible to monitor changes in bottomfast ice over time, with significant implications for erosion and spring flooding and the release of methane hydrates. We argue that the utility of InSAR and its potential applications also extend to maritime activities and shipping. With regards to the latter, vessel traffic typically does not traverse landfast ice. However, the assessment of landfast ice stability and spatiotemporal extent can potentially aid the management of conflicting ice uses such as in the case of the access route to the Voisey's Bay Mine in the Canadian Arctic, which cuts through the landfast ice that is part of a traditional Nunatsiavummiut use area (Bell et al., 2014).
With respect to ice users, sea ice navigation near or through landfast sea ice is presently predominately supported by sea ice charts used to map areas occupied by landfast ice. However, these charts do not provide information about the relative stability of the ice. The information provided here would likely be useful in the context of navigation and support of on-ice operations. The InSAR-based approach de-scribed here can potentially provide support by identifying the following stability-related features: 1. low-stability ice that may break off and drift into shipping lanes, 2. grounded ridges that may be problematic for ice navigation, but at the same time may provide added stability for on-ice operations, 3. stable areas to use for equipment staging by coastal community hunters and industry, 4. bottomfast ice for development of ice roads for transportation of heavy loads.
We further demonstrate the scientific and operational value of InSAR over sea ice through the examination of interferograms of ice arches. In this context, they can be considered part of an additional stability zone of quasi-landfast ice (i.e., temporarily stabilized pack ice). Preliminary analysis of the Nares Strait ice arch in 2017 suggests that interferograms may reveal precursors to the failure of ice arches. We further speculate that InSAR tools can be developed to inform stakeholders of changing landfast ice stability and ice movement. Such applications would have potential value for an early warning system designed to alert ice users of hazards related to ice movement and breakout events. The use of inverse modeling may further help derive the small-scale strain field from interferograms, which may improve our ability to predict their failure. We expect that InSAR can provide valuable information for stakeholders, enabling the tracking of ice dynamics and stability on seasonal timescales. The ability to provide stability information to stakeholders also opens up for the development of operational guidelines in terms of what stability zones should be prioritized or avoided. This work builds on previous applications of InSAR to study landfast ice and demonstrates what can be achieved over large areas with a standardized workflow. The year 2017 was the first in which Sentinel-1 covered the Arctic coast with the IW images necessary for this analysis. If this coverage continues, there will be considerable opportunity for development beyond what is presented here, including development of automated methods for mapping and classifying landfast ice suitable for incorporation into operational ice charts. Furthermore, through additional analysis of landfast ice and ice arches subject to different forcing conditions, we anticipate improving our understanding of stabilizing and destabilizing mechanisms. This would allow the improved prediction of formation and breakup. Not only will this enhance the operational sea ice information available to stakeholders, but it also allows us to better understand the response of coastal sea ice to a changing Arctic environment. Data availability. Sentinel-1 data from this analysis can be obtained free of charge from the Copernicus Open Access Hub (https://scihub.copernicus.eu/, European Space Agency, 2019) or the Alaska Satellite Facility Vertex interface (https://vertex.daac.asf. alaska.edu/, Alaska Satellite Facility, 2018). See the Supplement for a full list of images used. The shapefiles for bottomfast, stabilized, and nonstabilized ice used for the final map shown in Fig. 9 are included in the Supplement.
Author contributions. DOD conducted the interferometric processing and analysis and drafted the initial manuscript. LEBE and ARM provided critical guidance on all aspects of the analysis and manuscript. HE provided expertise relevant to sea ice processes in different Arctic regions. FJM provided expertise relevant to the interferometric processing and analysis. All co-authors also provided valuable recommendations and corrections resulting in the final manuscript.
Competing interests. The authors declare that they have no conflict of interest.