Intraplate Response to the Great 2004 Sumatra-Andaman Earthquake: A Study from the Andaman Segment

The 1300-km rupture of the 2004 interplate earthquake terminated at around 15 degrees N, in the northernmost segment of the Andaman-Nicobar subduction zone. This part of the plate boundary is noted for its generally lower level seismicity, compared with the southern segments. Based on the Global Centroid Moment Tensor (CMT) and National Earthquake Information Center (NEIC) data, most of the earthquakes of M-w >= 4.5 prior to 2004 were associated with the Andaman Spreading Ridge (ASR), and a few events were located within the forearc basin. The 2004 event was followed by an upward migration of hypocenters along the subducting plate, and the Andaman segment experienced a surge of aftershock activity. The continuing extensional faulting events, including the most recent earthquake (10 August 2009; M-w 7.5) in the northern end of the 2004 rupture, suggest the reduction of compressional strain associated with the interplate event. The style of faulting of the intraplate events before and after a great plate boundary earthquake reflects the relative influences of the plate-driving forces. Here we discuss the pattern of earthquakes in the Andaman segment before and after the 2004 event to appraise the spatial and temporal relation between large interplate thrust events and intraplate deformation. This study suggests that faulting mechanisms in the outer-ridge and outer-rise regions could be indicative of the maturity of interplate seismic cycles.

Critical earthquake input power spectral density function models for engineering structures

The problem of determining optimal power spectral density models for earthquake excitation which satisfy constraints on total average power, zero crossing rate and which produce the highest response variance in a given linear system is considered. The solution to this problem is obtained using linear programming methods. The resulting solutions are shown to display a highly deterministic structure and, therefore, fail to capture the stochastic nature of the input. A modification to the definition of critical excitation is proposed which takes into account the entropy rate as a measure of uncertainty in the earthquake loads. The resulting problem is solved using calculus of variations and also within linear programming framework. Illustrative examples on specifying seismic inputs for a nuclear power plant and a tall earth dam are considered and the resulting solutions are shown to be realistic.

Metamorphic-metasomatic fluids and Al, Si order/disorder of K-rich alkali feldspars from southern Indian high grade terrain

The structural state of K-feldspars in the quartzofeldspathic gneisses, charnockites, metapelites and pegmatites from the southern Kamataka, northern Tamil Nadu and southern Kerala high-grade regions of southern India has been characterized using petrographic and powder X-ray diffraction methods. The observed distribution pattern of structural state with a preponderance of disordered K-feldspar polymorphs in granulites compared to the ordered microclines in the amphibolite facies rocks is interpreted to reflect principally the varying H2O contents in the metamorphic-metasomatic fluids across metamorphic grade. The K-feldspars in the pegmatites of granitic derivation and in a pegmatite of inferred metamorphic origin also point to the important role of aqueous fluids in their structural state.

Revisiting the earthquake sources in the Himalaya: Perspectives on past seismicity

The similar to 2500 km-long Himalaya plate boundary experienced three great earthquakes during the past century, but none of them generated any surface rupture. The segments between the 1905-1934 and the 1897-1950 sources, known as the central and Assam seismic gaps respectively, have long been considered holding potential for future great earthquakes. This paper addresses two issues concerning earthquakes along the Himalaya plate boundary. One, the absence of surface rupture associated with the great earthquakes, vis-a-vis the purported large slip observed from paleoseismological investigations and two, the current understanding of the status of the seismic gaps in the Central Himalaya and Assam, in view of the paleoseismological and historical data being gathered. We suggest that the ruptures of earthquakes nucleating on the basal detachment are likely to be restricted by the crustal ramps and thus generate no surface ruptures, whereas those originating on the faults within the wedges promote upward propagation of rupture and displacement, as observed during the 2005 Kashmir earthquake, that showed a peak offset of 7 m. The occasional reactivation of these thrust systems within the duplex zone may also be responsible for the observed temporal and spatial clustering of earthquakes in the Himalaya. Observations presented in this paper suggest that the last major earthquake in the Central Himalaya occurred during AD 1119-1292, rather than in 1505, as suggested in some previous studies and thus the gap in the plate boundary events is real. As for the Northwestern Himalaya, seismically generated sedimentary features identified in the 1950 source region are generally younger than AD 1400 and evidence for older events is sketchy. The 1897 Shillong earthquake is not a decollement event and its predecessor is probably similar to 1000 years old. Compared to the Central Himalaya, the Assam Gap is a corridor of low seismicity between two tectonically independent seismogenic source zones that cannot be considered as a seismic gap in the conventional sense. (C) 2011 Elsevier B.V. All rights reserved.

Profiling of late Tertiary-early Quaternary surface in the lower reaches of Narmada valley using microtremors

In this paper, we propose the first approximation for thickness of Quaternary sediment and late Quaternary early Tertiary topography for the part of lower reaches of Narmada valley in a systematic way using the shallow seismic method, that records both horizontal and vertical components of the microtremor (ambient noise) caused by natural processes. The measurements of microtremors were carried out at 31 sites spaced at a grid interval of 5 km s using Lennartz seismometer (5 s period) and City shark-II data acquisition system. The signals recorded were analysed for horizontal to the vertical (H/V) spectral ratio using GEOPSY software. For the present study, we concentrate on frequency range between 0.2 Hz and 10 Hz. The thickness of unconsolidated sediments at various sites is calculated based on non-linear regression equations proposed by Ibs-von Seht and Wohlenberg (1999) and Parolai et al. (2002). The estimated thickness is used to plot digital elevation model and cross profiles correlating with geomorphology and geology of the study area. (C) 2011 Elsevier Ltd. All rights reserved.

Synthetic aperture radar: A focus on current problems

Synthetic aperture radar (SAR) is a powerful tool for mapping and remote sensing. The theory and operation of SAR have seen a period of intense activity in recent years. This paper attempts to review some of the more advanced topics studied in connection with modern SAR systems based on digital processing. Following a brief review of the principles involved in the operation of SAR, attention is focussed on special topics such as advanced SAR modelling and focussing techniques, in particular clutterlock and autofocus, Doppler centroid (DC) estimation methods involving seismic migration technique, moving target imaging, bistatic radar imaging, effects of system nonlinearities, etc.

Reliability Based Design Optimization for Internal Stability of Reinforced Soil Structures Using Pseudo‐Static Method

Seismic design of reinforced soil structures involves many uncertainties that arise from the backfill soil properties and tensile strength of the reinforcement which is not addressed in current design guidelines. This paper highlights the significance of variability in the internal stability assessment of reinforced soil structures. Reliability analysis is applied to estimate probability of failure and pseudo‐static approach has been used for the calculation of the tensile strength and length of the reinforcement needed to maintain the internal stability against tension and pullout failures. Logarithmic spiral failure surface has been considered in conjunction with the limit equilibrium method. Two modes of failure namely, tension failure and pullout failure have been considered. The influence of variations of the backfill soil friction angle, the tensile strength of reinforcement, horizontal seismic acceleration on the reliability index against tension failure and pullout failure of reinforced earth structure have been discussed.

Site Characterization of Bangalore using Geophysical Method

A number of geophysical methods have been proposed for near-surface site characterization and measurement of shear wave velocity by using a great variety of testing configurations, processing techniques,and inversion algorithms. In particular, two widely-used techniques are SASW (Spectral Analysis of SurfaceWaves) and MASW (Multichannel Analysis of SurfaceWaves). MASW is increasingly being applied to earthquake geotechnical engineering for the local site characterization, microzonation and site response studies.A MASW is a geophysical method, which generates a shear-wave velocity (Vs) profile (i.e., Vs versus depth)by analyzing Raleigh-type surface waves on a multichannel record. MASW system consisting of 24 channels Geode seismograph with 24 geophones of 4.5 Hz frequency have been used in this investigation. For the site characterization program, the MASW field experiments consisting of 58 one-dimensional shear wave velocity tests and 20 two-dimensional shear wave tests have been carried out. The survey points have been selected in such a way that the results supposedly represent the whole metropolitan Bangalore having an area of 220 km2.The average shear wave velocity of Bangalore soils have been evaluated for depths of 5m, 10m, 15m, 20m, 25m and 30 m. The subsoil site classification has been made for seismic local site effect evaluation based on average shear wave velocity of 30m depth (Vs30) of sites using National Earthquake Hazards Reduction Program (NEHRP) and International Building Code (IBC) classification. Soil average shearwave velocity estimated based on overburden thickness from the borehole information is also presented. Mapping clearly indicates that the depth of soil obtained from MASW is closely matching with the soil layers in bore logs. Among total 55 locations of MASW survey carried out, 34 locations were very close to the SPT borehole locations and these are used to generate correlation between Vs and corrected “N” values. The SPT field “N” values are corrected by applying the NEHRP recommended corrections.

Model track studies on fouled ballast using ground penetrating radar and multichannel analysis of surface wave

Ballast fouling is created by the breakdown of aggregates or outside contamination by coal dust from coal trains, or from soil intrusion beneath rail track. Due to ballast fouling, the conditions of rail track can be deteriorated considerably depending on the type of fouling material and the degree of fouling. So far there is no comprehensive guideline available to identify the critical degree of fouling for different types of fouling materials. This paper presents the identification of degree of fouling and types of fouling using non-destructive testing, namely seismic surface-wave and ground penetrating radar (GPR) survey. To understand this, a model rail track with different degree of fouling has been constructed in Civil engineering laboratory, University of Wollongong, Australia. Shear wave velocity obtained from seismic survey has been employed to identify the degree of fouling and types of fouling material. It is found that shear wave velocity of fouled ballast increases initially, reaches optimum fouling point (OFP), and decreases when the fouling increases. The degree of fouling corresponding after which the shear wave velocity of fouled ballast will be smaller than that of clean ballast is called the critical fouling point (CFP). Ground penetrating radar with four different ground coupled antennas (500 MHz, 800 MHz, 1.6 GHz and 2.3 GHz) was also used to identify the ballast fouling condition. It is found that the 800 MHz ground coupled antenna gives a better signal in assessing the ballast fouling condition. Seismic survey is relatively slow when compared to GPR survey however it gives quantifiable results. In contrast, GPR survey is faster and better in estimating the depth of fouling. (C) 2011 Elsevier B.V. All rights reserved.

Evaluation of spatial variation of peak horizontal acceleration and spectral acceleration for south India: a probabilistic approach

In this work, an attempt has been made to evaluate the spatial variation of peak horizontal acceleration (PHA) and spectral acceleration (SA) values at rock level for south India based on the probabilistic seismic hazard analysis (PSHA). These values were estimated by considering the uncertainties involved in magnitude, hypocentral distance and attenuation of seismic waves. Different models were used for the hazard evaluation, and they were combined together using a logic tree approach. For evaluating the seismic hazard, the study area was divided into small grids of size 0.1A degrees A xA 0.1A degrees, and the hazard parameters were calculated at the centre of each of these grid cells by considering all the seismic sources within a radius of 300 km. Rock level PHA values and SA at 1 s corresponding to 10% probability of exceedance in 50 years were evaluated for all the grid points. Maps showing the spatial variation of rock level PHA values and SA at 1 s for the entire south India are presented in this paper. To compare the seismic hazard for some of the important cities, the seismic hazard curves and the uniform hazard response spectrum (UHRS) at rock level with 10% probability of exceedance in 50 years are also presented in this work.

Frequency response analysis of layered ground in Ahmedabad during Bhuj earthquake

This paper presents the results of seismic response analysis of layered ground in Ahmedabad City during the earthquake in Bhuj on 26(th) January 2001. An attempt has been made to understand the reasons for the failure of multistoreyed buildings founded on soft alluvial deposits in Ahmedabad. Standard Penetration test at a site very close to the Sabarmati river belt was carried out for geotechnical investigations. The program SHAKE91, widely used in the field of earthquake engineering for computing the seismic response of horizontally layered soil deposits, was used to analyse the soil profile at the selected site considering the ground as one dimensional layered elastic system. The ground accelerations recorded at the ground floor of the Regional Passport Staff Quarters building, which is very close to the investigated site, was used as input motion. Also, Finite Element Analysis was carried out for different configurations of multistorey building frames for evaluating their natural frequencies and is compared with the predominant frequency of the layered soil system. The results reveal that the varying degree of damage to multistorey buildings in the close proximity of Sabarmati river area was essentially due to the large amplification of the ground and the near resonance condition.

Review of correlations between SPT N and shear modulus: A new correlation applicable to any region

A low strain shear modulus plays a fundamental role in earthquake geotechnical engineering to estimate the ground response parameters for seismic microzonation. A large number of site response studies are being carried out using the standard penetration test (SPT) data, considering the existing correlation between SPT N values and shear modulus. The purpose of this paper is to review the available empirical correlations between shear modulus and SPT N values and to generate a new correlation by combining the new data obtained by the author and the old available data. The review shows that only few authors have used measured density and shear wave velocity to estimate shear modulus, which were related to the SPT N values. Others have assumed a constant density for all the shear wave velocities to estimate the shear modulus. Many authors used the SPT N values of less than 1 and more than 100 to generate the correlation by extrapolation or assumption, but practically these N values have limited applications, as measuring of the SPT N values of less than 1 is not possible and more than 100 is not carried out. Most of the existing correlations were developed based on the studies carried out in Japan, where N values are measured with a hammer energy of 78%, which may not be directly applicable for other regions because of the variation in SPT hammer energy. A new correlation has been generated using the measured values in Japan and in India by eliminating the assumed and extrapolated data. This correlation has higher regression coefficient and lower standard error. Finally modification factors are suggested for other regions, where the hammer energy is different from 78%. Crown Copyright (C) 2012 Published by Elsevier Ltd. All rights reserved.

Recent microtremors near the Idukki Reservoir, Kerala, South India

The continuing low-level seismicity in the vicinity of the Idukki Reservoir, Kerala, is interesting from the perspective of hydrologically triggered earthquakes. While the frequency of triggered earthquakes in the vicinity of a reservoir usually reduces with time and the largest earthquake usually occurs within a few years on the initial filling, the triggered seismicity in the proximity of the Idukki Reservoir seems to be showing a second, delayed peak, as the 1977 (M 3.5) tremor was followed by a slightly larger event in 2011, 24 years after the first burst of activity. Quite unprecedented in the context of reservoir-triggered sequences, we consider this delayed sequence as the hydrologic response of a critically stressed hypocentral region, to monsoonal recharging. The sustained activity several decades after the impoundment and the temporal relation with the monsoon suggest that at least some parts of the reservoir region continue to retain the potential for low-level seismic activity in response to hydrologic cycles.

The hazard potential of the western segment of the Makran subduction zone, northern Arabian Sea

Evaluating the hazard potential of the Makran subduction zone requires understanding the previous records of the large earthquakes and tsunamis. We address this problem by searching for earthquake and tectonic proxies along the Makran Coast and linking those observations with the available constraints on historical seismicity and the tell-tale characteristics of sea floor morphology. The earthquake of Mw 8.1 of 1945 and the consequent tsunami that originated on the eastern part of the Makran are the only historically known hazardous events in this region. The seismic status of the western part of the subduction zone outside the rupture area of the 1945 earthquake remains an enigma. The near-shore shallow stratigraphy of the central part of Makran near Chabahar shows evidence of seismically induced liquefaction that we attribute to the distant effects of the 1945 earthquake. The coastal sites further westward around Jask are remarkable for the absence of liquefaction features, at least at the shallow level. Although a negative evidence, this possibly implies that the western part of Makran Coast region may not have been impacted by near-field large earthquakes in the recent past-a fact also supported by the analysis of historical data. On the other hand, the elevated marine terraces on the western Makran and their uplift rates are indicative of comparable degree of long-term tectonic activity, at least around Chabahar. The offshore data suggest occurrences of recently active submarine slumps on the eastern part of the Makran, reflective of shaking events, owing to the great 1945 earthquake. The ocean floor morphologic features on the western segment, on the contrary, are much subdued and the prograding delta lobes on the shelf edge also remain intact. The coast on the western Makran, in general, shows indications of progradation and uplift. The various lines of evidence thus suggest that although the western segment is potentially seismogenic, large earthquakes have not occurred there in the recent past, at least during the last 600 years. The recurrence period of earthquakes may range up to 1,000 years or more, an assessment based on the age of the youngest dated coastal ridge. The long elapsed time points to the fact that the western segment may have accumulated sufficient slip to produce a major earthquake.

On the Compatibility of a Flux Transport Dynamo with a Fast Tachocline Scenario

The compatibility of the fast-tachocline scenario with a flux-transport dynamo model is explored. We employ a flux-transport dynamo model coupled with simple feedback formulae relating the thickness of the tachocline to the amplitude of the magnetic field or to the Maxwell stress. The dynamo model is found to be robust against the nonlinearity introduced by this simplified fast-tachocline mechanism. Solar-like butterfly diagrams are found to persist and, even without any parameter fitting, the overall thickness of the tachocline is well within the range admitted by helioseismic constraints. In the most realistic case of a time-and latitude-dependent tachocline thickness linked to the value of the Maxwell stress, both the thickness and its latitudinal dependence are in excellent agreement with seismic results. In nonparametric models, cycle-related temporal variations in tachocline thickness are somewhat larger than admitted by helioseismic constraints; we find, however, that introducing a further parameter into our feedback formula readily allows further fine tuning of the thickness variations.