From a geological standpoint, northern Pakistan is one of the most active and unstable areas in the world. As a consequence, many massive landslides have occurred in the area in historical times that have destroyed infrastructure, blocked the Hunza River, and damaged the Karakoram Highway repeatedly. However, despite the high frequency of large magnitude landslide events, and the consequent damages, the entire area is largely understudied, mainly due to the difficult logistics and the large distances involved. This work is aimed at applying the potential use of Interferometric Synthetic Aperture Radar (InSAR) for landslide identification and investigation for the Hunza-Nagar Region. Sentinel-1 images covering a period of more than two years (February 2017–August 2019) were used and processed by adopting the small baseline subset (SBAS) method. The obtained deformation rate measured along the line of sight (VLOS) varies from −114 to 20 mm/year. The downslope velocity deformation rates (Vslope) range from 0 to −300 mm/year. The Vslope stability threshold for our study area was calculated to be −14 mm/year from the Vslope standard deviation. Four active landslides with Vslope exceeding 14 mm/year were recognizable and have been confirmed by field inspection. The identified landslides listed from the most active to least active are the Humarri, Mayoon, Khai, and Ghulmet landslides, respectively. VLOS exceeding 114 mm/year was observed in the Humarri landslide, which posed a threat of damming a lake on the Hispar River and was also a risk to the Humarri Village located below the landslide. The maximum mean deformation detected in the Ghulmet, and Mayoon landslide was in the order of 30 mm/year and 20 mm/year, respectively. More importantly, it was found that in places, the slope deformation time series showed a patchy correlation with precipitation and seismic events in the area. This may indicate a complex, and possibly uncoupled, relationship between the two controlling agents promoting the deformation. However, the collective impact of the two factors is evident in the form of a continuously descending deformation curve and clearly indicates the ground distortion. The results indicate a potentially critical situation related to the high deformation rates measured at the Humarri landslide. On this specific slope, conditions leading to a possible catastrophic failure cannot be ruled out and should be a priority for the application of mitigation measures.

No abstract

The Hunza Valley, in the northwestern Karakoram Mountains, North Pakistan, is a typical region with many towns and villages, and a dense population and is prone to landslides. The present study completed landslide identification, updating a comprehensive landslide inventory and analysis. First, the ground surface deformation was detected in the Hunza Valley by SBAS-InSAR from ascending and descending datasets, respectively. Then, the locations and boundaries were interpreted and delineated, and a comprehensive inventory of 118 landslides, including the 53 most recent InSAR identified active landslides and 65 landslides cited from the literature, was completed. This study firstly named all 118 landslides, considering the demand for globally intensive research and hazard mitigation. Finally, the deformation, spatial–topographic development, and distribution characteristics in the Hunza Valley scale and three large significant landslides were analyzed. Information on 72 reported landslides was used to construct an empirical power law relationship linking landslide area (AL) to volume (VL) (VL = 0.067 × AL1.351), and this formula predicted the volume of 118 landslides in this study. We discovered that the landslides from the literature, which were interpreted from optical images, had lower levels of velocity, area, elevation, and height. The SBAS-InSAR-detected active landslide was characterized by higher velocity, larger area, higher elevation, larger slope gradient, larger NDVI (normalized difference vegetation index), and greater height. The melting glacier water and rainfall infiltration from cracks on the landslide’s upper part may promote the action of a push from gravity on the upper part. Simultaneously, the coupling of actions from river erosion and active tectonics could have an impact on the stability of the slope toe. The up-to-date comprehensive identification and understanding of the characteristics and mechanism of landslide development in this study provide a reference for the next step in landslide disaster prevention and risk assessment.