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

Journal of Seismology -     

渤海地区地震危险性特征及对工程抗震设防的启示

本文分析了渤海海域地震构造环境及地震活动特征,总结了近年来对该地区潜在震源区划分的一些新认识,分析了地震危险性区划结果.在此基础上,引入地震危险性特征参数K对研究区进行地震危险性分区,通过不同场点的地震危险性曲线证明了地震环境造成的危险性差异,并以研究区内3个概率水准的地震动峰值加速度的比值统计证明了引入地震危险性分区的必要性.研究表明:(1)不同概率水准的地震动参数之间不存在固定的比例关系,受地震环境影响,地震危险性分区和地震动参数区划图结合使用可以更准确地反应地震危险性特征;(2)对于一般工程,根据地震动参数区划图所给的50年超越概率10%的地震动参数折算所得到的小震的地震作用不具有统一的概率水准,抗震设防目标的实现没有科学保证;海洋平台抗震设计中利用强度水平地震推算得到韧性水平地震作用也不具有统一的概率水准;(3)将地震危险性分区和地震动参数区划结合使用或编制统一抗倒塌水准的区划图是一般工程抗震设防的科学做法,对于海洋平台的抗震设防,宜明确韧性水平地震作用的概率水准.     

Seismic zoning of the southern slope of Greater Caucasus from the fractal parameters of the earthquakes, stress state, and GPS velocities

By complex analysis of GPS velocities, seismicity, fractal dimensions of the spatial distribution of seismic epicenters, focal mechanisms of the earthquakes, and stress state of the Earth’s crust, four seismic zones (Balaken-Zagatala, Sheki-Gabala, Shamakhy-Ismailly, and Absheron) are revealed within the southern slope of the Greater Caucasus. The suggested method can be used as a criterion in seismotectonic zoning; it could also be useful in the assessment of seismic hazards in the collision zones.     

Environmetric and GIS techniques for hydrochemical characterization of the Dal lake,Kashmir Himalaya,India

This research explains the background processes responsible for the spatial distribution of hydrochemical properties of the picturesque eutrophic Himalayan Lake, Dal, located in Kashmir valley, India. Univariate and multivariate statistical analyses were used to understand the spatiotemporal variability of 18 hydrochemical parameters comprising of 12,960 observations collected from 30 sampling sites well distributed within the lake at a grid spacing of 1 km² from March 2014 to February 2016. Hierarchical cluster analysis (HCA) grouped all the sampled data into three clusters based on the hydrochemical similarities, Discriminant analysis also revealed the same clusters and patterns in the data, validating the results of HCA. Wilk’s λ quotient distribution revealed the contribution of ions, nutrients, secchi disk transparency, dissolved oxygen and pH in the formation of clusters. The results are in consonance with the Principal Component Analysis of the whole lake data and individual clusters, which showed that the variance is maximally explained by the ionic component (46.82%) followed by dissolved oxygen and pH (9.36%), nitrates and phosphates (7.33%) and Secchi disk transparency (5.98%). Spatial variability of the hydrochemistry of the lake is due to the variations in water depth, lake water dynamics, flushing rate of water, organic matter decomposition, and anthropogenic pressures within and around the Dal lake ecosystem. Overall, the water quality of the lake is unfit for drinking due to the presence of coliform bacteria in the lake waters.     

Linking human-biophysical interactions with the trophic status of Dal Lake,Kashmir Himalaya,India

The study analyses the long-term biophysical and demographic changes in Dal lake, located in the heart of Srinagar city, Kashmir India, using a repository of historical, remote sensing, socio-economic and water quality data supported by the extensive field observations. The lake faces multiple pressures from the unplanned urbanization, high population growth, nutrient load from intensive agriculture and tourism. The data showed that the lake has shrunk from 31 km² in 1859–24 km² in 2013. Significant changes were observed in the land use and land cover (LULC) within the lake (1859–2013) and in the vicinity of the lake (1962–2013). Analysis of the demographic data indicates that the human population within the lake has shown more than double the national growth rate. Additionally, 7 important water quality parameters from 82 well distributed sites across the lake were analyzed and compared with the past data to determine the historical changes in the water quality from 1971 to 2014. The changes in the LULC and demography have adversely affected the pollution status of this pristine lake. Ortho-phosphate phosphorous concentration has increased from 16.75 μg L⁻¹ in 1977–45.78 μg L⁻¹ in 2014 and that of the nitrate-nitrogen from 365 μg L⁻¹ to 557 μg L⁻¹, indicating nutrient enrichment of the lake over the years. Built-up area within the lake has increased 40 times since 1859, which, together with the changes in the population and settlements, have led to the high discharge of untreated nutrient-rich sewage into the lake. Similarly the expansion of floating gardens within the lake and agriculture lands in the catchment has contributed to the increased nutrient load into the lake due to the increasing use of fertilizers. The information about the existing land cover, demography and water quality was integrated and analyzed in GIS environment to identify the trophic status of the lake. The analysis indicated that 32% of the lake falls under sever degradation, 48% under medium degradation while as 20% of the lake waters are relatively clean. It is believed that the results provide improved knowledge and insights about the lake health and causal factors of its degradation necessary for effectively restoring its ecological and hydrological functionality.     

Strong motion envelope modelling of the source of the Chamoli earthquake of March 28, 1999 in the Garhwal Himalaya, India

Garhwal Himalaya has been rocked by two major earthquakes in the span of just eight years, viz. Uttarkashi earthquake of 20th Oct, 1991 and Chamoli earthquake of 28th March, 1999. Chamoli earthquake of March 28, 1999 was recorded at 11 different stations of a strong motion array installed in the epicentral region. The maximum peak ground acceleration (353 cm/s²) was recorded at an accelerograph located at Gopeshwar. The data from eleven stations has been used for comparison with the simulated acceleration envelopes due to a model of the rupture responsible for this earthquake. For simulation of acceleration envelope the method of Midorikawa (1993) has been modified for its applicability to Himalayan region. This method has earlier been used by Joshi and Patel (1997) and Joshi (1999) for the studyof Uttarkashi earthquake of 20th Oct, 1991. The same method has been used for study of Chamoli earthquake. Layered earth crust has been introduced in place of homogeneous one in this method. The model of rupture is placed at a depth of 12 km below the Munsiari thrust for modelling Chamoli earthquake. Peak ground acceleration was calculated from simulated acceleration envelope using layered as well as homogeneous earth crust. For the rupture placed in a layered crust model peak ground acceleration of order 312 cm/s² was simulated at Gopeshwar which is quite close to actually recorded value. The comparison of peak ground acceleration values in terms of root mean square error at eleven stations suggests that the root mean square error is reduced by inclusion of a layered earth crust in place of homogeneous earth crust.     

提取视震源时间函数的PLD方法及其对2005年克什米尔Mw7．6地震的应用

通过重构用于确定视震源时间函数有效持续时间的判别函数,对提取视震源时间函数的PLD方法进行了改进;利用合成资料和实际资料,验证了改进后PLD方法的可行性和稳定性.将PLD方法应用于2005年克什米尔Mw7.6地震及其11个余震的1887条记录,在84个台站处获得了这次地震的视震源时间函数.分别平均从不同台站的P波、S波、Rayleigh波和Love波中得到的视震源时间函数,获取了主震的平均视震源时间函数.对视震源时间函数的分析表明,2005年克什米尔Mw7.6地震的持续时间大约为25 S,这是一次"急始型"地震,总体上表现为圆盘形破裂.但有迹象表明,破裂在初期有向西北方向发展的单侧传播趋势.     

Anthropogenic Eutrophication and Red Tide Outbreak in Lacustrine Systems of the Kashmir Himalaya

Red tide, a recurrent phenomenon has become conspicuous in several Kashmir lake ecosystems since 1991. The responsible organism (Euglena pedunculata), a rare flagellate rediscovered in the Kashmir Himalaya (Khan 1993) caused first and unprecedented red tide outbreak, constituting a maximum of 96% of resident numerical phytoplankton density in Dal Lake. At present, conflicting hypotheses exist on the generation of causal assemblage(s) imparting redness to waters: Jeeji Bai (1991) linked its origin to acid precipitation – a fallout of burning oil‐fields during the Gulf War – whilst Khan (1993) holds local factor(s) responsible. Field/experimental studies support the latter contention that the influx of untreated sewage, in unison with warm temperatures, high levels of PhAR, iron and interruption to hydrological flow‐pattern together with absence/or reduction in grazing activity created conducive environmental milieu for red tide outbreak. Dal Lake “red tide” drifted the bloom‐inoculum to other waters, including Lake Wular, where additional ecological niches were carved out, threatening the aesthetic value and biological diversity of Kashmir lakes. Ecological monitoring indicates frequent seasonal red tide occurrence in Dal Lake (including summer‐autumn event of 1998) which testifies its unabated eutrophication status. Further studies are needed on ecological adaptability and biogeographic distribution of this rare and unique red tide‐causing flagellate.     

Geological setting of the 8 October 2005 Kashmir earthquake

The source of the 8 October 2005 earthquake of M 7.6 was the northwest-striking Balakot–Bagh (B–B) fault, which had been mapped by the Geological Survey of Pakistan prior to the earthquake but had not been recognized as active except for a 16-km section near Muzaffarabad. The fault follows the Indus–Kohistan Seismic Zone (IKSZ); both cut across and locally offset the Hazara–Kashmir Syntaxis defined by the Main Boundary and Panjal thrusts. The fault has no expression in facies of the Miocene–Pleistocene Siwalik Group but does offset late Pleistocene terrace surfaces in Pakistan-administered Jammu-Kashmir. Two en-échelon anticlines near Muzaffarabad and Balakot expose Precambrian Muzaffarabad Limestone and are cut by the B–B fault on their southwest sides, suggesting that folding and exposure of Precambrian rocks by erosion accompanied Quaternary displacement along the fault. The B–B fault has reverse separation, northeast side up; uplift of the northeast side accompanied displacement, producing higher topography and steeper stream gradients northeast of the fault. No surface expression of the B–B fault has been found northwest of the syntaxis, although the IKSZ and steeper stream gradients continue at least as far as the Indus River, the site of the Pattan earthquake of M 6.2 in 1974. To the southeast, northwest-striking faults were mapped by the Geological Survey of Pakistan. One of these faults, the Riasi thrust, cuts across the southwest flank of an anticline exposing Precambrian limestone. Farther southeast, in Indian-administered territory, Holocene activity on the Riasi thrust has been described. In the Kangra reentrant still farther southeast, active faulting may follow the Soan thrust, along which Holocene and Pleistocene offsets have been described. The Soan thrust, rather than the south flank of the Janauri anticline, may represent the surface projection of the 1905 Kangra earthquake of M 7.8.     

Seismic Wave Attenuation in the Greater Cairo Region, Egypt

In the present study, a digital waveform dataset of 216 local earthquakes recorded by the Egyptian National Seismic Network (ENSN) was used to estimate the attenuation of seismic wave energy in the greater Cairo region. The quality factor and the frequency dependence for Coda waves and S-waves were estimated and clarified. The Coda waves (Q _c) and S-waves (Q _d) quality factor were estimated by applying the single scattering model and Coda Normalization method, respectively, to bandpass-filtered seismograms of frequency bands centering at 1.5, 3, 6, 12, 18 and 24?Hz. Lapse time dependence was also studied for the area, with the Coda waves analyzed through four lapse time windows (10, 20, 30 and 40?s). The average quality factor as function of frequency is found to be Q _c?=?35?±?9f ^0.9±0.02 and Q _d?=?10?±?2f ^0.9±0.02 for Coda and S-waves, respectively. This behavior is usually correlated with the degree of tectonic complexity and the presence of heterogeneities at several scales. The variation of Q _c with frequency and lapse time shows that the lithosphere becomes more homogeneous with depth. In fact, by using the Coda Normalization method we obtained low Q _d values as expected for a heterogeneous and active zone. The intrinsic quality factor (Q _i ^?1 ) was separated from the scattering quality factor (Q _s ^?1 ) by applying the Multiple Lapse Time Domain Window Analysis (MLTWA) method under the assumption of multiple isotropic scattering with uniform distribution of scatters. The obtained results suggest that the contribution of the intrinsic attenuation (Q _i ^?1 ) prevails on the scattering attenuation (Q _s ^?1 ) at frequencies higher than 3?Hz.     

Prediction of the maximum credible ground motion in Singapore due to a great Sumatran subduction earthquake: the worst‐case scenario

Although Singapore is located in a low‐seismicity region, huge but infrequent Sumatran subduction earthquakes might pose structural problems to medium‐ and high‐rise buildings in the city. Based on a series of ground motion simulations of potential earthquakes that may affect Singapore, the 1833 Sumatran subduction earthquake (M_w=9.0) has been identified to be the worst‐case scenario earthquake. Bedrock motions in Singapore due to the hypothesized earthquake are simulated using an extended reflectivity method, taking into account uncertainties in source rupture process. Random rupture models, considering the uncertainties in rupture directivity, slip distribution, presence of asperities, rupture velocity and dislocation rise time, are made based on a range of seismologically possible models. The simulated bedrock motions have a very long duration of about 250 s with a predominant period between 1.8 and 2.5 s, which coincides with the natural periods of medium‐ and high‐rise buildings widely found in Singapore. The 90‐percentile horizontal peak ground acceleration is estimated to be 33 gal and the 90‐percentile horizontal spectral acceleration with 5% damping ratio is 100 gal within the predominant period range. The 90‐percentile bedrock motion would generate base shear force higher than that required by the current design code, where seismic design has yet to be considered. This has not taken into account effects of local soil response that might further amplify the bedrock motion. Copyright © 2002 John Wiley & Sons, Ltd.     

提取视震源时间函数的PLD方法及其对2005年克什米尔MW7.6地震的应用

通过重构用于确定视震源时间函数有效持续时间的判别函数,对提取视震源时间函数的PLD方法进行了改进;利用合成资料和实际资料,验证了改进后PLD方法的可行性和稳定性.将PLD方法应用于2005年克什米尔M_W7.6地震及其11个余震的1887条记录,在84个台站处获得了这次地震的视震源时间函数.分别平均从不同台站的P波、S波、Rayleigh波和 Love波中得到的视震源时间函数,获取了主震的平均视震源时间函数.对视震源时间函数的分析表明,2005年克什米尔M_W7.6地震的持续时间大约为25 s,这是一次“急始型”地震,总体上表现为圆盘形破裂.但有迹象表明,破裂在初期有向西北方向发展的单侧传播趋势.     

云南耿马7.2级地震破裂与震源应力场

1 98 8年 1 1月 6日在云南省澜沧县和耿马县境内先后发生了 7.6级和 7.2级强烈地震。其中 7.2级地震在震中区小黑江两岸灰岩裸露区出现了具有构造意义的地裂缝带。采用断层滑动矢量法反演出形成上述地裂缝带的应力张量为 :σ1走向 2 2 3°,倾角 2 8°;σ2 走向 42°,倾角6 2°;σ3 走向 1 3 3°,倾角 0°;应力比 R=0 .5 7,总体显示出走滑兼压性特征 ,基本代表了震源机制解所反映的震源应力场特征。断面滑动矢量解析可做为一种判定地裂缝是否是构造成因的方法     

Study of seismicity in the NW Himalaya and adjoining regions using IMS network

The Reviewed Event Bulletin (REB) of the International Data Center (IDC) has been used in order to investigate the seismicity of the Northwest Himalaya and its neighboring region for the time period June 1999 to March 2015 within the geographical coordinates 25–40° N latitude and 65–85° E longitude. We have used a very precisely located earthquake dataset recorded by the International Monitoring System (IMS) Network containing 7,583 events with body wave magnitudes from 2.5 to 6.3. The study area has been subdivided into six regions based on the Flinn-Engdahl (F-E) seismic and geographical regionalization scheme, which was used as the region classifications of the International Data Center catalog. The examined region includes NW India, Pakistan, Nepal, Xizang, Kashmir, and Hindukush. For each region, Magnitudes of completeness (Mc) and Gutenberg-Richter (GR) recurrence parameters (a and b values) have been estimated. The Gutenberg-Richter analysis is preceded by an overview of the seismotectonics of the study area. The obtained Mc values vary from 3.5 to 3.9. The lower value of Mc was found mainly in Xizang region whereas the higher Mc threshold is evident in Pakistan region. However, the b values vary from 1.19 to 1.48. The lowest b value is recorded in Xizang region, which is mostly related to the Main Karakoram Thrust (MKT) fault, whereas the highest b values are recorded in NW India and Kashmir regions, which are mostly related to the Main Frontal Thrust (MFT) fault. The REB for the selected period has been compared to the most renowned bulletin of global seismicity, namely that issued by the National Earthquake Information Center (NEIC) of the United States Geological Survey (USGS). A study of 4,821 events recorded by USGS in the study region indicates that about 36 % of seismic events were missed and the catalog is considered as complete for events with magnitudes ≥4.0. However, both a and b values are obviously higher than those of IMS catalog. The a and b parameters in the Gutenberg-Richter magnitude–frequency relationship have been utilized to forecast the probability of future earthquakes of different magnitudes and returned periods (recurrence intervals).     

Late Quaternary tectono-geomorphological investigations of Zanskar Basin,NW Himalaya,India, using geospatial techniques

Geospatial techniques play a crucial role in geomorphic studies, particularly in the challenging terrains like mountainous regions, inaccessible areas and densely vegetated landscapes, where geomorphic features cannot be recorded easily. Tectono-geomorphologic observations provide important clues regarding the landscape evolution, morpho-dynamics and ongoing tectonism of the region. The present study has been carried out in the Zanskar Basin (ZB), located to the south of the Indus Tsangpo Suture Zone (ITSZ), in the hinterland of the NW Himalaya. This study has been carried out to assess and evaluate active tectonics by employing tectono-geomorphic analysis, dynamics in drainage networks, geomorphological field observations and the Geographic Information System (GIS) environment. High-resolution satellite images, topographic maps and the Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) were used to generate primary data sets, which were corroborated with field investigations for valid inferences. The geometry of the ZB suggests that continuous tectonic activity exerts first-order control on the overall shape, size and structure of the ZB. This first-order response is clearly reflected in the landforms modified by tectonic processes, namely, linear mountain fronts, elongated shape and tilting of the basin, braided and meandering river courses and lower stream length gradient index values in hard rock terrain. The ZB exhibits several eye-catching geomorphic features, such as well-defined triangular facets with wide base lengths and wine-glass valleys with small outlets along the footwall block of the Zanskar Shear Zone/South Tibetan Detachment System (ZSZ/STDS), as well as the presence of wind gaps, water gaps, bedrock incision, incised and entrenched valleys, narrow gorges and a high incision rate inferring active tectonics and recent uplift in the region. In addition, the existence of uplifted river terraces, as well as the stepped morphology of fans and strath terraces, suggests that the region is experiencing recent activity and ongoing tectonic uplift. These modified geomorphic characteristics suggest that the hinterland, which is part of the NW Himalaya, is tectonically quite active and has experienced a differential rate of tectonics during its evolution. The quantified geomorphic indices and their relations with the tectonics, climate and erosion activity infer that the basin geometry is mostly controlled by the ZSZ/STDS that dips 20°–70° NE, the south-dipping Zanskar Counter Thrust (ZCT) and other local tectonic elements like the Choksti Thrust (CT), Stondgey Thrust, Zangla Thrust and tectonic structures. The synergised results of quantified geomorphic indices and tectono-geomorphic evidence in the ZB strongly indicate that both the past and ongoing tectonism have significantly shaped and modified geomorphology of the ZB.     

GPS measurements of postseismic deformation due to October 8, 2005 Kashmir earthquake

Relaxation of the coseismic stresses following an earthquake causes postseismic crustal deformation, which can last for days to years. Continuous monitoring of postseismic deformation facilitates the understanding of the mechanism of deformation and postseismic relaxation and viscous rheology. After the October 8, 2005 Kashmir earthquake, global positioning system data for 8 months, starting from October, 2005 have been analyzed from three continuous sites located at Gulmarg, Amritsar, and Jaipur. The average velocity during the observation period at Gulmarg (8.6 cm/year) is significantly higher than the Indian plate velocity exhibiting postseismic crustal deformation. The velocity at Amritsar (5.9 cm/year) and Jaipur (5.1 cm/year) is comparable to the Indian plate velocity. At Gulmarg, the logarithmic function fits well to the north–south component of postseismic transients (~in the coseismic slip direction). The nature of decay in these transients suggests that the deformation is mainly due to an afterslip, and the second possible contribution may be from the viscous relaxation process. This paper presents the characteristics of postseismic transients and possible contributions from various postseismic mechanisms subsequent to the Kashmir earthquake.     

Variation of Seismic Coda Wave Attenuation in the Garhwal Region, Northwestern Himalaya

Seismic coda wave attenuation ( $ Q_{\text{c}}^{ - 1} $ ) characteristics in the Garhwal region, northwestern Himalaya is studied using 113 short-period, vertical component seismic observations from local events with hypocentral distance less than 250?km and magnitude range between 1.0 to 4.0. They are located mainly in the vicinity of the Main Boundary Thrust (MBT) and the Main Central Thrust (MCT), which are well-defined tectonic discontinuities in the Himalayas. Coda wave attenuation ( $ Q_{\text{c}}^{ - 1} $ ) is estimated using the single isotropic scattering method at central frequencies 1.5, 3, 5, 7, 9, 12, 16, 20, 24 and 28?Hz using several starting lapse times and coda window lengths for the analysis. Results show that the ( $ Q_{\text{c}}^{ - 1} $ ) values are frequency dependent in the considered frequency range, and they fit the frequency power law ( $ Q_{\text{c}}^{ - 1} \left( f \right) = Q_{0}^{ - 1} f^{ - n} $ ). The Q ₀ (Q _c at 1?Hz) estimates vary from about 50 for a 10?s lapse time and 10?s window length, to about 350 for a 60?s lapse time and 60?s window length combination. The exponent of the frequency dependence law, n ranges from 1.2 to 0.7; however, it is greater than 0.8, in general, which correlates well with the values obtained in other seismically and tectonically active and highly heterogeneous regions. The attenuation in the Garhwal region is found to be lower than the Q _c ^?1 values obtained for other seismically active regions of the world; however, it is comparable to other regions of India. The spatial variation of coda attenuation indicates that the level of heterogeneity decreases with increasing depth. The variation of coda attenuation has been estimated for different lapse time and window length combinations to observe the effect with depth and it indicates that the upper lithosphere is more active seismically as compared to the lower lithosphere and the heterogeneity decreases with increasing depth.     

Geochemical implications for the alteration of the Dras volcanic rocks from the ophiolites of Indus Suture Zone,Kashmir Himalaya

The Dras volcanic rocks form a part of the ophiolite belt along the Indus Suture Zone in the Kashmir Himalaya. These volcanic rocks have suffered alteration as in any other ophiolite zone. Three types of alterations spilitisation, submarine weathering and ridongitisation were suggested. The spilitic mineralogy appears to be secondary and must have developed due to the reaction of these rocks with hot sea water. Depletion of MgO and CaO and enrichment of K₂O of these rocks relate to the submarine weathering at lower temperatures. Rodingitisation effect is reflected in the chemistry of some rocks with enriched CaO and depleted SiO₂. The trace elements — Co, Cr, Ni, Cu, Pb and Rb do not show any considerable changes during alteration.     

Statistical analysis on yearly seismic moment release data to demarcate the source zone for an impending earthquake in the Himalaya

Tectonism in the Himalayan fold-thrust belt had generated great earthquakes in the past and will spawn more in the future. Sequential cumulative moment release data of macroearthquakes (Mb ≥ 4.5) over the years 1964–2006 in four zones of the Himalaya was analysed by nonparametric RUD method. The Z values of RUD analysis had neither rejected nor supported the null hypothesis of randomness. However, the Hurst analysis and plot, a statistical procedure to identify clustering of low and high values in a time series, brought out a pattern for earthquake prognostication. The pattern was a negative sloping segment representing a sluggish moment release over years, followed by a positive sloping segment indicating a sudden high moment release with occurrence of medium/large size earthquake(s). In recent past, such a negative sloping has been found in Zones B (1992–2006) and D (1998–2006), indicating an impending moderate/mega earthquake event in near future.