Critical response of structures to multicomponent earthquake excitation

This paper aims to develop an improved understanding of the critical response of structures to multicomponent seismic motion characterized by three uncorrelated components that are defined along its principal axes: two horizontal and the vertical component. An explicit formula, convenient for code applications, has been derived to calculate the critical value of structural response to the two principal horizontal components acting along any incident angle with respect to the structural axes, and the vertical component of ground motion. The critical response is defined as the largest value of response for all possible incident angles. The ratio r_cr/r_srss between the critical value of response and the SRSS response—corresponding to the principal components of ground acceleration applied along the structure axes—is shown to depend on three dimensionless parameters: the spectrum intensity ratio γ between the two principal components of horizontal ground motion characterized by design spectra A(T_n) and γA(T_n); the correlation coefficient α of responses r_x and r_y due to design spectrum A(T_n) applied in the x‐ and y‐directions, respectively; and β = r_y/r_x. It is demonstrated that the ratio r_cr/r_srss is bounded by 1 and . Thus the largest value of the ratio is , 1.26, 1.13 and 1.08 for γ = 0, 0.5, 0.75 and 0.85, respectively. This implies that the critical response never exceeds times the result of the SRSS analysis, and this ratio is about 1.13 for typical values of γ, say 0.75. The correlation coefficient α depends on the structural properties but is always bounded between −1 and 1. For a fixed value of γ, the ratio r_cr/r_srss is largest if β = 1 and α = ±1. The parametric variations presented for one‐storey buildings indicate that this condition can be satisfied by axial forces in columns of symmetric‐plan buildings or can be approximated by lateral displacements in resisting elements of unsymmetrical‐plan buildings. Copyright © 2000 John Wiley & Sons, Ltd.     

Application of MPA to estimate probability of collapse of structures

This paper develops a modal pushover analysis‐ (MPA) based approximate procedure to quantify the collapse potential of structural systems. The computationally demanding incremental dynamic analysis (IDA) of the structural system is avoided by MPA of the structure in conjunction with empirical equations for the collapse strength ratio for the first‐mode single‐degree‐of‐freedom (SDF) system; higher modes of vibration play essentially no role in estimating the ground motion intensity required to cause collapse of the structure. Presented are collapse fragility curves for 6‐, 9‐, and 20‐story regular special moment‐resisting teel frames computed by the exact and approximate procedures, demonstrating that the MPA‐based approximate procedure requires only a small fraction (1% in one example) of the computational effort inherent in exact IDA and still achieves highly accurate results. Copyright © 2010 John Wiley & Sons, Ltd.     

Modeling of strong ground motions for 1991 Uttarkashi, 1999 Chamoli earthquakes, and a hypothetical great earthquake in Garhwal–Kumaun Himalaya

In this study, stochastic finite fault modeling is used to simulate Uttarkashi (1991) and Chamoli (1999) earthquakes using all available source, path, and site parameters available for the region. These two moderate earthquakes are recorded at number of stations of a strong motion network. The predicted peak ground accelerations at these stations are compared with the observed data and the ground motion parameters are constrained. The stress drop of Uttarkashi and Chamoli earthquakes is constrained at 77 and 65?bars, respectively, whereas the quality factor Q _C is 112 $ f^{0.97} $ and 149 $ f^{0.95} $ for these two regions. The high-frequency attenuation parameter Kappa is in the range 0.04?C0.05. The constrained ground motion parameters are then used to simulate Mw 8.5 earthquake in central seismic gap region of Himalaya. Two scenarios are considered with epicenter of future great earthquake at locations of Uttarkashi and Chamoli earthquakes using above constrained parameters. The most vulnerable towns are the towns of Dehradun and Almora where expected PGA is in excess of 600?cm/s² at V_S30 520?m/s when the epicenter of the great earthquake is at the location of Uttarkashi (1991) earthquake. The towns of Shimla and Chandigarh can expect PGA close to 200?cm/s². Whereas when the epicenter of the great earthquake is at the location of Chamoli (1999) earthquake, the towns of Dehradun and Almora can expect PGA of around 500 and 400?cm/s², respectively, at V_S30 620?m/s. The National Capital Region, Delhi can expect accelerations of around 80?cm/s² in both the cases. The PGA contour maps obtained in this study can be used to assess the seismic hazard of the region and identify vulnerable areas in and around central Himalaya from a future great earthquake.     

Estimation of radioactive leakages into the Pacific Ocean due to Fukushima nuclear accident

High concentrations of several radionuclides were reported in the sea near the Fukushima Daiichi Nuclear Power Station (FDNPS) in Japan due to the nuclear accident that occurred on 11 March 2011. The main source of these concentrations was leakage of highly radioactive liquid effluent from a pit in the turbine building near the intake canal of Unit-2 of FDNPS through a crack in the concrete wall. In the immediate vicinity of the plant, seawater concentrations reached 68 MBq m^?3 for ¹³⁴Cs and ¹³⁷Cs, and exceeded 100 MBq m^?3 for ¹³¹I in early April 2011. These concentrations began to fall as of 11 April 2011 and, at the end of April, had reached a value close to 0.1 MBq m^?3 for ¹³⁷Cs. Following the nuclear accident, the Tokyo Electric Power Company (TEPCO) had initiated intense monitoring of the environment including the Pacific Ocean. Seawater samples were collected and the concentrations of few radionuclides were measured on a wide spatial and temporal scale. In this study, the measured concentrations of different radionuclides near the south discharge canal of the FDNPS were used to estimate their leakages into the Pacific Ocean. The method is based on estimating the release rates that reproduce the concentration of radionuclides in seawater at a chosen location using a two-dimensional advection–dispersion model in an iterative manner. The radioactive leakages were estimated as 5.68 PBq for ¹³¹I, 2.24 PBq for ¹³⁴Cs and 2.25 PBq for ¹³⁷Cs. Leakages were also estimated for ^99mTc, ¹³⁶Cs, ¹⁴⁰Ba and ¹⁴⁰La and they range between 0.02 PBq (^99mTc) and 0.53 PBq (¹⁴⁰Ba). It was estimated that about 11.28 PBq of radioactivity in total was leaked into the Pacific Ocean from the damaged FDNPS. Out of this, ¹³¹I constitutes 50.3 %; ¹³⁴Cs 20 %; ¹³⁷Cs 20 %; ¹⁴⁰Ba 4.6 %; ¹³⁶Cs 2.6 %; ¹⁴⁰La 2.3 % and ^99mTc 0.2 % of the total radioactive leakage. Such quantitative estimates of radioactive leakages are essential prerequisites for short-term and local-scale as well as long-term and large-scale radiological impact assessment of the nuclear accident.     

The 2007 Talala,Saurashtra, western India earthquake sequence: Tectonic implications and seismicity triggering

A damaging and widely felt moderate (M_w 5.0) earthquake occurred in the Talala region of Saurashtra, Gujarat (western India) on November 6, 2007. The highly productive sequence comprised about 1300 micro earthquakes (M > 0.5) out of which 325 of M ? 1.5 that occurred during November 6, 2007–January 10, 2008 were precisely located. The spatial aftershock distribution revealed a NE–SW striking fault in accordance with the centroid moment tensor solution, which in turn implies left-lateral motion. The orientation and sense of shear are consistent with similarly orientated geological fault identified in the area from satellite imagery and field investigation.The aftershocks temporal decay, b-value of frequency–magnitude distribution, spatial fractal dimension, D, and slip ratio (ratio of the slip occurred on the primary fault to the total slip) were examined with the purpose to identify the properties of the sequence. The high b-value (1.18 ± 0.01) may be attributed to the paucity of the larger (M ? 4.0) aftershocks and reveals crustal heterogeneity and low stress regime. The high p-value (1.10 ± 0.39), implying fast decay rate of aftershocks, evidences high surface heat flux. A value of the spatial fractal dimension (D) equal to 2.21 ± 0.02 indicates random spatial distribution and source in a two-dimensional plane that is being filled-up by fractures. A slip ratio of 0.42 reveals that more slip occurred on secondary fault systems.The static Coulomb stress changes due to the coseismic slip of the main shock, enhanced off fault aftershock occurrence. The occurrence of a moderate earthquake (M_w 4.3) on October 5, 2008 inside a region of positive Coulomb stress changes supports the postulation on aftershock triggering. When the stress changes were resolved on a cross section including the stronger (M4.8) foreshock plane that is positioned adjacent to the main fault, it became evident that the activity continued there due to stress transfer from the main rupture.     

Attenuation of High Frequency P and S Waves in the Gujarat Region,India

The local earthquake waveforms recorded on broadband seismograph network of Institute of Seismological Research in Gujarat, India have been analyzed to understand the attenuation of high frequency (2–25 Hz) P and S waves in the region. The frequency dependent relationships for quality factors for P (Q _P) and S (Q _S) waves have been obtained using the spectral ratio method for three regions namely, Kachchh, Saurashtra and Mainland Gujarat. The earthquakes recorded at nine stations of Kachchh, five stations of Saurashtra and one station in mainland Gujarat have been used for this analysis. The estimated relations for average Q _P and Q _S are: Q _P = (105 ± 2) f ^{0.82 ± 0.01}, Q _S = (74 ± 2) f ^{1.06 ± 0.01} for Kachchh region; Q _P = (148 ± 2) f ^{0.92 ± 0.01}, Q _S = (149 ± 14) f ^{1.43 ± 0.05} for Saurashtra region and Q _P = (163 ± 7) f ^{0.77 ± 0.03}, Q _S = (118 ± 34) f ^{0.65 ± 0.14} for mainland Gujarat region. The low Q (<200) and high exponent of f (>0.5) as obtained from present analysis indicate the predominant seismic activities in the region. The lowest Q values obtained for the Kachchh region implies that the area is relatively more attenuative and heterogeneous than other two regions. A comparison between Q _S estimated in this study and coda Q (Qc) previously reported by others for Kachchh region shows that Q _C > Q _S for the frequency range of interest showing the enrichment of coda waves and the importance of scattering attenuation to the attenuation of S waves in the Kachchh region infested with faults and fractures. The Q _S/Q _P ratio is found to be less than 1 for Kachchh and Mainland Gujarat regions and close to unity for Saurashtra region. This reflects the difference in the geological composition of rocks in the regions. The frequency dependent relations developed in this study could be used for the estimation of earthquake source parameters as well as for simulating the strong earthquake ground motions in the region.     

Evaluation of three‐dimensional modal pushover analysis for unsymmetric‐plan buildings subjected to two components of ground motion

The accuracy of the three‐dimensional modal pushover analysis (MPA) procedure in estimating seismic demands for unsymmetric‐plan buildings due to two horizontal components of ground motion, simultaneously, is evaluated. Eight low‐and medium‐rise structures were considered. Four intended to represent older buildings were designed according to the 1985 Uniform Building Code, whereas four other designs intended to represent newer buildings were based on the 2006 International Building Code. The median seismic demands for these buildings to 39 two‐component ground motions, scaled to two intensity levels, were computed by MPA and nonlinear response history analysis (RHA), and then compared. Even for these ground motions that deform the buildings significantly into the inelastic range, MPA offers sufficient degree of accuracy. It is demonstrated that PMPA, a variant of the MPA procedure, for nonlinear systems is almost as accurate as the well‐known standard response spectrum analysis procedure is for linear systems. Thus, for practical applications, the PMPA procedure offers an attractive alternative to nonlinear RHA, whereby seismic demands can be estimated directly from the (elastic) design spectrum. In contrast, the nonlinear static procedure specified in the ASCE/SEI 41‐06 Standard is demonstrated to grossly underestimate seismic demands for some of the unsymmetric‐plan buildings considered. Copyright © 2011 John Wiley & Sons, Ltd.     

Estimation of seismic hazard in Gujarat region, India

The seismic hazard in the Gujarat region has been evaluated. The scenario hazard maps showing the spatial distribution of various parameters like peak ground acceleration, characteristics site frequency and spectral acceleration for different periods have been presented. These parameters have been extracted from the simulated earthquake strong ground motions. The expected damage to buildings from future large earthquakes in Gujarat region has been estimated. It has been observed that the seismic hazard of Kachchh region is more in comparison with Saurashtra and mainland. All the cities of Kachchh can expect peak acceleration in excess of 500?cm/s² at surface in case of future large earthquakes from major faults in Kachchh region. The cities of Saurashtra can expect accelerations of less than 200?cm/s² at surface. The mainland Gujarat is having the lowest seismic hazard as compared with other two regions of Gujarat. The expected accelerations are less than 50?cm/s² at most of the places. The single- and double-story buildings in Kachchh region are at highest risk as they can expect large accelerations corresponding to natural periods of such small structures. Such structures are relatively safe in mainland region. The buildings of 3?C4 stories and tall structures that exist mostly in cities of Saurashtra and mainland can expect accelerations in excess of 100?cm/s² during a large earthquake in Kachchh region. It has been found that a total of 0.11 million buildings in Rajkot taluka of Saurashtra are vulnerable to total damage. In Kachchh region, 0.37 million buildings are vulnerable. Most vulnerable talukas are Bhuj, Anjar, Rapar, Bhachau, and Mandvi in Kachchh district and Rajkot, Junagadh, Jamnagar, Surendernagar and Porbandar in Saurashtra. In mainland region, buildings in Bharuch taluka are more vulnerable due to proximity to active Narmada-Son geo-fracture. The scenario hazard maps presented in this study for moderate as well as large earthquakes in the region may be used to augment the information available in the probabilistic seismic hazard maps of the region.