Remote Sensing Applications for Landslide Monitoring and Investigation in Western Canada

: Transportation infrastructure in mountainous terrain and through river valleys is exposed to a variety of landslide phenomena. This is particularly the case for highway and railway corridors in Western Canada that connect towns and industries through prairie valleys and the Canadian cordillera. The ﬂuidity of these corridors is important for the economy of the country and the safety of workers, and users of this infrastructure is paramount. Stabilization of all active slopes is ﬁnancially challenging given the extensive area where landslides are a possibility, and monitoring and minimization of slope failure consequences becomes an attractive risk management strategy. In this regard, remote sensing techniques provide a means for enhancing the monitoring toolbox of the geotechnical engineer. This includes an improved identiﬁcation of active landslides in large areas, robust complement to in-place instrumentation for enhanced landslide investigation, and an improved deﬁnition of landslide extents and deformation mechanisms. This paper builds upon the extensive literature on the application of remote sensing techniques and discusses practical insights gained from a suite of case studies from the authors’ experience in Western Canada. The review of the case studies presents a variety of landslide mechanisms and remote sensing technologies. The aim of the paper is to transfer some of the insights gained through these case studies to the reader.


Introduction
The integrity of the transportation network in Western Canada for securing uninterrupted fluidity of goods and people through the network are essential for the economy of the region. Transportation corridors through the Canadian Cordillera and through river valleys in the prairies are exposed to natural slope instabilities and steep cuts required to accommodate railroads and highways, where slope-related hazards are common. Landslide types in these regions include rock falls, rock slides and topples, earth flows, deep seated rotational and translational soil landslides, and compound landslides [1][2][3][4][5][6]. In this regard, a robust approach to manage the risks associated with landslides in natural and cut slopes becomes essential.
The risk management framework for landslides in natural and cut slopes has been previously discussed elsewhere [7][8][9][10][11]. After defining the scope of analysis, the framework requires characterizing the magnitude and likelihood of the slope failure, as well as the magnitude and likelihood of the potential consequences, which become the basis for estimating risk (combination of likelihood and consequence for each failure scenario analyzed). The level of risk calculated in this risk analysis (qualitative or quantitative) is evaluated against some adopted criteria to define its tolerability (risk assessment). The outcome of the risk assessment informs decision-making regarding the requirement for risk mitigation and control strategies (e.g., protection, stabilization, early warning), as well as the need for increased knowledge of the slope physical characteristics (e.g., geology, groundwater), failure mechanisms and potential triggers; and monitoring requirements for managing warning), as well as the need for increased knowledge of the slope physical characteristics (e.g., geology, groundwater), failure mechanisms and potential triggers; and monitoring requirements for managing residual risks. This is where remote sensing techniques become provide increased benefit to in-place instrumentation, particularly given the more frequent adoption of quantitative risk analyses since the late 1990s and early 2000's [12][13][14][15]. This risk management framework is illustrated in Figure 1 [11]. Figure 1. General risk management framework for landslides in natural and cut slopes. After [11].
Within the risk management framework in Figure 1, monitoring plays a fundamental role. Not only is robust monitoring the cornerstone of effective early warning, but also fundamental for increased understanding of the landslide mechanisms, kinematics and triggers. In this regard, remote sensing technologies have continued to develop to the point where they are becoming routine practice for landslide monitoring and for landslide investigation as part of the hazard analysis in Figure 1 [16]. These include ground-based and airborne (from manned and unmanned aircrafts) photogrammetry, laser scanning, and satellite and ground-based radar interferometry [17][18][19][20][21][22][23][24][25][26][27][28].
This paper presents a review of six case studies in Western Canada from the authors' experience; where remote sensing was used in combination with other site investigation and monitoring technologies to gain understanding of the slope deformation mechanisms and develop early warning systems. A brief description of the geologic context, landslide extents and deformation patterns are followed by a description of the instrumentation and remote sensing techniques deployed and the insights gained. It is noted that the work presented here reflects the methods that were considered most suited for the case studies discussed, which may vary in other contexts and for other jurisdictions. This work is not an exhaustive literature review of remote sensing techniques for landslide investigation; however, it presents valuable examples for how the information captured from remote sensing techniques enhance the implementation of the landslide risk management framework.

Location and Geology
The Chin Coulee landslide is located near the town of Taber in southern Alberta, Canada. The landslide is within a broad valley that was dammed to develop the Chin Coulee water reservoir. The landslide was likely triggered by a combination of the works required to relocate Highway 36 from near the toe of the slope to the crest of the slope (which required fill placement in the upper part of the slope) and filling of the reservoir Figure 1. General risk management framework for landslides in natural and cut slopes. After [11].
Within the risk management framework in Figure 1, monitoring plays a fundamental role. Not only is robust monitoring the cornerstone of effective early warning, but also fundamental for increased understanding of the landslide mechanisms, kinematics and triggers. In this regard, remote sensing technologies have continued to develop to the point where they are becoming routine practice for landslide monitoring and for landslide investigation as part of the hazard analysis in Figure 1 [16]. These include ground-based and airborne (from manned and unmanned aircrafts) photogrammetry, laser scanning, and satellite and ground-based radar interferometry [17][18][19][20][21][22][23][24][25][26][27][28].
This paper presents a review of six case studies in Western Canada from the authors' experience; where remote sensing was used in combination with other site investigation and monitoring technologies to gain understanding of the slope deformation mechanisms and develop early warning systems. A brief description of the geologic context, landslide extents and deformation patterns are followed by a description of the instrumentation and remote sensing techniques deployed and the insights gained. It is noted that the work presented here reflects the methods that were considered most suited for the case studies discussed, which may vary in other contexts and for other jurisdictions. This work is not an exhaustive literature review of remote sensing techniques for landslide investigation; however, it presents valuable examples for how the information captured from remote sensing techniques enhance the implementation of the landslide risk management framework.

Location and Geology
The Chin Coulee landslide is located near the town of Taber in southern Alberta, Canada. The landslide is within a broad valley that was dammed to develop the Chin Coulee water reservoir. The landslide was likely triggered by a combination of the works required to relocate Highway 36 from near the toe of the slope to the crest of the slope (which required fill placement in the upper part of the slope) and filling of the reservoir [29][30][31]. The location of the Chin Coulee landside and an aerial view are shown in Figure 2 (modified from [30]).
hard. The thickness of the till varies from approximately 20 m at the toe of the slope, to 35 m near the head scarp. Highly fractured shales and coal were encountered underneath the till [31]. The interpreted piezometric elevation is also shown in Figure 2. The landslide stratigraphy and piezometric elevation have been interpreted based on borehole logs (BH-1 to BH-7 in Figure 2), inclinometer readings (in BH-2, 3, 4 and 5), and piezometric readings (BH-1, 6 and 7) [30].

Landslide Extents and Deformation Patterns
The landslide extents were interpreted based on field and photographic mapping of the surface expression of the lateral failure surfaces and the head scarp. The landslide is The local stratigraphy of the landslide is shown in Figure 2. Medium plasticity clay fill material from the initial highway construction and subsequent realignments and regrading is exposed at the head scarp of the landslide. A silty clay till with traces of fine gravel, mantles the valley slope. The till has low to medium plasticity and is very stiff to hard. The thickness of the till varies from approximately 20 m at the toe of the slope, to 35 m near the head scarp. Highly fractured shales and coal were encountered underneath the till [31]. The interpreted piezometric elevation is also shown in Figure 2. The landslide stratigraphy and piezometric elevation have been interpreted based on borehole logs (BH-1 to BH-7 in Figure 2), inclinometer readings (in BH-2, 3, 4 and 5), and piezometric readings (BH-1, 6 and 7) [30].

Landslide Extents and Deformation Patterns
The landslide extents were interpreted based on field and photographic mapping of the surface expression of the lateral failure surfaces and the head scarp. The landslide is approximately 350 m wide and up to 45 m deep. The length of the landslide is approximately 200 m long with the toe of landslide within the reservoir. The failure surface has been inferred based on the location of the head scarp, the shear zones identified in the slope inclinometers, and the location of a low-strength coal seam within the shales and siltstones ( Figure 2). The total landslide volume is estimated at approximately 2 million cubic meters, based on measurement along the proposed failure plane and current landslide extents.
The Chin Coulee landslide has been interpreted as a translational retrogressive landslide seated within the highly fractured bedrock zone and sliding over the low-strength coal seam in Figure 2. Retrogression is inferred from surface tension cracks and concentrated deformations measured in the slope inclinometers near the head scarp. Displacement rates vary between 10 and 50 mm/year.

Instrumentation and Remote Sensing Techniques
Inclinometer readings are available for BH-2, 3, 4 and 5. Other instrumentation at the site includes a differential GPS system consisting of ten units (installed in July of 2018). This system is acquiring high frequency data, and details of its installation and interpretation can be found in [29]. Here, we will focus on the use of remote sensing techniques to obtain quantitative landslide information by leveraging historical aerial photography and state of practice remote sensing acquisition and processing.
Historic stereo pairs were available for 1945,1960,1970,1982,1993,1999 and 2012. The stereopairs in 1982 had the more detailed scale (1:2500), with five photographs overlapping the scene of the landslide. These photographs were used to develop a surface model of the landslide using the software Pix4D (www.pix4d.com), and validated against the location of large boulders outside the landslide area that were assumed to be stable. The surface model was converted (sampled) to a point cloud with a minimum density of 500 points/m 2 .
Laser scanning of the landslide was initiated in July of 2018 using a long-range scanner, ILRIS-LR, with a wavelength of 1064 nm, a pulse frequency of 10 kHz, and a beam divergence of 0.014324 • [32]. Scanning was done from 3 locations across the reservoir from the landslide [29], with the average scan distance on Chin Coulee at 1000 m. The system has an accuracy of ±7 mm, which can be improved to ±4 mm through shot repetition and an averaging option. The scanner rendered over 8 million points, after removing tall vegetation and artificial structures and equipment. Point density ranged between 1000 and 2000 points/m2. Quantitative information of the amount of slope deformation since the movement initiated was sought as a means to improve the understanding on the landslide kinematics for a better interpretation of the mechanisms behind the slope movement. Changes in the slope surface were calculated through comparing the point clouds from the 1982 aerial photographs and a laser scan from July 2019 (Change Detection). The software used for point cloud analysis and change detection was CloudCompare V2.10 (available through www.cloudcompare.org). The method used for change detection was M3C2 [33]. In the M3C2 method, change is quantified as the distance between point clouds measured parallel to the normal vector of the surface in the vicinity of each measurement point. The threshold for change detectability (limit of detection-LOD) was set at two standard deviations of measurements within areas considered stable (zero change). Changes calculated within ± two standard deviations were considered potential measurement randomness. The limit of detection between the 1982 photogrammetric model and the 2019 LiDAR model was 1 m, attributed to a lower accuracy of the 1982 model derived from historic photos and surficial changes between 1982 and 2019 due to erosional processes [30].

Results and Discussion
The results of the Change Detection between the 1982 photographs and the laser scanning in 2019 are shown in Figure 3 [30]. Negative values correspond to material loss and positive values correspond to movement towards the scanner (out of slope).
pushing a passive wedge outward on a sub-horizontal shear zone (seen in the calculated change as gain in material or movement towards the scanner). The sharp change to outward movement suggests the location of the contact between the driving and passive wedges ( Figure 3) [30].
This enhanced knowledge regarding the mechanisms of the Chin Coulee Landslide allows improved determination of the slope stability as the reservoir fluctuates and evaluate the potential for sudden acceleration under reservoir increases and draw down conditions, supporting landslide risk evaluation under diverse scenarios.  Loss of material in the upper half of the slope was measured at approximately 4 m between 1982 and 2019, which is validated by inspection of the landslide scarp when compared to a 2.5-m-high fence and a 1.5-m-long guardrail post ( Figure 3a). Although there are differences in resolution and accuracy between the models, this comparison increased the confidence in the magnitudes of the calculated change for interpretation of the landslide kinematics. Changes in the lower half of the landslide suggest 2 to 3 m of outward movement at the toe of the slope (and up to 4 m in some areas). Unlike the gradational change in the direction of deformation expected from rotational landslides, the transition is sharp between the areas showing material loss (interpreted as downward movement of the landslide mass) and the outward displacement in the lower half of the landslide. This deformation trend is consistent with compound landslides mechanisms, where an active wedge moves downward (seen in the change as a loss of material) pushing a passive wedge outward on a sub-horizontal shear zone (seen in the calculated change as gain in material or movement towards the scanner). The sharp change to outward movement suggests the location of the contact between the driving and passive wedges ( Figure 3) [30].

Location and Geology
This enhanced knowledge regarding the mechanisms of the Chin Coulee Landslide allows improved determination of the slope stability as the reservoir fluctuates and evaluate the potential for sudden acceleration under reservoir increases and draw down conditions, supporting landslide risk evaluation under diverse scenarios.

Location and Geology
Site C018 is a weathered rock slope adjacent to Highway 837 along the south bank of the Red Deer River, approximately 10 km northwest of the town of Drumheller, Alberta ( Figure 4). This slope is predominantly in a weak, poorly indurated sedimentary rock with layers of stronger sandstone, siltstone, and coal [34]. Bedrock is weak to extremely weak, with a UCS of less than 5 MPa and GSI of 20 to 25. The exposed bedding planes are visible in the valley slopes due to the variations of color and texture of the different materials. Bedding structures have an approximately sub-horizontal inclination of 10 • [35].
The slope at site C018 shows extensive signs of weathering processes including d gullies, sinkholes, weak runoff marks, and pipe erosion, where weathering processes ubiquitous and result in a soft soil-like material at surface, with the bedrock struct almost eradicated. Discontinuous blocks of more indurated bedrock remain as weathered rock blocks [34].

Landslide Extents and Deformation Patterns
Site C018 is a slope approximately 520 m long and 90 m high ( Figure 4). This site a long history of instabilities with fallen material reaching and blocking the highway. slope displays three modes of failure, which largely depend on the moisture condition the materials forming the slope: (1) earthflows associated with heavy rainfall; (2) rock f that result from differential weathering rates and undermining of various rock units; (3) slides of frozen slabs of highly weathered rock [34].

Remote Sensing Technique and Analyses
Site C018 presents an example of the use of inexpensive unmanned aerial veh (UAV) photogrammetry as the basis for monitoring surface changes and for develop virtual surface models for geomorphic analysis (in this example, preferential drain paths and rock fall trajectory models).
Three UAV photogrammetry surveys were conducted in December 2017, May 2 and November 2018, respectively. The UAVs used in the surveys had a 12-MP (m pixels) camera for the first and third flight and a 17-MP camera for the second flight. 12-MP camera had a sensor 1/2.3″ CMOS (6.3 mm width and 4.7 mm height), a 94° f of view, a 20 mm focal length, and aperture f 1:2.8. The 17-MP camera had a sensor 1/ CMOS (6.3 mm width and 4.7 mm height), an 82° field of view, a 44 mm focal length, aperture f 1:3.3. Minimum overlap between photographs was 60%, which required The slope at site C018 shows extensive signs of weathering processes including deep gullies, sinkholes, weak runoff marks, and pipe erosion, where weathering processes are ubiquitous and result in a soft soil-like material at surface, with the bedrock structure almost eradicated. Discontinuous blocks of more indurated bedrock remain as less weathered rock blocks [34].

Landslide Extents and Deformation Patterns
Site C018 is a slope approximately 520 m long and 90 m high ( Figure 4). This site has a long history of instabilities with fallen material reaching and blocking the highway. The slope displays three modes of failure, which largely depend on the moisture conditions of the materials forming the slope: (1) earthflows associated with heavy rainfall; (2) rock falls that result from differential weathering rates and undermining of various rock units; and (3) slides of frozen slabs of highly weathered rock [34].

Remote Sensing Technique and Analyses
Site C018 presents an example of the use of inexpensive unmanned aerial vehicle (UAV) photogrammetry as the basis for monitoring surface changes and for developing virtual surface models for geomorphic analysis (in this example, preferential drainage paths and rock fall trajectory models).
Three were developed using the software Pix4D and the surface model was converted (sampled) to a point cloud with a minimum density of 500 points/m 2 for Change Detection using the M3C2 method [35]. The surfaces were also used to develop a surface drainage network following the methods outlined in [36]; and for 2.5-dimensional rock fall trajectory simulations following the methods outlined in [37] and [38].

Results and Discussion
The results obtained by [35] are summarized in Figure  600 and 910 photographs for each of the three surveys in chronological order. Surface models were developed using the software Pix4D and the surface model was converted (sampled) to a point cloud with a minimum density of 500 points/m 2 for Change Detection using the M3C2 method [35]. The surfaces were also used to develop a surface drainage network following the methods outlined in [36]; and for 2.5-dimensional rock fall trajectory simulations following the methods outlined in [37] and [38].

Results and Discussion
The results obtained by [35] are summarized in Figure   The Change Detection results highlighted the state of activity of the areas of the slope that had recorded instabilities (Figure 5a-c). This analysis provided insight to prioritize area F1 for further hazard analysis, as some material was actively being eroded mid-slope and accumulating near the toe. This active erosion was more pronounced on weather- The Change Detection results highlighted the state of activity of the areas of the slope that had recorded instabilities (Figure 5a-c). This analysis provided insight to prioritize area F1 for further hazard analysis, as some material was actively being eroded mid-slope and accumulating near the toe. This active erosion was more pronounced on weathersusceptible rock units, undermining stronger rock above, which could lead to rock fall events with the potential to reach the highway. This information, in combination with the rock fall trajectory simulations in Figure 5c, allowed the development of hazard mitigation strategies at this site. The surface drainage models provided valuable information to evaluate the impact of surface runoff on the potential for earthflow occurrences. These results suggest surface water management strategies as one slope hazard and risk control option, diverting surface runoff away from the crest of the slope.

Location and Geology
The 10-Mile Slide is an active portion of a post-glacial earthflow, and is located in the Province of British Columbia, Canada, North of Vancouver [39] (Figure 6a). The landslide has been monitored for over 40 years, as its deformation has led to frequent maintenance of a highway and railway track that cross the continuously deforming mass. susceptible rock units, undermining stronger rock above, which could lead to rock fall events with the potential to reach the highway. This information, in combination with the rock fall trajectory simulations in Figure 5c, allowed the development of hazard mitigation strategies at this site. The surface drainage models provided valuable information to evaluate the impact of surface runoff on the potential for earthflow occurrences. These results suggest surface water management strategies as one slope hazard and risk control option, diverting surface runoff away from the crest of the slope.

Location and Geology
The 10-Mile Slide is an active portion of a post-glacial earthflow, and is located in the Province of British Columbia, Canada, North of Vancouver [39] (Figure 6a). The landslide has been monitored for over 40 years, as its deformation has led to frequent maintenance of a highway and railway track that cross the continuously deforming mass. The landslide is within the Fraser River valley, where glacial and post-glacial sediments are common along the valley slopes. The bedrock lithology in the area includes andesite to dacite volcanic rocks, sandstone, and shale. These rocks are overlain by quaternary deposits (including glacial drift blankets, colluvium, alluvium, and landslide The landslide is within the Fraser River valley, where glacial and post-glacial sediments are common along the valley slopes. The bedrock lithology in the area includes andesite to dacite volcanic rocks, sandstone, and shale. These rocks are overlain by quaternary deposits (including glacial drift blankets, colluvium, alluvium, and landslide deposits) [39,40]. This is reflected in the stratigraphy of the landslide, shown in Figure 6b, Remote Sens. 2021, 13, 366 9 of 23 which consists of a layer (up to 20 m thick) of mixed landslide deposits including medium to high plastic clays, silts, and the presence of zones of varying sand and gravel content. This layer overlies colluvium materials (sand, silt, and clay), which in turn overlie glacial deposits (stiff, gravelly clayey, and silty till with the presence of glaciofluvial well graded, rounded gravel, and sand deposits).

Landslide Extents and Deformation Patterns
The 10-Mile Slide is approximately 200 m wide and 140 m high, with an approximate volume of 750,000 m 3 [40]. Since the 1980s, displacement rates have varied spatially and temporally, with recorded rates up to 10 mm/day. The area of deformation was inferred from sharp scarps and tension cracks. The railway operator installed a 128-m-long retaining wall to control deformations, to an approximate depth of 16.5 m. and the highway operator has frequently realigned the highway as a response to the landslide movement.
The landslide is sliding on a continuous shear surface that extends upslope beyond the location of the railway tracks ( Figure 6). This continuous shear surface appears to have a dip of approximately 22 • , sub-parallel to the ground surface. Piezometers installed at different locations within the 10-Mile Slide have measured pore water pressure equal to or less than atmospheric pressure, suggesting that increases in pore pressures are not the driving mechanism for slope movements [39,40].

Instrumentation and Remote Sensing Techniques
Several slope inclinometers (SI) have been installed at the site. These were typically sheared within weeks of installation due to the fast deformation velocities [39]. Survey prisms (19 in total) were installed along the retaining wall and have been monitored since 2011. These were manually surveyed at different frequencies, with daily measurements in periods when the landslide showed the fastest deformation velocities. A differential GPS system consisting of 11 GPS units was installed in April of 2017 (same system to the one installed at the Chin Coulee Landslide). This system acquires high frequency data, and details of its installation and interpretation can be found in [41].
Ground-based (also terrestrial) laser scanning (TLS) was conducted in May and July of 2017 using a long-range scanner, ILRIS-LR, the same equipment used at Chin Coulee. Scanning was done from two locations across the river, with the average scan distance approximately 600m and a point density of approximately 1000 points/m 2 [41]. The point clouds from the TLS surveys were used for Change Detection with the software Cloud-Compare V2 and the M3C2 method (same procedures discussed before). The threshold for change detectability (limit of detection-LOD) was set at two standard deviations of measurements within areas considered stable; 5 cm in the case of the two scans.

Results and Discussion
The results of the Change Detection at the 10-Mile Slide between May and July of 2017 are shown in Figure 7. These results are shown together with the deformations measured by the GPS system between April and November of 2017. Figure 7a shows the location of the GPS units, together with the deformation vectors (to a scale of 1:500) and Figure 7b shows the cumulative deformations in selected units for the same period of time.
GPS monitoring showed that left and central areas of the landslide (North West side) have the highest displacement magnitudes (GPS units 46, 47 and 48). At the center of the landslide and near the highway cut, Unit 46 showed 400 mm displacement between April and November 2017. This is an order of magnitude larger than the units at the south-west end of the landslide. Similarly, Unit 43 showed displacement of approximately 30 mm during the same time frame.  Importantly, the areas showing increased deformation within the landslide are consistent between both monitoring techniques. TLS Change Detection allows for better definition of the extent of the most active areas of the landslide, which allows for better interpretation of the landslide kinematics. It has been interpreted that an active block at the left side of the landslide toe is moving out of slope at a faster rate, leading to a decrease in the stability of the slope mass behind, which tends to follow the movements of the lower, left block. TLS Change Detection complemented with the detailed temporal resolution of the GPS Units allows for better interpretation of the landslide mechanisms, and therefore improved modeling of potential stabilization strategies. The setup also allows for the near real-time monitoring frequency required for landslide early warning systems.

Location and Geology
Fourteen landslides have been identified along a 10 km stretch of the Thomson River Valley, south of Ashcroft, British Columbia, Canada [42]. The location of this section of valley is shown in Figure 8.
Eight soil units overlie andesitic bedrock in the study area [43]. These units were deposited by the advance and retreat of three glaciations (Figure 8). Unit 2 is comprised of highly plastic clay interbedded with glaciolacustrine silt (unit 2), heavily consolidated by glacial action. Unit 3 presents large quantities of silt, sand, gravel, and boulders near its base, also consolidated by glacial action. Unit 6 corresponds to glacial till (gravel size particles in a sandy, silty and clayey matrix). The stratigraphy of the valley is further complicated by an ancient Pleistocene river, which had cut down into unit 2, to the east of The results from Change Detection are compatible in magnitude with the GPS units (considering the two-month period between scans). At Unit 48, the cumulative displacement was 75 mm, with a range of 30 to 90 mm with Change Detection. At Unit 47, the displacement was 80 mm, with a range of 60 to 100 mm using Change Detection. At Unit 46, the cumulative displacement was 180 mm. However, Change Detection ranged between 50 and 125 mm. This difference near Unit 46 could be associated with the lower density of TLS points on this area due to the oblique angle of the LiDAR shot.
Importantly, the areas showing increased deformation within the landslide are consistent between both monitoring techniques. TLS Change Detection allows for better definition of the extent of the most active areas of the landslide, which allows for better interpretation of the landslide kinematics. It has been interpreted that an active block at the left side of the landslide toe is moving out of slope at a faster rate, leading to a decrease in the stability of the slope mass behind, which tends to follow the movements of the lower, left block. TLS Change Detection complemented with the detailed temporal resolution of the GPS Units allows for better interpretation of the landslide mechanisms, and therefore improved modeling of potential stabilization strategies. The setup also allows for the near real-time monitoring frequency required for landslide early warning systems.

Location and Geology
Fourteen landslides have been identified along a 10 km stretch of the Thomson River Valley, south of Ashcroft, British Columbia, Canada [42]. The location of this section of valley is shown in Figure 8. the current location of the Thompson River and deposits of boulders, gravels, and sands at depths below the interface of units 2 and 3 have been reported in borehole logs [42].

Landslide Extents and Deformation Patterns
The landslides in the 10 km section of the valley range in volume between 1 and 15 million m 3 and their velocities range from 4 to 30 mm/year to over 100 mm/year, although some have had faster displacement episodes of meters per day to meters per hour [42]. The areal extent of these landslides is shown in Figure 8b. These slides typically develop as compound slides, with an initial sliding mass that translates sub horizontally along the surfaces of rupture. The surfaces of rupture of the landslides through this section of the Thompson River Valley appear to occur along a common weak layer of clay located below the bottom of the river. After some movements, displaced masses become unstable and active wedges form at the main scarps. The weight of these wedges drives the mass towards the river [42].
One of these landslides, the Ripley Landslide, has been extensively investigated and shows an annual cycle of deformation accelerating in the Fall (September) and decelerating in the Spring (April), as a response to changes in river elevation. This landslide has shown annual deformations over 100 mm/year [42,44].

Instrumentation and Remote Sensing Techniques
Only a subset of these landslides have been instrumented, most with slope inclinometers that are manually read, annually. The Ripley landslide, however, has become a field laboratory for investigation and monitoring. Monitoring instrumentation at this landslide include grouted-in vibrating wire piezometers, inclinometer casings, Shape Accel Arrays (SAA), Differential GPS systems, and a wave-guide system that uses waves emitted as the landslide deforms and interprets them as relative displacement; Eight soil units overlie andesitic bedrock in the study area [43]. These units were deposited by the advance and retreat of three glaciations (Figure 8). Unit 2 is comprised of highly plastic clay interbedded with glaciolacustrine silt (unit 2), heavily consolidated by glacial action. Unit 3 presents large quantities of silt, sand, gravel, and boulders near its base, also consolidated by glacial action. Unit 6 corresponds to glacial till (gravel size particles in a sandy, silty and clayey matrix). The stratigraphy of the valley is further complicated by an ancient Pleistocene river, which had cut down into unit 2, to the east of the current location of the Thompson River and deposits of boulders, gravels, and sands at depths below the interface of units 2 and 3 have been reported in borehole logs [42].

Landslide Extents and Deformation Patterns
The landslides in the 10 km section of the valley range in volume between 1 and 15 million m 3 and their velocities range from 4 to 30 mm/year to over 100 mm/year, although some have had faster displacement episodes of meters per day to meters per hour [42]. The areal extent of these landslides is shown in Figure 8b. These slides typically develop as compound slides, with an initial sliding mass that translates sub horizontally along the surfaces of rupture. The surfaces of rupture of the landslides through this section of the Thompson River Valley appear to occur along a common weak layer of clay located below the bottom of the river. After some movements, displaced masses become unstable and active wedges form at the main scarps. The weight of these wedges drives the mass towards the river [42].
One of these landslides, the Ripley Landslide, has been extensively investigated and shows an annual cycle of deformation accelerating in the Fall (September) and decelerating in the Spring (April), as a response to changes in river elevation. This landslide has shown annual deformations over 100 mm/year [42,44].

Instrumentation and Remote Sensing Techniques
Only a subset of these landslides have been instrumented, most with slope inclinometers that are manually read, annually. The Ripley landslide, however, has become a field laboratory for investigation and monitoring. Monitoring instrumentation at this landslide include grouted-in vibrating wire piezometers, inclinometer casings, Shape Accel Arrays (SAA), Differential GPS systems, and a wave-guide system that uses waves emitted as the landslide deforms and interprets them as relative displacement; among others. Site investigation techniques deployed at the site include borehole logging, geophysical investigations (seismic, electric resistivity), real time electric resistivity to calculate changes in unsaturated water contents, TLS and aerial laser scanning (ALS); among others [44][45][46][47]. Some of the in-place instrumentation at the Ripley Landslide are shown in Figure 9. These instruments have made it possible to arrive at the interpretation of the landslide kinematics as described above.
Remote Sens. 2021, 13, x FOR PEER REVIEW 12 of 22 among others. Site investigation techniques deployed at the site include borehole logging, geophysical investigations (seismic, electric resistivity), real time electric resistivity to calculate changes in unsaturated water contents, TLS and aerial laser scanning (ALS); among others [44][45][46][47]. Some of the in-place instrumentation at the Ripley Landslide are shown in Figure 9. These instruments have made it possible to arrive at the interpretation of the landslide kinematics as described above. Large-scale Satellite Interferometric Synthetic Aperture Radar (InSAR) was used in this section of the Thomson River Valley to understand the state of activity of slopes that are not being monitored at frequent as the Ripley Landslide, and potentially detect other active slopes not being monitored. The images used were captured by the Canadian satellite RADARSAT-2. East-facing images (group F4N) were collected between 2011 and 2015, and west-facing images (groups U5 and U21) were collected between 2013 and 2015. Images were acquired every 24 days, with some exceptions when longer gaps occurred. SAR images were processed using Persistent Scatterer Interferometry (PSI), and displacements measured are the average displacements for areas up to 15 m × 15 m [48].

Results and Discussion
The LOS displacement rate measured by the image group U21 at active landslides is shown in Figure 9 [48]. The average precision of measurements was 1 mm considering the Standard Deviation of measurements at stable areas. Images were acquired every 24 days, with some exceptions when longer gaps occurred. SAR images were processed using Persistent Scatterer Interferometry (PSI), and displacements measured are the average displacements for areas up to 15 m × 15 m [48].

Results and Discussion
The LOS displacement rate measured by the image group U21 at active landslides is shown in Figure 9 [48]. The average precision of measurements was 1 mm considering the Standard Deviation of measurements at stable areas.
The work in [48] used information about the LOS orientation of the satellite image groups and assumptions of the likely landslide displacement vectors based on the landslide kinematics prevalent in the area. Their analysis identified six areas of concentrated displacements within the study area, all within or adjacent to the footprints of identified landslides (Figure 9). Average LOS displacement rates ranged between 11 mm/year (Barnard) and 39 mm/year (Ripley Landslide). The maximum downslope displacement rate (Projection of the LOS to expected landslide displacement vector) was calculated at 120 mm/year (toe of the Goddard Slide). Maximum downslope displacements were between 73 and 89 mm/year for the other landslides, except for the Barnard Slide, with a maximum downslope displacement of 50 mm/year. Use of InSAR allowed the identification of landslides showing an Active state of activity, including the Ripley Landslide and the identification of an active mass adjacent to the South landslide (South Extension in Figure 9a), which had not been previously identified. It also allowed the identification of landslides potentially at a Reactivated state (Goddard Toe, Barnard, North Toe, and Red Hill landslides). These insights provide valuable information for an enhanced understanding of the corridor's risk level associated with potential landslide events and for improved landslide stability analyses for risk assessment purposes [49,50].
The InSAR analysis also identified that most of the Thompson River Valley landslides exhibit seasonal changes in the rate of displacement (median displacements particularly at the Ripley Landslide, Red Hill toe, Goddard toe, and South landslides). These typically show increased displacements during the late fall and winter months versus relative stability during the spring and summer. This suggests that factors that affect the stability of these landslides are common across the study area, and that these landslides are marginally stable [48].

Location and Geology
The Checkerboard Creek rock slope is located within the Revelstoke Reservoir, near the town of Revelstoke, British Columbia, Canada. Figure 10 shows the location of this landslide and a front view from across the reservoir.
The typical cross section in Figure 10 presents a sketch of the materials and structure forming the slope. Massive to weakly foliated granodiorite overlie gneiss and schist of a major regional fault system (Columbia River Fault), which has developed a broad, regional, brittle deformation zone of altered and mechanically deformed rock. Shears zones and joints in the slope are frequent, and dip steeply into and out of the slope at angles between 60 • and 90 • from horizontal. The rock mass quality ranges from very strong, fresh, undisturbed and blocky rock to highly weathered and altered, weak and disturbed. The lower quality rock masses contain frequent sheared and crushed zones and are typically found within 60 m of the slope surface [11,51].
have shown steady rates of deformation [11,52,53]. The displacement rate is 0.5 mm/year at the boundaries of the actively deforming mass and up to 13 mm/year within the most active area. Deformation rates are greatest at the surface of the slope and decrease progressively, with depths to approximately 50 to 60 m below surface, where no deformations are detected. This depth coincides with the weathered and sheared volume of rock described above.

Instrumentation and Remote Sensing Techniques
The Checkerboard Creek rock slope has an extensive array of in-place instrumentation, mostly installed between 1984 and 2000. The instrumentation includes surface survey monuments (read annually), manual slope inclinometers (read annually) and in-place probes (read every 6 h), nested standpipe piezometers, multipoint borehole extensometers (read every 4 h), surface cable extensometers (read annually), strain meters, borehole and surface thermistors, and a weather station that records daily temperature and precipitation values [11,52]. A global navigation satellite system (GNSS) network has also been installed at the slope, which comprises of four sensors. Airborne laser scanning (ALS) was also carried out in September 2014 to acquire a high-resolution digital elevation model.
An IBIS-L GB-InSAR (ground-based, interferometric, synthetic aperture radar) manufactured by IDS Georadar was installed in October 2016 across the Revelstoke Reservoir from the Checkerboard Creek rock slope. The system has a radar head that acquires data as it travels on a 2-m-long rail to obtain a two-dimensional image from the area of interest (synthetic aperture). The system can target areas at distances of up to 4

Landslide Extents and Deformation Patterns
The Checkerboard Creek rock slope is 260 m high, from Highway 23 to the middle reach of the Checkerboard Creek ( Figure 10). The width of the slope is 600 m and the overall slope angle is 30 • , steeper at the toe (45 • ) and flatter in the upper area (25 • ) [11]. The active zone of deformation has an average slope angle of 45 • , and deformations have been detected down to 50 to 60 m deep. The volume of this active zone is estimated to be between 2 and 3 million m 3 .
Slope deformations are concentrated at the central part of the slope, near the crest of the slope cut required to accommodate the highway. Displacements measured by instruments near the surface have revealed an annual cycle dominated by an active period from early October to April/May (early autumn throughout late winter), followed by a relatively quiet period between May and September (spring and summer). These cycles have been interpreted as temperature-dependent, where expansion and contraction of the rock blocks near the surface cause changes in the normal stresses at the contact of joint surfaces, and therefore changing the available shear strengths and deformation characteristics. Other instruments and measurements deeper into the active rock mass have shown steady rates of deformation [11,52,53]. The displacement rate is 0.5 mm/year at the boundaries of the actively deforming mass and up to 13 mm/year within the most active area. Deformation rates are greatest at the surface of the slope and decrease progressively, with depths to approximately 50 to 60 m below surface, where no deformations are detected. This depth coincides with the weathered and sheared volume of rock described above.

Instrumentation and Remote Sensing Techniques
The Checkerboard Creek rock slope has an extensive array of in-place instrumentation, mostly installed between 1984 and 2000. The instrumentation includes surface survey monuments (read annually), manual slope inclinometers (read annually) and in-place probes (read every 6 h), nested standpipe piezometers, multipoint borehole extensometers (read every 4 h), surface cable extensometers (read annually), strain meters, borehole and surface thermistors, and a weather station that records daily temperature and precipitation values [11,52]. A global navigation satellite system (GNSS) network has also been installed at the slope, which comprises of four sensors. Airborne laser scanning (ALS) was also carried out in September 2014 to acquire a high-resolution digital elevation model.
An IBIS-L GB-InSAR (ground-based, interferometric, synthetic aperture radar) manufactured by IDS Georadar was installed in October 2016 across the Revelstoke Reservoir from the Checkerboard Creek rock slope. The system has a radar head that acquires data as it travels on a 2-m-long rail to obtain a two-dimensional image from the area of interest (synthetic aperture). The system can target areas at distances of up to 4 km, covering a horizontal angle of up to 80 • , and vertical angle of up to 40 • . The resolution of the image therefore depends on the distance between the slope and the radar.
Challenges for system deployment include the inaccessibility to the electric power grid, which requires a self-sustained power supply. This is compounded with the climate of the site, which includes cold and dark winters, which decrease the solar power potential in the area and decreases the storage of the battery system. Given the very slow deformation rates of the Checkerboard Creek rock slope, and the challenges of power supply at the site, GB-InSAR monitoring was set at to D-InSAR mode (discontinuous, short duration monitoring campaigns), with monitoring campaigns of 1 to 2 weeks long [53]. Within each monitoring campaign, images were acquired at an average frequency of 6 images per hour. Due to the dense vegetation cover, GB-InSAR images were restricted to the exposed rock face adjacent to the highway.

Results and Discussion
The map of displacements in the LOS between the GB-InSAR and the Checkerboard Creek rock slope is shown in Figure 11. The extensive in-place instrumentation at the site allowed interpretation of temporal trends and insights into the spatial deformation patters; however, the GB-InSAR results allowed for an enhanced definition of the active portions of the exposed rock face. This active area detected by the GB-InSAR is bounded by the crest of the exposed rock face and the dark dashed line in Figure 11a. The average deformations in this area are consistent with those measured by in-place instruments. The threshold values used for identifying the most active area was 2 mm/year [52,53].
It was also identified that there is a transition between faster deformation rates near the crest of the expose rock face and its toe. This is consistent with the interpreted deformation kinematics in which the rock blocks within a heavily disaggregated rock mass are showing toppling like behavior and relative displacements between one another [53].
A critical insight from the GB-InSAR monitoring campaign was the identification of a potentially unstable block that appears to be sliding on a ledge defined by a persistent discontinuity that dips 15 • to 20 • out of slope ( Figure 11). This newly identified potentially unstable block is within an area of the exposed rock face that was previously interpreted as less active. This block was most active between September 2016 and September 2017; however, it did not show increased activity when compared to the surrounding area after that period. This area of the exposed rock face was identified for continuous assessment as part of the slope's risk management system [53].
Creek rock slope is shown in Figure 11. The extensive in-place instrumentation at the site allowed interpretation of temporal trends and insights into the spatial deformation patters; however, the GB-InSAR results allowed for an enhanced definition of the active portions of the exposed rock face. This active area detected by the GB-InSAR is bounded by the crest of the exposed rock face and the dark dashed line in Figure 11a. The average deformations in this area are consistent with those measured by in-place instruments. The threshold values used for identifying the most active area was 2 mm/year [52,53]. Figure 11. GB-InSAR displacement rate in the LOS at the Checkerboard Creek rock slope; crest of exposed rock face is shown with solid dark line, boundary of active deformation is shown in dark dashed line (a). View of the radar system looking towards the slope (b). D-InSAR average LOS displacement results for the active crest of the exposed rock face compared against a survey point and an extensometer located near the crest of the exposed rock face (c). Adapted from [53].
It was also identified that there is a transition between faster deformation rates near the crest of the expose rock face and its toe. This is consistent with the interpreted deformation kinematics in which the rock blocks within a heavily disaggregated rock Figure 11. GB-InSAR displacement rate in the LOS at the Checkerboard Creek rock slope; crest of exposed rock face is shown with solid dark line, boundary of active deformation is shown in dark dashed line (a). View of the radar system looking towards the slope (b). D-InSAR average LOS displacement results for the active crest of the exposed rock face compared against a survey point and an extensometer located near the crest of the exposed rock face (c). Adapted from [53].

Location and Geology
Tornado Mountain is north of the town of Sparwood, at the boundary between the provinces of British Columbia and Alberta, Canada. The study site corresponds to a steep rock outcrop on the British Columbia side, which in reality is part of Mount Lyne, adjacent to Tornado Mountain ( Figure 12). The outcrop is 50 m high and near vertical, in strong, blocky limestone above a steep (approximately 45 • ) talus slope. The slope of the mountain side becomes gradually softer to an average of 22 • as it approaches a railway alignment. This lower section of the slope is covered by rock fragments, soil, and is sparsely vegetated. This site has witnessed large blocks (up to 1.5 m in equivalent diameter) rolling downslope to the vicinity of the railway tracks [54].
The S042 Site is located adjacent to Highway 742, approximately 5 km southwest of Canmore, Alberta, Canada ( Figure 12). This site consists of a near-vertical, 80-m-high rock slope on limestones, dolomitic limestones, dolostones and shale. Rock falls and erosional processes have created a talus slope of 40 • between the slope and the highway. Rock falls mostly are captured on the talus slope; however; blocks of less than 0.3 m in equivalent diameter are commonly found on the road [55].

Results and Discussion
Ground-based, aerial and UAV photogrammetry have become standard practice for developing high-resolution surface models of rock slopes for structural mapping, both in civil and mining projects. This information is then used as input to estimate rock mass characteristics (e.g., expected behavior and strength) and to develop geotechnical stability and hydrogeologic slope models. The analyses at Tornado Mountain and S042 Site are an example of how this information can also be used to calculate rock fall hazards in terms of volume-frequency curves. The method and results at Tornado Mountain and the S042 Site are detailed in [54] and [56]. A similar approach is described in [57] and tested for a rock fall site in Italy.
The results of the analyses at Tornado Mountain are shown in Figure 13. This figure shows the results of some steps of the analysis, including virtual mapping of discontinuities on the photogrammetric model (Figure 13a,b), as well as one iteration of the DFN developed based on the virtual mapping ( Figure 13c) and the resulting fragmented rock fall volume-frequency curve (Figure 13d). Similar results were obtained for the S042 site ( Figure 14). This figure shows the resulting discontinuity sets from virtual mapping in the photogrammetric model (Figure 14a), one iteration of the DFN developed based on the virtual mapping (Figure 14b), and the resulting fragmented rock fall volumefrequency curve (Figure 14c). Unfortunately, a view of the structures being mapped during the virtual mapping process was not available.
For the two cases analyzed, the rock fall volume-frequency relationships obtained with the DFN's were consistent with those from surveyed blocks. These results show that the approach could be used to estimate the fragmented rock fall volume-frequency relationships in the absence of rock fall records or block surveys along rock fall trajectories, as part of rock fall hazard and risk analyses. The approach presented in [54] and [56] was based on the hypothesis that fragmentation in strong rocks (particularly limestones characterized by Unconfined Compressive Strengths over 50 MPa) will tend to occur through weak planes and by extension of non-persistent discontinuities as the in situ blocks detach, fall, bounce, and roll. In this regard the method should not be applied to other, weaker rocks, without validation. Both sites have been the matter of rock fall hazard investigations as part of assessments into the requirements and dimensioning of rock fall protection structures, for the railway tracks and highway 742 at the Tornado Mountain and Site S042, respectively. As part of these investigations, fragmented rock fall volume-frequency relationships were calculated on the basis of rock block surveys and photogrammetric techniques for virtual structural mapping (fragmented rock falls refer to fallen block volumes after detachment and disaggregation along their trajectory).

Photogrammetry and Structural Mapping
Ground-based photogrammetry at the study areas were performed with digital singlelens reflex (DSLR) cameras. The sensor size used was 35.9 mm by 24 mm, with 36.3 million pixels. The lens used at the Tornado Mountain site had a 200 mm fixed focal length and the lens used for Site S042 had a 150 mm fixed focal length. The distance between the camera and the slope were between 500 m and 600 m at the Tornado Mountain site, and between 200 m and 600 m at Site S042. The setup described rendered ground pixel sizes between 0.7 cm and 2 cm [56].
The photographs at both sites were processed in ADAM Technology's photogrammetry software suite (https://www.adamtech.com.au), which builds a digital surface model using photogrammetric principles and allows automatic and manual detection and mapping of rock structures. Orientations were mapped manually within the software for both sites.
The dip, dip direction and spacing of the sets of discontinuities were fitted to probability distributions and used to develop multiple iterations of discrete fracture networks (DFN) of persistent discontinuities for both sites using Wolfram's Mathematica software package (https://www.wolfram.com). These DFNs were used to calculate the fragmented rock fall volume-frequency distributions through numerical integration [56].

Results and Discussion
Ground-based, aerial and UAV photogrammetry have become standard practice for developing high-resolution surface models of rock slopes for structural mapping, both in civil and mining projects. This information is then used as input to estimate rock mass characteristics (e.g., expected behavior and strength) and to develop geotechnical stability and hydrogeologic slope models. The analyses at Tornado Mountain and S042 Site are an example of how this information can also be used to calculate rock fall hazards in terms of volume-frequency curves. The method and results at Tornado Mountain and the S042 Site are detailed in [54,56]. A similar approach is described in [57] and tested for a rock fall site in Italy.
The results of the analyses at Tornado Mountain are shown in Figure 13. This figure shows the results of some steps of the analysis, including virtual mapping of discontinuities on the photogrammetric model (Figure 13a,b), as well as one iteration of the DFN developed based on the virtual mapping ( Figure 13c) and the resulting fragmented rock fall volume-frequency curve (Figure 13d). Similar results were obtained for the S042 site ( Figure 14). This figure shows the resulting discontinuity sets from virtual mapping in the photogrammetric model (Figure 14a), one iteration of the DFN developed based on the virtual mapping (Figure 14b), and the resulting fragmented rock fall volume-frequency curve (Figure 14c). Unfortunately, a view of the structures being mapped during the virtual mapping process was not available. For the two cases analyzed, the rock fall volume-frequency relationships obtained with the DFN's were consistent with those from surveyed blocks. These results show that the approach could be used to estimate the fragmented rock fall volume-frequency relationships in the absence of rock fall records or block surveys along rock fall trajectories, as part of rock fall hazard and risk analyses. The approach presented in [54,56] was based on the hypothesis that fragmentation in strong rocks (particularly limestones characterized by Unconfined Compressive Strengths over 50 MPa) will tend to occur through weak planes and by extension of non-persistent discontinuities as the in situ blocks detach, fall, bounce, and roll. In this regard the method should not be applied to other, weaker rocks, without validation. Figure 13. Virtual structural mapping (a) and resulting discontinuity sets (b) for the Tornado Mountain study site. One iteration of the DFN for the Tornado Mountain site (c). fragmented rock fall volume cumulative distribution based on the DFN realizations compared to block volumes surveyed along the rock fall trajectories, DFN volumes are truncated to the minimum range of surveyed block volumes for accurate comparison (lack of survey of small fraction of blocks) (d). Adapted from [55,56].

General Discussion and Conclusions
Monitoring plays a fundamental role in landslide risk management. Not only is robust monitoring the cornerstone of effective early warning, but also fundamental for increased understanding of the landslide mechanisms, kinematics and triggers. There is extensive experience in landslide site characterization and monitoring with the use of in situ testing, instrumentation, and remote sensing techniques; and technology keeps evolving such that instruments become more accurate and reliable. In this regard, remote sensing technologies have continued to develop to the point where they are becoming routine practice as part of landslide hazard and risk assessments and management. These include ground-based and airborne (from manned and unmanned aircrafts) photogrammetry, laser scanning, and satellite and ground-based radar interferometry [21][22][23][24][25][26][27][28]. This manuscript takes advantage of case studies from the authors' experience to illustrate the role of remote sensing approaches for implementing the systematic landslide risk management framework that has been widely adopted by industry [8][9][10]58].
Transportation corridors in Western Canada have to traverse the Canadian Cordillera and river valleys in the prairies, in which landslide hazards are common. The importance of these transportation corridors for the Canadian economy, and their vast extent, require cost effective techniques for landslide investigation and monitoring. Remote sensing techniques have been found to play a key role in achieving this. This paper synthesized six case studies in Western Canada from the authors' experience in which remote sensing was used in combination with other site investigation and monitoring technologies to gain understanding of the slope deformation mechanisms and develop early warning systems. These case studies cover different landslide types, including rock falls, compound soil and rock landslides, rock mass landslides, translational soil landslides and earth flows. Remote sensing techniques covered in these case studies include ground-based and satellite InSAR, ground-based and UAV photogrammetry, as well as photogrammetry on historical photos, and ground-based laser scanning. This illustrates the versatility of remote sensing techniques for landslide investigation showing diverse kinematics and deformation mechanisms, as well as the diversity of available practical remote sensing techniques.
The case studies presented here demonstrate the deployment of remote sensing techniques for landslide investigation and monitoring, and how the most gain is obtained when complementing these techniques with information from other, more traditional in-place instrumentation (e.g., slope inclinometers, GPS units, survey prisms). The insights gained in these six case studies, including better definition of landslide deformation trends, both temporally and spatially, improvements in the understanding of landslide mechanisms, insights into rock fall volume-frequency distributions, and identification of the state of activity of large landside at a regional extent, among others. The insights illustrated in this manuscript are based on projects in civil engineering works; however, these are consistent with the reliance on remote sensing techniques by the mining industry to manage potential open pit slope failures. In this regard, remote sensing techniques have already been identified as essential for the successful implementation of open pit slope risk management framework [59]. It is the intent of the authors that these case studies will provide useful examples of the use of remote sensing techniques from a practical perspective and in a variety of contexts. It is noted that the work presented here reflects the methods that were considered most suited for the case studies discussed, which may vary in other contexts and for other jurisdictions. This work is not an exhaustive literature review of remote sensing techniques for landslide investigation; however, they present valuable examples for how the information captured from remote sensing techniques enhance the implementation of the landslide risk management framework.