Earth has a complex structure and due to innumerable phenomena occurring on its surface and interior, disastrous events are likely to happen without any warning. Such adverse geological conditions, which are capable of causing damage or loss of property and life, are termed as geological hazard   . According to Coppola , hazard is the origin of the potential damage to community in which population, infrastructure, inhabitant’s property, business and environment are included, for example, it affects human life in the form of injuries, deaths; property damage and social-economic disturbances  . Whether it be an earthquake, volcanic activity, tsunami, landslide etc., it is almost difficult to stop such natural events due to their unpredictability . Due to advancement in the field of geology, we are able to understand such events, eventually enabling ourselves to predict them so that the damage can be reduced  . The phenomenon of movement of earth material from insignificant to catastrophic level either due to natural or artificial causes is termed as landslide   . Lundgren  described landslides as the movement of mass of earth, which initiated under the action of the gravitational forces. Singh and Goel  termed the advancement of mass of rock or soil in downward and outward directions i.e. the displacement of the mass on a down the hillside or on a sloping ground, as landslide. The nature and intensity of the landslides make them a threat to settlements, infrastructure and communication networks thereby increasing the concern of study and research.
Slope failure is the result of shock e.g. increase in destabilizing force or external load, earthquakes, undermining, cracks in rock with various types such as hair-line cracks, frost wedging, or by increased pressure of water in cracks, clay swelling, excavation and loss of the capillary tension can also become a cause of decrease in shear strength resulting in slope failure . Slope stability analysis can be performed by kinematic analysis, limit equilibrium and numerical modelling techniques. Kinematic analysis is a very useful method to understand types of failure by using discontinuities and joints orientations   . Problems related to slope stability can be reduced with refined assessment and detailed examination . The Northern Pakistan has abundance of various types of landslides, which usually occur along major thrusts and road cuts  . Many researchers studied landslides in detail to formulate the techniques to rectify and conciliate the issue of landslides       . The earthquake, which occurred on October 8th, 2005, caused numerous landslides in the Balakot area and this earthquake also activated various landslides in the nearby area   - .
2. Study Area
The area under study comprises of mountainous outcrops covering the city of Balakot, district Mansehra in the province of Khyber Pakhtunkhwa, Pakistan. Balakot city is situated on the banks of River Kunhar before it enters Kashmir Valley. Location map of the study area is given in Figure 1. The city is known to be a tropical area with a good vegetation and forestation cover along steep mountains. This study incorporates the area of Balakot city emphasizing upon
Figure 1. Location map of selected sites .
Based on visual observation and slope characteristics, two sites were selected for detailed studies and sites were marked as LS-1 and LS-2. The outcrops of site LS-1 lie between longitude and latitude 34˚31'39.5"N and 73˚20'49.1"E to 34˚31'30.3"N and 73˚20'43.6"E, whereas LS-2 lies between 34˚31'26.4"N and 73˚20'42.8"E to 34˚31'17.1"N and 73˚20'43"E.
Rock slope failures in this area mainly triggered after earthquake in 2005. Balakot city was destroyed and rock units seems unstable which caused causalities and damages to infrastructure   . The lithological units along the road cut composed of soft rocks as well as high precipitation in the area induce sliding of rock mass along the road.
The North-western Himalayas of Pakistan is formed by the collision of plate boundaries   . From north to south it is divided as Asian plate, shyok suture zone, Kohistan island arc and Indian plate     . Indian plate is divided into sub-Himalayas, lesser Himalayas and higher Himalayas bounded by major thrusts. The study area lies close to main boundary thrust (MBT) and Hazara Kashmir Syntaxis, which is a major syntaxial bend of Himalayas. Due to two active regional faults traversing right through the Balakot town, almost entire Balakot city falls in potential hazardous zone. The MBT is passing through North of the Balakot city as Balakot-Bagh thrust (Figure 2) .
Figure 2. Geological map of study area (reproduced after USGS  ).
The Balakot-Bagh thrust is the main identified feature in this zone . Formation observed in the emphasized area is Salkhala Formation. The formation is mainly termed as metamorphosed zone and comprised of Precambrian to Cambrian sequence of quartzite, slates, phyllite, schists and marble .
Salkhala Formation of Precambrian age consists of slates with graphitic beds interbedded with marble . It also contains quartzite, quartz schist, garnet-biotite and talc schist at different places. Near Balakot city, this formation thins out and remains to only few hundred meters  . It is further divided into two rock units; upper is called as Kundal unit and basal as Gandaf unit .
A systematic work methodology was developed in order to study the landslides in Balakot area by employing literature review of previous research especially after earthquake. An overview of the research methodology is giving in Figure 3 and it is explained briefly below. A reconnaissance study carried out to assess the nature of slopes and geology of study area. Geological map demarcating regional faults, rock and soil units was prepared (Figure 2).
Discontinuity Survey were carried out as per standard   to assess nature of jointed rock. The parameters such as orientation, spacing, persistence, groundwater conditions and infilling material were used to assess rock quality.
After gathering required data from the sites which seemed to be prone to sliding, the stability analyses of those slopes were conducted by using the empirical classification system, limit equilibrium analysis and kinematic analysis. The parameters which were used in this study are explained in few lines.
For in-situ evaluation of strength of rocks, an indirect tool, the Schmidt
Figure 3. Work methodology flow chart.
Rebound Hammer was employed in the field. The value of rebound number (Rn) was calculated and later it was correlated with the unit weight using correlation chart   . Geological strength index (GSI) was calculated by using the correlation given below:
The GSI value obtained from Equation (1) is one of the input parameters of the analysis using limit equilibrium method. Equation (2) was proposed by Palmström  that a correlation can be made for the indirect estimation of Rock Quality Designation (RQD) by using the joint volumetric count (Jv). The value of Jv represent the sum of reciprocal of all joints which are defined as follows:
Here, S1, S2 and represent the spacing of the joints depending on the spacing of joints in rock mass. NR represents the number of the random joints present in the unit and according to Palmström    spacing of 5 m should be used for the random joints.
Procedure of evaluating and predicting the possible direction of moment of rock is called the kinematic analysis . Kinematic analysis is a geometric method, which uses the angular interactions among the discontinuity planes in order to predict the possibility and failure of jointed rock mass unit . Kinematic analysis performed with DIPS v. 7.0 , a program of RocScience suite. The limit equilibrium analyses were performed to estimate global failure analysis, factor of safety and surface of failure. By employing the Slide 7.0 (RocScience), limit equilibrium analyses were performed in order to measure the factor of safety (FOS). As discuss earlier, the study area lies in seismically active zone, so to generate a realistic scenario as per field, twenty various analyses were performed in Slide. Analyses were performed in dry to wet condition, whereas the seismic loading of various peak ground acceleration (PGA) was also employed in the analyses. Four slices methods, Bishop , Janbu , Corps of Engineers , and Spencer  , were used in these analyses.
Geological map of Balakot city, and adjoining areas having scale 1:24,000 was prepared using satellite images from Google Earth and it was reproduced using ArcGIS 10.3. The major rock units exposed at selected slopes slates and quartzite of Salkhala Formation. The digitized map of Balakot city is given in Figure 2.
Along the road, slided mass was present which caused damaged to retaining structures. Sites that seem unstable and had a chance of landslide were marked and selected for detail studies (Figure 1 and Figure 2). Discontinuity data were conducted on selected sites. Samples were subjected to calculate density with average values ranging from 2.7 - 3.3 g/cm3. Schmidt Rebound (Rn) values ranges from 15 to 60 so the average strength of rock mass was calculated to be within the range of 84 - 134 MPa by using the correlation chart    . Three to four joint sets were observed with closed to moderately space with very low to low persistence and apertures were observed to be as closed to moderately wide.
Empirical classification of rock mass carried out by Rock Mass Rating (RMR) classification system    . Table 1 shows that rock mass rating is 62 and 63 for LS1 and LS-2 respectively that falls in good quality.
Kinematic analysis carried out by Dips 7.0 using the orientation data, the LS-1 is likely to have Planar sliding and direct toppling along joint J1 (Figure 4(a)). Flexural toppling is likely to happen along Random joints whereas the wedge sliding is not likely along this slope. On the other hand, LS-2 is likely to have planar sliding and direct toppling along joint J1 and there is probability of wedge failure along J1 and J2 (Figure 4(b)). Flexure toppling is likely to occur along the foliation at LS-2. The results of the kinematic analyses are summarized in Table 2. The results of kinematic analyses confirm the field observations as well.
The limit equilibrium analyses were performed for global failure analysis, factor of safety (FOS) and surface of failure. By employing the Slide 7.0 (RocScience) software, analyses were performed in dry to wet condition, whereas the seismic loading of 0.32 peak ground acceleration (PGA) and under extreme conditions using both wet and PGA.
For the LS-1, in existing conditions, Bishop, Janbu, Spencer and Corps of Engineer method calculated FOS is 5.932, 6.049, 6.125 and 6.098 respectively. Whereas in wet conditions with 10 m ground water level and 70% saturation,
Table 1. Rock mass rating of slopes in study area.
Figure 4. Kinematic analyses of slopes (a) Landslide-1; (b) Landslide-2.
FOS is 6.919, 6.949, 7.044 and 7.238. With maximum seismic loading where the PGA was set as 0.32, the results for aforementioned methods were, 4.205, 3.854, 4.569 and 3.691. In extreme conditions, where the PGA is 0.32 along with presence of water, Bishop, Janbu, Spencer and Corps of Engineer #1 methods calculated FOS as 4.69, 4.278, 5.079 and 4.458 respectively (Figure 5).
Table 2. Summary of results based on kinematic analyses.
Figure 5. Limit equilibrium analyses of LS-1 using Spencer slice method (a) under normal existing condition of slope (b) water table at 10 m with 70% saturation (c) seismic loading at 0.32 peak ground acceleration (PGA) (d) under extreme condition with 10 m ground water and 0.32 PGA.
At LS-2, in dry and existing conditions, FOS for Bishop, Janbu, Spencer and Corps of Engineer methods turned out to be 5.932, 6.049, 6.125 and 6.098 respectively. In wet conditions, FOS values were 5.932, 6.049, 6.125 and 6.098 for aforementioned methods. With the maximum seismic loading of 0.32 PGA, FOS values were 4.205, 3.854, 4.569 and 3.961 for Bishop, Janbu, Spencer and Corps of Engineer methods respectively. In the extreme conditions, the values of FOS for Bishop, Janbu, Spencer and Corps of Engineer methods were 4.205, 3.854, 4.569 and 3.961 respectively (Figure 6). FOS values obtained through all the analyses, despite extreme conditions, came out to be greater than 1 (Table 3) which means that the slopes are safe from failure.
The current study focuses on evaluation of slopes in static and pseudo-static condition, as Balakot lies in a tectonically active regime. During the field excursion,
Figure 6. Limit equilibrium analyses of LS-2 using Spencer slice method (a) under normal existing condition of slope (b) water table at 10 m with 70% saturation (c) seismic loading at 0.32-peak ground acceleration (PGA) (d) under extreme condition with 10.
Table 3. Results of factor of safety (FOS) using limit equilibrium slices method.
two sites were selected, and geological map was prepared. The rocks exposed on the sites were mainly Slates and phyllites of Salkhala Formation. Discontinuity surveys were conducted on the selected sites. The discontinuity data was used to evaluate the rock quality using Rock Mass Rating (RMR). RMR depicts that rock lies in good category.
Kinematic analysis was performed with computer code DIPS v. 7.0 to assess stability of slopes. At LS-1, plane failure and direct toppling are likely along the joints while there is no possibility of wedge failure. Flexural topping can also be expected at this site. While at LS-2, plane failure is likely along foliation planes and wedge failure can be expected at the intersection of joints of J1 and J2.
Limit equilibrium analysis was performed with computer code Slide v. 7.0 to assess FOS of slopes at static and pseudo-static conditions. At four various conditions such as dry, wet, seismic loading and extreme, it was revealed that FOS is >1 which concluded that there is no possibility of global failure along these slopes.
This study concludes that using empirical classification system, kinematic and limit equilibrium analyses simultaneously, stable slopes can be designed in seismically active regime, however, verification of the results is always emphasised. The results of study verify the study of Palmström and Stille , that is the designing of engineering projects should always be related to ongoing stage of the project with ongoing decisions. In the present study the decisions were made by the forward analyses using three different methods which were empirical classification schemes for the rock quality, kinematic analysis for the rock failures and limit equilibrium analyses for the stability of slopes in static and pseudo static conditions as the present study lies in active tectonic regime.
 Tilloy, A., Malamud, B.D., Winter, H. and Joly-Laugel, A. (2019) A Review of Quantification Methodologies for Multi-Hazard Interrelationships. Earth-Science Reviews, 196, Article ID: 102881.
 Van Westen, C.J., Castellanos, E. and Kuriakose, S.L. (2008) Spatial Data for Landslide Susceptibility, Hazard, and Vulnerability Assessment: An Overview. Engineering Geology, 102, 112-131.
 Akram, M.S., Zeeshan, M., Mirza, K., Ahmed, L., Noor, A. and Ali, M. (2019) Slope Stability Analyses Using Classification Systems and Numerical Methods: Case Study from Lower Dir, Khyber Pakhtunkhwa, Pakistan. Himalayan Geology, 40, 67-77.
 Lundgren, L. (1978) Studies of Soil and Vegetation Development on Fresh Landslide Scars in the Mgeta Valley, Western Uluguru Mountains, Tanzania. Geografiska Annaler: Series A, Physical Geography, 60, 91-127.
 Singh, B. and Goel, R.K. (2011) Engineering Rock Mass Classification: Tunneling, Foundations, and Landslides. Butterworth-Heinemann, Oxford.
 Gurocak, Z., Alemdag, S. and Zaman, M.M. (2008) Rock Slope Stability and Excavatability Assessment of Rocks at the Kapikaya Dam Site, Turkey. Engineering Geology, 96, 17-27.
 Qi, C., Wu, J., Liu, J. and Kanungo, D.P. (2016) Assessment of Complex Rock Slope Stability at Xiari, Southwestern China. Bulletin of Engineering Geology and the Environment, 75, 537-550.
 Ikenna Nkpadobi, J., Kuna Raj, J. and Fatt Ng, T. (2016) Classification of Cut Slopes in Weathered Meta-Sedimentary Bedrocks. Earth Sciences Research Journal, 20, 1-9.
 Aydan, O. (2006) Geological and Seismological Aspects of Kashmir Earthquake of October 8, 2005 and a Geotechnical Evaluation of Induced Failures of Natural and Cut Slopes. Journal of The School of Marine Science and Technology, 4, 25-44.
 Akram, M.S., Ahmed, L., Farooq, S., Ahad, M.A., Zaidi, S.M.H., Khan, M., et al. (2018) Geotechnical Evaluation of Rock Cut Slopes Using Basic Rock Mass Rating (RMR Basic), Slope Mass Rating (SMR) and Kinematic Analysis along Islamabad Muzaffarabad Dual Carriageway (IMDC), Pakistan. Journal of Biodiversity and Environmental Sciences, 13, 297-306.
 Owen, L.A., Kamp, U., Khattak, G.A., Harp, E.L., Keefer, D.K. and Bauer, M.A. (2008) Landslides Triggered by the 8 October 2005 Kashmir Earthquake. Geomorphology, 94, 1-9.
 Hussain, A. (2005) Geology and Tectonics of Northern Pakistan with Respect to October 8, 2005, Earthquake. In: Earthquake Rehabilitation Conference, Seismology, Structures and Codes, Islamabad, 18-19.
 Zaré, M. and Karimi-Paridari, S. (2009) An Investigation on Balakot, Muzaffarabad (Pakistan) Earthquake, 8 Oct. 2005, Mw 7.6; Geological Aspects and Intensity Distribution. Journal of Seismology, 13, 327-337.
 Aydan, O., Ohta, Y. and Hamada, M. (2009) Geotechnical Evaluation of Slope and Ground Failures during the 8 October 2005 Muzaffarabad Earthquake, Pakistan. Journal of Seismology, 13, 399-413.
 Treloar, P.J., Coward, M.P., Williams, M.P. and Khan, M.A. (1989) Basement-Cover Imbrication South of the Main Mantle Thrust, North Pakistan. Geological Society of America Special Papers, 232, 137-152.
 Anczkiewicz, R., Oberli, F., Burg, J.P., Meier, M., Dawood, H. and Hussain, S.S. (1998) Magmatism South of the Indus Suture, Lower Swat, Pakistan. 13th Himalaya-Karakoram-Tibet International Workshop, Peshawar, Pakistan, 20-22 April 1998, 7-9.
 Garrido, C.J., Bodinier, J.-L., Burg, J.-P., Zeilinger, G., Hussain, S.S., Dawood, H., et al. (2006) Petrogenesis of Maficgarnet Granulite in the Lower Crust of the Kohistan Paleo-Arc Complex (Northern Pakistan): Implications for Intra-Crustal Differentiation of Island Arcs and Generation of Continental Crust. Journal of Petrology, 47, 1873-914.
 Heuberger, S., Schaltegger, U., Burg, J.-P., Villa, I.M., Frank, M., Dawood, H., et al. (2007) Age and Isotopic Constraints on Magmatism along the Karakoram-Kohistan Suture Zone, NW Pakistan: Evidence for Subduction and Continued Convergence after India-Asia Collision. Swiss Journal of Geosciences, 100, 85-107.
 Calkins, J.A., Offield, T.W., Abdullah, S.K. and Ali, S.T. (1975) Geology of the Southern Himalaya in Hazara, Pakistan, and Adjacent Areas. U.S. Government Publishing Office, Washington DC.
 Calkins, J.A., Offield, T.W. and Ali, S.T. (1969) Geology and Mineral Resources of Southern Hazara District, West Pakistan and Part of Western Kashmir. United States Department of State and Government of Pakistan Project Report, 43, 92.
 ISRM (International Society for Rock Mechanics) (1981) Suggested Methods for the Quantitative Description of Discontinuities in Rock Masses. Rock Characterisation, Testing and Monitoring-ISRM Suggested Methods. Pergamon Press, Oxford, 319-368.
 Aydin, A. (2008) ISRM Suggested Method for Determination of the Schmidt Hammer Rebound Hardness: Revised Version. In: The ISRM Suggested Methods for Rock Characterization, Testing and Monitoring: 2007-2014, Springer, Amsterdam, 25-33.
 Palmstrom, A. (2005) Measurements of and Correlations between Block Size and Rock Quality Designation (RQD). Tunnelling and Underground Space Technology, 20, 362-377.
 Palmstrom, A. and Broch, E. (2006) Use and Misuse of Rock Mass Classification Systems with Particular Reference to the Q-System. Tunnelling and Underground Space Technology, 21, 575-593.
 Palmstrom, A. and Stille, H. (2007) Ground Behaviour and Rock Engineering Tools for Underground Excavations. Tunnelling and Underground Space Technology, 22, 363-376.