Correction to: Stress shadow,stress triggering,and recent earthquake activity in the Kashmir Himalaya,India

Journal of Seismology -     

Landslide hazard zonation using analytical hierarchy process along National Highway-3 in mid Himalayas of Himachal Pradesh,India

Landslide zonation studies emphasize on preparation of landslide hazard zonation maps considering major instability factors contributing to occurrence of landslides. This paper deals with geographic information system-based landslide hazard zonation in mid Himalayas of Himachal Pradesh from Mandi to Kullu by considering nine relevant instability factors to develop the hazard zonation map. Analytical hierarchy process was applied to assign relative weightages over all ranges of instability factors of the slopes in study area. To generate landslide hazard zonation map, layers in geographic information system were created corresponding to each instability factor. An inventory of existing major landslides in the study area was prepared and combined with the landslide hazard zonation map for validation purpose. The validation of the model was made using area under curve technique and reveals good agreement between the produced hazard map and previous landslide inventory with prediction accuracy of 79.08%. The landslide hazard zonation map was classified by natural break classifier into very low hazard, low hazard, moderate hazard, high hazard and very high landslide hazard classes in geographic information system depending upon the frequency of occurrence of landslides in each class. The resultant hazard zonation map shows that 14.30% of the area lies in very high hazard zone followed by 15.97% in high hazard zone. The proposed model provides the best-fit classification using hierarchical approach for the causative factors of landslides having complex structure. The developed hazard zonation map is useful for landslide preparedness, land-use planning, and social-economic and sustainable development of the region.     

Epidemic type aftershock sequence (ETAS) modeling of northeastern Himalayan seismicity

We analyzed the seismicity of northeastern Himalayan region of latitude (25 to 32°?N) and longitude (86–97°?E). The US Geological Survey catalogue is used in this study for a period from 1973 to June 2011. The seismicity of the region is modeled using epidemic type aftershock sequence (ETAS) model. The region is divided in three parts: (1) whole region, (2) subregion I, and (3) subregion II. The magnitude of completeness is found to be 4.6 for all the three regions. The ETAS parameters for all the regions are found same within the standard errors. There is no significant change observed in the seismicity since 1973 based on the ETAS modeling.     

Determination of Q β(f) in Different Parts of Kumaon Himalaya from the Inversion of Spectral Acceleration Data

This paper presents the results of a modified two-step inversion algorithm approach to find S wave quality factor Q _β(f) given by Joshi (Bull Seis Soc Am 96:2165–2180, 2006). Seismic moment is calculated from the source displacement spectra of the S wave using both horizontal components. Average value of seismic moment computed from two horizontal components recorded at several stations is used as an input to the first part of inversion together with the spectra of S phase in the acceleration record. Several values of the corner frequency have been selected iteratively and are used as inputs to the inversion algorithm. Solution corresponding to minimum root mean square error (RMSE) is used for obtaining the final estimate of Q _β(f) relation. The estimates of seismic moment, corner frequency and Q _β(f) from the first part of inversion are further used for obtaining the residual of theoretical and observed source spectra which are treated as site amplification terms. The acceleration record corrected for the site amplification term is used for determination of seismic moment from source spectra by using Q _β(f) obtained from first part of inversion. Corrected acceleration record and new estimate of seismic moment are used as inputs to the second part of the inversion scheme which is similar to the first part except for use of input data. The final outcome from this part of inversion is a new Q _β(f) relation together with known values of seismic moment and corner frequency of each input. The process of two-step inversion is repeated for this new estimate of seismic moment and goes on until minimum RMSE is obtained which gives final estimate of Q _β(f) at each station and corner frequency of input events. The Pithoragarh district in the state of Uttarakhand in India lies in the border region of India and Nepal and is part of the seismically active Kumaon Himalaya zone. A network of eight strong motion recorders has been installed in this region since March, 2006. In this study we have analyzed data from 18 local events recorded between March, 2006 and October, 2010 at various stations. These events have been located using HYPO71 and data has been used to obtain frequency-dependent shear-wave attenuation. The Q _β(f) at each station is calculated by using both the north-south (NS) and east-west (EW) components of acceleration records as inputs to the developed inversion algorithm. The average Q _β(f) values obtained from Q _β(f) values at different stations from both NS and EW components have been used to compute a regional average relationship for the Pithoragarh region of Kumaon Himalaya of form Q _β(f)?=?(29?±?1.2)f ^{(1.1 ± 0.06)}.     

Estimation of Source Parameters of M w 6.9 Sikkim Earthquake and Modeling of Ground Motions to Determine Causative Fault

In this study, source parameters of the September 18, 2011 M _w 6.9, Sikkim earthquake were determined using acceleration records. These parameters were then used to generate strong motion at a number of sites using the stochastic finite fault modeling technique to constrain the causative fault plane for this earthquake. The average values of corner frequency, seismic moment, stress drop and source radius were 0.12 Hz, 3.07 × 10²⁶ dyne-cm, 115 bars and 9.68 km, respectively. The fault plane solution showed strike-slip movement with two nodal planes oriented along two prominent lineaments in the region, the NE-oriented Kanchendzonga and NW-oriented Tista lineaments. The ground motions were estimated considering both the nodal planes as causative faults and the results in terms of the peak ground accelerations (PGA) and Fourier spectra were then compared with the actual recordings. We found that the NW–SE striking nodal plane along the Tista lineament may have been the causative fault for the Sikkim earthquake, as PGA estimates are comparable with the observed recordings. We also observed that the Fourier spectrum is not a good parameter in deciding the causative fault plane.     

Drainage Performance of Different Sizes Tire Chips used Alone and Mixed with Natural Aggregates as Leachate Drainage Layer Material

Results of compressibility and hydraulic conductivity testing performed under different levels of average vertical stress were used to determine the appropriate tire chips size (in particular aspect ratio) and tire chip—gravel random state mix combination. Proper sidewall treatment was provided to reduce side wall friction up to minimum level and constant head conditions were maintained with appropriate design flow rate. Results showed that tire chips (length 25 mm, aspect ratio 1) mixed with 40 mm gravel in the random state in mixing ratio of 3:1 (ds:dg) was observed to undergo minimum strain, maximum hydraulic conductivity and hence free drainage even under high vertical stress levels. By selecting tire chips of appropriate aspect ratio and optimum mix combination for the construction of drainage layer of MSW landfill leachate collection system lower strain response, higher hydraulic conductivity values and relatively free drainage even under anticipated higher vertical stress levels may be achieved.     

Earthquake hazard zonation of Sikkim Himalaya using a GIS platform

An earthquake hazard zonation map of Sikkim Himalaya is prepared using eight thematic layers namely Geology (GE), Soil Site Class (SO), Slope (SL), Landslide (LS), Rock Outcrop (RO), Frequency–Wavenumber (F–K) simulated Peak Ground Acceleration (PGA), Predominant Frequency (PF), and Site Response (SR) at predominant frequencies using Geographic Information System (GIS). This necessitates a large scale seismicity analysis for seismic source zone classification and estimation of maximum earthquake magnitude or maximum credible earthquake to be used as a scenario earthquake for a deterministic or quasi-probabilistic seismic scenario generation. The International Seismological Center (ISC) and Global Centroid Moment Tensor (GCMT) catalogues have been used in the present analysis. Combining b-value, fractal correlation dimension (Dc) of the epicenters and the underlying tectonic framework, four seismic source zones are classified in the northeast Indian region. Maximum Earthquake of M _W 8.3 is estimated for the Eastern Himalayan Zone (EHZ) and is used to generate the seismic scenario of the region. The Geohazard map is obtained through the integration of the geological and geomorphological themes namely GE, SO, SL, LS, and RO following a pair-wise comparison in an Analytical Hierarchy Process (AHP). Detail analysis of SR at all the recording stations by receiver function technique is performed using 80 significant events recorded by the Sikkim Strong Motion Array (SSMA). The ground motion synthesis is performed using F–K integration and the corresponding PGA has been estimated using random vibration theory (RVT). Testing for earthquakes of magnitude greater than M _W 5, a few cases presented here, establishes the efficacy and robustness of the F–K simulation algorithm. The geohazard coverage is overlaid and sequentially integrated with PGA, PF, and SR vector layers, in order to evolve the ultimate earthquake hazard microzonation coverage of the territory. Earthquake Hazard Index (EHI) quantitatively classifies the terrain into six hazard levels, while five classes could be identified following the Bureau of Indian Standards (BIS) PGA nomenclature for the seismic zonation of India. EHI is found to vary between 0.15 to 0.83 quantitatively classifying the terrain into six hazard levels as “Low” corresponding to BIS Zone II, “Moderate” corresponding to BIS Zone III, “Moderately High” belonging to BIS Zone IV, “High” corresponding to BIS Zone V(A), “Very High” and “Severe” with new BIS zones to Zone V(B) and V(C) respectively.     

Characteristics of seismic wave attenuation in the Kishtwar and its adjoining region of NW Himalaya

Journal of Seismology - We investigate the attenuation characteristics of high frequency seismic waves in the Kishtwar and its adjoining region of NW Himalaya using 161 local earthquakes (M...     

Largest earthquake in Himalaya: An appraisal

The largest earthquake (Mw 8.4 to 8.6) in Himalaya reported so far occurred in Assam syntaxial bend in 1950. However, some recent studies have suggested for earthquake of magnitude M_w 9 or more in the Himalayan region. In this paper, we present a detailed analysis of seismological data extending back to 1200 AD, and show that earthquake in Himalayan region may not be expected to be as large as those of subduction zones. Also, there appears to be a lateral variation in the earthquake magnitude, being lesser in the western syntaxial bend when compared close to the eastern syntaxial bend. This is attributed to the difference in the plate boundary scenario; dominance of strike-slip and thrusting along the western syntaxis as against thrusting and remnant subduction along the eastern syntaxis.