LESSON LEARNED FROM THE EFFECT OF RECENT FAR FIELD SUMATRA EARTHQUAKES TO PENINSULAR MALAYSIA

Similar documents
THE EFFECT OF THE LATEST SUMATRA EARTHQUAKE TO MALAYSIAN PENINSULAR

MICROZONATION MAPS FOR KUALA LUMPUR AND PUTRAJAYA

Usability of the Next Generation Attenuation Equations for Seismic Hazard Assessment in Malaysia

ATTENUATION FUNCTION RELATIONSHIP OF SUBDUCTION MECHANISM AND FAR FIELD EARTHQUAKE

Malaysian Journal of Civil Engineering 22(1) : (2010) Malaysia

SEISMIC HAZARD ANALYSIS. Instructional Material Complementing FEMA 451, Design Examples Seismic Hazard Analysis 5a - 1

Geo-Marine Letters Volume 36, 2016, electronic supplementary material

Part 2 - Engineering Characterization of Earthquakes and Seismic Hazard. Earthquake Environment

Investigation of long period amplifications in the Greater Bangkok basin by microtremor observations

Overview of Seismic PHSA Approaches with Emphasis on the Management of Uncertainties

Hazard and Vulnerability of Moderate Seismicity Regions

I. Locations of Earthquakes. Announcements. Earthquakes Ch. 5. video Northridge, California earthquake, lecture on Chapter 5 Earthquakes!

Modelling Strong Ground Motions for Subduction Events in the Wellington Region, New Zealand

VALIDATION AGAINST NGA EMPIRICAL MODEL OF SIMULATED MOTIONS FOR M7.8 RUPTURE OF SAN ANDREAS FAULT

STUDY ON THE BI-NORMALIZED EARTHQUAKE ACCELERATION RESPONSE SPECTRA

SITE CLASSIFICATION AND SEISMIC RESPONSE OF DHAKA CITY SOILS

Improvements to the Development of Acceleration Design Response Spectra. Nicholas E. Harman, M.S., P.E., SCDOT

A note on ground motion recorded during Mw 6.1 Mae Lao (Northern Thailand) earthquake on 5 May 2014

log 4 0.7m log m Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi. Module 1 Seismology Exercise Problems :

EVALUATION OF SEISMIC SITE EFFECTS FOR BANGKOK DEEP BASIN

Incorporating simulated Hikurangi subduction interface spectra into probabilistic hazard calculations for Wellington

Unique Site Conditions and Response Analysis Challenges in the Central and Eastern U.S.

Liquefaction Assessment using Site-Specific CSR

An Approach for Seismic Design in Malaysia following the Principles of Eurocode 8

SITE RESPONSE ANALYSIS FOR SEISMIC DESIGN OF A 48-STOREY TOWER BUILDING IN JAKARTA

GROUND RESPONSE ANALYSIS FOR SEISMIC DESIGN IN FRASER RIVER DELTA, BRITISH COLUMBIA

Micro Seismic Hazard Analysis

The Earthquake Cycle Chapter :: n/a

SHAKE TABLE STUDY OF SOIL STRUCTURE INTERACTION EFFECTS ON SEISMIC RESPONSE OF SINGLE AND ADJACENT BUILDINGS

Professor Terje Haukaas University of British Columbia, Vancouver terje.civil.ubc.ca. Earthquakes

Seismic Site Effects of Soil Amplifications in Bangkok

Occurrence of negative epsilon in seismic hazard analysis deaggregation, and its impact on target spectra computation

SEISMIC HAZARD ANALYSIS AND SEISMIC INPUT TO TSUNAMI MODELING FOR MICROZONATION OF MEURAXA DISTRICT CITY OF BANDA ACEH

Determination of the seismicity and peak ground acceleration for Lombok Island: An evaluation on tectonic setting

Mw 7.8, Southwest of Sumatra, Indonesia Wed, 2 March 2016 at 12:49:48 UTC M /03/03

Three Fs of earthquakes: forces, faults, and friction. Slow accumulation and rapid release of elastic energy.

Identifying the causes and effects of earthquakes

Modelling of Site Specific Response Spectrum for Buildings in Makassar

Development of U. S. National Seismic Hazard Maps and Implementation in the International Building Code

ACCOUNTING FOR SITE EFFECTS IN PROBABILISTIC SEISMIC HAZARD ANALYSIS: OVERVIEW OF THE SCEC PHASE III REPORT

Estimation of Peak Ground Acceleration for Delhi Region using Finsim, a Finite Fault Simulation Technique

SIMULATION OF NEAR FIELD STRONG GROUND MOTIONS AT TOMBAK SITE IN IRAN USING HYBRID METHOD

2C09 Design for seismic and climate changes

New Design Spectral Acceleration of Soft and Deep Deposits in Bangkok

COMPARISON BETWEEN 2D AND 3D ANALYSES OF SEISMIC STABILITY OF DETACHED BLOCKS IN AN ARCH DAM

FOCAL MECHANISMS OF SUBDUCTION ZONE EARTHQUAKES ALONG THE JAVA TRENCH: PRELIMINARY STUDY FOR THE PSHA FOR YOGYAKARTA REGION, INDONESIA

Effects of Surface Geology on Seismic Motion

Synthetic Earthquake Ground Motions for the Design of Long Structures

New Prediction Formula of Fourier Spectra Based on Separation Method of Source, Path, and Site Effects Applied to the Observed Data in Japan

A GEOTECHNICAL SEISMIC SITE RESPONSE EVALUATION PROCEDURE

Eleventh U.S. National Conference on Earthquake Engineering Integrating Science, Engineering & Policy June 25-29, 2018 Los Angeles, California

THE RESPONSE SPECTRUM

DIRECT HAZARD ANALYSIS OF INELASTIC RESPONSE SPECTRA

CHARACTERIZATION OF DIRECTIVITY EFFECTS OBSERVED DURING 1999 CHI-CHI, TAIWAN EARTHQUAKE

Updated Graizer-Kalkan GMPEs (GK13) Southwestern U.S. Ground Motion Characterization SSHAC Level 3 Workshop 2 Berkeley, CA October 23, 2013

Interpretive Map Series 24

Evaluating the Seismic Coefficient for Slope Stability Analyses

Ground-Motion Attenuation Relationships for Subduction- Zone Earthquakes in Northern Taiwan

Deterministic Seismic Hazard Assessment of Quetta, Pakistan

Seismic site response analysis in Perth Metropolitan area

Mechanics of Earthquakes and Faulting

Elastic rebound theory

Characteristics of a Force Loads on Structures. Dead Load. Load Types Dead Live Wind Snow Earthquake. Load Combinations ASD LRFD

GROUND MOTION SPECTRAL INTENSITY PREDICTION WITH STOCHASTIC GREEN S FUNCTION METHOD FOR HYPOTHETICAL GREAT EARTHQUAKES ALONG THE NANKAI TROUGH, JAPAN

Magnitude 6.9 GULF OF CALIFORNIA

NGA-Subduction: Development of the Largest Ground Motion Database for Subduction Events

Damping Scaling of Response Spectra for Shallow CCCCCCCCCrustalstallPaper Crustal Earthquakes in Active Tectonic Title Line Regions 1 e 2

I. What are Earthquakes?

STRONG GROUND MOTION ATTENUATION IN THE SEA OF JAPAN (OKHOTSK-AMUR PLATES BOUNDARY) REGION

Seismic Response of Sedimentary Basin Subjected to Obliquely Incident SH Waves

Development of acceleration time histories for Semarang, Indonesia, due to shallow crustal fault earthquakes

Estimation of Shear Wave Velocity Using Correlations

5. What is an earthquake 6. Indicate the approximate radius of the earth, inner core, and outer core.

Seismic Activity near the Sunda and Andaman Trenches in the Sumatra Subduction Zone

Once you have opened the website with the link provided choose a force: Earthquakes

RECIPE FOR PREDICTING STRONG GROUND MOTIONS FROM FUTURE LARGE INTRASLAB EARTHQUAKES

Centrifuge Evaluation of the Impact of Partial Saturation on the Amplification of Peak Ground Acceleration in Soil Layers

SURFACE WAVES AND SEISMIC RESPONSE OF LONG-PERIOD STRUCTURES

Seismic Hazard of Singapore and Malaysia

Important Concepts. Earthquake hazards can be categorized as:

Earthquakes How and Where Earthquakes Occur

Draft National Annex to Eurocode 8 for Malaysia and cost implication for residential buildings with thin size elements

SCALING LAWS FOR FOURIER ACCELERATION SPECTRA IN FORMER YUGOSLAVIA

PEAK GROUND MOTION CHARACTERISTICS OF 1995 KOBE EARTHQUAKE AND AN EXTRACTED SIMPLE EVALUATION METHOD

An entire branch of Earth science, called, is devoted to the study of earthquakes.

Deterministic Generation of Broadband Ground Motions! with Simulations of Dynamic Ruptures on Rough Faults! for Physics-Based Seismic Hazard Analysis

THE NATURE OF SITE RESPONSE DURING EARTHQUAKES. Mihailo D. Trifunac

GROUND MOTION TIME HISTORIES FOR THE VAN NUYS BUILDING

Simulation of Strong Ground Motions for a Shallow Crustal Earthquake in Japan Based on the Pseudo Point-Source Model

QUAKE/W ProShake Comparison

NORTH OF TEHRAN SITE EFFECT MICROZONATION

Seismic Analysis of an Instrumented Earth Dam in Terms of the Variability of the Response Design Spectra

Modelling Subduction Zone Seismogenic Hazards in Southeast Asia for Seismic Hazard Assessments

Seismic Hazard & Risk Assessment

A Visco-Elastic Model with Loading History Dependent Modulus and Damping for Seismic Response Analyses of Soils. Zhiliang Wang 1 and Fenggang Ma 2.

Magnitude 7.2 OAXACA, MEXICO

STUDY ON INFLUENCE OF LOCAL SOIL CONDITIONS ON GROUND MOTION AMPLIFICATION

Magnitude 7.0 NEW CALEDONIA

Dynamic modelling in slopes using finite difference program

Modifications to Risk-Targeted Seismic Design Maps for Subduction and Near-Fault Hazards

Transcription:

13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 416 LESSON LEARNED FROM THE EFFECT OF RECENT FAR FIELD SUMATRA EARTHQUAKES TO PENINSULAR MALAYSIA Azlan ADNAN 1, Aminaton MARTO 2, Hendriyawan 3 SUMMARY This paper presents the results of seismic hazard study of two large earthquakes occurred near Sumatra at the end of 2002 (M w =7.4) and in the early of 2003 (M w = 5.8). The objective of this study is to evaluate the effect of those Sumatra Earthquakes to Peninsular Malaysia. This paper discusses the following items: 1) tectonic setting of Sumatra, 2) location, mechanism, and size of recent earthquake, 3) analysis of ground acceleration at bedrock for Penang and Kuala Lumpur, and 4) analysis of local site effect for Penang and Kuala Lumpur. Seismic analysis was conducted in order to predict peak ground acceleration (PGA) and spectral acceleration at bedrock of Penang and Kuala Lumpur. Analysis was carried out using deterministic method and appropriate attenuation relationship. The selection of attenuation functions were based on the similarity of faulting mechanism between site region and that in which attenuation formulas were derived. The result showed that the PGA at bedrock for Penang and Kuala Lumpur range between 2 and 5 gals. Local site effect had also been analyzed for Penang and Kuala Lumpur. The analysis was carried out using 1- D shear wave propagation theory by using three-acceleration time histories records (i.e. Elcentro N-S 1940, Loma Prieta, and synthetic) as input data. Synthetic acceleration time histories were generated using random vibration theorem. The results of shear wave propagation analysis showed that strong motion data influenced amplification of motion and the shape of response spectrum. According to the results, acceleration at the bedrock had been amplified about 2 to 7 times at the surface of Penang and Kuala Lumpur. INTRODUCTION The researches regarding earthquake engineering in Malaysia are relatively behind compared to other engineering fields. This is due to the fact that Malaysia earthquake event in history is not so profound and the nearest distance of earthquake epicenter from Malaysia is approximately 350 km. Generally; earthquake can cause significant damages within 100-200 km radius from the fault or epicenter. At farther distance 1 2 3 Assoc. Prof. and Head of Structural Earthquake Engineering Research (SEER), University Teknologi Malaysia, azelan@fka.utm.my Assoc. Prof. Faculty of Civil Engineering, University Teknologi Malaysia, aminaton@utm.my Ph.D. Student, Faculty of Civil Engineering, University Teknologi Malaysia, h3ndri@excite.com

amplitudes of incoming seismic shear waves are generally small [1], however, the Bowl of Jelly phenomenon as what had happened to Mexico City in 1984 should be considered more seriously. The phenomena have shown that an earthquake could have a significant effect although at longer distance due to long period component of shear waves. Although Peninsula Malaysia is located in the stable Sunda Shelf with low to moderate seismic activity level, tremors due to Sumatra earthquake had been reported several times. For instance, there are two large earthquakes near Sumatra occurred at the end of 2002 (M w =7.4) and in the early of 2003 (M w = 5.8). Although no casualties or damages were reported due to those earthquakes, the tremors have been causing panic to a lot of people in several cities in Peninsula Malaysia. Cracks on a few buildings in Penang due to the earthquake on 2 November 2002 also have been reported as well. This research is proposed to analyze the effects of those two earthquakes to Wilayah Persekutuan Kuala Lumpur and Putrajaya cities. The considerations that are included in this study are geological and seismological conditions, attenuation of earthquake wave propagation in bedrock, specific acceleration time histories, and local soil condition. TECTONIC SETTING Malaysia is located on the southern edge of major tectonic plate of Eurasian. This position is relatively closed to the boundary between Eurasian Plate in the northern side and Australian Plate in the southern side. Based on that location, generally tectonic features that affect Peninsular Malaysia can be divided into two classifications. The first classification is subduction zone. All of those earthquakes that occurred near convergent boundaries where Indo-Australian plate is being subducted under Eurasian plate are classified into this zone. The Indo-Australian plate is sliding approximately northward beneath Sumatra and Java, where the direction of convergence is N20oE and the overall rate convergence is 7.7-cm/year [2]. Generally, this Subduction zone can be divided into three segments, i.e. Sumatra Segment, Sunda Strait segment, and Java Segment. The Sumatra subduction zone is a very active feature. The largest thrust-fault earthquakes in the Sumatra subduction zone in the last two centuries were in the year of 1833, with the magnitude of 8.8-9.2, and that of 1861, with the magnitude of 8.3-8.5. The second classification is transform zone. All of those earthquakes occurred due to strike slip movement along clearly defined fault in the frontal arc area such as Sumatra Fault are classified as transform fault. The Sumatra fault is about 1900 km long structure that accommodates right lateral strike slip associated with the oblique convergence along the plate margin. Several large earthquakes have occurred in this zone. These events were included the 1926 Padang Panjang (MS = 6.75), the 1933 Liwa (MS=7.5), the 1964 Aceh (mb = 6.7) and the 1993 Liwa (MS = 7.2) earthquakes. Tectonic setting around Peninsular Malaysia can be seen in Figure 1. According to USGS, the earthquake on 2 November 2002 (NS2002) occurred on Sumatra subduction zone whereas the earthquake on 22 January 2003 (NS2003) occurred as a result of strike slip mechanism of Sumatra fault. In order to verify the mechanism of that earthquake, spatial distribution of earthquakes nearby NS2002 and NS2003 were constructed from earthquake data occurred from the year 1900 to 2002. The distribution of epicenters is presented in Figure 2. The hypocentral profile revealed downward dipping zone of seismicity nearby those two earthquakes are depicted in Figure 3 and 4.

Sumatra Fault Subduction Zone Figure 1. Tectonic setting of Sumatra Figure 2. Location of epicenter of earthquake from 1900-2003 (USGS-NEIC)

Depth (km) 0 20 40 60 80 100 120 140 160 180 200 Distance (km) 10 20 30 40 50 60 Figure 3. Hypocentral profile around NS2002 4 < Mw < 5 5 < Mw < 6 6 < Mw<7 Mw > 7 Earthq. 2/11/02 0 Kuala Lumpur Melaka Penang Kuala Trengganu Distance (km) 0 100 200 300 400 500 600 700 800 900 Depth (km) 50 100 150 200 Epicenter 4<M<5 5<M<6 6<M<7 M>7 Epicenter (Mw =5.8) 250 Figure 4. Hypocentral profile around NS2003 Based on Figure 3 and Figure 4, it can be seen the NS2002 was located on subduction zone and the NS2003 occurred on transform zone of Sumatra. SEISMIC HAZARD ASSESSMENT Seismic hazard assessment is performed in order to obtain peak ground acceleration (PGA) at bedrock for Penang and Kuala Lumpur. Analysis is carried out using deterministic method and using several appropriate attenuation relationships. In the study, the attenuation function proposed by Atkinson and Boore [3] was used to calculate PGA for subduction earthquake event whilst Campbell [4] formula was used for strike-slip earthquake event. Based on the calculation, the PGA s at bedrock for Penang and Kuala Lumpur due to NS2002 event were 4.84 gal and 3.32 gal, respectively and the PGA s at bedrock for those two cities due to NS2003 event were 4.86 gal and 2.23 gal, respectively. The contours of iso-acceleration of those two events at bedrock of Peninsular Malaysia are presented in Figure 5 and 6.

2.6 179.3 88.6 52.3 16.8 34.0 23.4 9.3 12.3 5.4 3.3 7.0 4.2 Figure 5. PGA contour at bedrock due to NS2002 event using Atkinson and Boore [3] formula Figure 6. PGA contour at bedrock due to NS2003 event using Campbell [4] formula GENERATION OF SYNTHETIC GROUND MOTION Acceleration time-histories are required in the analysis of shear wave propagation in soil deposits. Selection of time-histories appropriates for specific geological and seismological conditions play an important role for obtaining accurate results because the characteristics of motion such as peak acceleration, duration and frequency content will affect the result of shear wave propagation analysis. Several procedures can be used to select earthquake acceleration time-histories at bedrock. These procedures include: (1) utilization of motions

previously detected near the site; (2) utilization of motions previously recorded at other locations during similar size earthquake and at distance comparable to those under consideration; (3) estimation of a target spectrum and then generation of a synthetic time history whose spectral ordinates provide a reasonable envelope to those of the target spectrum. Procedure number (1) is rare to be conducted in Malaysia because there are no representative strong motions that can be used for analysis. Procedure number (2) is difficult to utilize at most locations because the number of recorded motions is not extensively enough to cover a sufficiently wide range of possibilities. The most commonly method to cover this shortcoming is by modifying the actual ground motion records. By using this method, acceleration time-histories can be obtained from rescaling the actual ground motions record. The disadvantage of this procedure is the seismological and geological mechanism may not appropriate with the local condition. The ground motion time histories recorded from Elecntro, Pasadena, and Mexico have been widely used in analyzing dynamic response in soil deposits and structure in region in which a lack of recording data such as Malaysia. In this study, procedure (3) is performed to generate synthetic ground motions using random vibration theorem. Procedure for generating synthetic ground motion based on random vibration theorem has been described by Gasparini and Vanmarcke [5]. This procedure is based on the fact that any periodic function can be expanded into a series of sinusoidal waves: () t = A sin( ω t + φ ) x i i i (1) i Where A i is the amplitude and φ i is the phase angle of the i th contributing sinusoidal. The amplitude A i are related to the (one side) spectral density function G(ω) in the following way: ωi A i = 2 G( ωi ) dω (2) o A relationship between the response spectrum and the spectral density function of ground motion at site can be seen in the following equation. 1/ 2 2 1 n n Sv G( n ) ω ω ω = G( ω) dω (3) π r ω s,p 0 n 1 4 ζs where: ζ ζ s = (4a) 2ζ ωn t 1 e 2 rs,p = [ 2 log{ 2n[ 1 exp( δy (s) π log 2n ]}] 1/ (4b) 1/ 2 4ζ t δ y ( s) = (4c) π ω t 1 n n = (4d) 2π log 0.368 The response spectrum in the equation (3) is obtained from seismic hazard analysis by using attenuation relationships for response spectrum. The power of the motion produced by using equation (1) does not vary with time. To simulate the transient character of real earthquakes, the steady-state motions are multiplied by a deterministic envelope function I (t). The artificial motion X (t) becomes:

X (t) = a(t) = I(t) An sin( ωn t + φn ) (5) n There are three different intensity envelope functions available such as trapezoidal, exponential, and compound [5]. In this study, duration and envelope intensity function were calculated based on a procedure proposed by Kuda [6]. In this study, the generation of synthetic ground motions by random vibration theorem was performed using SIMQKE [5]. The results of analysis can be seen in Figure 7. a (gal) a (gal) 6.0 4.0 2.0-2.0-4.0-6.0 0 5 10 15 20 25 30 35 time (sec.) 3.5 2.5 1.5 0.5-0.5-1.5-2.5-3.5 (a) Penang NS2002 0 5 10 15 20 25 30 35 time (sec.) a (gal) a (gal) 6.0 4.0 2.0-2.0-4.0-6.0 0 1 2 3 4 5 6 7 8 9 10 11 time (sec.) (b) Penang NS2003 3.5 2.5 1.5 0.5-0.5-1.5-2.5-3.5 0 1 2 3 4 5 6 7 8 9 10 11 time (sec.) (c) KL NS2002 (d) KL NS2003 Figure 7. Synthetic Accelerations at Bedrock of Penang and Kuala Lumpur LOCAL SITE EFFECT Local site effect strongly influences the characteristic of ground surface motions, i.e. peak acceleration amplitudes and shapes of response spectra. Comparison of peak acceleration attenuation relationships for sites underlain by different types of soil profiles shows distinct trends in amplification behavior [7]. According to Idriss [7], peak accelerations at the surface of soil deposits are slightly greater than that on rock. At higher acceleration levels, however, the low stiffness and nonlinearity of soft soils often prevent them from developing peak ground accelerations as large as those observed on rock. Local site conditions also influence the frequency content of surface motions and hence the response spectra they produce. The average normalized response spectra from Seed et al. [8] also show that at longer periods, spectral amplification is much higher for soil sites than for rock sites. This effect can be very significant, particularly when long period structures such as bridges and tall buildings are founded on such deposits. In this study, site response analyses is carried out using 1-D shear wave propagation theory resulting in the estimate of ground motion parameter, i.e. peak surface acceleration (PSA) and surface spectral acceleration. The analysis is performed using program NERA [9], which stands for Nonlinear Earthquake Response Analysis. In this program, soil model proposed by Iwan [10] and Mroz [11] were used in order to model nonlinear stress-strain curves of soil. This model used series of n-mechanical elements, having different stiffness and sliding resistance.

Dynamic Soil Properties The dynamic soil properties that are required in a site response analysis are maximum shear modulus, G max or shear wave velocity, V s. Several researchers such as Seed [12] and Ohsaki and Iwasaki [13] have proposed correlations to convert the results from static soil test obtained from laboratory or in situ testing to obtain dynamic soil properties such as G max, and V s. Locations of study in this analysis are Sekolah Kebangsaan Sungai Penang, Penang and Bandar Baru Sentul, Kuala Lumpur. In this study, dynamic soil properties were obtained using empirical correlations proposed by Seed et al. [12] for non-cohesive soil and Ohsaki and Iwasaki [13] for cohesive soil. Shear Wave Propagation Analysis In this study, shear wave propagation analysis had been analyzed using three strong motion data i.e. Elcentro N-S 1940, Loma Prieta and synthetic. The average PGA obtained from the results of shear wave propagation analysis using NERA [9] is presented in Figure 8. The results can be summarized in Table 1. Table 1. The results of shear wave propagation analysis No. Event Location PGA (gal) PSA (gal) Amplification 1 2 nd Nov. 2002 Penang 4.84 25.3 5.2 2 Kuala Lumpur 3.32 7.3 2.2 3 22 nd Jan. 2003 Penang 4.86 20.8 4.3 4 Kuala Lumpur 2.23 11.6 5.2 PGA (gal) 5.0 1 15.0 2 25.0 3 PGA (gal) 0 5.00 10 15.00 5.0 1 5.0 Depth (km) 15.0 2 Depth (km) 1 25.0 15.0 3 EQ 2002 EQ 2003 2 EQ 2002 EQ 2003 35.0 4 25.0 (a) Penang (b) Kuala Lumpur Figure 8. Distribution of acceleration against depth The results of analysis show that the peak acceleration at bedrock increase about 2 to 5 times at the surface due to the effect of local soil condition and characteristic of ground motion. According to Modified Mercalli Intensity (MMI) Scale [14], generally the intensity at the surface due to those earthquakes in Penang and Kuala Lumpur approximately correlate with scale V and IV on MMI scale, respectively.

Figures 9 and 10 show the comparison of response spectrum between earthquake on 2 November 2002 and 22 January 2003. Response spectrum are obtained by averaging the results from three acceleration time histories, i.e. Elcentro N-S 1940, Loma Prieta and synthetic. These figures showed the strong motion data influence amplification of motion and the shape of response spectra. According to these figures, the locations of maximum spectral acceleration are between period of 0.1 sec and 1.0 sec. It means that practically, the maximum effect of the motion will occur on the 1 to 10- storey building in Penang and Kuala Lumpur. SA (gal) 20 18 16 14 12 10 8 6 4 2 EQ 2002 EQ 2003 1 0.10 1.00 10 Period (sec.) Figure 9. Response spectrum at the surface of Penang SA (gal) 8 7 6 5 4 3 2 1 EQ 2002 EQ 2003 1 0.10 1.00 10 Period (sec.) Figure 10. Response spectrum at the surface of Kuala Lumpur CONCLUSION Based on the results of the attenuation formula, the PGA s at bedrock for Penang and Kuala Lumpur due to earthquake on 2 November 2002 were 4.84 gal, and 3.32 gal, respectively and the PGA s at bedrock for Penang and Kuala Lumpur due to earthquake on 22 January 2003 were 4.86 gal and 2.23 gal, respectively. Acceleration time-histories are required in the analysis of shear wave propagation in soil deposits. Selection of time-histories appropriates for specific geological and seismological conditions play an important role for obtaining accurate results because the characteristics of motion such as peak acceleration, duration and

frequency content will affect the result of shear wave propagation analysis. In this study, synthetic ground motions have been generated in the frequency domain by using random vibration theorem. Local site effect strongly influences the characteristic of ground surface motions, i.e. peak acceleration amplitudes and shapes of response spectra. Therefore, site response analysis is needed in order to predict the effect of the motion at the surface. The analysis was carried out using 1-D shear wave propagation analysis theory at Sekolah Kebangsaan Sungai Penang and Bandar Baru Sentul The results showed that the peak accelerations at bedrock increase about 2 to 5 times at the surface due to the effect of local soil condition and characteristic of ground motions. It can be concluded that intensities at the surface due to those two earthquakes in Penang and Kuala Lumpur approximately correlate with scale V and IV on Modified Mercalli Intensity (MMI) scale, respectively. The effect of the earthquake to building depends on the natural period of the building. The maximum effect of the motion will occur on the 1 to 10- storey building in Penang and Kuala Lumpur. REFERENCES 1. Lee, S.L., Balendra, T., Tan, T.S. A Study of Earthquake Acceleration Response Spectra at Far Field. US-Asia Conference on Engineering for Mitigating Natural Hazards Damage, Bangkok, Thailand, 14 18 December 1987. 2. DeMets, C., R.G. Gordon, D.F. Argus, and S. Steia. Current Plate Motions. Geophysics Journal International, v.101. p. 425-478, 1990. 3. Atkinson, G.M. and Boore, D.M. Stochastic Point-Source Modeling of Ground Motions in the Cascadia Region. Seismological Research Letters, Vol. 68, 1997. 4. Campbell, K.W. Prediction of Strong Ground Motion Using the Hybrid Empirical Method: Example Application to ENA. Bulletin of the Seismological Society of America, 2002. 5. Gasparini D.A Vanmarcke, E.H. SIMQKE (User s manual and Documentation). Dept. of Civil Eng. Massachusetts Institute of Technology, 1976. 6. Newjec Inc. Topical Report on Seismology (Step3) for Muria, Engineering Report. November 1996. 7. Idriss, I.M. Response of Soft Soil Sites During Earthquake. in J.M. Duncan, ed., Proceeding, H. Bolton Seed Memorial Symposium, BiTech Publishers, Vancouver, British Columbia, Vol. 2, pp. 273-289, 1990. 8. Seed, H.B., Ugas,C., and Lysmer, J. Site Dependent Spectra for Earthquake Resistance Design. Bulletin of the Seismological Society of America, Vol. 66, pp. 221-243, 1976. 9. Bardet, J.P., Ichii, K. and Lin, C.H. NERA-A Computer Program for Nonlinear Earthquake Site Response Analyses of Layered Soil Deposits. Manual, Department of Civil Engineering University of Southern California, 2001. 10. Iwan, W. D. On A Class Of Models For The Yielding Behavior Of Continuous And Composite Systems. Journal of Applied Mechanics, ASME, Vol. 34, pp.612-617, 1967. 11. Mróz, Z. On the Description Of Anisotropic Workhardening. Journal of Mechanics and Physics of Solids, Vol.15, pp.163-175, 1967. 12. [Seed, H.B., I.M. Idriss, and I. Arango. Evaluation of Liquefaction Potential Using Field Performance Data. J. Geotech. Eng., ASCE, 109(3): 458-482.,1983. 13. Ohsaki, Y. and R. Iwasaki On Dynamic Shear Moduli and Poisson's Ratio of Soil Deposits. Soils and Foundations, 13(4): 61-73, 1973. 14. Richter, C.F. Elementary Seismology. W.H. Freeman & Co., San Francisco, 1958.