Ground-motion predictions in Shillong region, northeast India

We deliver ground-motion prediction equations for Shillong region, northeast India, based on a database generated by finite-fault stochastic simulations. An examination of the regional seismic source characteristics is carried out beforehand. Micro/minor earthquakes (M _W?<?5.0) nucleating at hypocentral depth <21?km in the region recorded at broadband seismic stations are observed to have Brune stress-drop ranging between 2.8 and 99.9?bars. Likewise, macroseismic intensity data for the 1897 Shillong Earthquake that nucleated at a hypocentral depth of ~35?km places the associated stress-drop at 100?200?bars. The apparent variation of the stress-drop parameter with depth is considered with two source zones namely lower-crust and upper-crust. Equations for the lower-crust predict higher ground-motion levels and exhibit affinity to those developed for stable continental region of Eastern North America. The ground-motion levels predicted by the equations for the upper-crust are relatively lower but are still higher compared with those predicted for tectonically active regions, viz., the Himalayas and Western North America.     

A regularized Wiener–Hopf filter for inverting models with magnetic susceptibility

For a magnetic target, the spatial magnetic signal can be expressed as a convolutional integral over Green's function of an assumed model with susceptibility as its parameter. A filter can be used to obtain the susceptibility by minimizing the mismatch between observed and the computed magnetic anomalies. In this perspective, we report the development of an advanced digital filter, which is efficient and can be used to map rock susceptibility from the acquired magnetic data. To design the new filter, we modified the space‐domain standard Wiener–Hopf filter by imposing two different constraints: (i) the filter energy constraint; and (ii) normalization of the filter coefficients. These constraints make it capable to characterize source bodies from their produced magnetic anomalies. We assume that the magnetic data are produced by induced magnetization only and interpretation can be as good as the subsurface model. Our technique is less sensitive to the data noise, which makes it efficient in enhancing the interpretation model. The modified filter demonstrates its applicability over the synthetic data with additive white Gaussian noise. In order to check the efficacy and adaptivity of this tool in a more realistic perspective, it is also tested on the real magnetic data acquired over a kimberlitic district adjoining to the western margin of the Cuddapah Basin in India to identify the source bodies from the anomalies. Our result shows that the modified Wiener–Hopf filter with the constraint for the magnetic data is more stable and efficient than the standard Wiener–Hopf filter.     

Microtremor survey in Talchir, India to ascertain its basin characteristics in terms of predominant frequency by Nakamura's ratio technique

The Talchir Basin, one of India's oldest basins, has been a subject of interest because of its rich coal deposits. The maximum thickness of the basin is about 1500 m. Beyond the basin is the hard metamorphic rocks of Precambrian age. The ambient noise survey data have been analyzed for the Talchir Basin using Nakamura's technique of horizontal–vertical-spectral-ratio (HVSR) to ascertain the basin structure in terms of the predominant frequency. The predominant frequency varies from 0.25 Hz to 7.8 Hz but a major portion of the basin comes under the range of 0.3 Hz–2.4 Hz while on the metamorphic rocks it is as high as 7.8 Hz. The variation in predominant frequency shows a good correlation with the sediment thickness of the basin. The results have been compared with the previous studies by other researchers and it shows consistency with the northerly dip of the basin. The present study has also been compared with the results of the synthetic seismogram that was performed for the Talchir Basin. The predominant frequency obtained from HVSR technique complements well with the frequency at which the peak response spectra ratio is observed. The present study of the predominant frequency identifies quite well the characteristics of Talchir Basin and is in good agreement with the synthetic ground motion modeling of the region.     

Influence of enzymatic lime on clay mineral behavior

Atterberg Limit Tests, unconfined compressive strength tests, and California Bearing Ratio Tests were conducted on the two minerals—kaolinite and bentonite—after subjecting them to stabilization with lime, enzyme, and enzymatic lime. The purpose of the study was to evaluate the effectiveness of enzymatic lime as a stabilizer over lime and enzyme. Results establish that enzymatic lime stabilization is the better stabilizer in terms of strength development. The unconfined compressive strength and California Bearing Ratio values of enzymatic lime stabilized samples have increased up to 5 and 30 times in the case of kaolinite mineral, but to only less than 1.5 times in the case of bentonite mineral. The study thus indicates that enzymatic lime stabilization is predominantly effective for clays containing kaolinite minerals.     

Exsolution textures in orthopyroxene in aluminous granulites as indicators of UHT metamorphism: New evidence from the Eastern Ghats Belt, India

Highly aluminous orthopyroxene, coexisting with sapphirine, cordierite, sillimanite, quartz and garnet in various combinations, constitute granoblastic mosaic peak metamorphic assemblages in aluminous granulites from three localities in the Eastern Ghats Belt, India. Orthopyroxene contains four types of intergrowths: (a) involving sapphirine with or without cordierite, (b) involving spinel, but without sapphirine, (c) involving cordierite, but without sapphirine and spinel, and (d) involving garnet, without sapphirine, spinel or cordierite. On the basis of textural and compositional data, origin of the intergrowths is ascribed to breakdown of Mg-Tschermak component, locally also involving Fe- and Ti-Tschermak. An attempt is made to compute the “pre-breakdown” compositions of orthopyroxene by image analysis, which shows maximum Al₂O₃ content of 13.4 wt.% in the pristine orthopyroxene. Geothermometry, phase equilibria consideration and application of existing experimental data on alumina solubility in orthopyroxene coexisting with sapphirine and quartz, collectively indicate extreme thermal conditions of metamorphism (> 1000 °C) for the studied assemblages. This re-affirms the notion that Al₂O₃ solubility in orthopyroxene is the most powerful indicator of UHT metamorphism (Harley, S.L., 2004. Extending our understanding of ultrahigh temperature crustal metamorphism. J. Mineral. Petrol. Sci. 99, 140–158). The intergrowths are considered to have formed due to cooling from the thermal peak spanning a temperature range of approximately 150 °C. Appearance of diverse types of intergrowths is probably related to subtle differences in bulk composition, particularly Fe:Mg ratios.     

Evolution of Mayurbhanj Granite Pluton, eastern Singhbhum, India: a case study of petrogenesis of an A-type granite in bimodal association

The A-type Mayurbhanj Granite Pluton (3.09 Ga), occurring along the eastern margin of the Singhbhum-Orissa Craton, eastern India, represents the final phase of acid plutonism in this crustal block of Archean age. The granite shows a bimodal association with a voluminous gabbroid body, exposed mainly along its western margin, and is associated with the Singhbhum Shear zone. The granite pluton is composed mainly of a coarse ferrohastingsite–biotite granite phase, with an early fine-grained granophyric microgranitic phase and a late biotite aplogranitic phase. Petrogenetic models of partial melting, fractional crystallisation and magma mixing have been advocated for the evolution of this pluton. New data, combined with earlier information, suggest that two igneous processes were responsible for the evolution of the Mayurbhanj Granite Pluton: partial melting of the Singhbhum Granite; followed by limited amount of mixing of acid and basic magmas in an anorogenic extensional setting. The necessary heat for partial melting was provided by the voluminous basaltic magma, now represented by the gabbroid body, emplaced at a shallow crustal level and showing a bimodal association with the Mayurbhanj Granite Pluton. The Singhbhum Shear Zone provided a possible channel way for the emplacement of the basic magma during crustal extension. It is concluded that all three phases of the Mayurbhanj Granite Pluton were derived from the same parent magma, generated by batch partial melting of the Singhbhum Granite at relatively high temperatures (980 °C) and low pressures (4 to <2 kbar) under anhydrous conditions. The coarse ferrohastingsite biotite granite phase shows evidence of limited and heterogeneous assimilation of country rock metasediments. However, the early microgranite phase and late aplogranite phase have not assimilated any metasediments. Compositional irregularities observed along the western margin of the Mayurbhanj Granite Pluton in contact with the gabbro body including a continuous fractionating sequence from quartz diorite to alkali-feldspar granite in the Notopahar area. Gradational contacts between the gabbro and the Mayurbhanj Granite Pluton in the Gorumahisani area etc., may be attributed to a limited amount of mixing between the gabbroid magma and the newly generated Mayurbhanj Granite magma. The mixing was mainly of liquid–liquid diffusive type, with a subordinate amount of mixing of solid–liquid type. Although A-type granites are commonly described as having high total REE (e.g. 270–400 ppm), studies on the late aplogranite phase of the Mayurbhanj Granite show that total REE values (100 ppm) are low. This low REE abundance may be attributed to the progressive residual nature of the Singhbhum Granite source during continued partial melting, when the magmas of the microgranite and coarse granite phases had already been removed from the source region.     

Transport and Biological Fate of Toluene in Low-Permeability Soils

The effect of simultaneous sorption, diffusion, and biodegradation on the fate and transport of toluene in low-permeability soil formations was examined. A transport model accounting for vapor and liquid sorption, vapor diffusions, and first-order biodegradation was developed to describe the movement of volatile solute in unsaturated soils. Modeling studies were followed with laboratory batch and column studies on fine-grained soil samples obtained from a gasoline-contaminated site. Batch experiments yielded the sorption and diffusion coefficients for generating theoretical solute transport profiles. Column studies were conducted to examine toluene sorption, diffusion, and biodegradation under aerobic and denitrifying conditions. Results from the column studies indicated that vapor sorption onto the soil was minimal due to the high moisture content of the soil. Comparison of model predictions with experimental results indicated that the SASK model, which is based on the resistivity theory, provided a more accurate prediction of the vapor phase tortuosity than the frequently used Millington-Quirk equation. Laboratory results of toluene concentration profiles matched well with the model predictions and yielded degradation rates comparable to those obtained in the field. Column studies, examining toluene biodegradation under aerobic and denitrifying conditions in low-permeability soils, indicated that the presence of excess nitrate in aerobic environments yielded higher solute degradation rates than those observed under exclusively aerobic systems.     

Earthquake hazard in Northeast India — A seismic microzonation approach with typical case studies from Sikkim Himalaya and Guwahati city

A comprehensive analytical as well as numerical treatment of seismological, geological, geomorphological and geotechnical concepts has been implemented through microzonation projects in the northeast Indian provinces of Sikkim Himalaya and Guwahati city, representing cases of contrasting geological backgrounds — a hilly terrain and a predominantly alluvial basin respectively. The estimated maximum earthquakes in the underlying seismic source zones, demarcated in the broad northeast Indian region, implicates scenario earthquakes of M _W 8.3 and 8.7 to the respective study regions for deterministic seismic hazard assessments. The microzonation approach as undertaken in the present analyses involves multi-criteria seismic hazard evaluation through thematic integration of contributing factors. The geomorphological themes for Sikkim Himalaya include surface geology, soil cover, slope, rock outcrop and landslide integrated to achieve geological hazard distribution. Seismological themes, namely surface consistent peak ground acceleration and predominant frequency were, thereafter, overlaid on and added with the geological hazard distribution to obtain the seismic hazard microzonation map of the Sikkim Himalaya. On the other hand, the microzonation study of Guwahati city accounts for eight themes — geological and geomorphological, basement or bedrock, landuse, landslide, factor of safety for soil stability, shear wave velocity, predominant frequency, and surface consistent peak ground acceleration. The five broad qualitative hazard classifications — ‘low’, ‘moderate’, ‘high’, ‘moderate high’ and ‘very high’ could be applied in both the cases, albeit with different implications to peak ground acceleration variations. These developed hazard maps offer better representation of the local specific seismic hazard variation in the terrain.     

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.