Introduction

Where should landslide researchers look for changes in slope stability induced by a warming climate? The most straightforward answer, which has been recognized since at least 1990 (e.g., Slaymaker 1990; Haeberli and Beniston 1998; Geertsema et al. 2006; Gruber and Haeberli 2007; Coe and Godt 2012; Huggel et al. 2012; Deline et al. 2015a; Gariano and Guzzetti 2016), is steep-cryospheric terrain in mountainous regions of the world because ice (either in permafrost or glaciers) in these areas is being degraded and reduced by warming temperatures (e.g., Paul et al. 2004; Harris et al. 2009; Huss and Hock 2015; Pastick et al. 2015; Kos et al. 2016; Zemp et al. 2019; Hock et al. 2020). An added benefit to these areas is that landslides are unlikely to be influenced by land use changes brought on by climate change. Land use changes by themselves can alter the frequency and magnitude of landslides (Sidle and Burt 2009; Petley 2010; Crozier 2010; Wasowski et al. 2010).

About 10% of the worldwide populace lives in mountainous regions that are susceptible to hazards resulting from a degrading cryosphere (Hock et al. 2020). The loss of ice in these regions impacts rock slopes in a variety of ways (e.g., Evans and Clague 1994; McColl 2012) including causing a reduction in cohesion and an increase in pore-water pressure in fractures in rocks (e.g., Gruber and Haeberli 2007) and debutressing (e.g., Deline et al. 2015a) and/or thermomechanically damaging (Grämiger et al. 2018) rock slopes when it occurs in the form of receding glaciers. Rock-slope failures in mountainous regions can trigger cascading hazard events, particularly when they enter lakes or fiords, where outburst floods, debris flows, and tsunamis can be generated (e.g., Evans and Delaney 2015; Haeberli et al. 2017; Higmann et al. 2018). The formation of new lakes and the expansion of fiords in areas of retreating glaciers increase the probability for cascading events that initiate as slope failure, and create hazards and risks in downslope areas where they may have not existed in the past (e.g., Haeberli et al. 2017). Risks to humans from these events increase as growing populations, tourism, transportation, and hydropower development expand into mountainous regions (Hock et al. 2020). A better understanding of changes in rock-slope stability caused by the loss of ice could lead to improved delineation of hazard zones, improved models for forecasting landslides and associated cascading hazards, and identification of areas and populations with the greatest exposure to these hazards (Haeberli et al. 2017; Hock et al. 2020).

Worldwide, only a few studies have built systematic, long-term (≥ 20 years) landslide catalogs that can be used to evaluate and compare the magnitude and frequency of landslide events through time in mountainous regions (Schlögel et al. 2011; Bessette-Kirton and Coe 2016; Coe et al. 2108; Hilbert et al. 2019), a problem that has been previously recognized by landslide researchers (e.g., Cox et al. 2015; Wood et al. 2015) and a recent Intergovernmental Panel on Climate Change (IPCC) special report (Hock et al. 2020, p. 175) that lamented a lack of studies to “systematically characterise and compare magnitude and extent of past impacts and future risk across high-mountain regions.”

Remotely sensed imagery and seismic data are well suited to systematic monitoring in mountainous regions because they can provide observations of areas that are usually very difficult to access on the ground (e.g., Kääb 2002; Hibert et al. 2019). However, the lack of systematic landslide catalogs is due in large part to a lack of scale-appropriate remote monitoring tools, either (1) long-term, readily accessible imagery with high (≤ 10 m) and consistent resolution (to avoid mapping bias due to changing resolution) acquired at short (< 30-day) acquisition intervals or (2) dense networks of permanent seismic monitoring stations (station spacing <~ 10 km(?)) in large (hundreds to thousands of km2) mountainous regions that are capable of routinely detecting a wide range of landslide sizes. Recent work by Dunning et al. (2015), who documented how a rock avalanche onto a glacier in New Zealand became masked by snow in remotely sensed imagery within 3 months after its occurrence, highlights the need for acquiring imagery at short time intervals.

The increasing number of Earth-observing satellites with high resolution, global coverage, and short acquisition periods, as well as the emerging field of environmental seismology (e.g., Kanamori and Given 1982; Faillettaz et al. 2008; Walter et al. 2008; and a special journal issue edited by Allstadt et al. 2019) applied to landslides (e.g., Berrocal et al. 1978; McSaveney and Downes 2002; Caplan-Auerbach and Huggel 2007; Helmstetter and Garambois 2010; Yamada et al. 2012; Allstadt 2013; Ekström and Stark 2013; Hibert et al. 2014; Allstadt et al. 2017), has the potential to revolutionize systematic remote monitoring of landslides. For example, in terms of new satellites, the joint U.S. National Aeronautics and Space Administration (NASA) and Indian Space Research Organisation (ISRO) Synthetic Aperture Radar Mission (NISAR) will launch in 2021 with the goal of repeat coverage of 70% Earth’s ice-covered surfaces every 12 days and a displacement measuring capability of 100 mm or less (see https://nisar.jpl.nasa.gov/). Additional examples include the European Space Agency (ESA)’s Sentinel-2 (https://sentinel.esa.int/web/sentinel/missions/sentinel-2) and Planet Labs (planet.com) satellites, both of which acquire moderate- to high-resolution (10–60 m for Sentinel-2 and ≤ 5 m for Planet Labs) optical imagery of the Earth’s landmass surface at short time intervals (< 10 days for Sentinel-2 and every day for Planet Labs). Given the abundance of these and other new imaging opportunities, as well as previously existing opportunities (e.g., Landsat [https://landsat.gsfc.nasa.gov/] and DigitalGlobe [https://www.digitalglobe.com]), focusing resources on a limited number of sites, rather than on the world as a whole, would seem to be the most productive way to make long-term systematic observations. This brings us back to the question of where to look for landslide changes.

This paper focuses on this question and proposes five bellwether sites for long-term slope stability monitoring using remote methods. The various definitions of bellwether include “an indicator or predictor of something,” “one that takes the lead or initiative,” and “an indicator of trends.” These definitions all capture the essence of why the term is appropriate for selecting landslide monitoring sites for studies of the impacts from climate change. That is, the sites will potentially provide indications of changes that will eventually occur over broader mountainous areas. The types of landslides that could be monitored at these sites include rock falls, rockslides, rock avalanches, and debris flows. Below, I describe the site selection methodology and site characteristics and discuss possible approaches for making systematic observations at the sites. I have written this paper for the 5th World Landslide Forum (WLF), Theme 7, “Frontiers in landslide science and innovative practices.” In addition to proposing bellwether sites, my intention in this paper is to provide the scientific justification for the work.

Methods

Selection of bellwether sites

To identify potential bellwether sites, I first searched the literature looking for mountainous areas with the likely presence of mountain permafrost, extensive glaciers that were receding or losing volume, past and future projected increases in air temperatures, and previous work and publications that documented landslide responses to warming temperatures and degrading permafrost or glaciers. After searching the literature, I spent several weeks looking through imagery in Google Earth and DigitalGlobe for mountain locations with a visual abundance of landslide deposits and/or sediment movement in the form of extensive glacial debris, recent glacial outburst flood deposits, or fresh landslide tracks on slopes. My philosophy in looking for recent landslide and sediment movement was that sites with existing moderate to high levels of activity would be most productive for gaining new insights into changes in landslide processes and hazards. By moderate to high levels of activity, I mean levels of activity that appear anomalous compared to their immediate surroundings, as well as compared to landslide activity in mountainous areas of other geographic regions. At the U.S. Geological Survey (USGS), we have had success using this approach to select in situ instrumental monitoring sites (e.g., the Chalk Cliffs debris-flow site in Colorado; Coe et al. 2008) and it appears that many other landslide researchers have used a similar approach to select their in situ monitoring sites (e.g., Hürlimann et al. 2019; Weber et al. 2019).

In addition to considering existing research and current levels of landslide activity, I also strove to identify potential bellwether sites with areas of about 5000 km2 and with a broad worldwide geographic distribution. I used 5000 km2 as target size because it was sufficiently large to characterize activity in a region, without being so large that it would be unmanageable. Additionally, that size had worked well for previous long-term systematic landslide observations in Alaska (Bessette-Kirton and Coe 2016; Coe et al. 2018).

After identifying potential sites, I used a tiered, semi-quantitative selection procedure with four levels (Fig. 1) to select bellwether sites. The criteria that I used for potential sites to advance to each successively higher level are shown in Fig. 1. Criteria became more restrictive as progressing from level 1 to level 4, with level 4 being five bellwether sites that I selected. At level 1, I began with the seven continents and a goal of selecting one site from each continent. To advance to level 2, sites had to have extensive (≥ 10) glaciers, mountain permafrost, and evidence for at least one landslide in the last 10 years. No sites from Antarctica or Africa met the level 2 criteria. To advance to level 3, sites needed to have met the previous criteria but also needed to have evidence for at least three landslides in the last 10 years. The criteria to advance to level 4 were more complicated and somewhat more subjective. Level 4 criteria included evidence for long-term records of landslides, existing landslide-monitoring infrastructure, and elevation and geographic position characteristics that were advantageous to achieve a broad range in bellwether site distribution characteristics. For example, the site that I selected in the Karakoram Range had a range in elevations unlike any of the other sites, which was advantageous in reaching level 4.

Fig. 1
figure 1

Bellwether site selection tree

For each bellwether site that I selected, I also compiled information on processes not related to climate change that could impact landslide magnitude and frequency. These processes included uplift related to plate tectonics or isostatic rebound, and earthquake occurrence. In an ideal situation, all landslide bellwether sites would not be impacted by any of these other processes, thus making the impact from warming temperatures easier to isolate and identify. Unfortunately, nearly all mountain ranges in cryospheric regions are impacted by at least one of these processes, and thus, they cannot be ignored.

Results

Selected bellwether sites

I have selected five bellwether sites (Fig. 2) that are geographically distributed in the northern and southern hemispheres and on five continents. All sites are in mid-latitude areas between 30 and 60° (Fig. 2, Table 1), are partially covered by glaciers, and contain permafrost as modeled by Gruber (2012a, 2012b). All of the sites are at active, or formerly active, plate boundaries (Table 1) and are therefore susceptible to earthquake-triggered landslides, as well as climatically triggered landslides. Since 1901, mean air temperatures at all of the sites have increased by at least 0.5 °C (IPCC 2013; Fig. 3) and, using the medium-low Representative Concentration Pathway (RCP) 4.5 emissions scenario that stabilizes radiative forcing in the year 2100 (Thomson et al. 2011), future mean air temperatures at all of the bellwether sites are projected to increase by at least another 1 °C by 2065 (van Oldenborgh et al. 2013; Fig. 4). Projections of future changes in precipitation at the sites are uncertain, with estimated changes falling within a standard deviation of present-day variability (van Oldenburgh et al. 2013, pp. 1320–1321). For this reason, in this paper, I focus on changes in temperature, whose effects on the cryosphere will probably dominate over those induced by changes in precipitation (e.g., Alifu et al. 2016).

Fig. 2
figure 2

Index map showing proposed bellwether sites in red font. Land masses are shaded light gray, tectonic plates are labeled, and boundaries are shown by black lines

Table 1 Characteristics of the five bellwether sites
Fig. 3
figure 3

Observed changes in mean surface air temperature from 1901 to 2012 (modified from IPCC 2013, p. 6). Bellwether sites are shown with white circle and are labeled. Observations in white grid cells were insufficient to estimate changes. Plus (+) signs show areas where observed trends are statistically significant

Fig. 4
figure 4

Projected changes in median values (i.e., the 50th percentile from a distribution of multi-model simulation results) of mean air temperature for the period 2046–2065 compared to 1986–2005 using the RCP 4.5 emissions scenario (modified from van Oldenborgh et al. 2013). Simulation results from the Coupled Model Intercomparison Project 5 (CMIP5). Seasonal periods shown are December–February (top) and June–August (bottom)

Below, I characterize each of the five bellwether sites in alphabetical order by country name. At the end of this section, I briefly describe other sites that were considered. I hope that readers of this article will have suggestions for additional sites. My lists of sites are not meant to be absolute and unalterable. Rather, they are simply intended as a starting point and discussion starter.

Eastern and northern flanks of the Northern Patagonia Ice Field, Chile

In South America, an overall warming trend has been observed over the entire continent between 1901 and 2012 (Fig. 3). Since 1975, the largest magnitudes of 1 °C per decade have been observed in central Brazil (e.g., de Barros Soares et al. 2017). For the same period, in the Andes Mountains of southern Chile (i.e., Patagonia), observed warming has not been statistically significant (e.g., Carrasco et al. 2008; de Barros Soares et al. 2017), but valley glaciers on the flanks of the two major icefields (the Northern Patagonia Ice Field (NPI) and Southern Patagonia Ice Field (SPI)) have been rapidly receding and thinning (e.g., Davies and Glasser 2012; Willis et al. 2012). Combined, the NPI and SPI make up the largest extent of ice in the Southern Hemisphere outside of Antarctica, and their ongoing melting has contributed a disproportionate amount (9%) to sea-level rise compared to larger icefields in other parts of the world (Ringot et al. 2003). In agreement with these trends, Wilson et al. (2018) noted that within the Central and Southern Andes between 1986 and 2016, Patagonia had the highest number of glacial lakes and reported glacial-lake outburst floods (GLOFs) (e.g., Jacquet et al. 2017), as well as an increased rate of glacial-lake emergence since 2000. For the 2046–2065 time period, projections of air temperatures using the mid-low-range RCP 4.5 emissions scenario indicate that average winter and summer temperatures will be 0.5 to 1.5 °C warmer than they were between 1986 and 2005 (Fig. 4; also see van Oldenborgh 2013, pp. 1346–1347, for higher-resolution details).

I selected a bellwether site along the eastern and northern flanks of the rapidly deglaciating NPI (Figs. 2 and 5) because glaciers are rapidly receding and thinning there (e.g., Davies and Glasser 2012; Willis et al. 2012) and because I observed about six large (areas > 0.2 km2), recent rock avalanches in glacier ablation areas in post-2010 satellite imagery (e.g., Fig. 5b). Some of these rock avalanches were previously identified as rock-fall deposits (Glasser et al. 2016). The bellwether site is about 220 km east of the triple junction of the Nazca, Antarctic, and South American tectonic plates (Fig. 2). This triple junction has subduction (compressive) and transform components, with the major Luquiñe-Ofqui right lateral fault zone running through the northwest part of the NPI (e.g., Georgieva et al. 2016). Georgieva et al. (2016) describe both strike-slip faults and thrust faults within the bellwether site and document exhumation rates of 0.1 to 0.9 mm/year. Dietrich et al. (2010) used campaign-style GPS between 2003 and 2006 to measure vertical crustal uplift rates (from combined isostatic and tectonic deformation) of about 16 mm/year south of the bellwether site. The western and central parts of the site (about 80% of the total area) are underlain by granites, granodiorites, and tonalites of the Cretaceous South Patagonia batholith (SERNAGEOMIN 2003). The eastern part of the site is underlain by Paleozoic metasediments (meta-sandstones and phyllites) and Jurassic monzodiorites, diorites, and granodiorites (SERNAGEOMIN 2003).

Fig. 5
figure 5

a Bellwether site (outlined with a solid red line) along the eastern and northern flanks of the Northern Patagonia Ice Field, southern Chile. The October 2017 rock avalanche was described by Fernandez and Coloma (2018). Imagery from multiple sources available in ArcMap is shown. b Two rock avalanches on the Glaciar Leones about 5 km southwest of Lago Leones on the east flank of the Northern Patagonia Ice Field shown in a DigitalGlobe image acquired on January 29, 2019. Images in the DigitalGlobe catalog revealed that the rock avalanche on the right occurred first, between July 25 and October 4, 2014, and the rock avalanche on the left occurred second, between November 8, 2016, and January 12, 2017. The first avalanche was noted in an AGU blog posting by Pelto (2015) and by Glasser et al. (2016). The rock avalanches are herein named the Leones West rock avalanche (left) and Leones East rock avalanche (right)

Chile is susceptible to landslides induced by climatic, volcanic, and seismic triggers (e.g., Francis et al. 1985; Hauser 2002; Nadim et al. 2006; Hermanns et al. 2012; Mather et al. 2014; Alfaro et al. 2017; Korup et al. 2019; Serey et al. 2019). However, few studies have been done within and near the bellwether site in southern Chile. Within the site, Fernandez and Coloma (2018) documented a 200,000-m3 rockslide that mobilized as a rock avalanche and closed a road (for about 4 months) used by the tourist industry in the Exploradores Valley on the northern flank of the icefield (Fig. 5a). The rockslide occurred near the edge of a retreating glacier on October 27, 2017 (spring in the Southern Hemisphere), during an extended period (less than 1 month, but the exact length is unknown) of higher-than-normal air temperatures (Fernandez and Coloma 2018). About 200 km north of the bellwether site, landslides in Aysén Fjord triggered by a 2007 earthquake have been recently researched (e.g., Sepúlveda et al. 2010; Glüer 2011; Oppikofer et al. 2012). Also, more recently, Duhart et al. (2019) documented a 2017 rockslide and debris flow that destroyed part of the village of Santa Lucía. This rockslide began in altered volcanic rocks in an area of a retreating glacier at an altitude of about 1400 m above sea level, slid onto a glacier and ice-cored moraine, and then mobilized as a debris flow that traveled a total distance of 12 km. Duhart et al. (2019) attribute the initiation of the landslide to a 2-week period when maximum daily temperatures exceeded 22 °C, as well as 123 mm of rainfall in the 24-h period before to the rockslide. Work by Oppikofer et al. (2012) near Aysén Fjord suggests that landslides within the bellwether site likely have a maximum source area volume of about 22 × 106 m3.

Western European Alps, France, Switzerland, Italy

In the European Alps, one of the most densely populated mountainous regions in the world (Haeberli and Beniston 1998), an overall warming trend has been observed between 1901 and 2012, with temperature changes of about 1.8 °C in the Southwestern Alps in the border region between France, Switzerland, and Italy (Fig. 3). Observations in Switzerland indicate that mean warming between 1975 and 2004 was 0.57 °C per decade, compared to 0.14 °C per decade between 1901 and 2000 (e.g., Rebetez and Reinhard 2008). During the late twentieth century, from 1985 to 1999, glaciers throughout the Alps decreased in size at a rate that was 7 times greater than the rate between 1850 and 1973 (Paul et al. 2004). Statistical modeling by Zemp et al. (2006) indicated that glaciers would be reduced another 80% (from their size in 2000) with a 3 °C increase in summer air temperature, even with a 10% increase in annual precipitation. Furthermore, Zemp et al. (2006) pointed out that this scenario could take place by the mid-to-late twenty-first century. More recent modeling that incorporated locally specific ice dynamics (Zekollari et al. 2019) confirmed these earlier results and indicated glacier losses of 79 ± 9% by 2100 using the RCP 4.5 emissions scenario. For the 2046–2065 time period, projections of air temperatures using the RCP 4.5 emissions scenario indicate that average winter and summer air temperatures will be 1 to 3 °C warmer than they were between 1986 and 2005 (Fig. 4; also see van Oldenborgh 2013, pp. 1350–1351, for higher-resolution details).

I selected a bellwether site in the Western European Alps (Figs. 2 and 6) because glaciers are rapidly shrinking there (e.g., Zemp et al. 2006; Salzmann et al. 2012) and permafrost is likely degrading (e.g., Harris et al. 2009) and because of an abundance of well-documented recent rock fall, rock avalanche, and debris flow activity that has been attributed to glacier withdrawal and/or degrading permafrost (e.g., Fig. 6b; Gruber et al. 2004; Chiarle et al. 2007; Deline et al. 2011; Fischer et al. 2006, 2012). The Western Alps formed during the Oligocene and Miocene epochs from convergence and collision of the Adriatic micro tectonic plate with the Eurasian Plate (e.g., Nocquet et al. 2016). This collision resulted in crustal thickening, development of extensive thrust sheets, and exhumation of basement rocks (e.g., Nocquet et al. 2016). Currently, no horizontal tectonic convergence occurs in the Western Alps (Fox et al. 2015; Nocquet et al. 2016), but Nocquet et al. (2016) estimate uplift rates of ~ 0.25 mm/year, which they attribute to glacial erosion, isostatic rebound, and deep crustal dynamics. The bedrock geology underlying the bellwether site is complex and includes Carboniferous granites, Mesozoic marine limestones, marls and shales, and Cretaceous flysches, ophiolites, amphibolites, gneisses, and schists (e.g., Lacassin et al. 1990; Leloup et al. 2005; Bousquet et al. 2012).

Fig. 6
figure 6

a Bellwether site (outlined with a solid red line) in the Western European Alps of France, Switzerland, and Italy. Imagery from multiple sources available in ArcMap is shown. b December 24, 2008, ~ 0.5 Mm3 Mont Crammont rock avalanche in the Western European Alps (described by Deline et al. 2011) shown in a DigitalGlobe image from August 28, 2009. The source area was in the steep rock cliff in the lower left corner of the image

Landslides linked to climate change in the Western Alps are the best documented in the world. Evidence for this statement comes from a multitude of papers on rock falls, rockslides, rock avalanches, and debris flows published in the twentieth and twenty-first centuries, which are too abundant to cite here. Instead, below, I try to mention the most relevant publications needed to provide support for selection of the bellwether site. Three of the best-known landslide study areas in the Alps, Mont Blanc, Monte Rosa, and the Matterhorn are located within the bellwether site (Fig. 6). At Mont Blanc, work on recent rock falls and rock avalanches has been particularly prominent (e.g., Ravanel et al. 2010; Ravanel and Deline 2011; Deline et al. 2011; Magnin et al. 2017). This work has shown that (1) mountain permafrost is highly sensitive to air temperatures; (2) the permafrost limit has been increasing in altitude, especially since the late 1990s; and (3) slope failures in rock occurred in response to air temperature increases and permafrost degradation. At Monte Rosa, research on the east face (e.g., Haeberli et al. 2002; Fischer et al. 2006, 2013) has documented enhanced ice loss in the late twentieth and early twenty-first centuries and showed that recent rock falls, debris flows, and rock-ice avalanches have initiated from these areas. At the Matterhorn, in situ mountain permafrost monitoring campaigns began on the Swiss (e.g., Hasler et al. 2012; Weber et al. 2019) and Italian sides of the peak (e.g., Chiarle et al. 2015; Coviello et al. 2015) following the summer of 2003 extreme heat wave that triggered rock falls throughout the Alps (Gruber et al. 2004). Thus far, results from monitoring along a highly fractured ridgeline on the Swiss side have shown that for most of the year, fracture aperture was inversely related to temperature, except for a period during summer months when fracture aperture increased and shearing took place. This time period began with sustained melting and ended with the first freeze in the fall, suggesting the water reduced the strength in fractures containing perennial ice (Hasler et al. 2012; Weber et al. 2017). Another well-known landslide within the bellwether site is the 1991 Randa rockslide on the flank of a formerly glaciated valley that lacks permafrost (e.g., Gischig et al. 2011a, b). The studies at Randa suggest that landslides within the bellwether site have a maximum source area volume of about 22 × 106 m3.

Eastern Karakoram Range, India and Pakistan

In the Himalayan Mountains of Asia, an overall warming trend has been observed in the twentieth century, with temperature changes of 1.6 to 1.8 °C in the Karakoram Range in Pakistan and India (Bhutiyani et al. 2007; Fig. 3). Recent observations from weather stations in northern Pakistan indicate that mean warming between 1986 and 2015 was 0.25 °C per decade (Waqas and Athar 2019). Most glaciers in the Himalayas have receded during this time period (Kääb et al. 2012). In contrast, glaciers in the central Karakoram Range began expanding in the 1990s (Hewitt 2005). Hewitt (2005) named this pattern the “Karakoram anomaly” and underscored the need for better monitoring of changing glaciers and climate in the area. More recent work indicates that the growth anomaly is still ongoing (e.g., Kääb et al. 2012; Farinotti et al. 2020; Paul 2020) and that it may be caused by summer cooling in the region caused by large-scale atmospheric circulation conditions (Forsythe et al. 2017). For the 2046–2065 time period, projections of air temperatures using the mid-low-range RCP 4.5 emissions scenario indicate that average winter and summer air temperatures will be 1.5 to 3 °C warmer than they were between 1986 and 2005 (Fig. 4; also see van Oldenborgh 2013, pp. 1366–1367, for higher-resolution details).

I selected a bellwether site in the eastern Karakoram Range along the eastern edge of the Karakoram anomaly, roughly centered on Siachen and Baltoro glaciers (Figs. 2 and 7), with lengths of 70 km and 67 km, respectively. The site contains some of the highest peaks on Earth including K2 and Gasherbrum I, the 2nd and 11th highest at 8611 m and 8080 m above sea level, respectively. I selected this site because it provides an extreme elevation, end-member site where glaciers are currently not rapidly receding; monitoring in the area is sparse (Hewitt 2005); and rock avalanches have recently occurred (Fig. 7b; Petley 2012; Schneider et al. 2013; Bhutiyani and Mahto 2018; Berthier and Brun 2019).

Fig. 7
figure 7

a Bellwether site (outlined with a solid red line) in the Karakoram Range, Himalayan Mountains of India and Pakistan. Imagery from multiple sources available in ArcMap is shown. The approximate position of the Karakoram fault (from Searle and Phillips 2007) is shown by the thick dashed line. The 2012 ice-rock avalanche was described by Petley (2012). b September 2010 rock avalanche on the east side of the Siachen Glacier in the eastern Karakoram Range with an estimated volume of ~ 109 Mm3 (described by Berthier and Brun 2019). DigitalGlobe image acquired on October 9, 2010

The bellwether site is at the boundary between the Indian and Eurasian tectonic plates (Fig. 2). At this location, the plate boundary has compressive and transform components, with the major Karakoram right lateral strike-slip fault running through the central part of the bellwether site (e.g., Nakata 1989; Taylor and Yin 2009). The slip rate for the Karakoram fault, as well as vertical uplift in the Karakoram region, is poorly constrained (Taylor and Yin 2009; Wallis et al. 2016). Wallis et al. (2016) modeled average erosion exhumation rates of about 0.67 + 0.27/− 0.18 mm/year for the last 7.4 Ma and attributed variations in these rates to variations in vertical-rock uplift driven by tectonic-driven crustal thickening. The site is underlain by highly deformed plutonic and metamorphic rocks (Searle 1991; Fraser et al. 2001; Searle et al. 2010) including Tertiary Baltoro granites and K2 and Muztagh Tower gneisses, diorites, and schists.

The Karakoram Range is susceptible to landslides induced by climate and earthquake triggers (e.g., Khattak et al. 2010; Berthier and Brun 2019). Ken Hewitt of Wilfrid Laurier University has been studying rock avalanches and their impact on glaciers in the central Karakoram (west of the bellwether site) for more than 30 years (e.g., Hewitt 1988, 2009) and has repeatedly stressed the need for additional research and monitoring of mass movements in the area. I chose a bellwether site mostly to the east of Hewitt’s study area in the central Karakoram because it is generally at a higher elevation and has experienced ice loss between 2008 and 2010 and between 2014 and 2016 (Berthier and Brun 2019) and three large recent rock-avalanche events have been documented there: one in April 2012 (Petley 2012; Schneider et al. 2013), one in September 2010 on the Siachen Glacier (Fig. 7b; Ekström and Stark 2013; Berthier and Brun 2019), and another in April or May 2000 on the North Terong Glacier (Bhutiyani and Mahto 2018). The rock avalanche on April 7, 2012, was an ice-rock avalanche that killed ~ 139 people in a Pakistani military camp on a debris fan (Petley 2012; Schnieder et al. 2013). The trigger for the event is unknown. The rock avalanche on the Siachen Glacier was a series of seven rock-slope failures from September 9 to 12 (Ekström and Stark 2013). The cumulative landslide source area volume for all of the failures was estimated at 107 Mm3 by Berthier and Brun (2019). Additionally, Berthier and Brun (2019) point out that the cumulative 107 Mm3 volume is larger than all rock avalanches in the worldwide catalog of rock avalanches onto glaciers (Table 9.1 in Deline et al. 2015b). The rock-avalanche deposit enhanced melting of the Siachen Glacier, contributing to a 6% loss from 2010 to 2016 (Berthier and Brun 2019). The rock avalanche on the North Terong Glacier had an estimated volume of 7.6 Mm3 and may have occurred during a period of anomalously warm temperatures from May 3 to May 4, 2000 (Bhutiyani and Mahto 2018)

Aoraki/Mount Cook area, Southern Alps, New Zealand

In New Zealand, observations indicate warming air temperatures in the twentieth and twenty-first centuries. An analysis of air temperature data from seven low-altitude stations on the North and South Islands found a 0.9 to 1 °C increase for the century between 1909 and 2009 (Mullan et al. 2010; Mullan 2012). A similar estimate from the IPCC (2013) indicates an overall warming trend between 1901 and 2012 of about 1.3 °C (Fig. 3). On the west side of the Southern Alps, Salinger et al. (1995) documented warming of 0.7 °C between 1920 and 1990. For most of the twentieth and twenty-first centuries, glaciers in the Southern Alps have lost mass and retreated (e.g., Blair 1994; Chinn 1999, 2001; Hoelzle et al. 2007), with overall volume loss from ~ 1850 to 2007 estimated at about 61% (Hoelzle et al. 2007). However, for an approximately two-decade-long period from the late 1970s to the late 1990s, the volume of many large glaciers increased, while others, particularly debris-covered glaciers, continued to recede (Chinn 1999). Hoelzle et al. (2007) note that glaciers in the Southern Alps will probably be highly sensitive to future temperature increases because of their relatively low equilibrium-line altitudes (ELAs) and an increased likelihood of rain, as opposed to snow, in the future. For the 2046–2065 time period, projections of air temperatures using the mid-low-range RCP 4.5 emissions scenario indicate that average winter and summer air temperatures will be 0.5 to 1.5 °C warmer than they were between 1986 and 2005 (Fig. 4; also see van Oldenborgh 2013, pp. 1382–1383, for higher-resolution details).

I selected a bellwether site in the Southern Alps, centered on Aoraki/Mount Cook and encompassing the high-altitude, east and west sides of the north-south trending main divide (Figs. 2 and 8). I chose the site because glaciers are receding there (e.g., Hoelzle et al. 2007) and because historical slope failures in the area have been numerous, with indications that the frequency and size of events is increasing (McSaveney et al. 2014), although interpretations as to why this might be happening differ (McSaveney et al. 2014, 2015). The western edge of the bellwether site corresponds approximately with the oblique-reverse Alpine Fault (e.g., Cox and Barrell 2007; Little et al. 2005), which is the transpressive boundary between the Australian and Pacific tectonic plates (e.g., Walcott 1998; Fig. 2). Uplift occurs in the bellwether site on the east side of the fault. The estimated average fault slip rates for dextral and vertical components are ~ 27 mm/year and ~ 8–12 mm/year, respectively (Little et al. 2005). The site is underlain by variably metamorphosed and deformed Permian and Triassic sandstone with interbedded siltstone and mudstone (Cox and Barrell 2007). The metamorphic grade and internal deformation increase from east to west, with rocks along the eastern edge of the site being non-schistose cleaved sandstone and rocks on the western edge of the site at the Alpine Fault being highly foliated schist. On the west side of the bellwether site, west of the main divide, high rainfall rates result in abundant slope failures and erosion rates that are roughly equal to the rate of uplift (Koons 1990; Hovius et al. 1997; Little et al. 2005).

Fig. 8
figure 8

a Bellwether site (outlined with a solid red line) in the Aoraki/Mt. Cook region of the Southern Alps, New Zealand. Imagery from multiple sources available in ArcMap is shown. Approximate position of the main divide shown with thick, black-dashed line. Approximate position of the Alpine Fault (from Cox and Barrell 2007) is shown with thick, white-dashed line, with teeth shown on upthrust side. b January 21, 2013, rock avalanche from a ridge connecting Mount Haast and Mount Dixon in the Southern Alps of New Zealand. Estimated volume is ~ 2 Mm3 (Dunning et al. 2015). DigitalGlobe image acquired on February 27, 2013

Landslides have been extensively researched within the bellwether site (e.g., Whitehouse and Griffiths 1983; McSaveney 2002; Cox and Allen 2009; Allen et al. 2011; Cox et al. 2015; Cody et al. 2020), with a relatively long documented history of rock falls and rock avalanches, particularly from areas along and near the main divide (e.g., McSaveney 2002). In the early 1980s, Whitehouse and Griffiths (1983) used an inventory of 42 rock avalanches from a 1700-year period of time to estimate a frequency for avalanches > 1 Mm3 at 1 per 94 years, with the largest in any century estimated at 56 Mm3. Based in part on a temporal cluster of rock avalanches in the 1990s, McSaveney (2002) suggested a frequency for rock avalanches > 1 Mm3 at 1 per 20–30 years. More recent work (McSaveney et al. 2014, 2015) suggests that the rate of rock avalanche occurrence since about 2004 is 2 per year. The reasons for this increase are unclear, with Allen et al. (2009, 2011) and Allen and Huggel (2013) tending to favor a mountain permafrost degradation explanation and McSaveney et al. (2014, 2015) favoring rock deformation from a global increase in seismic moment release and/or accumulating elastic strain since the last great earthquake on the Alpine Fault in ~ 1717 AD. All researchers seem to agree that a fundamental problem is that many rock avalanches in the remote region are not recorded, and therefore, the existing record is too incomplete to definitively determine the impact of climate change on the distribution and frequency of recent landslides (e.g., Cox and Allen 2009; Allen et al. 2011; Cox et al. 2015). Work by McSaveney (2002) suggests that landslides within the bellwether site have a maximum total volume (including entrained ice) of about 80 × 106 m3.

Fairweather Range, Southeast Alaska, USA

In the USA, the greatest increases in air temperature and precipitation from climate change (e.g., van Oldenborgh et al. 2013; Walsh et al. 2014; Wuebbles et al. 2014; USGCRP 2018), and thus the most dramatic adverse impacts to slope stability, are expected in Alaska (Coe 2016). Mean air temperatures have increased in Alaska by at least 1.25 °C since 1901 (Fig. 3). In 2014, Walsh et al. (2014) estimated an increase of 1.5 °C compared to the period between 1951 and 1980, with the most dramatic increases occurring during winter months. During a period from 2014 to 2016, mean air temperatures were the warmest on record (NOAA 2017; Walsh et al. 2017), and a cluster of exceptionally large rock avalanches occurred in the Fairweather Range in Glacier Bay National Park and Preserve (GBNPP) in Southeast Alaska (Coe et al. 2018), including the 70 Mm3 Lamplugh rock avalanche (Bessette-Kirton et al. 2018).

Since 1995, glaciers in GBNPP have lost mass overall, but the mass balance is complicated locally, with some glaciers losing mass and others gaining mass (e.g., Johnson et al. 2013; Alifu et al. 2016). For the 2046–2065 time period, projections of air temperatures using the mid-low-range RCP 4.5 emissions scenario indicate that average temperatures in Southeast Alaska will be 1.5 to 3 °C warmer than they were between 1986 and 2005 (Fig. 4; also see van Oldenborgh 2013, pp. 1330–1331, for higher-resolution details).

I selected a bellwether site in the Fairweather Range in GBNPP (Figs. 2 and 9) because of its recent documented history of climatically induced rock avalanches (Fig. 9b; Bessette-Kirton and Coe 2016; Coe et al. 2018; Bessette-Kirton et al. 2018). The site is centered on a glaciated peninsula along the transform boundary of the Pacific and North American tectonic plates (Fig. 2). The plate boundary is marked by the strike-slip Fairweather Fault (Pflaker et al. 1978) which accommodates a large portion of the relative motion between the plates through a dextral slip of about ~ 43 mm/year (Elliott et al. 2010). Elevations in the area range from sea level up to 4670 m above sea level at Mount Fairweather at the U.S.-Canadian border. Much of the study area is covered by ice, either in the form of valley glaciers or from the Brady icefield. In the last 250 years, the area has been undergoing rapid uplift (16 to 30 mm/year) caused primarily by viscoelastic rebound from post-Little Ice Age (LIA) deglaciation (Larsen et al. 2005). Glacier Bay proper (a fjord, Fig. 9a) is a result of post-LIA deglaciation (e.g., Connor et al. 2009).

Fig. 9
figure 9

a Bellwether site (outlined with a solid red line) in the Fairweather Range, Glacier Bay National Park and Preserve, Southeast Alaska, USA. Imagery from multiple sources available in ArcMap is shown. Approximate position of the Fairweather Fault (from Pflaker et al. 1978) is shown with thick, black-dashed line. b June 28, 2016, Lamplugh rock avalanche in the Fairweather Range. Measured volume from Bessette-Kirton et al. (2018) is 69.9 Mm3 (also described by Bessette-Kirton (2017) and Dufresne et al. (2019)). DigitalGlobe image acquired on July 16, 2016

Bedrock in GBNPP is exposed at mountain peaks, flanks, and ridge lines and in coastal zones. The bedrock geology is complex and consists of Paleozoic and Mesozoic accretionary terranes containing Tertiary sedimentary, intrusive, and volcanic rocks (Brew et al. 1978; Wilson et al. 2015). Bedrock units that surround Glacier Bay are Tertiary, Oligocene, and Eocene granitic rocks; the Cretaceous Chugach accretionary complex consisting of flysch, graywacke, and basalts; Cretaceous and Jurassic quartz monzodiorite; Cretaceous foliated granitic rocks; and Devonian to Ordovician shale, chert, and argillite (Wilson et al. 2015).

Glacier Bay National Park and Preserve has a history of both earthquake and climate-induced subaerial rockslides and rock avalanches (Miller 1960; Post 1967; Evans and Clague 1999; Geertsema 2012; Bessette-Kirton and Coe 2016; Bessette-Kirton et al. 2018; Coe et al. 2018; Dufresne et al. 2019). Because of the site’s proximity adjacent to fjords, when these landslides rapidly enter the water, they can generate tsunamis with extraordinary destructive power. This destructive power was demonstrated by both the 1958 Lituya Bay earthquake-triggered rockslide (Miller 1960), which generated the largest tsunami runup height in history (524 m; Higmann et al. 2018). GBNPP’s close proximity to the Fairweather Fault guarantees that large earthquakes will impact the area and trigger landslides in the future. However, since the M 7.8 earthquake that triggered the 1958 Lituya Bay rockslide and tsunami, at least 90 M > 4 earthquakes have occurred within 100 km of GBNPP (the largest were two M 6.1 earthquakes: one in January 2000 and the other in July 2014). Coe et al. (2019) did not find any evidence that these earthquakes triggered rock-slope failures in GBNPP. Instead, all recent rockslides and rock avalanches in GBNPP have been associated with climatic triggers (i.e., most likely warm air temperatures; Coe et al. 2018). Between 1984 and 2016, at least 24 climatically induced rock avalanches occurred in GBNPP (Bessette-Kirton and Coe 2016) with a mean recurrence interval of 505 days (Coe et al. 2018). Available evidence indicates that these avalanches have increased in size and travel distance with time. Work by Coe et al. (2019) suggests that landslides within the bellwether site have a maximum source area volume of about 80 × 106 m3.

Other sites considered as bellwether sites

The five sites that I have selected are not the only sites where long-term observations could take place. These five sites are only meant as a starting point. Other sites that I considered as bellwether sites (Table 2) were the Coast Mountains of southwest British Columbia, Canada; the central west coast of Greenland; the Tröllaskagi Peninsula in northern Iceland; the Caucus mountain range at the border between Georgia and Russia; the St. Elias Mountains in southern Alaska and the Yukon Territory; the Cordillera Blanca in Peru; and many sites in mountainous areas of Asia. Recent publications from each of these areas suggest that slope failures may be occurring in response to climate change (Table 2).

Table 2 Other sites considered as bellwether sites, from level 2 of Fig. 1

Discussion

Possible approaches for making systematic observations

As described above, landslide activity levels at the five bellwether sites are variable, with the Southern Alps, Fairweather Range, and the European Alps being most active; Northern Patagonia Ice Field being moderately active; and the Karakoram being the least active. Approaches for making long-term observations will need to account for differences in landslide sizes at the sites and ideally to establish minimum size detection thresholds early in the overall observation process.

A guide for how to proceed with long-term landslide observations could be taken from studies of glaciers. Glaciers are sensitive to climate change, and in most parts of the world, glaciers have receded since about 1985 (e.g., Paul et al. 2004). This period correlates well with the beginning of worldwide 30-m coverage by Landsat (the launch of Landsat 4 in 1982). As a first step, at each of the sites, landslides that are visible in Landsat imagery back to 1982 could be compiled. In Alaska, Landsat imagery provided a consistent medium that allowed Bessette-Kirton and Coe (2016) to systematically detect and map avalanches that covered areas larger than about 10 ha (100,000 m2, 0.1 km2) from 1984 to 2016. This approach would also be viable in the Southern Alps, the Northern Patagonia Ice Field, and the Karakoram, but possibly not in the European Alps, where most recent documented landslides have generally been smaller than the other four sites. In the European Alps, higher resolution imagery and/or redundant systems such as a dense seismic network would probably be needed to capture the relatively small landslides that have recently been documented (e.g., Ravanel et al. 2010).

For seismic monitoring to be used routinely for detection of small- and medium-sized landslide events, several things would need to happen. First, denser networks of seismometers would need to be installed in remote, difficult to access places, and second, better automatic detection algorithms would need to be developed to better distinguish landslide signals from signals generated from other types of events, such as glacier movement (e.g., Caplan-Auerbach and Huggel 2007; Faillettaz et al. 2008; Walter et al. 2008). Of course, the installation and maintenance of a high density of seismometers in rugged mountainous regions is probably not economically feasible or practical in areas of 5000 km2. However, recent work has shown that existing regional seismic networks can be very useful for establishing the precise timing and evolution of large landslide events (e.g., Guthrie et al. 2012; Allstadt 2013; Ekström and Stark 2013; Iverson et al. 2015; Coe et al. 2016; Dufresne et al. 2019). This approach could be helpful at the bellwether sites, particularly during times with extensive cloud cover for long periods of time. But, at least for the foreseeable future, satellite imagery will probably be the primary tool that is routinely available for systematic landslide detection and mapping (for a wide range of landslide sizes) in the bellwether sites.

Landslide inventory mapping from repeat imagery with a consistent resolution and acquired over long (> 30-year) time frames is the essential data collection element at each bellwether site. The consistent resolution is essential to avoid biasing collected data toward more landslides being mapped from higher-resolution imagery. Inevitably, as the resolution of imagery continues to improve, the lower size threshold of landslides that can be detected will become smaller (for example, when changing from 30-m Landsat images to 5-m Planet Labs images). However, changing the resolution of imagery used for mapping should not occur in a haphazard manner. Rather, a consistent resolution should be used for as long as possible. Then, when the imagery used for mapping changes, distinct breaks should be noted in the record of each bellwether site so that the change can be acknowledged and accounted for when analyzing for changes in the frequency and magnitude in the time series of mapped landslides. At a minimum, data collected for each mapped landslide should include a date of occurrence or range in possible dates of occurrence, time of occurrence (if available), type of landslide trigger, landslide polygon, landslide area, slope direction (aspect) of source area, curvilinear travel distance from headscarp to toe, headscarp elevation, toe elevation, and a height/length calculation. Bessette-Kirton and Coe (2016) provide a possible template for items that should be recorded. Determining a specific type of trigger may not be possible, but at a minimum, there should be a binary answer (yes or no) as to whether an earthquake was involved, as well as measured values of air temperature and precipitation for antecedent time periods of 1 day, 2 days, 1 week, 1 month, 3 months, 6 months, and possibly, 12 months. Earthquake records from regional or local seismic networks will need to be searched for possible earthquakes at the time that landslides occurred. One possible search method is given in Coe et al. (2018). Volume estimates would also be helpful, but they are typically difficult to accurately measure. If accurate (± 20%) volumes are not possible, they should not be reported.

Overall, research at the bellwether sites should focus on interconnected processes responsible for slope-stability changes in degrading cryospheric mountainous terrain. For establishing links between landslide characteristics (e.g., time of occurrence and size) and warming temperatures, as well as for gaining a better understanding of the physical processes controlling observed slope failures, long-term, in situ instrumental monitoring at high elevations in each bellwether will be required. An archetypical example of the type of in situ monitoring ideally needed is provided in a recent paper by Weber et al. (2019). However, at a minimum, in situ monitoring should consist of automated instrumentation to measure air temperature and precipitation.

A couple reasonable questions for landslide researchers to ask about this proposed work are who should do the work; how should long-term landslide observations at bellwether sites be funded, especially if commercial imagery needs to be purchased; and what group, if any, should be responsible for coordinating the efforts? In writing this paper, I have tried to provide a thought-provoking proposal for Theme 7, “Frontiers in landslide science and innovative practices,” of the 5th WLF in Kyoto, Japan (Sassa 2019a). I have tried to tailor the paper to fit well within Theme 7 by presenting a possible landslide research frontier, but with the full realization, those answers to the questions listed above still need to be determined. My intent with the paper was to provide the scientific basis and justification for bellwether sites, followed by a discussion of the proposal with other members of the international landslide community at the 5th WLF. My hope is that these discussions could lead to a clear path forward in terms of work distribution, coordination, and funding opportunities. To partly address these questions, in March 2020, I consulted with Duncan Quincey and Marta Chiarle, two of the officers of GAPHAZ (https://www.gaphaz.org/), an international scientific standing committee concerned with glacier and permafrost hazards in mountains, about the level of interest in the work, as well as possible collaborators and funding sources. After this consultation, my current opinion is that work at bellwether sites would best be done by landslide scientists and students with vested interest and expertise in regions where the proposed bellwether sites are located. The level of overall project coordination that will be needed, and what group (if any) should have that responsibility, remains to be determined and could be a topic of discussion at the 5th WLF. However, some possibilities include having the bellwether sites function in a manner similar to the way that U.S. National Science Foundation Critical Zone Observatories (CZOs) do for watersheds within the USA (https://criticalzone.org/national/infrastructure/observatories-1national/), or perhaps the sites could become part of the network of GEO Geohazard Supersites (https://geo-gsnl.org/) that operate throughout the world. Additionally, the proposed, systematic, long-term approach to understanding interconnected processes responsible for slope-stability changes at bellwether sites fits well with the goals of the Mountain Research Initiative (https://www.mountainresearchinitiative.org/) and their Group on Earth Observatories Global Network for Observations and Information in Mountain Environments (GEO GNOME) (https://mountainresearchinitiative.org/activities/projects/geo-gnome) program.

Furthermore, in terms of funding, one possible source could be directly related to the 5th WLF, a meeting with a focus on strengthening worldwide landslide collaboration and partnerships, particularly between the International Consortium on Landslides (Sassa 2019b), the United Nations University (OKI 2019), and the United Nations Educational, Scientific and Cultural Organization (UNESCO) (Nair-Bedouelle 2020). A study of climate change and its impact on slope stability in cryospheric mountainous regions throughout the world seems ideally suited for funding by such groups. Partnerships and collaborations between governments would also be beneficial to securing funding, establishing common mapping approaches, and integrating results.

Conclusions

The diversity of possible landslide responses to climate change (e.g., Gariano and Guzzetti 2016) places the landslide community in somewhat of a quandary as to the best way to proceed with research on the topic. In this paper, I have presented a proposal that calls on the international community of landslide researchers to use high-resolution satellite imagery, as well as seismic networks where possible, to systematically monitor five 5000-km2 cryospheric mountainous regions throughout the world. The five sites are the northern and eastern flanks of the Northern Patagonia Ice Field in southern Chile, the Western European Alps, the eastern Karakoram Range in the Himalayan Mountains, the Southern Alps of New Zealand, and the Fairweather Range in Southeast Alaska, USA. Some of these sites will be readily recognizable to landslide researchers who are already doing work related to climate change in these areas. I suggest that successful long-term observations at these “bellwether” sites will address a concern expressed in a recent Intergovernmental Panel on Climate Change special report (Hock et al. 2020) regarding a lack of studies to characterize and compare the magnitude of landslides in mountainous regions, and provide an indication of slope-stability changes at similar sites in other parts of the world.