Seismic hazards assessment of Kumaun Himalaya and adjacent region

The ongoing continent?Ccontinent collision between Indian and Eurasian plates houses a seismic gap in the geologically complex and tectonically active central Himalaya. The seismic gap is characterized by unevenly distributed seismicity. The highly complex geology with equally intricate structural elements of Himalaya offers an almost insurmountable challenge to estimating seismogenic hazard using conventional methods of Physics. Here, we apply integrated unconventional hazard mapping approach of the fractal analysis for the past earthquakes and the box counting fractal dimension of structural elements in order to understand the seismogenesis of the region properly. The study area extends from latitude 28°N?C33°N and longitude 76°E?C81°E has been divided into twenty-five blocks, and the capacity fractal dimension (D ₀) of each block has been calculated using the fractal box counting technique. The study of entire blocks reveal that four blocks are having very low value of D ₀ (0.536, 0.550, 0.619 and 0.678). Among these four blocks two are characterized by intense clustering of earthquakes indicated by low value of correlation fractal dimension (D _c ) (0.245, 0.836 and 0.946). Further, these two blocks are categorized as highly stressed zones and the remaining two are characterized by intense clustering of structural elements in the study area. Based on the above observations, integrated analysis of the D _c of earthquakes and D ₀ of structural elements has led to the identification of diagnostic seismic hazard pattern for the four blocks.     

Estimation of strain distribution using GPS measurements in the Kumaun region of Lesser Himalaya

The continuous process of continent–continent collision between the Indian and the Eurasian plates has led to the formation of the Himalayan range and continuously caused earthquakes in the region. Large earthquakes with magnitudes of 8 and above occur in this region infrequently, releasing the elastic strain accumulated over years around the plate boundary. Geodetic measurements can help estimate the strain distribution along the fault system. These measurements provide information on active deformations and associated potential seismic hazards along the Himalayan arc. In order to understand the present deformation around the plate boundary, we collected GPS data during three campaigns in the years of 2005–2007 at 16 sites in the Kumaun region of the Lesser Himalaya. Horizontal velocity vectors estimated in ITRF2000 are found to be in the range of 41–50 mm/yr with an uncertainty level of the order of 1 mm/yr. The velocity field indicates that the present convergence of around 15 mm/yr takes place in the Kumaun Himalaya. Further, we estimate the strain components in the study area for understanding the currently active tectonic process in the region. The estimated dilatational strain indicates that the northern part near the Main Central Thrust (MCT) is more compressional than the southern part. Maximum shear strain is mostly accommodated in the northern part too. The maximum shear and dilatational strain rates are about 1.0 and 0.5 μstrain/yr. It is seen that the distribution of high shear strain spatially correlates with seismicity. The maximum of extensional and compressional strains due to the force acting along the Main Central Thrust (MCT) in the NW–SE direction are found to be 0.4 and 0.1 μstrain/yr, respectively. The maximum shear strain in the northern part of the Himalaya appears to be associated with the convergence of the region proposed by other geophysical studies.     

Structural set up in northern part of Kumaun Himalaya

The Main Central Thrust (MCT) and the Main Boundary Thrust (MBT) are the two major thrusts in Kumaun, the MCT forming the boundary between highly sheared, deformed and mylonitized rocks of the Great Himalayan Central Crystallines and the Lesser Himalayan metasedimentaries. While in the Central Crystallines four-folding episodes are observed of which two are of the Precambrian age, the Lesser Himalayan rocks show only two phases of folding. MCT has its own distinctive structural history and the crystalline mass comprises an integral part of peninsular India.     

Seismicity and source parameters of local earthquakes in Bilaspur region of Himachal Lesser Himalaya

The pattern of local seismicity (110 events) and the source parameters of 26 local events (1.0?≤?Mw?≤?2.5) that occurred during May 2008 to April 2009 in Bilaspur region of Himachal Lesser Himalaya were determined. The digital records available from one station have been used to compute the source parameters and f _max based on the Brune source model (1970) and a high-frequency diminution factor (Boore 1983) above f _max. The epicentral distribution of events within 30 km of local network is broadly divided into three clusters of seismic activity: (1) a cluster located to the south of the Jamthal (JAMT) station and falls to the north of the Main Boundary Thrust (MBT) which seems to reflect the contemporary local seismicity of the segment of the MBT, (2) an elongated zone of local seismicity NE–SW trending, delineated NE of JAMT station that falls in the Lesser Himalaya between the MBT and the Main Central Thrust, and (3) NE–SW trending zone of local seismic activity located at about 10 km east of NHRI station and about 15 km northeast of NERI station and extending over a distance of about 20 km. Majority of events occur at shallow depths up to 20 km, and the maximum number of events occurs in the focal depth range between 10 and 15 km. The entire seismic activity is confined to the crust between 5 and 45 km. The average values of these source parameters range from 3.29?×?10¹⁷ to 3.73?×?10¹⁹?dyne-cm for seismic moment, 0.1 to 9.7 bars for stress drops, and 111.78 to 558.92 m for source radii. The average value of f _max for these events varies from 7 to 18 Hz and seems to be source dependent.     

Morphotectonic analysis of Kosi River basin in Kumaun Lesser Himalaya: an evidence of neotectonics

The Himalayan arc is one of the prominent sites on Earth, for ongoing research on active tectonics because of the frequent occurrence of earthquakes of low-moderate intensity that occur in various sectors of the region. The present study is an attempt to decipher the active uplift and relative tectonic activity in Kosi River basin, a part of the southern Kumaun Himalayas. Several morphotectonic parameters such as asymmetry factor (AF), hypsometric integral (HI), mountain front sinuosity (S_mf), channel sinuosity (S), and basin elongation ratio (R_e) have been calculated with an objective to compare different sub-basins in the Kosi watershed that may prove useful in deciphering of relative tectonic activity. The watershed delineation of the Kosi River basin as well as its sub-basins, and detailed drainage network has been accomplished by using the CARTOSAT-1 DEM with the help of ArcGIS 10.3 software, using TauDEM tool and Global Mapper 18. Analysis of these morphotectonic parameters reveals that although the whole of the Kosi River basin lies in the seismically active zone, but the northern part along North Almora Thrust (NAT), central part around South Almora Thrust (SAT), Ramgarh Thrust (RT), and southern part along Himalayan Frontal Thrust (HFT) are tectonically more active and undergoing neotectonic rejuvenation. The information derived would prove beneficial in identification of hazard prone areas and in planning of socio-economic development in mountainous terrain.     

Siwalik sedimentation in a part of the Kumaun Himalaya,India

The Siwalik Group in a part of the Kumaun Himalaya has been studied with respect to its sedimentologic properties. Size-based environmental data indicate a fluviatile environment for the Middle and Upper Siwalik sediments. The Lower Siwalik samples indicate a border-line environment, possibly a fluvial-deltaic complex. Petrologically, the Siwalik samples are essentially sublitharenites and litharenites. Variation in petrological character in successive Siwalik units is not very marked, although the heavy-mineral assemblages serve the purpose of stratigraphic identification.Sedimentary structures, though not profuse, show a well-developed cyclic development corresponding to the idealised fining-upward sequence of alluvial sediments. They indicate deposition by laterally shifting braided streams. A major portion of the Siwalik detritus may be considered to have its provenance in the Himalayan metamorphic areas.     

Sediment and hydro biogeochemistry of Lake Nainital, Kumaun Himalaya, India

The picturesque Nainital Lake, in the Uttarakhand state of India, is one of the major tourist attractions in the northern part of India. The increasing tourism and population around these lakes are a major concern for the ecology and good sustenance of the lakes. The present study is aimed to understand the behaviour of nutrients and metals in the sediment and their association with chemical forms in the lake. The study was accomplished by studying the water, interstitial water and sediments for major oxides, nutrients and metals in the lake. The different chemical forms of phosphorus and metals in the sediments were done using sequential extraction procedures. The water chemistry (Ca + Mg:Na + K) and the sediment chemistry (CIA and Al₂O₃/K₂O) show that the rocks in the catchment area play an important role in the geochemistry of the lake. The metals in the water also show that the Tallital basin is more polluted than the Mallital basin, may be due to the influence of Bus station. The high concentration of chloride, NH₄, SO₄ and metals in the sediment water interface and the interstitial water shows denitrification, sulfidisation and sulfide oxidation in the anoxic bottom water. The sediment composition shows that the phosphorus in the water is sequestered as carbonate flour apatite, and the metals precipitate as carbonate. The geo- accumulation index shows that the metals zinc, cobalt and nickel show moderate polluted nature than other metals. In general, the lake is less affected by anthropogenic activities. The chemical processes undergoing within the lake, like sulfidisation and sulfide oxidation, oxide dissolution and denitrification and organic matter degradation play an important role in the remobilization of the metals from the lake sediments.     

Petrogenesis of the granitoid rocks from Askot crystallines,Kumaun Himalaya

The Askot crystallines form a doubly plunging synformal belt and occurs as a detached crystalline belt or klippen in the vast sedimentary terrain lying between Central crystallines towards north and the Almora crystallines to the south. It is dominated by granite gneiss and augen gneiss, and also comprise of metapelites, migmatites and basic intrusives. In this paper, the geochemical studies of the granite gneiss and augen gneiss from the Askot crystallines, Kumaun Himalaya were carried out in order to understand their origin and evolution. The granite gneiss is generally foliated, with less foliated and porphyritic variety seen in the core part. The K-feldspar shows Carlsbad twinning, while plagioclases show complex twinning. They show euhedral zircon and apatite along with titanite as accessory minerals. The granite gneiss is moderately evolved (Mg# ∼50) and has granodiorite composition with metaluminous, calc-alkaline trends. They show higher concentration of Ti, Ca, Mg and low abundance of ∑REE (∼165 ppm) in comparison to augen gneiss. They show volcanic arc signatures and compare well with Lateorogenic granites of Proterozoic times distributed world wide. These calc-alkaline granites appear derived from a Paleoproterozoic mafic/intermediate lower-crust reservoir probably involving arc magma underplating. Granite gneiss is also peraluminous with molar A/CNK>1.1, and the heterogeneity of granite gneiss can be explained with the precursor melts, experiencing assimilation during up-rise through crust or contamination of source itself involving sediments from the subduction zone.     

Fault rupture directivity of large earthquakes in Himalaya

The paper examines the predominant fault rupture directivity during large earthquakes in different sectors of the Himalaya which influences strong ground motion and damage scenario. The nature of the faulting of earthquakes vis-à-vis their rupture directivity has been discussed. It is found that the rupture directivity near the Indo-Eurasian plate boundary varies from place to place i.e. either along the strike direction of the faults or at right angles to it. The secondary meizoseismal areas as observed for 1505 Dharchula, 1803 Uttarakhand, 1905 Kangra earthquakes in the Himalaya and 2001 Bhuj earthquake in stable continental region suggest that they are a fairly good indicator of predominant rupture directivity since the latter accentuates the site response up to a longer distance. The resulting larger ground motions, therefore, need to be incorporated in the design of engineering structures by suitable modifications in the BIS code.     

Back-thrusting in Lesser Himalaya: Evidences from magnetic fabric studies in parts of Almora crystalline zone,Kumaun Lesser Himalaya

The present study aims to understand evolution of the Lesser Himalaya, which consists of (meta) sedimentary and crystalline rocks. Field studies, microscopic and rock magnetic investigations have been carried out on the rocks near the South Almora Thrust (SAT) and the North Almora Thrust (NAT), which separates the Almora Crystalline Zone (ACZ) from the Lesser Himalayan sequences (LHS). The results show that along the South Almora Thrust, the deformation is persistent; however, near the NAT deformation pattern is complex and implies overprinting of original shear sense by a younger deformational event. We attribute this overprinting to late stage back-thrusting along NAT, active after the emplacement of ACZ. During this late stage back-thrusting, rocks of the ACZ and LHS were coupled. Back-thrusts originated below the Lesser Himalayan rocks, probably from the Main Boundary Thrust, and propagated across the sedimentary and crystalline rocks. This study provides new results from multiple investigations, and enhances our understanding of the evolution of the ACZ.     

Petrology of thermally metamorphosed pelitic rocks in the Champawat Area,Kumaun Himalaya

In the Champawat area, Kumaun Himalaya, greenschist facies regionally metamorphosed rocksviz chlorite-phyllite and schist have been subjected to thermal metamorphism due to emplacement of batholithic granite/granodiorite body. As a consequence, biotite, garnet, andalusite, fibrolite, sillimanite and perthite minerals have formed in the contact rocks. The conspicuous absence of cordierite and staurolite reported from such aureole rocks is due to higher FeO/MgO ratio of the bulk rock composition in the former while the absence of staurolite is due to low Al₂O₃/FeO+MgO ratio in the schists. AFM diagram demonstrates that in muscovite-bearing schist, the bulk composition of chlorite- and cordierite-bearing rocks are restricted to low FeO/MgO side and thus the restricted occurrence of former and the absence of latter in the contact rocks of the area. This is further evident from the common occurrence of almandine-rich garnet in the rocks.     

Seismicity and source parameters of moderate earthquakes in Sikkim Himalaya

In this study, we accurately relocate 360 earthquakes in the Sikkim Himalaya through the application of the double-difference algorithm to 4?years of data accrued from a eleven-station broadband seismic network. The analysis brings out two major clusters of seismicity??one located in between the main central thrust (MCT) and the main boundary thrust (MBT) and the other in the northwest region of Sikkim that is site to the devastating Mw6.9 earthquake of September 18, 2011. Keeping in view the limitations imposed by the Nyquist frequency of our data (10?Hz), we select 9 moderate size earthquakes (5.3????Ml????4) for the estimation of source parameters. Analysis of shear wave spectra of these earthquakes yields seismic moments in the range of 7.95?×?10²¹ dyne-cm to 6.31?×?10²³ dyne-cm and corner frequencies in the range of 1.8?C6.25?Hz. Smaller seismic moments obtained in Sikkim when compared with the rest of the Himalaya vindicates the lower seismicity levels in the region. Interestingly, it is observed that most of the events having larger seismic moment occur between MBT and MCT lending credence to our observation that this is the most active portion of Sikkim Himalaya. The estimates of stress drop and source radius range from 48 to 389?bar and 0.225 to 0.781?km, respectively. Stress drops do not seem to correlate with the scalar seismic moments affirming the view that stress drop is independent over a wide moment range. While the continental collision scenario can be invoked as a reason to explain a predominance of low stress drops in the Himalayan region, those with relatively higher stress drops in Sikkim Himalaya could be attributed to their affinity with strike-slip source mechanisms. Least square regression of the scalar seismic moment (M ₀) and local magnitude (Ml) results in a relation LogM ₀?=?(1.56?±?0.05)Ml?+?(8.55?±?0.12) while that between moment magnitude (M _w ) and local magnitude as M _w ?=?(0.92?±?0.04)Ml?+?(0.14?±?0.06). These relations could serve as useful inputs for the assessment of earthquake hazard in this seismically active region of Himalaya.     

Probabilities of occurrence of great earthquakes in the Himalaya

Long-term conditional probabilities of occurrence of great earthquakes along the Himalaya plate boundary seismic zone have been estimated. The chance of occurrence of at least one great earthquake along this seismic zone over a period of 100 years (beginning the year 1999) is estimated to be about 0.89. The 100-year probability of such an earthquake occurring in the Kashmir seismic gap is about 0.27, in the central seismic gap about 0.52 and in the Assam gap about 0.21. The 25-year probabilities of their occurrence in these gaps are 0.07, 0.17, and 0.05 respectively. These probability estimates may be used profitably to assess the seismic hazard in the Himalaya and the adjoining Ganga plains.     

Geomorphic positioning and depositional dynamics of river systems in Lower Siwalik basin,Kumaun Himalaya

Qualitative and quantitative analysis of river systems in the Lower Siwalik sequence has enabled characterization of channel patterns, river metamorphosis and resulting sandstone body evolution in time and space. Processes related six lithofacies repeat to generate 8–10 m thick multistoried sandstone complexes deposited in perennial channel belts. Based on lateral mapping of the sandstone bodies, the surfaces of genetic significance ranging from 3^rd, 4^th and 5^th order, suggest presence of meandering, braided and anastomosing river patterns that were responsible for the Lower Siwalik sedimentation. Variation in local base-level in response to allogenic factors including climate and tectonics forced river systems to acquire different patterns. Eustasy seems to control large-scale basin level changes. Quantitatively reconstructed morphological parameters and their comparison with modern and ancient analogues, supported by other independent evidences such as stratigraphical position of sandstone bodies in vertically measured columns and mineralogical characteristics of channel sandstones, enabled to decipher the geomorphic positioning of the Lower Siwalik channels in distal parts of megafan and interfluve areas within the foreland basin setting.     

Recession of Milam Glacier,Kumaun Himalaya,Observed via Lichenometric Dating of Moraines

Glaciers being very sensitive to climate change have been identified as one of the best indicators of climate change and evidences have proved that most of the Himalayan glaciers have receded with an increased rate during the recent past under the influence of global warming. Lichenometric study was carried out on the moraines of Milam glacier (located in Pithoragarh district of Uttarakhand) with the help of lichen species Dimelaena oreina having an average annual growth rate of 1.31 mm. The study revealed that Milam glacier has receded 1450 m in last 69.37 years with an average recession rate of 20.90 m/year. Since lichenometric studies are cost effective and ecofriendly in comparison to carbon dating, satellite and remote sensing based studies and also reliable, hence, it should be promoted in Himalaya which is an abode of glaciers.     

Tectonochemistry and p-t Conditions of Ramgarh and Almora Gneisses from Askot Klippe,Kumaun Lesser Himalaya

The chemical and petrological correlation of metamorphic nappes and klippes overlying the Proterozoic sedimentary units in the Kumaun Himalaya is still debated. The Ramgarh and Almora gneisses, not previously distinguished in the Askot Klippe, show distinct field, petrological and chemical signatures markedly similar to the tectonostratigraphic disposition of the Almora Nappe. A negative Eu anomaly in the Ramgarh granitic gneisses indicates lesser plagioclase fractionation while the Eu anomaly in the Almora pelitic gneisses is likely to have been controlled by feldspar crystallization in restites. During the anatexis at 776°C temperature and 6.6 kbar pressure, the melt moved slightly away to its crystallization sites. The Rb/Sr ratio ?0.54 and Nb ?10 ppm is consistent with the granodioritic composition. The negative Sr anomaly in the underlying Ramgarh granitic gneisses indicates a distinct mantle derived source/plagioclase fractionation with a notable correspondence to other late orogenic granites, particularly the basement Ulleri gneisses from the Nepal Himalaya. Ramgarh gneisses plot in the late-and post-COLG field. The Askot ensemble is likely to be the tectonometamorphically reworked basement, viz. the Ramgarh Group along with its metapelitic cover o f the Almora Group, together comprising southward thrust remnants of the leading edge of the Indian Plate that collided with Tibet during the Tertiary Himalayan orogeny.     

Sediment accumulation rate in the lakes of Kumaun Himalaya,India using210Pb and226Ra

The rate of sedimentation and the source of sediments in the lake basins of Nainital region, Kumaun Himalaya, have been estimated employing²¹⁰Pb and²¹⁰Ra methods. This has yielded a rate of sedimentation of 11.5, 4.70, 3.72, and 3.00 mm/yr in Nainital, Bhimtal, Naukuchiyatal, and Sattal lakes, respectively. The higher rate of sedimentation in Nainital lake, compared to other lakes, is related to faster erosion in the catchment aided by greater anthropogenic activity, while the slowest rate in Sattal lake is due to less erosion and more input of soil-derived material involving a slow rate of accumulation.     

Petrology and geochemistry of the mafic dyke rocks from precambrian almora crystallines of Kumaun Lesser Himalaya

Mafic dykes of Almora region intrude the Precambrian crystalline rocks of Kumaun Lesser Himalaya. Mafic dykes exhibit fine grained margin and medium to coarse grained core, melanocratic, low to highly ferromagnetic (MS=0.85?38.58×10^?3SI) in nature commonly showing subophitic to ophitic textures with ol-pl-cpx-hbl-bt-mt-ap-sp assemblage, and modally correspond to leucogabbro and olivinegabbro (sensu stricto). Olivine (Fo₆₁-Fo₃₃), clinopyroxene (Wo₄₆-En₄₂-Fs₂₂ to Wo₄₀-En₃₆-Fs₁₅) and plagioclase (An₅₈-An₁₂) have crystallized in the temperature range of ca1400–980°C at pressure <2 kbar in an olivine tholeiitic basalt parent. Low acmite (Na_pfu=0.033?0.025), (Mg^#=0.64–0.82), Ti-Al contents of clinopyroxenes and their evolution along enstatite-ferrosilite join (i.e. Mg?Fe substitution) strongly suggest tholeiitic nature of mafic dyke melt with changing activities of alumina and silica. Clinopyroxene compositions of mafic dykes differ markedly as compared to those observed for adjoining Bhimtal volcanics but closely resemble to that crystallized in tholeiitic melts of Deccan province. Observed Cr vs Mg^# variation, enriched LILE (Sr, Ba)-LREE and positive Eu-anomaly of the studied mafic dykes are indicative of fractional crystallization of olivine-clinopyroxene -plagioclase from a crustally-contaminated tholeiitic basalt magma derived from enriched mantle source. The mafic dykes of Almora are geochemically identical to mafic dykes of Nainital, but are unrelated to Precambrian mafic volcanic flow and dykes of NW Himalaya and dykes of Salma and Rajmahal regions.     

Neotectonic Landform Development and Associated Mass Movements along Eastern Ramganga Valley in the Kumaun Himalaya,India

Geotectonics - The present study focuses on the neotectonics activity and related mass movement along the Eastern Ramganga River Basin (ERRB), Kumaun Himalaya, India. To understand the...