Back analysis of landslide susceptibility zonation mapping for the 2005 Kashmir earthquake: an assessment of the reliability of susceptibility zoning maps

The October 2005 earthquake triggered several thousand landslides in the Lesser Himalaya of Kashmir in northern Pakistan and India. Analyses of ASTER satellite imagery from 2001 were compared with a study undertaken in 2005; the results show the extent and nature of pre- and co-/post-seismic landsliding. Within a designated study area of ~2,250 km², the number of landslides increased from 369 in 2001 to 2,252 in October 2005. Assuming a balanced baseline landsliding frequency over the 4 years, most of the new landslides were likely triggered by the 2005 earthquake and its aftershocks. These landslides mainly happened in specific geologic formations, along faults, rivers and roads, and in shrubland/grassland and agricultural land. Preliminary results from repeat photographs from 2005 and 2006 after the snowmelt season reveal that much of the ongoing landsliding occurred along rivers and roads, and the extensive earthquake-induced fissuring. Although the susceptibility zoning success rate for 2001 was low, many of the co-/post-seismic landsliding in 2005 occurred in areas that had been defined as being potentially dangerous on the 2001 map. While most of the area in 2001 was (very) highly susceptible to future landsliding, most of the area in 2005 was only moderate to low susceptible, that is, most of the landsliding in 2005 actually occurred in the potentially dangerous areas on the 2001 map. This study supports the view that although susceptibility zoning maps represent a powerful tool in natural hazard management, caution is needed when developing and using such maps.     

国内外地震滑坡研究：现状、问题与展望

地震滑坡一直是地学界关注的研究热点之一,特别是2005年10月8日克什米尔M_s 7.8地震、2008年5月12日汶川M_s 8.0地震、2010年玉树M_s 7.1地震、2014年云南鲁甸M_s 6.5地震、2015年4月25日尼泊尔M_s 8.1地震以及2017年8月8日九寨沟M_s 7.0地震等一系列强震活动,都不同程度地触发了大量地震滑坡,并伴生了许多次生灾害,从而促进地震滑坡研究进入了新的阶段。本文通过总结国内外地震滑坡研究现状,归纳了地震滑坡的特征和分布规律、地震动参数与地震滑坡关系、边坡地质条件与地震滑坡关系、地震滑坡的动力响应特征、地震滑坡预测以及地震滑坡危险性与风险区划等6方面的研究现状,在此基础上进一步提出了未来地震滑坡应关注的主要研究方向和重点。     

同震滑坡发生概率研究——新一代地震滑坡危险性模型

地震滑坡发生真实概率研究基本空白。本研究创新性的利用贝叶斯概率方法与机器模型开展了中国地震滑坡危险性真实概率研究,制作了第一代中国地震滑坡危险性概率图。基于9个地震案例开展研究,包括1999年台湾集集、2005年克什米尔、2008年汶川、2010年玉树、2013年芦山、2013岷县、2014鲁甸、2015尼泊尔、2017九寨沟地震,这9次地震中7次发生在中国,2005年克什米尔与2015尼泊尔地震均发生在中国邻区,可以更好的控制模型预测精度。这些地震事件均有详细完整的,利用面要素标识的地震滑坡数据,包括306 435处真实的地震滑坡记录。考虑到真实的地震滑坡发生区域,滑坡面积规模的差别,滑坡与不滑样本的比例等因素,共选取了5 117 000个模型训练样本。选择绝对高程、相对高差、坡度、坡向、斜坡曲率、坡位、地形湿度指数、土地覆盖类型、植被覆盖度、与断层距离、地层、年均降水量、地震动峰值加速度共13个地震滑坡影响因子。采用贝叶斯概率方法与机器学习模型相结合,建立地震滑坡发生的多因素影响模型,得到各个连续因子的权重与分类因子的各个分类的权重。再将模型应用到整个中国研究区,地震动峰值加速度因子为触发因子。分别考虑研究区在经历不同地震动峰值加速度（0.1~1 g,每0.1 g一个结果,共10个结果）下的地震滑坡发生真实概率。此外,还结合中国地震动峰值加速度分布图,得到了中国地震动峰值加速度背景下的地震滑坡发生真实概率分布。     

Earthquake-induced landslide hazard and vegetation recovery assessment using remotely sensed data and a neural network-based classifier: a case study in central Taiwan

A catastrophic earthquake with a Richter magnitude of 7.3 occurred in the Chi-Chi area of Nantou County on 21 September 1999. Large-scale landslides were generated in the Chiufenershan area of Nantou County in central Taiwan. This study used a neural network-based classifier and the proposed NDVI-based quantitative index coupled with multitemporal SPOT images and digital elevation models (DEMs) for the assessment of long-term landscape changes and vegetation recovery conditions at the sites of these landslides. The analyzed results indicate that high accuracy of landslide mapping can be extracted using a neural network-based classifier, and the areas affected by these landslides have gradually been restored from 211.52 ha on 27 September 1999 to 113.71 ha on 11 March 2006, a reduction of 46.24%, after six and a half years of assessment. In accordance with topographic analysis at the sites of the landslides, the collapsed and deposited areas of the landslide were 100.54 and 110.98 ha, with corresponding debris volumes of 31,983,800 and 39,339,500 m³. Under natural vegetation succession, average vegetation recovery rate at the sites of the landslides reached 36.68% on 11 March 2006. The vegetation recovery conditions at the collapsed area (29.17%) are shown to be worse than at the deposited area (57.13%) due to topsoil removal and the steep slope, which can be verified based on the field survey. From 1999 to 2006, even though the landslide areas frequently suffered from the interference of typhoon strikes, the vegetation succession process at the sites of the landslides was still ongoing, which indicates that nature, itself, has the capability for strong vegetation recovery for the denudation sites. The analyzed results provide very useful information for decision-making and policy-planning in the landslide area.     

Characteristics of moderate- to large-scale landslides triggered by the <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> 7.8 2015 Gorkha earthquake and its aftershocks

The M _w 7.8 2015 Gorkha earthquake and its aftershocks significantly impacted the lives and economy of Nepal. The consequences of landslides included fatalities, property losses, blockades of river flow, and damage to infrastructural systems. Co-seismic landslides triggered by this earthquake were significantly widespread and pose a major geodisaster. There were tens of thousands of landslides triggered by the earthquake, majority of which were distributed in between the epicenter of the main shock and the M _w 7.3 aftershock. Although 14,670 landslides triggered by this earthquake were identified, only approximately 23% of them were of moderate to large scale with areas greater than 100 m². Of the moderate- to large-scale landslides identified, just over 90% were triggered by the main shock and smaller aftershocks prior to the major (M _w 7.3) aftershock, while nearly 10% were triggered by the ground shaking induced by the major aftershock. Moreover, the number of landslides triggered by the 2015 Gorkha earthquake, specifically by the main shock, was slightly more than the expected number of landslides for the recorded maximum peak ground acceleration (PGA) in comparison to the co-seismic landslides triggered by 26 earthquakes. Over 90% of those moderate- to large-scale landslides were concentrated within the estimated fault rupture surface. Majority of these moderate- to large-scale landslides were disrupted failures with over 96% of which were classified as earth falls. However, the majority of small-scale landslides were rock or boulder falls. The most number of moderate- to large-scale landslides were triggered in the slate, shale, siltstone, phyllite, and schist of the Lesser Himalayan formation followed by an equally significant number in both schist, gneiss, etc. of the Higher Himalayan formation and the phyllite, metasandstone, schist, etc. of the Lesser Himalayan formation. The sizes (i.e., areas) of the landslides were lognormally distributed, with a mode area of 322.0 m². Slope inclinations of the moderate- to large-scale landslides followed a normal distribution with a mean slope inclination of 32.6° and standard deviation of 13.5°. There exists a strong correlation between the number of landslides and the peak ground acceleration within the study area, specific for different geological formations.     

Landslides triggered by slipping-fault-generated earthquake on a plateau: an example of the 14 April 2010, Ms 7.1, Yushu, China earthquake

On 14 April 2010 at 07:49 (Beijing time), a catastrophic earthquake with Ms 7.1 struck Yushu County, Qinghai Province, China. A total of 2,036 landslides were interpreted from aerial photographs and satellite images, verified by selected field checking. These landslides cover about a total area of 1.194 km². The characteristics and failure mechanisms of these landslides are presented in this paper. The spatial distribution of the landslides is evidently strongly controlled by the locations of the main co-seismic surface fault ruptures. The landslides commonly occurred close together. Most of the landslides are small; there were only 275 individual landslide (13.5 % of the total number) surface areas larger than 1,000 m². The landslides are of various types. They are mainly shallow, disrupted landslides, but also include rock falls, deep-seated landslides, liquefaction-induced landslides, and compound landslides. Four types of factors are identified as contributing to failure along with the strong ground shaking: natural excavation of the toes of slopes, which mean erosion of the base of the slope, surface water infiltration into slopes, co-seismic fault slipping at landslide sites, and delayed occurrence of landslides due to snow melt or rainfall infiltration at sites where slopes were weakened by the co-seismic ground shaking. To analyze the spatial distribution of the landslides, the landslide area percentage (LAP) and landslide number density (LND) were compared with peak ground acceleration (PGA), distance from co-seismic main surface fault ruptures, elevation, slope gradient, slope aspect, and lithology. The results show landslide occurrence is strongly controlled by proximity to the main surface fault ruptures, with most landslides occurring within 2.5 km of such ruptures. There is no evident correlation between landslide occurrences and PGA. Both LAP and LND have strongly positive correlations with slope gradient, and additionally, sites at elevations between 3,800 and 4,000 m are relatively susceptible to landslide occurrence; as are slopes with northeast, east, and southeast slope aspects. Q₄ al-pl, N, and T₃ kn ¹ have more concentrated landslide activity than others. This paper provides a detailed inventory map of landslides triggered by the 2010 Yushu earthquake for future seismic landslide hazard analysis and also provides a study case of characteristics, failure mechanisms, and spatial distribution of landslides triggered by slipping-fault generated earthquake on a plateau.     

Climate extremes: an observation and projection of its impacts on food production in ASEAN

Multifractal behaviour of interevent time sequences is investigated for the earthquake events in the NW Himalaya, which is one of the most seismically active zones of India and experienced moderate to large damaging earthquakes in the past. In the present study, the multifractal detrended fluctuation analysis (MF-DFA) is used to understand the multifractal behaviour of the earthquake data. For this purpose, a complete and homogeneous earthquake catalogue of the period 1965–2013 with a magnitude of completeness M _w 4.3 is used. The analysis revealed the presence of multifractal behaviour and sharp changes near the occurrence of three earthquakes of magnitude (M _w ) greater than 6.6 including the October 2005, Muzaffarabad–Kashmir earthquake. The multifractal spectrum and related parameters are explored to understand the time dynamics and clustering of the events.     

Controlling parameter analyses and hazard mapping for earthquake-triggered landslides: an example from a square region in Beichuan County, Sichuan Province, China

On May 12, 2008, at 1428 hours (Beijing time), a catastrophic earthquake, with a magnitude of Ms 8.0, struck the Sichuan Province, China. About 200,000 landslides, as a secondary geological hazard associated with the earthquake, were triggered over a broad area. These landslides were of almost all types such as shallow, disrupted landslides, rock falls, deep-seated landslides, and rock avalanches. Some of these landslides damaged and destroyed large part of some towns, blocked roads, dammed rivers, and caused other serious damages. The purpose of this study is to detect correlations between landslide occurrence and the surface rupture plane, ground shaking conditions (measured by peak ground acceleration, PGA), lithology, slope gradient, slope aspect, topographic position, and distance from drainages by using two indices, landslide area percentage (LAP) and the landslide number density (LND), based on geographic information system (GIS) technology and statistical analysis method in a square region (study area) of Beichuan County, Sichuan Province, China. There were 5,096 landslides related with the earthquake which were delineated by visual interpretation and selected field checking throughout the study area. The total area (horizontal projection) of the 5,096 landslides is about 41.103 km². The LAP, which is defined as the percentage of the plane area affected by landslides, was 10.276 %, and the LND, means the number of landslides per square kilometers, was 12.74 landslides/km². Statistical analysis results show that both LAP and LND have a positive correlation with slope gradient and a negative correlation with distance from the surface rupture. However, the correlation between the occurrence of landslides with PGA, topographic position, and distance from drainages are uncertain, or has just a little positive correlation. The correlation between landslide and slope aspect also shows the effect of the directivity of the seismic wave. The Zbq formation had the most concentrated landslide activity with the LND value of 21.78 landslides/km _, ² and the ∈₁ q Gr. geological units had the highest LAP value. Furthermore, weight index (W _i) model is performed with a GIS platform to derive landslide hazard index map. The success rate of the model was 71.615 % and, thus, it was valid. In addition, comparison of five landslide controlling parameters’ influence on landslide occurrences was also carried out.     

Comparing the Hawkes and Trigger Process Models for Aftershock Sequences Following the 2005 Kashmir Earthquake

In an earlier study (Van Lieshout and Stein in Math Gesoci 44(3):309–326, 2012) we postulated the existence of two major earthquakes in the 2005 Kashmir disaster instead of a single one, based upon the pattern of aftershocks. In this study, we explore this hypothesis further by fitting several spatial point pattern models. In particular, we discuss the Hawkes and the trigger process models for earthquake aftershock sequences following the Kashmir catastrophe in 2005. The minimum contrast method is used for estimation of the parameters. The study shows that the trigger model fits better than the Hawkes model. The most likely number of main shocks is rounded to 2 generating the almost 200 aftershocks, whereas the Hawkes model would estimate a parent process of approximately 18 parents with on average about 10 descendants. We conclude that the spatial pattern of aftershocks can best be understood as a mixture of two bivariate normal distributions centered around two major shocks and estimate the parameters.     

An open dataset for landslides triggered by the 2016 Mw 7.8 Kaikōura earthquake,New Zealand

On November 14, 2016, the northeastern South Island of New Zealand was hit by the magnitude Mw 7.8 Kaikōura earthquake, which is characterized by the most complex rupturing mechanism ever recorded. The widespread landslides triggered by the earthquake make this event a great case study to revisit our current knowledge of earthquake-triggered landslides in terms of factors controlling the spatial distribution of landslides and the rapid assessment of geographic areas affected by widespread landsliding. Although the spatial and size distributions of landslides have already been investigated in the literature, a polygon-based co-seismic landslide inventory with landslide size information is still not available as of June 2021. To address this issue and leverage this large landslide event, we mapped 14,233 landslides over a total area of approximately 14,000 km². We also identified 101 landslide dams and shared them all via an open-access repository. We examined the spatial distribution of co-seismic landslides in relation to lithologic units and seismic and morphometric characteristics. We analyzed the size statistics of these landslides in a comparative manner, by using the five largest co-seismic landslide inventories ever mapped (i.e., Chi-Chi, Denali, Wenchuan, Haiti, and Gorkha). We compared our inventory with respect to these five ones to answer the question of whether the landslides triggered by the 2016 Kaikōura earthquake are less numerous and/or share size characteristics similar to those of other strong co-seismic landslide events. Our findings show that the spatial distribution of the Kaikōura landslide event is not significantly different from those belonging to other extreme landslide events, but the average landslide size generated by the Kaikōura earthquake is relatively larger compared to some other large earthquakes (i.e., Wenchuan and Gorkha).

     

Topographic and geologic controls on landslides induced by the 2015 Gorkha earthquake and its aftershocks: an example from the Trishuli Valley,central Nepal

The devastating Gorkha earthquake (M _w 7.8) on April 25, 2015 and its aftershocks triggered numerous landslides across the Lesser and Higher Himalayas of central Nepal. This study aims to characterize these landslides, based on the local topography and geology, and to develop data for landslide hazard zoning. This study focused on a mountainous catchment of the Trishuli River, where a digital elevation model was used to examine hilllslope and river profiles, aerial photos were used to identify 155 coherent landslides, and satellite images were used to map 912 earthquake-induced landslides. The topography of this area is mainly characterized by incised V-shaped inner gorges and steep (> 35°) SW-facing scarp slopes. Although most of the coherent landslides were not reactivated by the earthquakes, the Gogane landslide was affected by the earthquake and partly failed. A majority of the earthquake-induced landslides (91%) were new landslides, while the others were enlarged old landslides. The earthquake-induced landslides occurred mainly on the steep slopes of V-shaped inner gorges and scarp slopes, in gneiss and quartzite strata of the Lesser Himalayas, and they were primarily associated with fractured rock masses. This analysis provides a framework for zoning areas vulnerable to earthquake-induced landslides.     

Interpretation of landslide distribution triggered by the 2005 Northern Pakistan earthquake using SPOT 5 imagery

The 2005 northern Pakistan earthquake (magnitude 7.6) of 8 October 2005 occurred in the northwestern part of the Himalayas. We interpreted landslides triggered by the earthquake using black-and-white 2.5-m-resolution System Pour l’Observation de la Terre 5 (SPOT 5) stereo images. As a result, the counts of 2,424 landslides were identified in the study area of 55 by 51 km. About 79% or 1,925 of the landslides were small (less than 0.5 ha in area), whereas 207 of the landslides (about 9%) were large (1 ha and more in area). Judging from our field survey, most of the small landslides are shallow rock falls and slides. However, the resolution and whitish image in the photos prevented interpreting the movement type and geomorphologic features of the landslide sites in detail. It is known that this earthquake took place along preexisting active reverse faults. The landslide distribution was mapped and superimposed on the crustal deformation detected by the environmental satellite/synthetic aperture radar (SAR) data, active faults map, geological map, and shuttle radar topography mission data. The landslide distribution showed the following characteristics: (1) Most of the landslides occurred on the hanging-wall side of the Balakot–Garhi fault; (2) greater than one third of the landslides occurred within 1 km from the active fault; (3) the greatest number of landslides (1,147 counts), landslide density (3.2 counts/km²), and landslide area ratio (2.3 ha/km²) was found within Miocene sandstone and siltstone, Precambrian schist and quartzite, and Eocene and Paleocene limestone and shale, respectively; (4) there was a slight trend that large landslides occurred on vertically convex slopes rather than on concave slopes; furthermore, large landslides occurred on steeper (30° and more) slopes than on gentler slopes; (5) many large landslides occurred on slopes facing S and SW directions, which is consistent with SAR-detected horizontal dominant direction of crustal deformation on the hanging wall.     

Landslides triggered by the 2016 Mj 7.3 Kumamoto,Japan, earthquake

The aim of this study is to establish a detailed and complete inventory of the landslides triggered by the Mj 7.3 (Mw 7.0) Kumamoto, Japan, earthquake sequence of 15 April 2016 (16 April in JST). Based on high-resolution (0.5–2 m) optical satellite images, we delineated 3,467 individual landslides triggered by the earthquake, occupying an area of about 6.9 km². Then they were validated by aerial photographs with very high-resolution (better than 0.5 m) and oblique field photos. Of them, 3,460 landslides are distributed in an elliptical area about 6000 km², with a NE-SW directed 120-km-long long axis and a 60-km-long NW-SE trending short axis. Most of the landslides are shallow, disrupted falls and slides, with a few flow-type slides and rock and soil avalanches. The analysis of correlation between the landslides and several control factors shows the areas of elevation 1000–1200 m, stratum of Q₃-Hvf, seismic intensity VIII and VIII+, and peak ground acceleration (PGA) 0.4–0.6 g register the highest landslide abundance. This study also discussed the relationship between the spatial pattern of the landslides and the seismotectonic structure featured by a strike-slip fault with a normal component and the volcanism in the study area.     

The Hattian Bala rock avalanche and associated landslides triggered by the Kashmir Earthquake of 8 October 2005

The Kashmir Earthquake of the 8 October killed an estimated 87 350 people, 25 500 through co-seismic landslides. The largest landslide associated with the earthquake was the 68 × 10⁶ m³ Hattian Bala rock avalanche that destroyed a village and killed around 1000 people. The deposit blocks the valley to a depth of 130 m impounding a lake that reached the dam-crest in April 2007. An outburst flood now threatens a major settlement 3 km downstream. A series of space images reveals landslide clusters in the rock avalanche source area prior to the earthquake. The images also reveal a large slow-moving landslide with its toe in the lake, failure of this landslide may induce dam failure through overtopping and scour. Eighty five landslides in the Hattian Bala catchment predate the shaking of 8 October 2005, a further 73 are co-seismic with the main shock, and 21 postdate it in the period up to October 2006. Landslide magnitude–frequency distribution plots derived from satellite images allow an assessment of the contribution of seismically triggered events as compared to background rates of activity.     

Landslides induced by the 2017 Mw7.3 Sarpol Zahab earthquake (Iran)

Landslides are the main secondary effects of earthquakes in mountainous areas. The spatial distribution of these landslides is controlled by the local seismic ground motion and the local slope stability. While gravitational instabilities in arid and semi-arid environments are understudied, we document the landslides triggered by the Sarpol Zahab earthquake (November 12, 2017, Mw7.3, Iran/Iraq border), the largest event ever recorded in the semi-arid Zagros Mountains. An original earthquake-induced landslide inventory was derived, encompassing landslides of various sizes and velocities (from rapid disrupted rockfalls to slow-moving coherent landslides). This inventory confirms the low level of triggered landslides in semi-arid environments. It also displays clear differences in the spatial and volumetric distributions of earthquake-induced landslides, having 386 rockfalls of limited size triggered around the epicenter, and 9 giant (areas of ca. 10⁶ m²) active and ancient deep-seated landslides coseismically accelerated at locations up to 180 km from the epicenter. This unusual distant triggering is discussed and interpreted as an interaction between the earthquake source properties and the local geological conditions, emphasizing the key role of seismic ground motion variability at short spatial scales in triggering landslides. Finally, the study documents the kinematics of slow-moving ancient landslides accelerated by earthquakes, and opens up new perspectives for studying landslide triggering over short (~?1–10 years) and long-time (~?1000–10,000 years) periods.

     

Multifractality in seismic sequences of NW Himalaya

Multifractal behaviour of interevent time sequences is investigated for the earthquake events in the NW Himalaya, which is one of the most seismically active zones of India and experienced moderate to large damaging earthquakes in the past. In the present study, the multifractal detrended fluctuation analysis (MF-DFA) is used to understand the multifractal behaviour of the earthquake data. For this purpose, a complete and homogeneous earthquake catalogue of the period 1965–2013 with a magnitude of completeness M _w 4.3 is used. The analysis revealed the presence of multifractal behaviour and sharp changes near the occurrence of three earthquakes of magnitude (M _w ) greater than 6.6 including the October 2005, Muzaffarabad–Kashmir earthquake. The multifractal spectrum and related parameters are explored to understand the time dynamics and clustering of the events.

     

Landslides induced by a low magnitude seismic sequence at continental interiors: a case study of the Santa Rosa Canyon,northeastern Mexico

Landslides caused by a low magnitude earthquake swarm (2.8?≤?M?≤?3.6) in 2012 were documented at the Santa Rosa Canyon in northeastern Mexico. Disrupted landslides from falls and slides, in both rocks and soils, were identified based on fieldwork and high-resolution satellite imagery along stream banks from natural cliffs and along the road cut in the epicentral area. Most of the landslides occurred on slopes greater than 40°, where geological features played a key role in triggering slope instabilities. The maximum distance limit for disrupted slides from the epicentral area was 7 km. The area affected by landslides during the early stage of the seismic sequence (July through August 2012) was 90 km². Landslide identification was limited in some areas by the resolution of the satellite imagery and dense cloud coverage. Both the epicentral distance and the area affected by landslides are above the global bounds reported in literature. The final landslide inventory is the first documented case of earthquake-induced landslides in northeastern Mexico.     

Seismic landslides hazard zoning based on the modified Newmark model: a case study from the Lushan earthquake,China

Seismically triggered landslides can cause great damage to the road construction in mountainous areas. The permanent displacement analysis based on Newmark sliding-block model can evaluate risk of these landslides from the perspective of deformation damage and overall failure probability of slopes. However, the sliding-block model does not consider the attenuation effect of the shear strength on the sliding surface during earthquake, causing the calculated value of Jibson method to be less than the actual value. Therefore, the Newmark model was modified by adding attenuation coefficients to the effective internal friction angle and the effective cohesion of geologic units. The landslide areal density was proposed for hazard zoning with the Wenchuan earthquake data. The results showed that the predicted values agreed well with the real distribution of the landslides triggered by the Lushan earthquake. The proposed hazard zoning method in this paper can predict the severity of seismic landslides in consideration of the environmental changes in mountainous regions after the earthquake and provide support for the site selection in highly seismic areas.

     

Integrating intensity–duration-based rainfall threshold and antecedent rainfall-based probability estimate towards generating early warning for rainfall-induced landslides in parts of the Garhwal Himalaya,India

In order to generate early warning for landslides, it is necessary to address the spatial and temporal aspects of slope failure. The present study deals with the temporal dimension of slope failures taking into account the most widespread and frequent triggering factor, i.e. rainfall, along the National Highway-58 from Rishikesh to Mana in the Garhwal Himalaya, India. Using the post-processed three-hourly rainfall intensity and duration values from the Tropical Rainfall Measuring Mission-based Multi-satellite Precipitation Analysis and the time-tagged landslide records along this route, an intensity–duration (I–D)-based threshold has been derived as I?=?58.7D ^?1.12 for the rainfall-triggered landslides. The validation of the I–D threshold has shown 81.6 % accuracy for landslides which occurred in 2005 and 2006. From this result, it can be inferred that landslides in the study area can be initiated by continuous rainfall of over 12 h with about 4-mm/h intensity. Using the mean annual precipitation, a normalized intensity–duration relation of NI?=?0.0612D ^?1.17 has also been derived. In order to account for the influence of the antecedent rainfall in slope failure initiation, the daily, 3-day cumulative, and 15- and 30-day antecedent rainfall values associated with landslides had been subjected to binary logistic regression using landslide as the dichotomous dependent variable. The logistic regression retained the daily, 3-day cumulative and 30-day antecedent rainfall values as significant predictors influencing slope failure. This model has been validated through receiver operating characteristic curve analysis using a set of samples which had not been used in the model building; an accuracy of 95.1 % has been obtained. Cross-validation of I–D-based thresholding and antecedent rainfall-based probability estimation with slope failure initiation shows 81.9 % conformity between the two in correctly predicting slope stability. Using the I–D-based threshold and the antecedent rainfall-based regression model, early warning can be generated for moderate to high landslide-susceptible areas (which can be delineated using spatial integration of preconditioning factors). Temporal predictions where both the methods converge indicate higher chances of slope failures for areas predisposed to instability due to unfavourable geo-environmental and topographic parameters and qualify for enhanced slope failure warning. This method can be verified for further rainfall seasons and can also be refined progressively with finer resolutions (spatial and temporal) of rainfall intensity and multiple rain gauge stations covering a larger spatial extent.     

Seismic hazard and probability assessment of Kashmir valley,northwest Himalaya,India

Seismic hazard analysis of the northwest Himalayan belt was carried out by using extreme value theory (EVT). The rate of seismicity (a value) and recurrence intervals with the given earthquake magnitude (b value) was calculated from the observed data using Gutenberg–Richter Law. The statistical evaluation of 12,125 events from 1902 to 2017 shows the increasing trend in their inter-arrival times. The frequency–magnitude relation exhibits a linear downslope trend with negative slope of 0.8277 and positive intercept of 4.6977. The empirical results showed that the annual risk probability of high magnitude earthquake M?≥?7.7 in 50 years is 88% with recurrence period of 47 years, probability of M?≤?7.5 in 50 years is 97% with recurrence period of 27 years, and probability of M?≤?6.5 in 50 years is 100% with recurrence period of 4 years. Kashmir valley, located in the NW Himalaya, encompasses a peculiar tectonic and structural setup. The patterns of the present and historical seismicity records of the valley suggest a long-term strain accumulation along NNW and SSE extensions with the decline in the seismic gap, posing a potential threat of earthquakes in the future. The Kashmir valley is characterized by the typical lithological, tectono-geomorphic, geotechnical, hydrogeological and socioeconomic settings that augment the earthquake vulnerability associated with the seismicity of the region. The cumulative impact of the various influencing parameters therefore exacerbates the seismic hazard risk of the valley to future earthquake events.