Sentinel-1 and RADARSAT Constellation Mission InSAR Assessment of Slope Movements in the Southern Interior of British Columbia, Canada

: Landslides are the most common natural hazard in British Columbia. The province has recorded the largest number of historical landslide fatalities in Canada, and damage to infrastructure comes at a great cost. In order to understand the potential impacts of landslides, radar remote sensing has become a cost-effective method for detecting downslope movements. This study investigates downslope movements in the Southern Interior of British Columbia, Canada, with Sentinel-1 and RADARSAT Constellation Mission (RCM) interferometric synthetic aperture radar (InSAR) data. The 2-dimensional time-series analysis with Sentinel-1 ascending and descending InSAR pairs from October 2017 to June 2021 observed distinct earthﬂow movements of up to ~15 cm/year in the east– west direction. The Grinder Creek, Red Mountain, Yalakom River, and Retaskit Creek earthﬂows previously documented are still active, with east–west movements of ~30 cm over the past four years. New RCM data acquired from June 2020 to September 2020 with a 4-day revisit capability were compared to 12-day Sentinel-1 InSAR pairs. The 4-day RCM InSAR pairs at higher spatial resolution showed better performance by detecting relatively small-sized slope movements within a few hundred meters, which were not clearly observed by Sentinel-1. The temporal variabilities observed from the RCM InSAR showed great potential for observing detailed slope movements within a narrower time window.


Introduction
Historically, landslides have had catastrophic impacts, in Canada, in terms of loss of life and damage to infrastructure. Since 1771, there have been at least 786 fatalities due to landslides [1]. Damage to infrastructure has cost millions of Canadian dollars annually in repair and mitigation [2,3]. Furthermore, Porter et al. (2019) estimated that, in the Western Canadian Basin alone, costs related to landslide damage or prevention exceed CAD 281-450 million per year [4]. In British Columbia (BC)'s mountainous regions, landslides pose an ever-present risk [1, 5,6]. Since 1880, there have been at least 356 landslide fatalities in BC caused by a variety of landslide types, such as rockfalls, rockslides, rock avalanches, debris flows, etc. [1]. Recently, in November 2020, a massive landslide occurred along the steep mountainous slopes over Elliot Creek in BC, which triggered a tsunami in a glacial lake and a subsequent outburst flood carrying~7.7 million cubic meters of mud and rock debris into Bute Inlet, reshaping the coastal landscape [7]. In January 2021, mud and rocks slid into Harrison Lake from a steep slope of Mount Breakenridge in BC over two days. A bathymetric survey of Harrison Lake reported that~4 million cubic meters of landslide debris were deposited on the lake floor, which implies the potential The study area of roughly 50,000 km 2 within British Columbia's Southern Interior (Figure 1) belongs to the northwest-southeast trending morphogeological belts of the Coast and Intermontane belts [29]. The Coast belt is mainly formed of granitic plutons ranging from the Jurassic to Cenozoic in age. Other rock types are volcanic and sedimentary in origin, dating from the late Proterozoic, to the Paleozoic and Holocene. Some were highly metamorphosed from mid-Cretaceous to early Ceneozoic time. The Intermontane rocks consist of a variety of volcanic, sedimentary, and granitic rocks from the Devonian to the Jurassic, to early Cenozoic [29].
The Fraser Glaciation in southern British Columbia began in Late-Wisconsinan time, about 29,000 years BP, culminated about 14,000 years BP, with deglaciation almost complete around 11,500 years BP [30]. The highest elevation in the study area is roughly 2600 m and the lowest is about 350 m (Figure 1).
The study area encompasses at least three biogeoclimatic zones: Montane spruce, Engelman Spruce-Subalpine Fir, Interior Douglas Fir [31]. The Montane Spruce zone is a narrow mid-elevation band (300-400 m) of vegetation. However, it can be found at high-elevation plateaus between 1100-1650 m. It reflects a cool dry climate with average temperatures above 10°C for 2-4 months a year and below 0°C for 5 months a year. Precipitation varies between 300 and 900 mm per year. The Engelman Spruce-Subalpine Fir zone occupies the highest elevations in British Columbia. It reflects cold long snowy winters and short cool summers. Average annual temperatures range from −2 to 2°C, and average annual precipitation is highly variable from 400-500 mm in drier areas up to 2200 mm in wetter areas. The Interior Douglas Fir zone occupies low-mid elevations where summers are short, warm and dry with mild winters. There is often a shortage of moisture during the growing season. Average annual temperatures vary from 1.6 to 9.5°C and annual precipitation ranges from 300 to 750 mm, but can exceed 1000 mm in wetter areas [31]. An historical photogrammetric survey of earthflows within the study area by Bovis and Jones (1992) reported that these were "glacier-like masses of soil and rock, typically 0.5-1 km wide and 2-6 km long, which move at velocities of 0.5 m/year" (p. 1746) [32]. The study revealed, through dendrochronology and tephrochronology, that many responded to the same hydroclimatic changes over time. Reactivation of the earthflows reflected cooler and wetter conditions, which produced an increase in groundwater recharge and pore water pressure at different times during the Holocene [32]. Some of these earthflows are further discussed in the Results (Section 4) and Discussion (Section 5).

Sentinel-1 and RCM Data
A total of 102 ascending (path: 64, frame: 163, April 2017-June 2021) and 91 descending (path: 13, frame: 421, October 2017-June 2021) Sentinel-1 Terrain Observation with Progressive Scan (TOPS) SAR data of VV polarization were processed (Table 1 and Figure  2). The Sentinel-1 TOPS SAR single look complex (SLC) data were acquired at a spatial resolution of 2.3 m by 13.9 m in range and azimuth, respectively, at incidence angles of ~30-45˚. The Sentinel-1 TOPS SAR interferometric wide (IW) mode data, with a ~250 km wide swath, consist of 27 and 28 burst images for ascending and descending, respectively. They were acquired with a 12-day repeat cycle (cf. 6-day revisit available with the Sentinel-1A and Sentinel-1B constellation) [33]. A total of 8 descending RCM 5M18 mode SLC images and 14 descending RCM 5M23 mode SLC images of HH polarization acquired from June 2020 to September 2020 were processed. The RCM SLC data were acquired at a spatial resolution of 2.7 m by 2.5 m in range and azimuth at incidence angles of 47˚ (5M18 mode) and 53˚ (5M23 mode), respectively. The 5 m resolution RCM data provide a ~30 km swath and a repeat cycle up to 4-days. An historical photogrammetric survey of earthflows within the study area by Bovis and Jones (1992) reported that these were "glacier-like masses of soil and rock, typically 0.5-1 km wide and 2-6 km long, which move at velocities of 0.5 m/year" (p. 1746) [32]. The study revealed, through dendrochronology and tephrochronology, that many responded to the same hydroclimatic changes over time. Reactivation of the earthflows reflected cooler and wetter conditions, which produced an increase in groundwater recharge and pore water pressure at different times during the Holocene [32]. Some of these earthflows are further discussed in the Results (Section 4) and Discussion (Section 5).

Sentinel-1 and RCM Data
A total of 102 ascending (path: 64, frame: 163, April 2017-June 2021) and 91 descending (path: 13, frame: 421, October 2017-June 2021) Sentinel-1 Terrain Observation with Progressive Scan (TOPS) SAR data of VV polarization were processed (Table 1 and Figure 2). The Sentinel-1 TOPS SAR single look complex (SLC) data were acquired at a spatial resolution of 2.3 m by 13.9 m in range and azimuth, respectively, at incidence angles of~30-45 • . The Sentinel-1 TOPS SAR interferometric wide (IW) mode data, with a~250 km wide swath, consist of 27 and 28 burst images for ascending and descending, respectively. They were acquired with a 12-day repeat cycle (cf. 6-day revisit available with the Sentinel-1A and Sentinel-1B constellation) [33]. A total of 8 descending RCM 5M18 mode SLC images and 14 descending RCM 5M23 mode SLC images of HH polarization acquired from June 2020 to September 2020 were processed. The RCM SLC data were acquired at a spatial resolution of 2.7 m by 2.5 m in range and azimuth at incidence angles of 47 • (5M18 mode) and 53 • (5M23 mode), respectively. The 5 m resolution RCM data provide a~30 km swath and a repeat cycle up to 4-days.  . Perpendicular baselines of Sentinel-1 (S1 ASC: ascending, S1 DES: descending) and RCM (RCM 5M18: descending 5M18 mode, RCM 5M23: descending 5M23 mode) InSAR pairs (red: Sentinel-1 ascending, green: Sentinel-1 descending, blue: RCM descending 5M18 mode, yellow: RCM descending 5M23 mode).

InSAR Processing
InSAR measures the phase difference in line-of-sight (LOS) between two SAR acquisitions. The phase difference modulated in 2π can be unwrapped into the displacements at a cm-scale wavelength. The InSAR processing was performed with an automated In-SAR processor based on the GAMMA software of the NRCan's Earth Observation Data Management System (EODMS) [34,35]. Based on small perpendicular and temporal baselines, multi-master InSAR pairs were selected and precisely co-registered for each pair. The interferograms were generated by multilooking, with a window of 4 by 1 pixels in range and azimuth (i.e., corresponding to 9.2 m by 13.9 m) for Sentinel-1 and a window of 2 by 4 pixels in range and azimuth (i.e., corresponding to 5.4 m by 10 m) for RCM. Note that a window of 2 by 1 pixels in range and azimuth (i.e., corresponding to 4.6 m by 13.9 m) was applied for Sentinel-1 processing for the comparison with RCM acquisitions from June 2020 to September 2020 to minimize the difference of the range resolution. The topographic phase was removed with the 30 m resolution Canadian Digital Elevation Model (CDEM). The multilooked DInSAR interferograms were filtered by the Goldstein adaptive filter [36] and unwrapped by the minimum cost flow (MCF) algorithm [37]. One-dimensional (1D, LOS) and two-dimensional (2D, east-west and vertical) time-series analyses were performed by applying the multidimensional small baseline subset (MSBAS) technique [38][39][40]. The 2D east-west and vertical movements were reconstructed by using the following Equation (1), Figure 2. Perpendicular baselines of Sentinel-1 (S1 ASC: ascending, S1 DES: descending) and RCM (RCM 5M18: descending 5M18 mode, RCM 5M23: descending 5M23 mode) InSAR pairs (red: Sentinel-1 ascending, green: Sentinel-1 descending, blue: RCM descending 5M18 mode, yellow: RCM descending 5M23 mode).

InSAR Processing
InSAR measures the phase difference in line-of-sight (LOS) between two SAR acquisitions. The phase difference modulated in 2π can be unwrapped into the displacements at a cm-scale wavelength. The InSAR processing was performed with an automated InSAR processor based on the GAMMA software of the NRCan's Earth Observation Data Management System (EODMS) [34,35]. Based on small perpendicular and temporal baselines, multi-master InSAR pairs were selected and precisely co-registered for each pair. The interferograms were generated by multilooking, with a window of 4 by 1 pixels in range and azimuth (i.e., corresponding to 9.2 m by 13.9 m) for Sentinel-1 and a window of 2 by 4 pixels in range and azimuth (i.e., corresponding to 5.4 m by 10 m) for RCM. Note that a window of 2 by 1 pixels in range and azimuth (i.e., corresponding to 4.6 m by 13.9 m) was applied for Sentinel-1 processing for the comparison with RCM acquisitions from June 2020 to September 2020 to minimize the difference of the range resolution. The topographic phase was removed with the 30 m resolution Canadian Digital Elevation Model (CDEM). The multilooked DInSAR interferograms were filtered by the Goldstein adaptive filter [36] and unwrapped by the minimum cost flow (MCF) algorithm [37]. One-dimensional (1D, LOS) and two-dimensional (2D, east-west and vertical) time-series analyses were performed by applying the multidimensional small baseline subset (MSBAS) technique [38][39][40]. The 2D east-west and vertical movements were reconstructed by using the following Equation (1), Remote Sens. 2021, 13, 3999 A is a matrix constructed from the time intervals of ascending and descending SAR acquisitions, ϕ is the azimuth angle, θ is the incidence angle, and λ and L are a regularization parameter and the Tikhonov regularization matrix, respectively [38][39][40]. V E and V V are vectors of the unknown east-west and vertical velocities (cf. the north-south component, V N , is neglected as V N V E, V ),Φ is a vector of the observed DInSAR deformation values, andÂ + is the pseudo-inverse of matrixÂ [38][39][40]. Figure 3 shows the linear deformation rate maps of east-west and vertical movements, which were reconstructed from Sentinel-1 ascending and descending time-series InSAR pairs. Slope movements are distinctly characterized by east-west components, which range from~-12 cm/year to~15 cm/year (i.e., positive values mean eastward movements). The east-west movements correspond to the downward directions of slopes and are dominantly observed from shallow to moderate slopes in valleys, including river channels. Vertical movements are relatively small and not as distinct as the east-west movements, as shown in Sites 1, 2, and 9 ( Figure 3b). Vertical movements are rarely observed from slopes shallower than~10 • . Most of the significant vertical movements are observed around mountain tops above~1700 m, which indicate downward movements of~-10 cm/year.  Figure 3). Additional slope movements other than the earthflows documented by [32,41] are observed on the east-west movement map (Figure 3a), which include Sites 1, 3, 6, 8, and 9 showing east-west movements of up to~50 cm over four years (Figures 4-6). Figure 7 shows the RCM interferograms of Sites 5, 6, and 9, compared to the Sentinel-1 interferograms. More detailed fringe patterns from small slope movements between tens of meters (e.g., fringes close to P1 and P1 in Figure 7a and P3 in Figure 7c) and~300 m in length are better captured by the RCM, which are not clearly differentiated from noise with Sentinel-1. In particular, the RCM captured a full DInSAR fringe from very small slope movements within~300 m in length in Site 9 (Figure 7e). The time-series profile analysis from June 30 to September 10, 2020, shows a similar range of displacements in both Sentinel-1 and RCM (Figures 8-10), but the RCM provided better spatial and temporal resolutions. In particular, significant movements are observed between mid-July and mid-August, based on the RCM results (e.g., 12 July-13 August in Figures 8 and 9, 15 July-16 August in Figure 10). For slope movements within~300 m in length in Site 9, the Sentinel-1 results (~1 cm) are significantly fewer than the RCM results (~2 cm) ( Figure 10). the Sentinel-1 results (~1 cm) are significantly fewer than the RCM results (~2 cm) ( Figure  10).

Discussion
The extensive Sentinel-1 TOPS SAR global coverage archives provide a great advantage for monitoring time-series slope movements from a regional perspective. Eastwest and vertical movements are reconstructed from ascending and descending LOS displacements. LOS measurements are insensitive to north-south movements due to the imaging geometry of near-polar orbiting sensors [38]. The observed slope movements in the area of interest are slow-moving earthflows and are better measurable in the east-west direction (< ~15 cm/year). The east-west movement rates estimated from the Sentinel-1 DInSAR correspond to the lower range of the earthflow movement rates at ~0.1 to 5 m/year documented in [32]. Vertical movements are small, relative to east-west movements and rarely observed from shallow slopes. Earthflows at shallow slope gradients show little vertical movements, though they show significant east-west movements (see Sites 1, 2, 9 in Figures 3a and 3b). The significant vertical movements around the mountain tops are observed only during winter months (i.e., ~December-March), which are likely to be related to ice and/or snow dynamics at high elevations.

Discussion
The extensive Sentinel-1 TOPS SAR global coverage archives provide a great advantage for monitoring time-series slope movements from a regional perspective. East-west and vertical movements are reconstructed from ascending and descending LOS displacements. LOS measurements are insensitive to north-south movements due to the imaging geometry of near-polar orbiting sensors [38]. The observed slope movements in the area of interest are slow-moving earthflows and are better measurable in the east-west direction (<~15 cm/year). The east-west movement rates estimated from the Sentinel-1 DInSAR correspond to the lower range of the earthflow movement rates at~0.1 to 5 m/year documented in [32]. Vertical movements are small, relative to east-west movements and rarely observed from shallow slopes. Earthflows at shallow slope gradients show little vertical movements, though they show significant east-west movements (see Sites 1, 2, 9 in Figure 3a,b). The significant vertical movements around the mountain tops are observed only during winter months (i.e.,~December-March), which are likely to be related to ice and/or snow dynamics at high elevations.  Figure 5a represents the documented earthflow at Red Mountain and Figure 5b represents the largest documented earthflow at Yalakom River [32].
Overall, the LOS displacements estimated from Sentinel-1 and RCM show similar patterns and amounts, but more detailed spatial and temporal variabilities observed from the RCM are not achieved by Sentinel-1. This is because the azimuth resolution of Sentinel-1 (13.9 m) is much coarser than that of the RCM (2.5 m), while the range resolution of Sentinel-1 (2.3 m) is comparable to that of the RCM (2.7 m). Though the range spacing of Sentinel-1 is achieved close to that of the RCM during the processing, the details of slope movements, particularly from earthflows within ~300 m in length, can be lost due to the azimuth spacing of more than 5 times, possibly resulting in saturation or underestimation of displacement measurements. There are differences in polarization and incidence angles between the Sentinel-1 and RCM acquisitions, but these have little impact on displacement measurements, given the similarities in the pattern and amounts of displacements observed from large-sized slope movements over a few hundred meters (~1 km 2 ).
Given the limited number of 4-day revisit RCM acquisitions and the slow movement rates of < ~15 cm/year in this region, it is not enough to assess the performance of 4-day repeat pass InSAR for monitoring slope movements. However, temporal variabilities observed from RCM 4-day (e.g., between 2020/09/02, 2020/09/06, and 2020/09/10, in Figure  8a) and 8-day (e.g., between 2020/08/13 and 2020/08/21, in Figure 8c) repeat passes show that the RCM 4-day repeat pass InSAR has great potential for observing more detailed slope movements with a better chance to obtain coherence, though the maximum detectable deformation gradient is limited to ~1.4 ꓫ 10 -3 (i.e., for a multilooked range pixel of ~10 m) [42].  Figure 5a represents the documented earthflow at Red Mountain and Figure 5b represents the largest documented earthflow at Yalakom River [32].
Overall, the LOS displacements estimated from Sentinel-1 and RCM show similar patterns and amounts, but more detailed spatial and temporal variabilities observed from the RCM are not achieved by Sentinel-1. This is because the azimuth resolution of Sentinel-1 (13.9 m) is much coarser than that of the RCM (2.5 m), while the range resolution of Sentinel-1 (2.3 m) is comparable to that of the RCM (2.7 m). Though the range spacing of Sentinel-1 is achieved close to that of the RCM during the processing, the details of slope movements, particularly from earthflows within~300 m in length, can be lost due to the azimuth spacing of more than 5 times, possibly resulting in saturation or underestimation of displacement measurements. There are differences in polarization and incidence angles between the Sentinel-1 and RCM acquisitions, but these have little impact on displacement measurements, given the similarities in the pattern and amounts of displacements observed from large-sized slope movements over a few hundred meters (~1 km 2 ).
Given the limited number of 4-day revisit RCM acquisitions and the slow movement rates of <~15 cm/year in this region, it is not enough to assess the performance of 4-day repeat pass InSAR for monitoring slope movements. However, temporal variabilities observed from RCM 4-day (e.g., between 2 September 2020, 6 September 2020, and 10 September 2020, in Figure 8a) and 8-day (e.g., between 13 August 2020 and 21 August 2020, in Figure 8c) repeat passes show that the RCM 4-day repeat pass InSAR has great potential for observing more detailed slope movements with a better chance to obtain coherence, though the maximum detectable deformation gradient is limited to~1.4 × 10 −3 (i.e., for a multilooked range pixel of~10 m) [42]. Figure 6. Close-ups of Sites 7-9 (a-c, respectively) of Google Earth imagery (left) and the Sentinel-1 east-west movement map overlaid on the hillshade DEM (middle) with Sentinel-1 2D time-series plots (right; cumulative east-west (blue) and vertical (purple) movements extracted from 5-by-5 pixels marked in black squares in the middle panels). The color-coded contours and dashed black arrows in the left panels represent the east-west movement rates estimated by Sentinel-1 and the direction of earthflows, respectively. The error bars in the time-series plots represent the standard deviations of the 5by-5 pixel measurements. Figure 7a represents the documented earthflow at Retaskit Creek [41].

Implications for Landslide Hazards and Risks
Both the Sentinel-1 and RCM DInSAR results demonstrate that certain types of landslides in the Southern Interior of British Columbia, especially earthflows, can be detected. Cumulative displacement amounts and displacement rates over a long period of time can be measured with time-series InSAR pairs. At a first glance, earthflows may not represent a high hazard, or risk, but slow and steady earthflows should be taken into consideration when building infrastructure, such as pipelines, roads, power lines, mines, etc. For example, Figure 11 shows close-up views of Google Earth imagery and the east-west movement rate of the largest earthflow (~2.4 km by 0.35 km) and other surrounding smaller earthflows at Yalakom River (i.e., Site 5 shown in Figures 3, 5b and 7a). The Sentinel-1 time-series results show that the average of the cumulative east-west movements is estimated at ~20 cm over four years for the largest earth flows at Yalakom River (see Figure  5b). This particular earthflow, triggered in weathered bedrock consisting of harzburgite, dunite, orthopyroxenite, and serpentinite mélange from the Shulaps Ultramafic complex [43], seems to have diverted the drainage of Yalakom River. This is slowly moving towards the road, causing riverbank erosion on the east side of the river channel. Hence, in this case, mitigation measures will have to be considered for future development projects. Figure 6. Close-ups of Sites 7-9 (a-c, respectively) of Google Earth imagery (left) and the Sentinel-1 east-west movement map overlaid on the hillshade DEM (middle) with Sentinel-1 2D time-series plots (right; cumulative east-west (blue) and vertical (purple) movements extracted from 5-by-5 pixels marked in black squares in the middle panels). The color-coded contours and dashed black arrows in the left panels represent the east-west movement rates estimated by Sentinel-1 and the direction of earthflows, respectively. The error bars in the time-series plots represent the standard deviations of the 5-by-5 pixel measurements. Figure 7a represents the documented earthflow at Retaskit Creek [41].

Implications for Landslide Hazards and Risks
Both the Sentinel-1 and RCM DInSAR results demonstrate that certain types of landslides in the Southern Interior of British Columbia, especially earthflows, can be detected. Cumulative displacement amounts and displacement rates over a long period of time can be measured with time-series InSAR pairs. At a first glance, earthflows may not represent a high hazard, or risk, but slow and steady earthflows should be taken into consideration when building infrastructure, such as pipelines, roads, power lines, mines, etc. For example, Figure 11 shows close-up views of Google Earth imagery and the east-west movement rate of the largest earthflow (~2.4 km by 0.35 km) and other surrounding smaller earthflows at Yalakom River (i.e., Site 5 shown in Figures 3, 5b and 7a). The Sentinel-1 time-series results show that the average of the cumulative east-west movements is estimated at~20 cm over four years for the largest earth flows at Yalakom River (see Figure 5b). This particular earthflow, triggered in weathered bedrock consisting of harzburgite, dunite, orthopyroxenite, and serpentinite mélange from the Shulaps Ultramafic complex [43], seems to have diverted the drainage of Yalakom River. This is slowly moving towards the road, causing riverbank erosion on the east side of the river channel. Hence, in this case, mitigation measures will have to be considered for future development projects. Remote Sens. 2021, 13, x FOR PEER REVIEW 10 of 17      . Earthflows at Yalakom River on Google Earth imagery (a) and its close-up view of the toe of the largest earthflow (b). The color-coded contours are the east-west movement rates estimated by Sentinel-1. The white rectangle represents the extent of Figure 11(b). The largest earthflow (solid white arrow) and surrounding smaller earthflows (dashed white arrows) are also represented in Figures 3, 5b, 7a, and 8. Sediments at the toe of the largest earthflow have diverted the river channel, which is causing riverbank erosion on the east side of the river channel (solid red arrows), next to the road.

Conclusions
Slow-moving earthflows in the southern interior of British Columbia were observed by the Sentinel-1 and RCM DInSAR analyses. The Sentinel-1 TOPS SAR provides both ascending and descending time-series data collected from 2017 with a broad coverage of ~250 km over this area. In addition to previously documented earthflows, numerous other slow-moving landslides were detected by the Sentinel-1 2D MSBAS analysis spanning four years. Earthflows are moving up to ~15 cm/year in the east-west direction, which correspond to the lower range of movement rates of documented earthflows. The Grinder Creek, Red Mountain, Yalakom River, and Retaskit Creek earthflows are very active, with Figure 11. Earthflows at Yalakom River on Google Earth imagery (a) and its close-up view of the toe of the largest earthflow (b). The color-coded contours are the east-west movement rates estimated by Sentinel-1. The white rectangle represents the extent of Figure 11b. The largest earthflow (solid white arrow) and surrounding smaller earthflows (dashed white arrows) are also represented in Figure 3, Figure 5b, Figure 7a, and Figure 8. Sediments at the toe of the largest earthflow have diverted the river channel, which is causing riverbank erosion on the east side of the river channel (solid red arrows), next to the road.

Conclusions
Slow-moving earthflows in the southern interior of British Columbia were observed by the Sentinel-1 and RCM DInSAR analyses. The Sentinel-1 TOPS SAR provides both ascending and descending time-series data collected from 2017 with a broad coverage of 250 km over this area. In addition to previously documented earthflows, numerous other slow-moving landslides were detected by the Sentinel-1 2D MSBAS analysis spanning four years. Earthflows are moving up to~15 cm/year in the east-west direction, which correspond to the lower range of movement rates of documented earthflows. The Grinder Creek, Red Mountain, Yalakom River, and Retaskit Creek earthflows are very active, with east-west movements of~30 cm over the past four years. The RCM provides much higher spatial resolution (i.e., 1 m spotlight, 3 m very-high-resolution and 5 m high-resolution modes) with a 4-day revisit capability, though the swath is limited to <~30 km with the high-resolution modes. Thus, it is recommended to initially observe slope movements at a regional scale by applying 2D time-series analysis with Sentinel-1 and then focus on slope movements of interest at high spatial and temporal resolution with the RCM. From a landslide hazard and risk perspective, slow-moving earthflows should be taken into consideration in development projects, as they can divert river drainage, which can, in turn, cause erosion of riverbanks and can affect infrastructure, as observed at Yalakom River.