Estimation of peak ground acceleration and spectral acceleration for South India with local site effects: probabilistic approach

In this work an attempt has been made to evaluate the seismic hazard of South India (8.0 degrees N-20 degrees N; 72 degrees E-88 degrees E) based on the probabilistic seismic hazard analysis (PSHA). The earthquake data obtained from different sources were declustered to remove the dependent events. A total of 598 earthquakes of moment magnitude 4 and above were obtained from the study area after declustering, and were considered for further hazard analysis. The seismotectonic map of the study area was prepared by considering the faults, lineaments and the shear zones in the study area which are associated with earthquakes of magnitude 4 and above. For assessing theseismic hazard, the study area was divided into small grids of size 0.1 degrees x0.1 degrees, and the hazard parameters were calculated at the centre of each of these grid cells by considering all the seismic sources with in a radius of 300 km. Rock level peak horizontal acceleration (PHA) and spectral acceleration (SA) values at 1 corresponding to 10% and 2% probability of exceedance in 50 years have been calculated for all the grid points. The contour maps showing the spatial variation of these values are presented here. Uniform hazard response spectrum (UHRS) at rock level for 5% damping and 10% and 2% probability of exceedance in 50 years were also developed for all the grid points. The peak ground acceleration (PGA) at surface level was calculated for the entire South India for four different site classes. These values can be used to find the PGA values at any site in South India based on site class at that location. Thus, this method can be viewed as a simplified method to evaluate the PGA values at any site in the study area.

Testing and Modeling of MR Damper and Its Application to SDOF Systems Using Integral Backstepping Technique

Magnetorheological dampers are intrinsically nonlinear devices, which make the modeling and design of a suitable control algorithm an interesting and challenging task. To evaluate the potential of magnetorheological (MR) dampers in control applications and to take full advantages of its unique features, a mathematical model to accurately reproduce its dynamic behavior has to be developed and then a proper control strategy has to be taken that is implementable and can fully utilize their capabilities as a semi-active control device. The present paper focuses on both the aspects. First, the paper reports the testing of a magnetorheological damper with an universal testing machine, for a set of frequency, amplitude, and current. A modified Bouc-Wen model considering the amplitude and input current dependence of the damper parameters has been proposed. It has been shown that the damper response can be satisfactorily predicted with this model. Second, a backstepping based nonlinear current monitoring of magnetorheological dampers for semi-active control of structures under earthquakes has been developed. It provides a stable nonlinear magnetorheological damper current monitoring directly based on system feedback such that current change in magnetorheological damper is gradual. Unlike other MR damper control techniques available in literature, the main advantage of the proposed technique lies in its current input prediction directly based on system feedback and smooth update of input current. Furthermore, while developing the proposed semi-active algorithm, the dynamics of the supplied and commanded current to the damper has been considered. The efficiency of the proposed technique has been shown taking a base isolated three story building under a set of seismic excitation. Comparison with widely used clipped-optimal strategy has also been shown.

P-T evolution of Glenelg eclogites, NW Scotland: Did they experience ultrahigh-pressure metamorphism?

Eclogites and their retrogressed equivalents from the eastern unit of the Glenelg-Attadale Inlier in NW Scotland preserve much microstructural evidence that indicates that very high-pressure/temperature eclogite facies conditions were reached, and followed by decompression and hydration during exhumation. Rutile exsolution in garnet and quartz exsolution in omphacite and titanite formed through mineral reactions during high P-T peak metamorphism. Isochemical phase diagrams modeled for samples from three different locations indicate that the outer part of the eastern unit preserves a peak metamorphic condition of c. 850-1000 degrees C at 18-25 kbar, whereas the central part has a similar pressure (c. 23 kbar), but a lower temperature (c. 670 degrees C). Due to the limitations in the phase diagram calculations the estimated P-T conditions represent the minimum conditions attained by the peak metamorphic assemblage, and the pre-exsoived peak assemblage probably stabilized at a higher pressure. This observation is strongly supported by the presence of exsolution microstructures. The present results demonstrate that the eastern unit experienced very high P-T conditions during peak metamorphism and a tight clockwise P-T trajectory and provide the first indication of possible ultrahigh-pressure metamorphism in the Glenelg eclogites. (C) 2009 Elsevier B.V. All rights reserved.

Mapping of Average Shear Wave Velocity for Bangalore Region: A Case Study

Mapping the shear wave velocity profile is an important part in seismic hazard and microzonation studies. The shear wave velocity of soil in the city of Bangalore was mapped using the Multichannel Analysis of Surface Wave (MASW) technique. An empirical relationship was found between the Standard Penetration Test (SPT) corrected N value ((N1)60cs) and measured shear wave velocity (Vs). The survey points were selected in such a way that the results represent the entire Bangalore region, covering an area of 220 km2. Fifty-eight 1-D and 20 2-D MASW surveys were performed and their velocity profiles determined. The average shear wave velocity of Bangalore soils was evaluated for depths of 5 m, 10 m, 15 m, 20 m, 25 m and 30 m. The sub-soil classification was made for seismic local site effect evaluation based on average shear wave velocity of 30-m depth (Vs30) of sites using the National Earthquake Hazards Reduction Program (NEHRP) and International Building Code (IBC) classification. Mapping clearly indicates that the depth of soil obtained from MASW closely matches with the soil layers identified in SPT bore holes. Estimation of local site effects for an earthquake requires knowledge of the dynamic properties of soil, which is usually expressed in terms of shear wave velocity. Hence, to make use of abundant SPT data available on many geotechnical projects in Bangalore, an attempt was made to develop a relationship between Vs (m/s) and (N1)60cs. The measured shear wave velocity at 38 locations close to SPT boreholes was used to generate the correlation between the corrected N values and shear wave velocity. A power fit model correlation was developed with a regression coefficient (R2) of 0.84. This relationship between shear wave velocity and corrected SPT N values correlates well with the Japan Road Association equations.

Reliability assessment of internal stability of reinforced soil structures: A pseudo-dynamic approach

A methodology for reliability based optimum design of reinforced soil structures subjected to horizontal and vertical sinusoidal excitation based on pseudo-dynamic approach is presented. The tensile strength of reinforcement required to maintain the stability is computed using logarithmic spiral failure mechanism. The backfill soil properties, geometric and strength properties of reinforcement are treated as random variables. Effects of parameters like soil friction angle, horizontal and vertical seismic accelerations, shear and primary wave velocities, amplification factors for seismic acceleration on the component and system probability of failures in relation to tension and pullout capacities of reinforcement have been discussed. In order to evaluate the validity of the present formulation, static and seismic reinforcement force coefficients computed by the present method are compared with those given by other authors. The importance of the shear wave velocity in the estimation of the reliability of the structure is highlighted. The Ditlevsen's bounds of system probability of failure are also computed by taking into account the correlations between three failure modes, which is evaluated using the direction cosines of the tangent planes at the most probable points of failure. (c) 2009 Elsevier Ltd. All rights reserved.

Reliability Analysis of Earth Dams

In the present study, results of reliability analyses of four selected rehabilitated earth dam sections, i.e., Chang, Tapar, Rudramata, and Kaswati, under pseudostatic loading conditions, are presented. Using the response surface methodology, in combination with first order reliability method and numerical analysis, the reliability index (beta) values are obtained and results are interpreted in conjunction with conventional factor of safety values. The influence of considering variability in the input soil shear strength parameters, horizontal seismic coefficient (alpha(h)), and location of reservoir full level on the stability assessment of the earth dam sections is discussed in the probabilistic framework. A comparison of results with those obtained from other method of reliability analysis, viz., Monte Carlo simulations combined with limit equilibrium approach, provided a basis for discussing the stability of earth dams in probabilistic terms, and the results of the analysis suggest that the considered earth dam sections are reliable and are expected to perform satisfactorily.

Relationship between Low Strain Shear Modulus and Standard Penetration Test N Values

A low strain shear modulus plays a fundamental role in the estimation of site response parameters In this study an attempt has been made to develop the relationships between standard penetration test (SPT) N values with the low strain shear modulus (G(max)) For this purpose, field experiments SPT and multichannel analysis of surface wave data from 38 locations in Bangalore, India, have been used, which were also used for seismic microzonation project The in situ density of soil layer was evaluated using undisturbed soil samples from the boreholes Shear wave velocity (V-s) profiles with depth were obtained for the same locations or close to the boreholes The values for low strain shear modulus have been calculated using measured V-s and soil density About 215 pairs of SPT N and G(max) values are used for regression analysis The differences between fitted regression relations using measured and corrected values were analyzed It is found that an uncorrected value of N and modulus gives the best fit with a high regression coefficient when compared to corrected N and corrected modulus values This study shows better correlation between measured values of N and G(max) when compared to overburden stress corrected values of N and G(max)

Geochemistry of fluoride rich groundwater in Kolar and Tumkur Districts of Karnataka

Groundwater is a significant water resource in India for domestic, irrigation, and industrial needs. By far the most serious natural groundwater-quality problem in India, in terms of public health, derives from high fluoride, arsenic, and iron concentrations. Hydrogeochemical investigation of fluoride contaminated groundwater samples from Kolar and Tumkur Districts in Karnataka are undertaken to understand the quality and potability of groundwater from the study area, the level of fluoride contamination, the origin and geochemical mechanisms driving the fluoride enrichment. Majority of the groundwater samples did not meet the potable water criteria as they contained excess (>1.5 mg/L) fluoride, dissolved salts (>500 mg/L) and total hardness (75-924 mg/L). Hydrogeochemical facies of the groundwater samples suggest that rock weathering and evaporation-crystallization control the groundwater composition in the study area with 50-67% of samples belonging to the Ca-HCO3 type and the remaining falling into the mixed Ca-Na-HCO3 or Ca-Mg-Cl type. The saturation index values indicated that the groundwater in the study area is oversaturated with respect to calcite and under-saturated with respect to fluorite. The deficiency of calcium ion concentration in the groundwater from calcite precipitation favors fluorite dissolution leading to excess fluoride concentration.

Multiple support excitations of open-plane frames by a filtered white noise and soil-structure interaction

Seismic structural design is essentially the estimation of structural response to a forced motion, which may be deterministic or stochastic, imposed on the ground. The assumption that the same ground motion acts at every point of the base of the structure (or at every support) is not always justifiable; particularly in case of very large structures when considerable spatial variability in ground motion can exist over significant distances example long span bridges. This variability is partly due to the delay in arrival of the excitation at different supports (which is called the wave passage effect) and due to heterogeneity in the ground medium which results in incoherency and local effects. The current study examines the influence of the wave passage effect (in terms of delay in arrival of horizontal ground excitation at different supports and neglecting transmission through the structure) on the response of a few open-plane frame building structures with soil-structure interaction. The ground acceleration has been modeled by a suitably filtered white noise. As a special case, the ground excitation at different supports has also been treated as statistically independent to model the extreme case of incoherence due to local effects and due to modifications to the ground motion resulting from wave reflections and refractions in heterogeneous soil media. The results indicate that, even for relatively short spanned building frames, wave passage effect can be significant. In the absence of soil-structure interaction, it can significantly increase the root mean square (rms) value of the shear in extreme end columns for the stiffer frames but has negligible effect on the flexible frames when total displacements are considered. It is seen that pseudo-static displacements increasingly contribute to the rms value of column shear as the time delay increases both for the stiffer and for the more flexible frames. When soil-structure interaction is considered, wave passage effect (in terms of total displacements) is significant only for low soil shear modulus, G. values (where soil-structure interaction significantly lowers the fundamental frequency) and for stiff frames. The contribution of pseudo-static displacement to these rms values is found to decrease with increase in G. In general, wave passage effect for most interactive frames is insignificant compared to the attenuating effect a decrease in G, has on the response of the interactive structure to uniform support excitation. When the excitations at different supports are statistically independent, it is seen that for both the stiff and flexible frames, the rms value of the column shear in extreme end columns is several times larger (more for the stiffer frames) than the value corresponding to uniform base excitation with the pseudo-static displacements contributing over 99% of the rms value of column shear. Soil-structure interaction has an attenuating effect on the rms value of the column shear, the effect decreasing with increase in G,. Here too, the pseudo-static displacements contribute very largely to the column shear. The influence of the wave passage effect on the response of three 2-bay frames with and without soil-structure interaction to a recorded horizontal accelerogram is also examined. (C) 2010 Elsevier Ltd. All rights reserved.

Total risk rating and stability analysis of embankment dams in the Kachchh Region, Gujarat, India

A method for total risk analysis of embankment dams under earthquake conditions is discussed and applied to the selected embankment dams, i.e., Chang, Tapar, Rudramata, and Kaswati located in the Kachchh region of Gujarat, India, to obtain the seismic hazard rating of the dam site and the risk rating of the structures. Based on the results of the total risk analysis of the dams, coupled non-linear dynamic numerical analyses of the dam sections are performed using acceleration time history record of the Bhuj (India) earthquake as well as five other major earthquakes recorded worldwide. The objective of doing so is to perform the numerical analysis of the dams for the range of amplitude, frequency content and time duration of input motions. The deformations calculated from the numerical analyses are also compared with other approaches available in literature, viz, Makdisi and Seed (1978) approach, Jansen's approach (1990), Swaisgood's method (1995), Bureau's method (1997). Singh et al. approach (2007), and Saygili and Rathje approach (2008) and the results are utilized to foresee the stability of dams in future earthquake scenario. (C) 2010 Elsevier B.V. All rights reserved.

Magnetic fabrics indicating Late Quaternary seismicity in the Himalayan foothills

The anisotropy of magnetic susceptibility (AMS) study was performed on soft sediment samples from a trenched fault zone across the Himalayan frontal thrust (HFT), western Himalaya. AMS orientation of K-min axes in the trench sediments is consistent with lateral shortening revealed by geometry of deformed regional structures and recent earthquakes. Well-defined vertical magnetic foliation parallel to the flexure cleavage in which a vertical magnetic lineation is developed, high anisotropy, and triaxial ellipsoids suggest large overprinting of earth-quake- related fabrics. The AMS data suggest a gradual variation from layer parallel shortening (LPS) at a distance from the fault trace to a simple shear fabric close to the fault trace. An abrupt change in the shortening direction (K-min) from NE-SW to E-W suggests a juxtaposition of pre-existing layer parallel shortening fabric, and bending-related flexure associated with an earthquake. Hence the orientation pattern of magnetic susceptibility axes helps in identifying co-seismic structures in Late Holocene surface sediments.

Channel Volatiles Of South Indian Cordierites As Indicators Of Metamorphic Fluid Composition

The channel volatiles in cordierites of the Precambrian high-grade metapelites from southern and eastern Karnataka northern Tamil Nadu and southern Kerala were analyzed in an attempt to use them as metamorphic fluid fugacity indicators. Infrared powder absorption spectra, used to characterize the channel volatiles, showed that all the 21 analyzed cordierites have H2O and CO2 as the channel volatiles, indicating the predominantly H2O-CO2 composition of the metamorphic fluids. The H2O fraction in the metamorphic fluid was computed using a published thermodynamic method in conjunction with gravimetrically determined cordierite channel H2O content, available P - T estimates and an appropriate equation of state for the H2O - CO2 fluids. The IR data and these calculated X(H2O) values indicate an overall correlation between the variation in the relative proportion of H2O and CO2 in the fluids and the metamorphic grade. The average computed X(H2O) values are: 0.78 for the amphibolite facies eastern Karnataka pelites, 0.36 for the amphibolite facies southern Karnataka pelites, 0.19 for the southern Karnataka transitional zone rocks and 0.13 for the northern Tamil Nadu granulites. Consistently low X(H2O) values, at about 0.2, were obtained for the orthopyroxene-bearing assemblages.

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.