Updated Understanding of the Ripley Landslide Kinematics Using Satellite InSAR

: The Thompson River valley hosts 14 landslides along a 10 km section, which threaten the two major railroads connecting the Port of Vancouver and the interior provinces in Canada. The Ripley landslide is one of the active landslides in this section of the valley. Previous research at this site included an analysis of landslide deformations using satellite radar interferometry focusing on deformations measured in the line of sight between the satellite and the slopes, and average downslope displacement (deformations projected in the average downslope direction). Since then, further stratigraphic interpretation has provided an enhanced understanding of the Ripley landslide. In this update, the new stratigraphic interpretation is supplemented with satellite InSAR data from May 2015 to May 2017 to enhance the current understanding of the landslide kinematics. The results indicate that the Ripley landslide has been moving at a rate between 2 and 82 mm per year, corresponding to a very slow to slow landslide. It is also observed that the movements tend to be near-horizontal on areas closer to the toe of the landslide, while the vertical component of deformation increases near the scarp of the landslide. This, together with the interpreted stratigraphy, indicates the kinematics corresponds to a compound landslide. This is consistent with interpreted landslide kinematics of older, more mature landslides in the area that have shown episodes of retrogression and suggests the possibility of a similar future behaviour of the Ripley landslide.


Introduction
Geohazards threaten many infrastructures, as well as lives, all around the world.Historically, there are many reported cases of loss of life and damage to infrastructure in Canada as a consequence of landslides [1][2][3][4].Since 1771 there have been hundreds of fatalities and damage to infrastructure caused by landslides in Canada, which have cost millions of dollars annually [5].A recent reminder of the harmful potential of landslides in Canada was the four fatalities confirmed in a recent landslide on British Columbia's Highway 99 in November 2021 (Simon Little, Global News, 20 November 2021).In addition, climate change and accelerated glacier melting and retreat can increase landslide activity in northern regions in Canada [3]; therefore, continued monitoring and understanding of landslide phenomena becomes important for resilience against climate change.
Identification of landslides and monitoring the terrain's movement is one of the most effective risk management approaches for areas susceptible to landslides, particularly when stabilization or avoidance are technically or economically challenging.Monitoring can help to adapt to these landslide phenomena and inform adequate maintenance scheduling of infrastructures.Monitoring can also provide robust early warning systems [6,7].Particularly, remote monitoring information can be used to increase our understanding of landslide kinematics and triggers, therefore enhancing landslide risk management strategies [8][9][10][11][12][13].This paper presents an updated understanding of the landslide kinematics using satellite InSAR for the Ripley landslide in the province of British Columbia, Canada.The Ripley landslide has been studied by a number of researchers [14][15][16][17].Previous research included an analysis of landslide deformations using satellite radar interferometry focusing on deformations measured in the line of sight between the satellite and the slopes and average downslope displacement (deformations projected in the average downslope direction) [16].Since then, further stratigraphic interpretation has provided an enhanced understanding of the geometry of the shear surfaces in the Ripley [6,14,[17][18][19].In this update, the new stratigraphic interpretation is supplemented with satellite InSAR data from May 2015 to May 2017 (Sentinel 1) to enhance the current understanding of the landslide kinematics.

Landslide Remote Sensing and Satellite InSAR
New technologies facilitate geotechnical engineers' real-time observation of a terrain's movement.Remote sensing, GPS monitoring, geophysical imaging and geotechnical investigation are commonly applied for landslides characterization [20].Remote sensing methods such as interferometric synthetic aperture radar (InSAR), light detection and ranging (LiDAR), and unmanned aerial vehicle (UAV) photogrammetry are applied to monitor cut slopes and natural slopes in current routine geohazard monitoring programs.These technologies help engineers to acquire data over large areas with high spatial resolution in comparison to in-place instruments such as slope inclinometers or survey monuments [7].
In this method, radar satellites carry SAR sensors that utilize wavelength information to calculate distance to the surface and that can be used for calculating relative change [24].
InSAR has become a reliable, accurate and near-real-time method to monitor potential landslide events and is particularly useful in remote areas where in-place ground monitoring instruments are technically or economically challenging.Although InSAR has been extensively used for monitoring landslides in the last three decades, particularly slow moving landslides, it is still challenging to find adequate data to generate InSAR calculations for rapid landslide acceleration.Other challenges include atmospheric noise, vegetation cover and seasonal effects, such as rainfall and snow cover [5].Notwithstanding these limitations, the ability of InSAR to provide ground deformation measurements with high accuracy makes this method an effective technique in many applications focused on deformation phenomena that occur on the Earth's surface [33].Satellite InSAR has been proven to be an effective method for identifying and monitoring ground displacement of slow to extremely slow moving slopes.InSAR information can be used for enhanced understanding of landslide mechanisms when complemented with other monitoring information such as precipitation and ground water fluctuation.The resolution of InSAR monitoring (sub cm/year) and its significant ability to capture massive areas make it very useful for widespread surface monitoring of landslides.It is important to note that satellite revisiting times (now between 1 day and 1 week) limit the temporal resolution of InSAR monitoring.Ground-based InSAR can provide adequate high-frequency data for forecasting the rapid acceleration that develops in a window of hours.No sites will have more frequent access than a one-day revisiting time.In this regard, the information is valuable for identifying the initiation of ground displacements, changes in displacement trends and seasonal displacements; however, it would not provide adequate warning for failure modes that lead to rapid acceleration followed by collapse that develop in just a few days as the variability of monitoring points could mask reliable identification of acceleration [34].

The Thompson River Valley Landslides South of Ashcroft and the Ripley Landslide
The Thompson River valley south of Ashcroft hosts one of the most important transportation corridors in Canada.Both the Canadian Pacific Railway (CP) and Canadian National Railway (CN) use this corridor to transport goods between the Port of Vancouver and the rest of Canada (Clague and Evans, n.d.).The Thompson River valley became a monitoring priority for Natural Resources Canada (NRCAN), the Geological Survey of Canada (GSC) and Transport Canada (TC) because of its economic importance and effects on transportation services in Canada and it is considered essential to understand and manage the risks associated with the landslides that threaten the corridor [34].
The Thompson River valley was formed by a series of glaciations and later exposed to river erosion and incision.The valley is infilled with a complex sequence of a wide range of deposits, from poorly sorted sand and gravel to rhythmically interbedded silt and clay.This corresponds to multiple glacial advance and retreat intervals in the Pleistocene [16,35,36].The extensive terrace that hosts several landslides was formed during postglacial times when the southward-flowing Thompson River down-cut 150 m of deposited sediments gradually along the valley [16,37].In these landslides, the existence of steep internal shears and basal through-going shear surfaces, as well as the presence of multiple graben and horst features, confirms the presence of a weak silt and clay layer that is considered responsible for the instability of many of these landslides [16,35,38].This weak layer contains presheared surfaces with residual shear strengths characterized by no cohesion and friction angles between 9 and 16 degrees, according to landslide back-analyses [14].
Fourteen landslides have been identified in the Thompson River valley, 12 of them traversed by one or both railway companies.The first reported landslide dates to the 1860s and periodic reactivation and slope movements have occurred since then.The location and extents of the landslides are shown in Figure 1 [14].
Landslides in this valley have volumes between 0.75 million cubic meters (Ripley landslide-Figure 1) and 15 million cubic meters (north slide-Figure 1), and their movement velocity range between 10 mm/year (extremely slow) and reported accelerated episodes of up to 50 mm/s historically (very rapid) [16].
The Ripley landslide is approximately 200 m long, 40 m high and 300 m wide.A front view of the landslide is shown in Figure 2 [39].The Ripley landslide was recognized by Charles Ripley in the 19th century by observation of an offset in a fence on the uphill side of CP railway tracks.Since then, the landslide has been inactive or moved slow enough that regular maintenance of track allows for continued operations of the CN and CP railways, which traverse this landslide.Shortly after construction of a section of track in 2005, the landslide showed some deformation and this movement was followed by the development of cracks and a scarp in 2007 [14].Since then, the railway track and ballast has required adjustments at higher frequencies than for other sections of track [17].The Ripley landslide is the only landslide currently being monitored at almost real-time, as it has become a field laboratory for novel investigation technology [6,[14][15][16][17][18]40].The other landslides in this area are being monitored with periodic slope inclinometer readings.
The Ripley landslide's estimated volume of 750,000 m 3 is relatively small in comparison to the other landslides in the Thompson River valley; however, it is one of the most active landslides with velocities between 25 and 180 mm/year (very slow, according to the classification by Cruden and Varnes [41].Landslides in this valley have volumes between 0.75 million cubic meters (Rip landslide-Figure 1) and 15 million cubic meters (north slide-Figure 1), and their mo ment velocity range between 10 mm/year (extremely slow) and reported accelerated e sodes of up to 50 mm/s historically (very rapid) [16].
The Ripley landslide is approximately 200 m long, 40 m high and 300 m wide.A fr slides in this area are being monitored with periodic slope inclinometer readings.
The Ripley landslide's estimated volume of 750,000 m 3 is relatively small in comparison to the other landslides in the Thompson River valley; however, it is one of the most active landslides with velocities between 25 and 180 mm/year (very slow, according to the classification by Cruden and Varnes [41].

Published Displacement Monitoring Data at the Ripley Landslide
The GPS system monitoring data used in this paper was acquired through an installation by CP and analyzed by the University of Alberta under the Canadian Railway Hazard Research Program.The system consists of three GPS monitoring stations on the Ripley landslide, near the railway tracks and one reference station located on stable bedrock outside the unstable area.GPS stations are Leica single-phase receivers in a differential GPS mode, able to detect 12.5 mm of cumulative ground movement with a variability of ±1 mm [15].Lateral and vertical displacements are processed by the Leica GeoMOS software and information on the system's latitude, longitude and vertical coordinates is provided at an hourly frequency.Hourly readings are used to calculate daily positioning of the GPS stations.The GPS displacement data in this paper corresponds to that published by Macciotta et al. [15] for the period between April 2008 and October 2014.(All three GPS started recording data in April 2008.GPS-1 and GPS-3 recorded data until 21 October 2014, while GPS-2 recorded data only until 29 April 2014.All three systems showed some data gaps associated with repairs and maintenance of the system.)These data are used to validate the landslide deformation trends (directionality) obtained by InSAR, assuming the landslide has not suffered changes in the direction of movement.Unfortunately, no displacement vector datasets coincide with the InSAR processing period in this paper.Positive values for horizontal movement show displacements in the downslope direction and negative vertical values show downward movement [15].Figure 3 shows the proximate locations of each installed GPS system.

Published Displacement Monitoring Data at the Ripley Landslide
The GPS system monitoring data used in this paper was acquired through an installation by CP and analyzed by the University of Alberta under the Canadian Railway Hazard Research Program.The system consists of three GPS monitoring stations on the Ripley landslide, near the railway tracks and one reference station located on stable bedrock outside the unstable area.GPS stations are Leica single-phase receivers in a differential GPS mode, able to detect 12.5 mm of cumulative ground movement with a variability of ±1 mm [15].Lateral and vertical displacements are processed by the Leica GeoMOS software and information on the system's latitude, longitude and vertical coordinates is provided at an hourly frequency.Hourly readings are used to calculate daily positioning of the GPS stations.The GPS displacement data in this paper corresponds to that published by Macciotta et al. [15] for the period between April 2008 and October 2014.(All three GPS started recording data in April 2008.GPS-1 and GPS-3 recorded data until 21 October 2014, while GPS-2 recorded data only until 29 April 2014.All three systems showed some data gaps associated with repairs and maintenance of the system.)These data are used to validate the landslide deformation trends (directionality) obtained by InSAR, assuming the landslide has not suffered changes in the direction of movement.Unfortunately, no displacement vector datasets coincide with the InSAR processing period in this paper.Positive values for horizontal movement show displacements in the downslope direction and negative vertical values show downward movement [15].Figure 3 shows the proximate locations of each installed GPS system.

Satellite InSAR Displacement Data
Satellite InSAR displacement data used in this paper were collected from Sentinel 1 and processed by TRE ALTAMIRA CLS Group Company.Processed information included displacements in the line-of-sight (LOS), which is the visual line between the satellite sensor and the points on the land being monitored between the end of 2014 to the Spring of 2018.The data with minimum gaps between May 2015 to May 2017 were selected to investigate the ground displacement of the Ripley landslide.The radar orbit is called ascending when a satellite moves from south to north and captures images, in our case, towards the east.The orbit is called descending when the satellite travels from north to south and captures images, in our case, to the west.There are two important angles for each satellite LOS to interpret the calculated displacements.These are the angle between the LOS and the vertical plane, or theta (θ), and the angle between each satellite's orbit and the geographic north, or delta (δ). Figure 4 shows the typical angle layout for both ascending and descending satellite orbits.

Satellite InSAR Displacement Data
Satellite InSAR displacement data used in this paper were collected from Sentinel correspond to non-point features that are multiple pixels in size.
Table 1 illustrates a summary of properties of the data that were collected from each Sentinel orbit.
Both PS and DS data were provided by TRE using their SqueeSAR method, which incorporates PSInSAR processing methods.Scatterers refer to data points that provide information for measuring the ground movement by LOS displacements.Permanent scatterers (PS) consist of some buildings, fences and other man-made structures and also some natural features such as rocks or exposed ground that are likely stable.Distributed scatterers (DS) correspond to large areas (up to hundreds of square meters) and were identified from exposed ground or fields.Although each DS presents an exact point for clarity of presentation, these points actually correspond to non-point features that are multiple pixels in size.
Table 1 illustrates a summary of properties of the data that were collected from each Sentinel orbit.Both PS and DS data were provided by TRE using their SqueeSAR method, which incorporates PSInSAR processing methods.

Sentinel 1
Sentinel 1 was launched by the European Space Agency (ESA).Satellite images used in this study that cover the Ripley landslide were collected between 3 November 2014 and 17 March 2018 for the ascending orbit, and between 6 November 2014 and 1 April 2018 for the descending orbit.LOS angles are shown in Table 2.In this paper, we propose a method for interpreting the geometry of the deformation based on the satellite InSAR, noting this should be complemented with a classical evaluation of the line-of-sight (LOS) deformations and velocities.
The LOS InSAR displacement is a projection of the true landslide surface displacement on its axis.The procedure followed in the paper aims at calculating a close approximation to these real displacements, although a few assumptions are required.For simplification, the calculated displacements are referred to as the "Real" displacement vectors (R), where "bold" represents a vector.Given the orbits of the satellites are sub-parallel to the north-south direction, resolving the north-south displacement component is not reliable.Therefore, a common assumption for landslide InSAR interpretation is to consider landslide motion in the downslope direction (azimuth and inclination) [16].The work presented in this paper does not assume the directionality of movement in the vertical plane; however, it assumes movement in the horizontal plane with an azimuth parallel to the slope azimuth in the area of measurement.Therefore, knowing the topography of the landslide provides the assumed unit vector of the horizontal component of R.This can be accomplished by calculating the unit vector of R in the horizontal plane, considering the vertical component as an unknown, and projecting R to the LOS unit vectors (based on the satellite geometry).Knowing the LOS displacement magnitude, the vertical component of R and its magnitude can be solved.
The unit vectors of LOS are denoted by S a for Sentinel ascending and S d for Sentinel descending.The projections of R on these directions would be the measured LOS displacements and are denoted by vectors P a and P d , respectively, each with scalar magnitudes M a and M d .Therefore, Equation (1) presents the expression for P i , where i represents the ascending or descending orbits.
The components of R can be denoted as (x,y,z), and the components of the unit vectors of the LOS as S i = (e, n, u).To find the LOS unit vectors, it is essential to use the angles between the LOS, the vertical plane and azimuth from Table 2 (Equations ( 2)-( 4)).
The coordinate system used in this paper corresponds to (east, north, up).The sign of u is positive and sign of n is negative for both orbits' LOS, while the e sign is negative for the ascending orbit and it is positive for the descending orbit.The unit vectors are then calculated as: S a = (−0.613,−0.123, 0.781) Replacing Equations ( 5) and ( 6) into Equation (1) and solving for M a and M d provides the measured movements in the LOS in terms of the components of R (Equations ( 7) and ( 8)).
The third equation to resolve the components of R follows the assumption that the total horizontal vector component of R is parallel to the slope's azimuth in the area of measurement (α) as follows: x y = tan α (9) Equations ( 7)-(9) provide a system of three equations, three unknowns for estimating R. Scatterer locations for the different satellites and for different orbits are not the same.Therefore, this required selecting areas of the landslide where information from two orbits is available, and averaging all scatterers for each orbit within that area to calculate representative LOS displacements (Figure 5).The measurements are also not synchronous between orientations.Data from Sentinel ascending are gathered from November 2014 to March 2018 every 25 days while reported data is from November 2014 to April 2018 in 25-day intervals for Sentinel descending.This required the selection of relatively long periods of time for calculating average R (annual basis), such that differences of a few days would not have a meaningful impact on the results.The slope displacement rate of the Ripley landslide had been reported to be up to 150 mm/year (or 0.4 mm/day), suggesting that asynchronous measurements by less than two weeks for a total period of analysis of two years could represent approximately a 2% error.

of 24
Figure 5. InSAR data coverage for both Sentinel ascending and descending orbits on the Ripley landslide and areas selected for analysis (containing LOS changes during the selected time frame).

Satellite InSAR Results at the Ripley Landslide
The LOS vectors in the horizontal plane for each area of analysis are shown in Figure 6.The magnitude and direction of LOS vectors for each area are located on the center of the area while it is representative of the average values for all captured points within the area.LOS directions for the same orbit are all parallel, given that LOS direction is defined by the orbit geometry.This figure also shows the direction of the average slope azimuth.It is observed that LOS magnitudes are large and of similar magnitude near the river, and significantly decrease near the back scarp of the landslide.Figure 7 shows the cumulative LOS deformations for each area for both Sentinel ascending and Sentinel descending orbits.Acceleration periods tend to start in the fall and deacceleration tends to start in the spring.These show clear episodes of acceleration and deceleration, consistent with previous observations [14][15][16].In this published study the acceleration corresponds to the drawdown effect of the Thompson River on the Ripley landslide.

Satellite InSAR Results at the Ripley Landslide
The LOS vectors in the horizontal plane for each area of analysis are shown in Figure 6.The magnitude and direction of LOS vectors for each area are located on the center of the area while it is representative of the average values for all captured points within the area.LOS directions for the same orbit are all parallel, given that LOS direction is defined by the orbit geometry.This figure also shows the direction of the average slope azimuth.It is observed that LOS magnitudes are large and of similar magnitude near the river, and significantly decrease near the back scarp of the landslide.Figure 7 shows the cumulative LOS deformations for each area for both Sentinel ascending and Sentinel descending orbits.Acceleration periods tend to start in the fall and deacceleration tends to start in the spring.These show clear episodes of acceleration and deceleration, consistent with previous observations [14][15][16].In this published study the acceleration corresponds to the drawdown effect of the Thompson River on the Ripley landslide.Figure 8 shows the plan view of the total horizontal component of the calculated R for the selected areas on the Ripley landslide.The satellite data used for calculating R correspond to the period between 26 May 2015 and 21 May 2017 for ascending orbit data and between 10 June 2015 and 12 May 2017 for descending orbit data to estimate annual displacement.The results are shown as average velocity of the real displacement vectors for each selected area in mm/year.Velocities are less than 3 mm/year behind the known back scarp and ranged between 12 and 53 mm/year at the back scarp within the landslide Figure 8 shows the plan view of the total horizontal component of the calculated R for the selected areas on the Ripley landslide.The satellite data used for calculating R correspond to the period between 26 May 2015 and 21 May 2017 for ascending orbit data and between 10 June 2015 and 12 May 2017 for descending orbit data to estimate annual displacement.The results are shown as average velocity of the real displacement vectors for each selected area in mm/year.Velocities are less than 3 mm/year behind the known back scarp and ranged between 12 and 53 mm/year at the back scarp within the landslide body, and up to 81 mm/year close to the river, which shows very slow to slow ground movement according to the velocity classification by [41,42].These results are in agreement with the reported results by Huntely et al. [43], both in magnitudes and directions, which report displacements between 68 and 82 mm/year for different locations within the Ripley landslide body with the maximum ground movement occurring in the middle of the toe slide [43].
12, 298 14 of 24 body, and up to 81 mm/year close to the river, which shows very slow to slow ground movement according to the velocity classification by [41,42].These results are in agreement with the reported results by Huntely et al. [43], both in magnitudes and directions, which report displacements between 68 and 82 mm/year for different locations within the Ripley landslide body with the maximum ground movement occurring in the middle of the toe slide [43].Figure 9 illustrates the ground movement velocities in section view on different cross-sections defined on the Ripley landslide.These sketched section views were inspired by previous research [14,18].GPS average annual velocity vectors are also shown in the cross-sections intersecting the location of the GPS units.Table 3 presents the magnitudes   9 illustrates the ground movement velocities in section view on different crosssections defined on the Ripley landslide.These sketched section views were inspired by previous research [14,18].GPS average annual velocity vectors are also shown in the cross-sections intersecting the location of the GPS units.Table 3 presents the magnitudes of R, azimuth and horizontal angles (angle between R and the horizontal plane).

Comparison with in-Place Displacement Monitoring Instrumentation
Published GPS displacement results at the Ripley landslide were used to evaluate the validity of the estimated R.
GPS data were selected between the installation day on 11 April 2008 and 11 April 2014.Annual results for the 6 years are averaged to get an estimated characteristic annual ground movement.In this regard, Macciotta et al. [15] showed that, in general, the acceleration and deceleration trend of the Ripley landslide measured by the GPS system was consistent for this period.The time series of GPS cumulative displacements in this period are presented as vertical and total horizontal components in Figure 10. Figure 11 shows each GPS horizontal displacement in plan-view together with the calculated R vectors.The summary of recorded displacement by each GPS is shown in Table 4.The GPS velocities in Table 4 are consistent with the calculated velocities for the R displacement vectors near the toe of the landslide.

Discussion
Figure 5 illustrated the scarcity of PS on the Ripley landslide from Sentinel ascending and descending orbits.Despite this low scatterer density, selected regions on the landslide did include scatterers for both satellite orientations so that an average R vector could be calculated.LOS results show faster deformations at the landslide's toe and slowing towards the upper areas of the landslide.The LOS displacement for both satellite orientations is significantly lower outside the landslide extents to the north and south.
The fastest velocity was calculated in region 5 with 82 mm/year.This region is located near railway tracks on the center of the Ripley landslide as it is shown in Figure 8.The Ripley landslide is a very slow to slow moving landslide based on the velocity classification by [41,42].Figure 7 indicates there are episodes in which the landslide is showing an accelerated rate of deformation as opposed to episodes in which the landslide moves slower.Although it is not as explicit as the seasonal trend of cumulative displacement of GPS data, which is shown in Figure 10, the cumulative LOS changes for both Sentinel ascending and descending satellites in Figure 7 introduces a faint trend in LOS changes that shows two acceleration seasons between October 2015 to May 2016 and November 2016 to the end of data, while they show deceleration between May 2016 and November 2016.This agrees with previous research that showed a correlation between ground movement and river water table fluctuation.Reductions in river flows (and drop in water head) has been associated with episodes on landslide acceleration [14,34].The seasonal movement is more obvious for areas near the river, which showed more movement than other areas near the back scarp on the landslide body and also outside landslide boundaries.The inability of data to perform a strong trend can be referred to the relatively longer time intervals between InSAR readings for both satellites and seasonal impacts on InSAR results (e.g., vegetation growth, weather changes and ground moisture).
Calculated R vectors are in general good agreement with displacement measurements from GPS-1 and GPS-2 (Figure 10).The velocity for R4 is 58 mm/year with an azimuth of 302 degrees (Table 3), while the average velocity and azimuth angle are reported as 67 mm/year and 293 degrees, respectively, for GPS-1 in Table 4.For R5, the azimuth angle is 307 degrees (Table 3), which is consistent with GPS-2 in Table 4 (298 degrees), and the recorded ground displacement velocity by GPS-2 is 74 mm/year, compared to 82 mm/year for R5.
However, GPS-3 shows a faster landslide velocity and a larger vertical component as opposed to other locations near the track.Although there are no InSAR regions coincidental with GPS-3, this unit shows different displacement vector characteristics than the InSAR regions in the vicinity.It has been interpreted that the vertical movement of GPS-3 may be exaggerated because it is installed on a retaining structure that was built as part of track construction in 2005, and the weight of the structure over materials loosened by landslide deformations could lead to increased ground settlement [39].
Figure 8 shows that the magnitude of landslide displacements near the railway tracks and towards the river is significantly larger than deformations upslope, closer to the back scarp.It is illustrated in the cross sections in Figure 9 that the magnitude of R vectors tends to decrease with distance from the river; however, the vertical components of ground movement are greater at locations near the back scarp.Displacements near the railway tracks and closer to the river are predominantly sub-horizontal while the vertical component of calculated R vectors near the back scarp slightly increases in proportion to the total R magnitude.The magnitude of R vector velocities in Table 3 is between 17 and 82 mm/year and the angle with vertical plane is from 75 to 86 degrees near the landslide's toe, while this angle is less than 60 degrees for the section near the back scarp.This would be consistent with the interpretation of the kinematics corresponding to a compound landslide, where a driving wedge in the upper section of the slope moves downslope pushing a passive wedge sliding over sub-horizontal shear surfaces.The vertical components of GPS displacements are also more notable than the vertical components of calculated R located on the landslide's toe.The differences between displacement direction in section-view for GPS and InSAR results can be assigned to their different distances from the landslide's back scarp.
Calculated R velocities were very small behind the back scarp of the landslide (between 2 and 5 mm/year).These values could suggest some initiation of movement behind the interpreted back scarp of the Ripley landslide; however, these small magnitudes could be within the limits of detection of the approach used, particularly when considering the assumptions required for the calculation of R. Any initiation of movement upslope from the known extents of the landslide could represent landslide retrogression, which is a common feature of most other landslides in the area.However, further monitoring would be required upslope from the known active area of the landslide to identify if retrogression of the landslide has initiated.

Conclusions
The Ripley landslide is a very slow moving landslide moving on a sub-horizontal weak clay layer.The landslide is located on an important transportation corridor in Canada along the Thompson River valley, and the two largest railway main lines traverse this landslide.This paper presents an update on the Ripley landslide displacement trends and kinematics using InSAR monitoring data.To obtain an improved understanding of InSAR displacement measurements, the landslide is divided into 12 areas that correspond to coverage from both Sentinel ascending and Sentinel descending orbits.A method is utilized where the horizontal component of movement is considered parallel to the slope azimuth to calculate the three-dimensional ground displacement vector.Ground displacement monitored with three GPS units on the landslide generally support the results from the method adopted.
Total InSAR LOS displacement velocities ranged between 2 mm/year and 82 mm/year, with faster sections near railroad tracks and the toe of the landslide, and slower sections at the back scarp of the landslide and upslope from the known active area.It is observed that landslide movement has a larger relative vertical component near the back scarp and becomes closer to horizontal direction near the railway tracks and near the toe.This would be consistent with the interpretation of the kinematics corresponding to a compound landslide, where a driving wedge in the upper section of the slope moves downslope pushing a passive wedge sliding over sub-horizontal shear surfaces.
Calculated ground velocities were very small behind the back scarp of the landslide (between 2 and 5 mm/year).These values could suggest some initiation of movement behind the interpreted back scarp of the Ripley landslide; however, these small magnitudes could be within the limits of detection of this approach, and further monitoring with in-situ instruments such as GPS, robotic total station and prisms to complement remote sensing techniques, e.g., satellite or ground-based InSAR, would be required upslope from the known active area of the landslide for identifying any potential landslide retrogression.

Figure 1 .
Figure 1.Thompson River valley location (a,b) plan-view and active landslides across the valley (c).Coordinates in UTM (units in meters for scale purposes).BC and AB in (a) correspond to th provinces of British Columbia and Alberta, respectively.

Figure 1 .
Figure 1.Thompson River valley location (a,b) plan-view and active landslides across the valley (c).Coordinates in UTM (units in meters for scale purposes).BC and AB in (a) correspond to the provinces of British Columbia and Alberta, respectively.

Figure 2 .
Figure 2. Ripley slide geometry including CN and CP railways and installed GPS on the site.

Figure 2 .
Figure 2. Ripley slide geometry including CN and CP railways and installed GPS on the site.

Figure 3 .
Figure 3. Location of installed GPS systems on Ripley landslide.

Figure 3 .
Figure 3. Location of installed GPS systems on Ripley landslide.

Figure 4 .
Figure 4. Geometry of the image acquisitions along the ascending (a) and descending (b) orbits (inspired by ground deformation InSAR analysis over the Thompson Canyon, British Columbia, Technical Details, August 2018, by TRE Group).

Figure 4 .
Figure 4. Geometry of the image acquisitions along the ascending (a) and descending (b) orbits (inspired by ground deformation InSAR analysis over the Thompson Canyon, British Columbia, Technical Details, August 2018, by TRE Group).

Figure 5 .
Figure 5. InSAR data coverage for both Sentinel ascending and descending orbits on the Ripley landslide and areas selected for analysis (containing LOS changes during the selected time frame).

Figure 6 .
Figure 6.LOS vectors projected on the horizontal plane for both Sentinel ascending and Sentinel descending orbits.The average direction of the slope azimuth is also shown.

Figure 6 .
Figure 6.LOS vectors projected on the horizontal plane for both Sentinel ascending and Sentinel descending orbits.The average direction of the slope azimuth is also shown.

Figure 8 .
Figure 8. Calculated total horizontal component of deformation (R) in plan-view on Ripley landslide.

Figure 8 .
Figure 8. Calculated total horizontal component of deformation (R) in plan-view on Ripley landslide.

Figure
Figure9illustrates the ground movement velocities in section view on different crosssections defined on the Ripley landslide.These sketched section views were inspired by previous research[14,18].GPS average annual velocity vectors are also shown in the cross-sections intersecting the location of the GPS units.Table3presents the magnitudes of R, azimuth and horizontal angles (angle between R and the horizontal plane).

Figure 9 .
Figure 9. Real vectors (ground movement) in section view on Ripley landslide (sliding surfaces are drawn based on the position of tension tracks and interpreted shear surface by Hendry et al. [14] based on previous BH and SAA installation).

Figure 9 .
Figure 9. Real vectors (ground movement) in section view on Ripley landslide (sliding surfaces are drawn based on the position of tension tracks and interpreted shear surface by Hendry et al. [14] based on previous BH and SAA installation).

Figure 10 .
Figure 10.Cumulative horizontal and vertical displacements of installed GPS systems on Ripley landslide.

Figure 11 .
Figure 11.R vectors and GPS average velocity vectors (mm/y) in plan-view.Figure 11.R vectors and GPS average velocity vectors (mm/y) in plan-view.

Figure 11 .
Figure 11.R vectors and GPS average velocity vectors (mm/y) in plan-view.Figure 11.R vectors and GPS average velocity vectors (mm/y) in plan-view.

Table 1 .
Details of the processed satellite InSAR data.

Table 2 .
Satellite viewing (LOS) angles for the Sentinel and Radarsat-2 imagery.

Table 3 .
Magnitudes of R and direction.

Table 3 .
Magnitudes of R and direction.

Table 4 .
Movement Records of Installed GPS on Ripley Slide.