State of tectonic stress in Shillong Plateau of northeast India

Tectonic stress regime in the Shillong plateau, northeast region of India, is examined by stress tensor inversion. Some 97 reliable fault plane solutions are used for stress inversion by the Michael and Gauss methods. Although an overall NNW-SSE compressional stress is observed in the area, the stress regime varies from western part to eastern part of the plateau. The eastern part of the plateau is dominated by NNE-SSW compression and the western part by NNW-SSE compression. The NNW-SSE compression in the western part may be due to the tectonic loading induced by the Himalayan orogeny in the north, and the NNE-SSW compression in the eastern part may be attributed to the influence of oblique convergence of the Indian plate beneath the Indo-Burma ranges. Further, Gravitational Potential Energy (GPE) derived stress also indicates a variation from west to east.     

Trends in reference evapotranspiration in the humid region of northeast India

In the present study, the trends in the reference evapotranspiration (ET_O) estimated through the Penman‐Monteith method were investigated over the humid region of northeast (NE) India by using the Mann‐Kendall (MK) test after removing the effect of significant lag‐1 serial correlation from the time series of ET_O by pre‐whitening. During the last 22 years, ET_O has been found to decrease significantly at annual and seasonal time scales for 6 sites in NE India and NE India as a whole. The seasonal decreases in ET_O have, however, been more significant in the pre‐monsoon season, indicating the presence of an element of a seasonal cycle. The decreases in ET_O are mainly attributed to the net radiation and wind speed, which are also corroborated by the observed trends in these two parameters at almost all the times scales over most of the sites in NE India. The steady decrease in wind speed and decline in net radiation not only balanced the impact of the temperature increases on ET_O, but may have actually caused the decreases in ET_O over the humid region of northeast India. Copyright © 2011 John Wiley & Sons, Ltd.     

Seismic source characterization for the Shillong Plateau in Northeast India

The present study aims at understanding the seismotectonic province of the Shillong Plateau (SP) by identifying the potential seismic source zones within a radius of 500 km from the centre of the SP. From existing literature and earthquake (EQ) data, the seismotectonic region is found to vary in terms of seismicity, tectonic features, geology, thickness of overburden, rupture characteristics and rate of movement. Thus, entire 500-km-radius seismotectonic region is divided into four seismic source zones: namely (1) the Shillong Plateau–Assam Valley Zone (SP-AVZ), (2) the Indo-Burma Ranges Zone (IBRZ), (3) the Bengal Basin Zone (BBZ) and (4) the Eastern Himalaya Zone (EHZ). EQ catalogues for each source zone is analysed for completeness of magnitude and time. Seismic parameter b estimated using a maximum likelihood method is found to be 0.91 ± 0.03, 0.94 ± 0.02, 0.80 ± 0.03 and 0.89 ± 0.03 for the SP-AVZ, IBRZ, BBZ and EHZ, respectively. In addition, the maximum likelihood method is used to estimate the mean annual activity rate, maximum possible magnitude (m _max), return period and probability of exceedance for the four zones. The b values estimated suggest that the BBZ is seismically more active; however, the rate of occurrence of EQs is highest in the IBRZ. Findings from this study are an indication of the relative contribution from each of the four seismic source zones towards a seismic hazard of the SP.     

Coupling Between Seismic Activity and Hydrogeochemistry at the Shillong Plateau, Northeastern India

Transient hydrogeochemical anomalies were detected in a granite-hosted aquifer, which is located at a depth of 110 m, north of the Shillong Plateau, Assam, India, where groundwater chemistry is mainly buffered by feldspar alteration to kaolinite. Their onsets preceded moderate earthquakes on December 9, 2004 (M_W = 5.3) and February 15, 2005 (M_W = 5.0), respectively, 206 and 213 km from the aquifer. The ratios [Na+K]/Si, Na/K and [Na+K]/Ca, conductivity, alkalinity and chloride concentration began increasing 3–5 weeks before the M_W = 5.3 earthquake. By comparison with field, experimental and theoretical studies, we interpret a transient switchover between source aquifers, which induced an influx of groundwater from a second aquifer, where groundwater chemistry was dominantly buffered by the alteration of feldspar to smectite. This could have occurred in response to fracturing of a hydrological barrier. The ratio Ba/Sr began decreasing 3–6 days before the M_W = 5.0 earthquake. We interpret a transient switchover to anorthite dissolution caused by exposure of fresh plagioclase to groundwater interaction. This could have been induced by microfracturing, locally within the main aquifer. By comparison with experimental studies of feldspar dissolution, we interpret that hydrogeochemical recovery was facilitated by groundwater interaction and clay mineralization, which could have been coupled with fracture sealing. The coincidence in timing of these two hydrogeochemical events with the only two M_W ≥ 5 earthquakes in the study area argues in favor of cause-and-effect seismic-hydrogeochemical coupling. However, reasons for ambiguity include the lack of similar hydrogeochemical anomalies coupled with smaller seismic events near the monitoring station, the >200 km length scale of inferred seismic-hydrogeochemical coupling, and the potential for far-field effects related to the Great Sumatra–Andaman Islands Earthquake of December 26, 2004.     

Algal Dynamics in Relation to Physico-Chemical Factors in a High Altitude Pond of Shillong,India

A one year phycological study of a high altitude pond has been made with a view to studying the algal dynamics as affected by physico-chemical characteristics of the water. This has revealed the presence of 34 genera of algae belonging to Chlorophyta, Bacillariophyta, Cyanophyta, and Euglenophyta. Chlorophytes were the most dominant and diatoms the most frequent representatives of algae. The community diversity along with water chemistry are indicative of the oligotrophic nature of this pond.     

Structural information derived from ambient noise tomography over a hydrocarbon-producing region in the Cachar fold belt,lower Assam,northeast India

Ambient noise tomography is a powerful tool that has found increasing application in reservoir analysis and imaging. The Cachar fold belt in lower Assam, northeast India encompasses several wells under active hydrocarbon production, along with several dry wells. To overcome the lack of active seismic data over the entire fold belt, a passive seismic study was carried out to image the concealed three-dimensional sub-surface structures. The data were recorded from February to November 2011 by a network of 65 wideband seismometers spanning an area of about 40 × 60 km². The data are crosscorrelated in the 2–5 s band, followed by phase-weighted stacking to estimate noise correlation functions with surface wave signatures. The traveltimes picked from the frequency-time analysis are utilized in a tomographic inversion for Rayleigh wave group velocities. The group velocity anomalies have a lateral resolution of ~ 3.5 × 5.5 km² and variations of up to $\pm 20\%$ for each period. The group velocities are in turn inverted for S-wave velocity distribution as a function of depth. The three-dimensional S-wave velocity tomograms reveal the tight anticlines and broad synclines, with high- and low-velocity zones corresponding to structural highs and lows, respectively. The structural interpretation is supported for the part of the region with producing wells and covered by active seismic data, wherein the post-stack time migrated seismic section shows anticlinal and synclinal features similar to those obtained from ambient noise tomography. The structures revealed by ambient noise tomography can help identify zones of interest to be targeted by active seismic surveys in the Cachar fold belt.     

Earthquake prediction in northeast India—A review

A multidimensional geophysical approach is made to study earthquake precursors in northeast India. Resistivity, gravity and seismological precursor anomalies for the two recent damaging earthquakes, one on December 30, 1984, magnitude 5.8, and the other on August 6, 1988, magnitude 7.5, are very encouraging to intensify these studies in the region. Improved methods of magnetic and geohydrological precursor studies are suggested.     

The Effects of Attenuation and Site on the Spectra of Microearthquakes in the Shillong Region of Northeast India

Microearthquake spectra from the Shillong region are analyzed to observe the effect of attenuation and site on these spectra. The spectral ratio method is utilized to estimate the Q values for both P- and S-waves in the subsurface layer, wherein the ratio of spectral amplitudes at lower and higher frequencies are taken into consideration for three stations at varying epicentral distances. Average estimates of Q _P and Q _S are 178 and 195. The ratio of Q _S to Q _P is estimated to be greater than 1 in major parts of the Shillong area, which can be related to the dry crust prevailing in the Shillong region. Typically, the variation in corner frequencies for these spectra is inferred to be characteristic of the site. Simultaneously, observations from spectral content of local earthquakes recorded at two different stations with respect to the reference site yield greater amplification of incoming seismic signals in the frequency range of 2–5 Hz, which is found to be well supported by the existing local lithology pertinent to that region.     

Ground motion parameters of Shillong plateau: One of the most seismically active zones of northeastern India

Strong ground motion parameters for Shillong plateau of northeastern India are examined. Empirical relations are obtained for main parameters of ground motions as a function of earthquake magnitude, fault type, source depth, velocity characterization of medium and distance. Correlation between ground motion parameters and characteristics of seismogenic zones are established. A new attenuation relation for peak ground acceleration is developed, which predicts higher expected PGA in the region. Parameters of ...     

Ground-motion intensity measure correlations observed in Italian strong-motion records

Ground-motion models (GMMs) are widely used in probabilistic seismic hazard analysis (PSHA) to estimate the probability distributions of earthquake-induced ground-motion intensity measures (IMs) at a site, given an earthquake of a certain magnitude occurring at a nearby location. Accounting for spatial and cross-IM correlations in earthquake-induced ground motions has important implications on probabilistic seismic hazard and loss estimates. This study first develops a new Italian GMM with spatial correlation for 31 amplitude-related IMs, including peak ground acceleration (PGA), peak ground velocity (PGV), and 5%-damped elastic pseudo-spectral accelerations (PSAs) at 29 periods ranging from 0.01 to 4 seconds. The model estimation is performed through a recently developed one-stage nonlinear regression algorithm proposed by the authors, known as the Scoring estimation approach. In fact, current state-of-practice approaches estimate spatial correlation separately from the GMM estimation, resulting in inconsistent and statistically inefficient estimators of interevent and intraevent variances and parameters in the spatial correlation model. We test whether this affects the subsequent cross-IM correlation analysis. To this aim, based on the newly developed GMM, the empirical correlation coefficients from interevent and intraevent residuals are investigated. Finally, a set of analytical correlation models between the selected IMs are proposed. This is of special interest as several correlation models between different IMs have been calibrated and validated based on advanced GMMs and global datasets, lacking earthquakes in extensional regions; however, modeling the correlation between different IM types has not been adequately addressed by current, state-of-the-art GMMs and recent ground-motion records for Italy.     

Scaling relations of moment magnitude,local magnitude,and duration magnitude for earthquakes originated in northeast India

In this study, we aim to improve the scaling between the moment magnitude (M _W), local magnitude (M _L), and the duration magnitude (M _D) for 162 earthquakes in Shillong-Mikir plateau and its adjoining region of northeast India by extending the M _W estimates to lower magnitude earthquakes using spectral analysis of P-waves from vertical component seismograms. The M _W-M _L and M _W-M _D relationships are determined by linear regression analysis. It is found that, M _W values can be considered consistent with M _L and M _D, within 0.1 and 0.2 magnitude units respectively, in 90 % of the cases. The scaling relationships investigated comply well with similar relationships in other regions in the world and in other seismogenic areas in the northeast India region.     

CrustalP-wave velocity and velocity-ratio study in northeast India by a microearthquake survey

The upper crustal (20 km)P-wave velocity beneath the Shillong Plateau and Nowgong area has been studied by the time-distance plot method. TheP-arrival data of the shallow (20 km) microearthquakes from three temporary networks are used, and the average velocity is found to be 5.55 km/s. The velocity ratio (V _p /V _s ) for the upper crust (0–20 km) as well as for the lower crust (21–40 km) are determined by the Wadati-plot method and station-by-station method. The average value obtained by the two methods is compatible; theV _p /V _s ranges between 1.74 to 1.76. A generalized seismic velocity model of the area is suggested by this study, which has been very useful for microearthquake location.     

Attenuation of P,S, and coda waves in Koyna region,India

The attenuation properties of the crust in the Koyna region of the Indian shield have been investigated using 164 seismograms from 37 local earthquakes that occurred in the region. The extended coda normalization method has been used to estimate the quality factors for P waves and S waves , and the single back-scattering model has been used to determine the quality factor for coda waves (Q _c). The earthquakes used in the present study have the focal depth in the range of 1–9 km, and the epicentral distance vary from 11 to 55 km. The values of and Q _c show a dependence on frequency in the Koyna region. The average frequency dependent relationships (Q = Q ₀ f ⁿ) estimated for the region are , and . The ratio is found to be greater than one for the frequency range considered here (1.5–18 Hz). This ratio, along with the frequency dependence of quality factors, indicates that scattering is an important factor contributing to the attenuation of body waves in the region. A comparison of Q _c and in the present study shows that for frequencies below 4 Hz and for the frequencies greater than 4 Hz. This may be due to the multiple scattering effect of the medium. The outcome of this study is expected to be useful for the estimation of source parameters and near-source simulation of earthquake ground motion, which in turn are required in the seismic hazard assessment of a region.     

Peak ground motion predictions in India: an appraisal for rock sites

Proper selection and ranking of Ground Motion Prediction Equations (GMPEs) is critical for successful logic-tree implementation in probabilistic seismic hazard analysis. The present study explores this issue in predicting peak ground accelerations at the rock sites in India. Macroseismic intensity data complemented with limited strong ground-motion recordings are used for the purpose. The findings corroborate the possible conformity between the GMPEs developed for tectonically active shallow crust across the globe. On the other hand, the relevant GMPEs in the intraplate regions cluster into two different groups with the equations of lower ranks catering to higher ground motions. The earthquakes in the subduction zones have significant regional implications. However, affinity in the ground-motion attenuations between the major interface events (M _W > 7.4) in Andaman-Nicobar, Japan and Cascadia, respectively, is noted. This can be also observed for the intraslab events in the Hindukush and Taiwan respectively. Overall, we do not observe any significant advantage with the equations developed using the regional data. These findings are expected to be useful in probabilistic seismic hazard analysis across the study region.     

Observations on the Abiotic Factors and Planktonic Periodicity in a Shallow Pond of the Highlands of Shillong (India)

The pond, situated at an altitude of 4862 m, shows an annual variation of temperature between 15 and 27 °C in a close correlation with the air temperature. In the case of a weakly acidic reaction of water an inverse relation exists between pH and conductivity and between free CO₂ and alkalinity. The nutrient content of 0.05 mg/l PO₄-P and NO₃-N each is low, the phosphorus concentration being maintained by re-solution from the sediment. Chlorophyceae and diatoms are dominant in the phytoplankton alternately, and rotatoria and copepoda are dominant in the zoo-plankton in winter and in summer, respectively. The pronounced changes in the composition of plankton do not show any definite dependence on the physico-chemical factors, particularly since a significant share of the nutrient metabolism and the oxygen balance in the pond investigated is determined by macrophytes.     

Ground-motion characterization for the probabilistic seismic hazard assessment in Turkey

This study describes the methodology implemented to establish the ground-motion logic-tree for national probabilistic seismic hazard map of Turkey for shallow active crustal regions. The presented procedure provides quantitative information to guide the hazard experts while establishing the logic tree to capture the epistemic uncertainty in ground-motion characterization. It uses non-data-driven and data-driven testing methods to identify and rank candidate ground-motion prediction equations (GMPEs) under a specific ground-motion database. The candidate GMPEs are subjected to visual inspection and are classified into center, body and range (CBR) spectral estimates for a proper consideration of epistemic uncertainty. The GMPEs classified into CBR are then used in a suite of seismic hazard sensitivity analysis to establish the most suitable GMPE logic-tree whose spectral estimates are not biased by any one of the GMPEs in the logic-tree structure. The sensitivity analysis considers normalized spectral ordinates and is not manipulated by the spectral amplitudes. The proposed procedure is inherited from the relevant studies of the Earthquake Model of the Middle East (EMME; www.efehr.org:8080/jetspeed/portal/emme.psml) regional seismic hazard project. This paper also highlights the similarities and differences in ground-motion characterization between EMME and our approach.     

Recent seismicity in Northeast India and its adjoining region

Recent seismicity in the northeast India and its adjoining region exhibits different earthquake mechanisms – predominantly thrust faulting on the eastern boundary, normal faulting in the upper Himalaya, and strike slip in the remaining areas. A homogenized catalogue in moment magnitude, M _W, covering a period from 1906 to 2006 is derived from International Seismological Center (ISC) catalogue, and Global Centroid Moment Tensor (GCMT) database. Owing to significant and stable earthquake recordings as seen from 1964 onwards, the seismicity in the region is analyzed for the period with spatial distribution of magnitude of completeness m _t, b value, a value, and correlation fractal dimension D _C. The estimated value of m _t is found to vary between 4.0 and 4.8. The a value is seen to vary from 4.47 to 8.59 while b value ranges from 0.61 to 1.36. Thrust zones are seen to exhibit predominantly lower b value distribution while strike-slip and normal faulting regimes are associated with moderate to higher b value distribution. D _C is found to vary from 0.70 to 1.66. Although the correlation between spatial distribution of b value and D _C is seen predominantly negative, positive correlations can also be observed in some parts of this territory. A major observation is the strikingly negative correlation with low b value in the eastern boundary thrust region implying a possible case of extending asperity. Incidentally, application of box counting method on fault segments of the study region indicates comparatively higher fractal dimension, D, suggesting an inclination towards a planar geometrical coverage in the 2D spatial extent. Finally, four broad seismic source zones are demarcated based on the estimated spatial seismicity patterns in collaboration with the underlying active fault networks. The present work appraises the seismicity scenario in fulfillment of a basic groundwork for seismic hazard assessment in this earthquake province of the country.     

On the geodetic observations during 1965 to 1969 in the Koyna region, India

Limited reports of repeat triangulation and levelling surveys carried out in the Koyna region by Survey of India, before and after the 1967 Koyna earthquake, are available. We analyse these observations of horizontal and vertical displacements to investigate whether analysis of these observations can provide additional constraints on the source parameters of the 1967 Koyna earthquake. We calculate surface displacements due to coseismic slip on the subsurface rupture of the 1967 Koyna earthquake and also due to the Koyna reservoir water load. We find that the reported displacements are too large to be attributed to the coseismic displacements and/or to the displacements induced by the reservoir water load. We conclude that these observations contain large random and systematic errors.     

Probabilistic assessment of earthquake recurrence in the January 26, 2001 earthquake region of Gujrat, India

Kutch region of Gujrat is one of the most seismic prone regions of India. Recently, it has been rocked by a large earthquake (M _w = 7.7) on January 26, 2001. The probabilities of occurrence of large earthquake (M≥6.0 and M≥5.0) in a specified interval of time for different elapsed times have been estimated on the basis of observed time-intervals between the large earthquakes (M≥6.0 and M≥5.0) using three probabilistic models, namely, Weibull, Gamma and Lognormal. The earthquakes of magnitude ≥5.0 covering about 180 years have been used for this analysis. However, the method of maximum likelihood estimation (MLE) has been applied for computation of earthquake hazard parameters. The mean interval of occurrence of earthquakes and standard deviation are estimated as 20.18 and 8.40 years for M≥5.0 and 36.32 and 12.49 years, for M≥6.0, respectively, for this region. For the earthquakes M≥5.0, the estimated cumulative probability reaches 0.8 after about 27 years for Lognormal and Gamma models and about 28 years for Weibull model while it reaches 0.9 after about 32 years for all the models. However, for the earthquakes M≥6.0, the estimated cumulative probability reaches 0.8 after about 47 years for all the models while it reaches 0.9 after about 53, 54 and 55 years for Weibull, Gamma and Lognormal model, respectively. The conditional probability also reaches about 0.8 to 0.9 for the time period of 28 to 40 years and 50 to 60 years for M≥5.0 and M≥6.0, respectively, for all the models. The probability of occurrence of an earthquake is very high between 28 to 42 years for the magnitudes ≥5.0 and between 47 to 55 years for the magnitudes ≥6.0, respectively, past from the last earthquake (2001).