ARMY TM 5-809-10-1 NAVY NAVFAC P-355.1 AIR FORCE AFM 88-3, CHAP. 13, SEC A
a TECHNICAL MANUAL
SEISMIC DESIGN GUIDELINES FOR
ESSENTIAL BUILDINGS
This copy is a reprint which includes current I pages from Changes I 1 DEPARTMENTS OF THE ARMY, THE NAVY, AND THE AIR FORCE 27 FEBRUARY 1986
89:121402368391212 FDR WASTE WM-1 I PDC TM 5-809-10-1 /NAVFAC P-355.1 /AFM 88-3. Chapter 13, Section A C1 CHANGE DEPARTMENTS OF THE ARMY, THE NAVY AND THE AIR FORCE No. 1 WASHINGTON, DC, 15 December 1986 TECHNICAL MANUAL SEISMIC DESIGN GUIDELINES FOR ESSENTIAL BUILDINGS TM 5-809-10-1/NAVFAC P-355.1IAFM 88-3, Chapter 13, Section A, 27 Feb- ruary 1986, is changed as follows: 1. Remove and insert pages below. New or changed text material is indicated by a vertical bar in the margin. Remove pages Insertpages Cover 1 and Cover 2...... Cover 1 and Cover 2 2. This transmittal sheet should be filed in the front of the publication for refer- ence purposes. The proponent agency of this publication Isthe Office of the Chief of Engi- neers, United States Army. Users are Invited to send comments and sug- gested Improvements on DA Form 2028 (Recommended Changes to Publi- cations and Blank Forms) direct to HQDA (DAEN-ECE-D). WASH DC 20314-1000.
By Order of the Secretaries of the Army, the Air Force, and the Navy:
Official: MILDRED E. HEDBERG JOHN A. WICKHAM, JR. BrigadierGeneral, United States Army General, United States Army The Adjutant General Chief of Staff CHARLES G. GABRIEL General, USAF Chief of Staff Official: NORMAND G. LEZY J. P. JONES, JR. Colonel, USAF RearAdmiral, CEC, U.S. Navy Directorof Administration Commander, Naval Facilities EngineeringCommand Distribution: Army: To be distributed in accordance with DA Form 12-34B, Requirements for Seismic Design for Buildings. AirForce: F Navy: FOREWORD
The seismic design guidelines manual was developed to meet one of the objectives for earthquake hazards reduction measures as promulgated by the Earthquake Hazards Reduction Act of 1977 (Public Law 95-124). The objective is the development and implementation of a technologi- cally and economically feasible, improved design and construction meth- ods and practices in areas of seismic risk to provide earthquake resistant structures which are especially needed in time of disaster.
This guideline manual provides the latest seismic design concepts for earthquake resistant structures by utilizing the dynamic analysis ap- proach. The concept is for essential buildings but also includes design provisions for high risk and irregular buildings. This manual also pro- vides methodologies and procedures to determine site-dependent earth- quake ground motions for sites anywhere in the United States. Two levels of earthquake motion are considered. At the first level, the struc- ture will be designed to remain elastic for damage control at a moderate earthquake and at the second level, the criterion requires that the struc- ture remains functional after a major earthquake. Also, commentary and design examples are included to provide a comprehensive applica- tions of the design methodologies for earthquake resistant facilities. The general direction and detailed development of this manual was under the supervision and guidance of the Office of the Chief of Engi- neers, Headquarters, Department of the Army, Washington, DC and necessary coordination was maintained with the Naval Facilities En- gineering Command, Headquarters, Department of the Navy, Washing- ton, DC and Directorate of Engineering and Services, Headquarters, Department of the Air Force, Washington, DC. TM 5-809-10-1 NAVFAC P.355,1 AFM 88-3, Chapter 13, Section A
TECHNICAL MANUAL DEPARTMENTS OF THE ARMY, THE NAVY No. 5-809-10-1 AND THE AIR FORCE NAVY MANUAL WASHINGTON, DC, 27 February1986 NAVFAC P-355.1 AIR FORCE MANUAL No. 88-3, CHAPTER 13, SECTION A SEISMIC DESIGN GUIDELINES FOR ESSENTIAL BUILDINGS Paragraph Page CHAPTER 1. GENERAL * Purpose and scope ...... 1-1 1-1 Backgrounda. 1-2 1-1 Preparation of project documents. 1-3 1-2 References and bibliography. 1-4 1-3 * CHAPTER 2. INTRODUCTION TO SEISMIC ANALYSIS Introduction. 2-1 2-1 General. 2-2 2-1 Ground motion caused by earthquakes . 2-3 2-1 Site effects . 2-4 2-5 Dynamic analysis of structure. 2-5 2-6 Nonstructural elements. 2-6 2-14 CHAPTER 3. SPECIFICATIONS OF GROUND MOTION Section I. Basic steps for specification of ground motion Introduction. 3-1 3-1 Definition of terms, glossary, and symbols. 3-2 3-3 General overview of seismic hazard analysis. 3-3 3-3 Section II. Procedure for site specific ground motion Determination of source seismicity. 3-4 3-11 Selection of the attenuation relation. 3-5 3-30 Site specific response spectra. 3-6 3-40 Interpretation and summary. 3-7 3-54 Section III. The ATC-3-06 method The ATC--06 method . 3-8 3-57 CHAPTER 4. CRITERIA FOR STRUCTURAL ANALYSIS Introduction. 4-1 4-1 General requirements...... 4-2 4-1 Elastic design provisions. 4-3 4-2 Post-yield analysis provisions. 4-4 4-7 CHAPTER 5. STRUCTURAL DESIGN PROCEDURE Introduction. 5-1 5-1 Preliminary design considerations. 5-2 5-1 General design procedures. 5-3 5-1 Designing for EQ-I . 5-4 5-12 Designing for EQ-I . 5-5 5-19 CHAPTER 6. NONSTRUCTURAL ELEMENTS Introduction. 6-1 6-1 General requirements...... 6-2 6-1 EQ-I provisions. 6-3 6-1 EQ-I provisions. 6-4 6-6 Architectural elements ...... 6-5 6-6 Mechanical and electrical elements ...... 6-6 6-6 Essential systems. 6-7 6-7 CHAPTER 7. STRUCTURES OTHER THAN BUILDINGS Introduction. 7-1 7-1 General requirements...... 7-2 7-1 Elevated tanks and other inverted pendulum structures. 7-3 7-1 Vertical tanks (on ground). 7-4 7-2 Horizontal tanks (on ground). 7-5 7-2 Retaining walls . 7-6 7-2 Buried structures . 7-7 7-2 APPENDIX A. SYMBOLS AND NOTATIONS . A-I APPENDIX B. REFERENCES. B-1 APPENDIX C. GROUND MOTION BACKGROUND DATA ...... C-1
I TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986
SEISMIC DESIGN GUIDELINES FOR ESSENTIAL BUILDINGS Paragraph Page APPENDIX D. DESIGN EXAMPLES-GROUND MOTION ...... D-1 APPENDIX E. DESIGN EXAMPLES-STRUCTURES...... E-1 APPENDIX F. DESIGN EXAMPLES-EQUIPMENT IN BUILDINGS ...... F-1 BIBLIOGRAPHY ...... Bibliography 1 GLOSSARY ...... Glossary 1 INDEX ...... Index I
LIST OF FIGURES Figure 2-1. Recorded acceleration at ground level for three components of motion 2-2. Ground acceleration and integrated ground velocity and displacement curves 2-3. Description of acceleration response spectrum 2-4. Response spectra from recorded ground acceleration shown in Figure 2-1, transverse (north) 2-5. Single-degree-of-freedom system 2-6. Multi-degree-of-freedom system 2-7. Multi-mass system represented by a single-mass system 2-8. Design response spectra for examples in figures 2-9 and 2-10 2-9. Sample modal analysis of a 30-story building 2-10. Sample modal analysis of a 7-story building 2-11. Response of flexibility mounted equipment in buildings 3-1. Selection procedures 3-2. General flow diagram selection chart 3-3. General flow chart 3-4. Flow diagram for the Western United States 3-5. Hazard evaluation of Western United States 3-6. Flow diagram for the Eastern United States 3-7. Hazard evaluation of Eastern United States 3-8. Regional differences 3-9. Flow chart for Step 1, source identification and modeling for the Western United States 3-10. Point, line, and area sources 3-11. Dipping plane source 3-12. Flow chart for Step 1, source identification and modeling for the Eastern United States 3-13. Seismic sources after Algumissen and Perkins 3-14. Seismic sources after Hadley and Divine 3-15. Seismic sources after Tera 3-16. Flow chart for Step 11, seismicity and recurrence relationships for Western United States and Eastern United States 3-17. Linear recurrence relationship 3-18. Bilinear recurrence relationship 3-19. Recurrence relation for North San Andreas 3-20. Flow chart for Step 111, seismic forecasting model 3-21. Step IV, attenuation of ground motion from source to site 3-22. Attenuation distances 3-23. OASES attenuation 3-24. Attenuation relations 3-25. Comparison of ground motion models for Mb = 5.5 3-26. Description of sets of M and R required for a given PGA 3-27. Step V, site specific response spectra 3-28. Newmark-Hall spectrum 3-29. Statistical averaging of normalized spectra 3-30. Average acceleration spectra for different site conditions 3-31. Eighty-four percentile acceleration spectra for different site conditions 3-32. Predominant periods in rock, earthquake magnitude 7 3-33. Predominant periods for maximum acceleration in rock 3-34. Comparison of DAF from Kiremidjian and Shah to DAF from Seed et al, soil class = 0, damping = 5% 3-35. Comparison of DAF from Kiremidjian and Shah to DAF from Seed et al, soil class = 1, damping = 5% 3-36. Comparison of DAF from Kiremidjian and Shah to DAF from Seed et al, soil class = 2, damping = 5% 3-37. Factors affecting spectral shape 3-38. Envelope quality of the DAF shape 3-39. Hazard curve for site PGA with exposure time of 50 years 3-40. Contour map for effective peak acceleration ii 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A
Figure 3-41. Contour map for effective peak acceleration 3-42. Contour map for effective peak velocity-related acceleration coefficient 3-43. Contour map for effective peak velocity-related acceleration coefficient 3-44. Schematic representation showing how effective peak acceleration and effective peak velocity are obtained from a response spectrum 3-45. Annual risk of exceeding various peak accelerations for locations on the indicated contours of A. and A, 3-46. Tripartite representation of EQ-1i 3-47. EQ-Il spectra for A. = A, = 0.40 and i3 = 5 percent 3-48. Effective spectra] envelope 3-49. Regional shape differences 3-50. Las Vegas, Nevada, site spectra for soil type S2 3-51. Emeryville, California, site spectra for soil type S3 4-1. Definition of inelastic demand ratios for flexural members 4-2. Ductility check of steel columns 4-3. Ductility check for concrete columns 5-1. Gravity/seismic load relationships 5-2. Dynamic structural characteristics 5-3. Nonproportional relationship between peak ground acceleration and spectral acceleration 5-4. Sample EQ-I spectrum and ZICS curve 5-5. Force-displacement capacity curve 5-6. Capacity spectrum method 6-1. Design M.F. vs. period ratio 6-2. Sample roof response spectrum 6-3. Post-yield M.F. curve C-I. Earthquake source model C-2. Types of fault slips C-3. The Richter scale CG4. Relation between earthquake magnitude and intensity C-5. McCann and Shah relationship C-6. The PGA-MMI relationship with the intervals associated with each intensity C-7. Single-degree-of-freedom system C-8. Maximum dynamic load factor for sinusoidal load C-9. Judgmental averaging of empirical and analytical site spectra C-10. Relative degree of fault activity D-i. Source models and records for sources I and 2 D-2. Recurrence relation for source I D-3. Recurrence data plot for source I D-4. Source location and element properties D-5. Probability calculations for event combinations giving the hazard P (PGA > 0.20g) DU-. Site hazard curve and scaled site spectrum for EQ-I D-7. Scheme of present seismic hazard methodology U-8. General flow chart for seismic hazard analysis D-9. Earthquake listing for example 2 D-10. Output for recurrence relationship, example 2 D-lI. Recurrence relationship for example 2 D-12. Output for bilinear recurrence relationship, example 2 0-13. Bilinear recurrence relationship for example 2 D-14. Seismic sources for region of example 3 D-15. Earthquake listing for sources in example 3 D-16. Output for recurrence relationships and site PGA probability distribution for example 3 D-17. Complementary cumulative distribution function for example 3 D-18. Acceleration zone graph (AZG) for CITY 2 E-l. Sample modal analysis E-2. Building with a box system E-3. Building with steel moment-resisting frames and steel braced frames E-4. Seven-story ductile concrete frame building F-I. Cooling tower in building F-2. Unit heater-flexible brace F-3. Unit heater-rigid support '-4. Tank on a building
iii -
TM 5-809-1 0-1 /NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 LIST OF TABLES Table 3-1. Return period as a function of exposure time and probability of exceedance 3-2. OASES attenuation constants for median PGA values 3-3. Spectrum amplification factors for horizontal elastic response 3-4. Map contour and ground motion levels 3-5. Site soil profile types 3-6. Site profile coefficient 3-7. Damping adjustment factors 4-1. Damping values for structural systems 4-2. Inelastic demand ratios 5-. Seismic design procedures 5-la. Seismic design of essential facilities 5-lb. Seismic design of high-risk buildings 5-ic. Seismic design for other buildings 5-2. General modal relationships 5-3. Seven-story building-transverse direction-summary of modal analysis 5-4. Conversion of V and BN to S. and T 6-1. Example of a response amplification curve for the building's fundamental mode of vibration. 6-2. Data for the floor (roof) response spectrum example of figure 6-2 6-3. Essential nonstructural systems 7-1. Damping values for structures other than buildings C-l. Magnitude and seismic moment C-2. The Modified Mercalli intensity scale C-3. The Rossi-Forel scale CA-. Relationship between Modified Mercalli intensity (MM) and Rossi-Forel intensity (RF) C-5. Relationship between MMI and PGA C-6. Magnitude-displacement relationship C-7. Displacement-fault length relationship CG-. Magnitude-fault length relationship C-9. Magnitude-length times displacement relationship C-10. Magnitude-length times squared displacement relationship C-li. Degree of fault activity D-1. Return period vs PGA for CITY 2 E-l. Design examples-structural F-I. Design examples-equipment in buildings
iv 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A CHAPTER 1 GENERAL
1-1. Purpose and scope. are based on DOD standards; however, the risk a. Purpose. This manual prescribes criteria levels may be revised, as warranted, by approval and furnishes guidelines for the design of es- authorities. sential buildings, high-risk buildings, and other d. Classificationof structures. structures that may require analytical proce- (1) Hazardouscritical facilities. These fa- dures that are beyond the scope of TM 5-809- cilities (e.g., nuclear power plants, dams, and 10/NAVFAC P-355/AFM 88-3, chapter 13, "Seis- LNG facilities) are not included within the scope mic Design for Buildings." Methodologies and of this manual, but are covered by other publi- procedures are given for determining site-de- cations or regulatory agencies. For any facilities pendent ground motion and for the dynamic housing hazardous items not covered by criteria, analysis of buildings. These criteria apply to all advice should be sought from DAEN-ECE-D elements responsible for design of military con- (Army), NAVFAC Code 04BA (Navy), or HQ struction located in seismic regions. This man- USAF/LEEE (Air Force). ual is a supplement to TM 5-809-1O/NAVFAC (2) Essential facilities. These are struc- P-355/AFM 88-3, chapter 13, referred to herein tures housing facilities that are necessary for as the Basic Design Manual. post-disaster recovery and require continuous operation during and after an earthquake. This b. Scope. Approval from DAEN-ECE-D includes facilities where damage from an earth- (Army), NAVFAC Code 04BA (Navy), or HQ quake may cause significant loss of strategic and USAF/LEEE (Air Force) is required for the use general communications and disaster response of this manual as an alternative requirement to capability. This category also includes facilities applicable provisions of the Basic Design Man- serving an essential military function that must ual. This manual is for guidance in the design not be disrupted. Typical examples are listed in of buildings and other structures housing es- the Basic Design Manual, paragraph 3-Sa. sential mission-oriented facilities and those that (3) High-risk. This classification includes are vitally needed for post-disaster recovery that those structures where primary occupancy is for require continuous operation during and after assembly of a large number of people; where the an earthquake. This manual may also be used primary use is for people that are confined; or for guidance in the design of buildings that are where services are provided to a large area or classified in a high-risk category; buildings that large number of other buildings. Buildings in are irregular in shape, size, and configuration this classification may suffer limited damage in that require consideration of the dynamic char- a large earthquake, but are recognized as war- acteristics of the structure; and all other build- ranting a higher level of safety than the average ings as an alternative to the equivalent lateral building. Typical examples are listed in the Basic static force procedure for determination and Design Manual, paragraph 3-5b. distribution of seismic forces. These guidelines (4) All others. The provisions of this man- encompass: (1) assessment of the seismic haz- ual may be used for irregular buildings or as an ard at the site; and (2) seismic design of the option for all other buildings not covered by the structural and nonstructural systems for new above paragraphs only with the consent of the buildings and other structures. The problems approval authority. relating to earthquake-induced ground failure (e.g., liquefaction) are already stated in Basic Design Manual paragraph 2-7 and will not be 1-2. Background. covered in this manual. Alterations or evalua- a. Expectations. Current seismic design cri- tions of existing structures are not specifically teria, such as prescribed by the Basic Design covered by this manual; however, the principles Manual, consist of specified equivalent lateral and guidelines contained herein may be adapted static forces that are resisted by the designed for such use. structural systems. Structures designed in con- c. Seismic hazardrisk levels. Seismic ground formance with such provisions and principles motion input for two risk levels is specified in are expected to be able to: (1) resist minor chapter 3 for the prescribed structural perform- earthquakes without damage; (2) resist mod- ance criteria in chapter 4. The selected risk lev- erate earthquakes without structural dam- els of the two earthquakes (EQ-I and EQ-Il) age, but with some nonstructural damage; and 1-1 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 (3) resist major or severe earthquakes without occurring during the life of the building (e.g., major failure of the building or its component 10-percent chance of being exceeded in 100 years). members and equipment, and to maintain life In the first phase of the procedure, the building safety. For most structures, even in a major is structurally designed to resist the lower level earthquake, structural damage should be lim- earthquake within prescribed bounds of elastic- ited to repairable damage. It is also recognized linear procedures. In the second phase of the that for certain critical facilities, particularly procedure, the building is analyzed for its re- those essential to public safety and well-being sponse to the higher level earthquake by means in case of emergency, criteria should be avail- of procedures that account for inelastic behav- able to the designer that will permit design of ior, ductility demands, potential instability, and a facility that will remain operational during damage control. These guidelines are intended and after an earthquake. to insure that essential facilities will be capable b. Lessons learned. Recent earthquakes have of resisting the two levels of earthquake ground demonstrated that the existing seismic design motion as follows: (1) for ground motion as- requirements, as they have been implemented, sociated with the maximum probable earth- are not necessarily adequate to insure contin- quake, only minor damage, if any, will occur and ued operation of critical facilities vitally needed the facilities will not have any loss of function; after a major earthquake, such as hospitals, fire and (2) for ground motion associated with the stations, and communications centers. There- maximum theoretical earthquake, no cata- fore, there is a need for a more realistic ap- strophic failures will occur, damage will be re- proach to seismic-resistant design for buildings pairable, and essential facilities will remain that must remain continuously functional after functional. The definitions and the methodology a major earthquake. for determining these earthquakes are covered c. Recent developments. Earthquake engi- in chapter 3. The criteria and procedures for neering research and data collected from ground design are covered in chapters 4 and 5. motion instrumentations and earthquake-caused building responses during the last two decades 1-3. Preparation of project documents. have greatly increased knowledge in geotechni- a. Design analysis. A design analysis con- cal fields and have presented a clearer under- forming to agency standards will be provided standing of the performance of materials and with final plans. This design analysis will include ) structural elements. Therefore, practicing en- seismic design computations for the determi- gineers are able to become familiar with meth- nation of ground motion charateristics, for the ods of dynamic analysis as they are exposed to determination of dynamic characteristics of the new design procedures by means of technical structure, for the stresses in the lateral-force- publications, conferences, and continuing edu- resisting elements and their connections, and cation programs. for the resulting lateral deflections and inters- d. Design philosophy. One way of attempt- tory drifts. The first portion of the Design Anal- ing to reduce the risk of earthquake damage to ysis, called the Basis of Design, will contain the buildings is by imposing a higher design force following specific information: coefficient, such as an I-f actor of 1.5, for essen- (1) A statement on the methodology used tial facilities. This is not always a sufficient or for determining the ground motion criteria, and satisfactory approach to seismic design. In- a description of the response spectra for which creasing the design forces by 50 percent may be the structure will be designed. insignificant if a major earthquake results in (2) A description of the structural system demands several times the design capacity. On selected for resisting lateral forces and a dis- the basis of current knowledge, it appears that cussion of the reasons for its selection. A sym- a two-level (or two-phase) approach to design metrically configured lateral resisting framing will give better insight to postulated behavior system, without vertical irregularities, will be of structures. In this procedure, geotechnical data required. However, if irregular conditions are and probabilistic techniques are used to postu- unavoidable, a statement describing special late the motion for two earthquakes: (1) the analytical procedures to account for the irreg- maximum probable earthquake, which is likely ularities will be submitted for review and ap- to occur one or more times during the life of the proval by the approval authority. building (e.g., an earthquake with a 50-percent (3) A statement regarding compliance with chance of being exceeded in 50 years); and (2) this manual, including a list of the values se- the maximum theoretical earthquake that can lected for damping and maximum inelastic de- occur at the site, but has a low probability of mand ratios for critical structural elements. 1-2 27 February 1986 TM 5-809-lO-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A (4) Any possible assumed future~expansion riods of vibration and equivalent design lateral for which provisions are made. forces and other factors. b. Drawings. Preparation of drawings will (c) Assumptions made for future exten- conform to agency standards for ordinary con- sions or additions. struction, with the following additional specific (3) Site adaption of standard drawings will requirements for seismic construction: include design revisions for the seismic area as (1) Preliminary drawings will contain a required. statement that seismic design will be incorpo- rated in accordance with this manual. The Basis of Design will comply with paragraph a 1-4. References and bibliography. above. Publications that may be required to supple- (2) Construction drawings for seismic areas ment the provisions of this manual are listed in will include the following additional special in- appendix B, References. Publications that may formation: be useful as back-up material and are presented (a) A statement on the seismic ground as suggested reading are included in the bibli- motion criteria including the design peak ground ography. When pertinent to the subject, some accelerations and related response spectra. publications in the bibliography are noted in the (b) A statement on the lateral-force de- text by the bibliography number, in parenthesis. sign criteria including a tabulation of the pe- 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A CHAPTER 2 INTRODUCTION TO SEISMIC ANALYSIS
2-1. Introduction. 2-3. Ground motion caused by This chapter provides an introduction to the basic earthquakes. concepts of dynamic analysis for buildings re- A general introduction to earthquake ground sponding to the ground motions caused by motion is presented in the Basic Design Manual. earthquakes. General guidance is given in the The relationship of a ground motion to the site selection and use of various procedures for the and an introduction to time-history and re- design of structural systems. sponse spectra are presented herein. A detailed methodology for determining site-specific ground 2-2. General. motion characteristics is covered by chapter 3 An earthquake causes vibratory ground mo- of this manual. tions at the base of a structure and the structure a. General. actively responds to these motions. Seismic de- (1) Ground motion is generally strongest in sign involves two distinct steps: (1) determining the vicinity of its source (e.g., a rupturing fault), or estimating the forces that will act on the with the severity of shaking diminishing with structure; and (2) designing the structure to an increase in distance. resist these forces and to keep deflections within (2) The predominant periods of ground mo- prescribed limits. tion vibration generally lengthen as distance in- a. Determination of forces. There are two creases from the source (para 3-6f). general approaches to determining seismic (3) Deep deposits of soft soils tend to pro- forces: (1) an equivalent static force procedure, duce ground surface motions having predomi- such as presented in the Basic Design Manual; nantly long period characteristics. and (2) a dynamic analysis procedure. This man- (4) Deposits of stiff soils or rock result in ual illustrates the dynamic analysis procedure. ground motions having predominantly short pe- Seismic forces are determined from data derived riod characteristics. from the specification of ground motion. These b. Time history. The basic measurement of ground motion data will generally be given in earthquake ground motion is the accelerogram terms of a response spectrum; however, in some record taken by seismometers. When these in- cases the data may be described in terms of a strument records are properly corrected for digitized time history. elimination of recording noise and for base line b. Design of the structure. Structures are adjustment, a primary data base for seismic load generally designed to resist applied forces well specifications is provided. Data banks of past within the elastic capacity of their structural earthquake records from all parts of the world members. This is accomplished either by pre- are readily accessible from earthquake research scribing maximum allowable working stresses centers. A typical seismometer station provides for materials, or by using a strength design con- records of two orthogonal horizontal motions cept with prescribed load factors. However, for and one vertical motion, as illustrated in figure exceptional loading conditions, such as caused 2-1. The corresponding processed accelero- by major earthquakes, structures may be re- grams are intended to be the best representa- quired to resist deformations that exceed the tion of the actual ground acceleration at the elastic capacities of the structural elements. In recording site. For a given component, the time conventional methods of seismic design, it is as- derivative relations between ground displace- sumed that the design criteria will provide ad- ment, x(t); velocity x(t); and acceleration, x(t), equate safety by means of load factors and special allow the presentation of each of these motion details that provide the necessary ductility to histories, as shown in figure 2-2. The maximum resist major earthquake deformations. In the or peak values of displacement (PGD), velocity methods presented in this manual, the design (PGV), and acceleration (PGA) provide the most procedures will give a better insight as to the elementary and popular measures of an earth- performance of a structure when subjected to quake's severity. Duration (or bracketed dura- the exceptional loading conditions of a major tion) of strong motion is also an important earthquake. This method is generally referred measure, but it is not explicitly used in design to as a two-level approach to structural design. criteria at the present time. TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986
RECORDED ACCELERATION AT GROU11D LEVEL DURING THE 1971 SAN FERNANDO EARTHQUAKE
Transverse (North)
0.25^- !1 M A
-0.25 , V V. Y .
*0 0.25 Longitudinal (West) 0.2 CU OA R mI .x%' e.~r<-111
-q a)j -0.25
U
0.25 ) 0 -0.25
* 8 * U M81 Al 81 8& 8
0 3 6 9 12 15 18 21 24 27 30
Time (seconds)
Reprinted from "The San Fernando, Cali- fornia, Earthquake of February 9, 1971," U.S. Covernment Printing Office, 1971. Figure 2-1. Recorded acceleration at ground level for three components of motion.
)
2-2 27 February 1986 TM 5-809-10-1/NAFAC P-355.1/AFM 88-3, Chapter 13, Section A
IMPERIAL VALLEY EARTHOURKE MAY 18. 1940 - 2037 PST IRO01D40.001.0 EL CENTRO SITE IMPERIAL , VALLEY IMRIATION DISTRICT COMP SODE c PEAK VALUES: ACCEL: 341.7 cw(C/C VELOCITY = 33.4 CK/SEc OISPL = 10.9 cm -500
U TSJ V Li0 aJ w Cc "i IC x ki
500 -'0
0- 0 _ K LL. 4J
40 -20
z r
U.)U -J 02 20 I _ __ _ 0 30 MO so TIME - SECONDS
Reprinted from "United States Earthquakes." U.S. Coast and Geodetic Surveys, U.S. Govern- ment Printing Office, 1940. Figure 2-2. Ground accelerationand integratedground velocity and displacement curves. TM 5-809-10-1/NAVFAC P-3S5.1/AFM 88-3, Chapter 13, Section A 27 February 1986 c. Response spectra. For design purposes, it several records can be normalized, averaged, and would be ideal to forecast the acceleration time then scaled according to seismicity to predict history of a future earthquake having a given future ground motion at a given site. The phys- hazard of occurrence. However, the complex ical definition of an acceleration response spec- random nature of an accelerogram makes it nec- trum is shown in figure 2-3. A set of linear elastic essary to employ a more general characteriza- single-degree-of-freedom (SDOF) systems hav- tion of ground motion. Specifically, the most ing a common damping ratio, I, but each having practical representation is the earthquake re- different harmonic periods over the range 0, T1, sponse spectrum. This spectrum is used not only T2, etc. is subjected to a given ground motion to describe the intensity and vibration fre- accelerogram. The entire time history of accel- quency content of accelerograms, but also the eration response is found for each system, and most important advantage is that spectra from the corresponding maximum value, Sa, is plotted
JrJTEM RAspoJAVJX
______GOUIN/ ____PA CCCE(.eATION LINeA urn OF G/VCEN bAMP/NC _~~/iP? WITH AtANGC- OP NA7UIAL 'C.1M~ It7T.. ACCE4ERA7IOm )
ACCELI~eATION dtEJR0 ONJe ~~~~~c e ~l ,MAK.
GAOUN;O ACCEi.eACIGAM 7.
TIME
US Army Corps of Engineers Figure2-3. Description of acceleration response spectrum. 2-4 27 February 1986 TM 5-809-10-1/NAVFAC P-355.I/AFM 88-3, Chapter 13, Section A on the period axis for each system period. The curve provides the maximum response value for curve connecting these Sa values is the accel- any given system period, T. eration response spectrum for the given acce- lerogram and damping ratio. Actual spectra for 2-4. Site effects. the transverse (north) accelerogram of figure a. Response spectrum shape. Response spec- 2-1 are shown for several damping ratios in fig- tra shapes are determined largely by empirical ure 2-4. A smoothed individual spectrum (fig 2- data. Time history records of past earthquakes 4b), or averages of multiple record spectra, is are used to construct response spectra. As the employed as the seismic load input for the dy- data bank increases, average trends can be ob- namic analysis of structures. Note that the Sa served with respect to the general shape of re-
2.0- SAN FERNANDO EQ 2/9/71 VAN NUYS HOLIDAY INN 1ST FLOOR NORTH
Damping - 0.01, 0.02, 0.05, 0.10
1.2-
I-J .8- -,, = 0.01 Li .0 = 0.02 0 = 0.05 _B = 0.10 A.- 1^ .4 U - 0 .3 .6 .9 1.2 1.5 1.8 2.1 2.4 2.7 3.0 PERIOD, SEC (T) a. RESPONSE SPECTRA: FOUR VALUES OF DAMPING ($)
1 .2 4 B - 0.05 SPECTRUM (from a. above.)
.8 SMOOTH SPECTRUM ,
hExample: For T=0.3 sec. S_=O.7 Q 0 0.3 0.6 0.9 1.2 1.5 1.8 7.1 2.11 2.7 3.0
PERIOD, SEC. (T)
b. SMOOTH RESPONSE SPECTRUM: S - 0;05
US Anny Corps of Engineers Figure 2-4. Response spectra from recorded ground accelerationshown in figure 2-1. 2-5 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 sponse spectra curves. When these data are system is used in textbooks to illustrate prin- catalogued in terms of tectonic region, event in- ciples of dynamics. It represents two kinds of tensity, distance, and site characteristics, spe- real buildings: (1) a single-column structure with cific response spectra shapes can then be a relatively large mass at its top; and (2) a sin- developed that represent the conditions of par- gle-story frame structure with flexible columns ticular sites. Procedures for developing re- and rigid roof system. In the idealized system, sponse spectra are covered in chapter 3, and the mass (M) represents the weight (W) of the illustrative examples are included in appendix D. system divided by the acceleration of gravity (g) b. Soil column. Site soil characteristics can (M = W/g). The pole or columns represent the be used to develop a mathematical model of a stiffness (K) of the system, which is a ratio equal soil column at a building site. For a postulated to a horizontal force (F) applied to the mass bedrock earthquake, analytical procedures can divided by the displacement (8) resulting from be used to calculate the soil column's effect on that force (K = Fi8). If the mass is deflected and the ground motion at the surface or the base of then quickly released, it will freely vibrate at a a structure. These results can be used either to certain frequency, which is called its natural calculate the shape of the response spectrum of freuency of vibration. The period of vibration these particular conditions, or used directly for (T), which is the inverse of the frequency of time history analysis of the structure. vibration, is the time taken for the mass to move c. Foundation design. All inertia forces through one complete cycle (i.e., from one side originating from the masses on the structure to the other and back again (part b of fig 2-5). must be transmitted to and from the lateral- The period is equal to 2irV M/K. In an ideal sys- force-resisting elements, to the base of the tem having no damping (I = 0), the displaced structure, and into the ground. Foundations must system described above would vibrate forever. be designed to provide stability for response due In a real system where there is some damping, to maximum seismic ground motion. It should the amplitude of motion will decrease for each also be noted that the type, size, and depth of a cycle until the structure stops oscillating and foundation system can have an effect on a struc- comes to rest (part c of fig 2-5). The greater the ture's response to seismic motion and that the damping, the sooner the structure comes to rest. actual seismic input is a series of reversing load The amount of damping is defined in terms of a cycles. ratio, or percentage, of critical damping. If the structure has damping equal to 100 percent of 2-5. Dynamic analysis of structures. critical damping (I = 1.0), the displaced struc- Structures that are keyed into the ground and ture will come to rest without crossing the ini- extend vertically some distance above the ground tial point of zero displacement. If oscillating act either as simple or complex oscillators when motion is applied to the base of the system, the subjected to earthquake-caused ground motion. SDOF system will be forced to vibrate. If the Simple oscillators are represented by single-de- oscillating motion at the base is at a period equal, gree-of-freedom (SDOF) systems, and complex or nearly equal, to the period of the SDOF sys- oscillators are represented by multi-degree-of- tem, the motion of the mass will amplify until freedom (MDOF) systems. When a structure's it is substantially greater than the motion at the base is suddenly moved by earthquake ground base. This condition is called resonance. The motion, the upper part of the structure will not lower the value of A, the higher the amplifica- respond instantaneously, but will lag behind be- tion. cause of the structure's inertial resistance and b. Multi-degree-of-freedom systems. Multi- flexibility. This concept is illustrated in the Basic story buildings are analyzed as MDOF systems Design Manual, paragraph 2-4. As time pro- as shown in figure 2-6. They can be represented gresses during an earthquake, the structure's by lumped masses attached at intervals along various natural modes of vibration will be ex- the length of a vertically cantilevered pole (part cited to peak amplitudes of motion as described a of fig 2-6). Each mass can be deflected in one by the response spectrum (para 2-3c). direction or another; for example, all masses a. Single-degree-of-freed6m system. One may simultaneously deflect in the same direc- fundamental system that is investigated by dy- tion (the fundamental mode of vibration), or namic analysis is the simple oscillator or SDOF some masses may go to the left while others are system, as shown in figure 2-5. Represented by going to the right (higher modes of vibration). a single lump of mass on the upper end of a An idealized system, such as shown in part a of vertically cantilevered pole or by a mass sup- figure 2-6, has a number of modes equal to the ported by two columns (part a of fig 2-5), this number of masses. Each mode has its own nat- 2-6 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A
z FOR~CE M F : WEIGHT F rl r. MAJJ I--j rmi,10ACc-MEN P O*,K= F II K : JTIPFNE.hff 4 F GP&4V1TY
a. ID~EALIZED JINGLE- LUMMP-MASJ JYJ7GMJ
I 1PJL~kr PERII OF VAIA 4ArI OH: T: 2 Tr .F~w
I -I
/ ~04
b. FREE V/P.RAOTION ( No nAMPI NG) -
VhJ5RATE4 a 0-~'
N 1P.- t
- -
D~AMP'ER c, DAMPED FREAf VIQRA7TON (nAMIsE ATA RATIO 70 CRI7ICAL D4AMPIN& E(QUAL 70 a )
US Army Corps of Engineers Figure2-5. Single-degree-of-freedom system.
2-7 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986
4I m7 k 7 *m 6
l
4
min mL l
I k
I- - rn-_ -- Multi-Mass Fundamental Second Tiii rd System - Mode Miode Mode
a. Idealized Lumped b. Mode Shapes Mass System
US Army Corps of Engineers
Figure 2-6. Multi-degree-of-freedom system. ural modal period of vibration with a unique mode system. For many buildings, the partici- I/ ) mode shape being formed by a line connecting pation of the higher modes is negligible in re- the deflected masses (part b of fig 2-6). When lation to the participation of the fundamental oscillating motion is applied to the base of the modes of vibration. However, for tall, long-pe- multi-mass system, these masses move. The de- riod, and irregular buildings, the second, third, flected shape is a combination of all the mode and, possibly, higher modes may have a sub- shapes; but modes having periods that are near, stantial effect. The amount of higher mode par- or equal to, predominant periods of the base ticipation depends on both the building's modal motion will be amplified more than the other characteristics and the amplitude-period char- modes. Illustrative examples of MDOF systems acteristics of the response spectrum. Assuming are included in appendix E. that several modes are significant, one must se- c. Multi-mode response to ground mo- lect an appropriate method of combining the re- tion. Each mode of an MDOF system can be sults of the several modes. One method is simply represented by an equivalent SDOF system hav- to add up the effects of each mode (absolute ing a normalized mass (M*) and stiffness (K*) sum). This is an overly conservative approach where the period equals 21 7V'M-7IK (M* and K* because the response spectrum gives the peak are functions of mode shapes, mass, and stiff- response of each mode, and different modes reach ness). This concept, as shown in figure 2-7, pro- their peak amplitudes at different times during vides the computational basis for using site the earthquake. Since the spectrum gives only specific earthquake response spectra based on the maximum values and the time of occurrence SDOF systems for analyzing multi-storied build- is unknown, some approximate method of mode ings. With the period, mode shape, mass distri- combination must be used. The method most bution, and response spectrum, one can compute commonly employed is to combine the modes by the deflected shape, story accelerations, forces, the square-root-of-the-sum-of-the-squares and overturning moments. Using the response (SRSS) of the peak response of each mode (this spectrum method on MDOF systems requires is analogous to a vector sum). This offers a rea- analyzing each predominant mode separately. sonable value between the upper bound as the Results of each individual modal analysis must absolute sum of the modes and the lower bound then be combined in order to analyze the multi- as the maximum value of a single mode. To il- 2-8 27 February 1986 TM 5-809-10-1I/NVFAC P-355.1IAFM 883, Chapter 13, Section A
6 roof I
m7. t -- "- f7 = m 7a7
k7 M4 ---% f6 m6aG kb
m5 I9 -- f5 = m5as F = M*Sa
k5
m4 I ~~~f4 -M~a 4
m3 -'--f 3 m 3a3
k3 K* m12 f2 m2a2
f ml aI
First Mode of a Equivalent Multi-Mass System Single-mass System
M* and K: are normalized values of mass and stiff- ness that represent the equivalent combined effects of the story masses (M) in(I stiffnesses (K)
US Army Corps of Engineers Figure2-7. Multi-mass system representedby a single-mass system.
2-9 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 lustrate the multi-mode analysis of multi-sto- for modal analysis examples of a 30-story build- ried buildings, two examples are given. Figure ing and a 7-story building. 2-8 shows design response spectra that are used
0.8
0.7 . _ __ __
0.6
0.5 _B=
0.4
2%20
O.1 _...B=19
1.0 2.0 3.0
PERIOD, T(sec)
, ~~~~~SPECTRAL ACCELERATION, Sa~w
i._00.10.48 .50 0.80 1.0 1.25 1.75! 2.0 2.25: 2.5 3.0
2 1 0.64 0.64 0.59 0.37 0.30 0.24 0.201 .171 0.151 0.131 0.121 0.10
5961 % 0.50 0.50 0.48 0.301 0.24 0.192 0.1610-1371 0.12 0.107;0.096 ':0.08 0.50 __ ---. 0.480.24 IJ----- 0.1892j00.9'O0 -- - 1------
7' 0.44 0.44 0.44 0.28 0.22 0.18 0.l5' 0.13 0.11 0.10 0.09 0.07
I . . .0 0.08_..._1... _ --b---- 10% 0.38 0.38 0.38 0.25 0.20 0.16 0A13 0.11. l i .08 0.0661 . .. 27 -7 - 0I2 _ I 0 _08 _ ._.0_
2 0.27 0.27 0.27 0.20 0.16 0.12 0! 000 0-°ioo.081fi0.07 ,0.06 0.05
US Army Corps of Engineers Figure2-8. Design response spectra for examples in figures 2-9 and 2-10. 2-10 27 February 1986 TM 5-809-1O0-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A (1) Thirty-story building. The example in seconds. The periods of the second and third figure 2-9 summarizes the results of a modal modes of vibration are 1.00 seconds and 0.56 sec- analysis of a structural framing system that onds, respectively. From the response spectrum represents one principal axis of a 30-story build- curve in figure 2-8, which represents 5 percent ing. The fundamental period of vibration is 3.0 of critical damping (A = 0.05), it is determined
1ST MODE 2ND MODE 3RD MODE T,=0.56 T1 =3.00 sec T2 =1.00 sec sec 29 ' .29 25 -. 02 21 -. 28 17 13 9 C Groun~d
U- (K) 25 5790 .104 .104 .094 .175 -' - - 25 584t .097 .037 .113 21 5841 .02 7 [-0_ - o - -- 17 58L1 L_' _- .D08h 8( 13 5841 .1321*1!5 036 -;097 .083 _ 9 5841 .020 -. 065 6211 'F, .11_5 Groun c 0 _C_CE E RT4 __g_ (b) MODAL STORY ACCELERATIO~NS (gl!S.)
N: 2C 250 * 601 604 -. E 543 _ __ .. - - 25 215 , 568 7 335 3_ _ 21 180 502 145 * 425 17 /* a 110 * 326 _ 13 I 212 I j -569 484 9 75 V- 40 121 405 57(, - 5 0 Grou~nd 0 (c) MODAL STORY FORCES !1.ip-)
th 29 35 601 604 543 __ 1010 __ 25 35 1169 511 ~~~~1584L 21 35 1677 905 -5 1906 1 17 35* 2)02 ! ;~~~~ 521 T ~ ~ 592225 13 35* 2428 y 4_) J -6~~~~4 5'-2e 2486 t 9 35I 2640 -633 -44 2715 5 40* 2671 -1038 2997 Ground (d) MODAL STORY SHEARS kips)
29 0 0 0 0 25 21035 21140 K 19005 335363 21 61950 4nODS 36890 90084 I8 17 120645 i5640 36715 152461 13 194215 _ 3ID3915 18900 2210771 5 279195 101675 420 297133 5 371595 79520 -1120 3800)0 " Ground 482035 - -- '20160 -___ - -*t3e5_ (e) MODAL STORY OVERTURNING MOMENTS (kip-ft) *Story 29 represents the roof, floors 29 and 28, and one-half of floor 27. Other story designations represent the reference story plus one-and-one-half stories above and below.
US Army Corps of Engineers Figure 2-9. Sample modal analysis of a 30-story building. 2-11 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 that the second mode spectral acceleration higher modes become somewhat insignificant (0.240g) is triple that of the first mode spectral because of the reversal of force directions. The acceleration (0.080g), and that the third mode SRSS curve is essentially equal to the first mode spectral acceleration (0.45g) is over 5 times that curve at the lower stories of the building. of the first mode spectral acceleration. On the (2) Seven-story building. The example in basis of mode shapes and modal participation figure 2-10 summarizes the results of a modal factors (chap 5), modal story displacements, ac- analysis of a structural framing system that celerations, forces, shears, and overturning mo- represents one principal axis of a 7-story build- ments can be determined. For ease of comparison ing. Back-up data for this example are included to the 7-story example (para (2). below), the 30- in appendix E (design example E-1). The pe- story building is compacted to seven lumped riods of vibration are roughly 30 percent of the masses, each representing four stories. Back-up periods of the 30-story building (fig 2-9); pe- data for this example are included in appendix riods of the first, second, and third modes being E (design example E-1). The modal analysis 0.880 seconds, 0.288 seconds, and 0.164 seconds, procedure is covered in chapter 5. respectively. From the 5-percent damped re- (a) Diagram (a) of figure 2-9 shows the sponse spectrum (P = 0.05) of figure 2-8, both modal displacements. Note that the funda- the second and third mode spectral accelera- mental mode (first mode) predominates, while tions (0.500g) are 80 percent greater than the second and third mode displacements are rela- first mode spectral acceleration (0.276g). tively insignificant. The SRSS combination does (3) Comparisons. By comparing figures 2- not differ greatly from the fundamental mode. 9 and 2-10, it can be seen that the influences of (b) Diagram (b) shows story accelera- the second and third modes in relation to the tions. In this form, the second and third modes first mode are larger for the 30-story building do play a significant role in the structure's max- than for the 7-story building. For taller build- imum response. While the shape of an individual ings with longer periods of vibration, the influ- mode is the same for displacements and accel- ences of the higher modes may become larger, erations, accelerations are proportional to dis- and participation of additional modes of vibra- placements divided by the squared value of the tion (e.g., fourth and fifth modes) may become modal period, which accounts for the greater significant. accelerations from the higher modes. The shape d. Response of irregular buildings. When of the SRSS combination of the accelerations is buildings are eccentric or have areas of discon- substantially different from shapes of any of tinuity or other irregularities, the behavioral the individual modes because it accounts for the characteristic are very complex; whereas build- predominance of the various modes at different ings with symmetrical shape, stiffness, and mass story levels. Note that the maximum accelera- distribution and with vertical continuity and tions on stories 5 through 25 do not vary by more uniformity behave in a fairly predictable man- than 10 percent from the mean value, indicating ner. In addition to the single axis of response that the maximum acceleration felt at most floor shown in figures 2-9 and 2-10, the torsional re- levels is fairly constant. However, these maxi- sponse (twisting about a vertical axis) as well mum values would not occur simultaneously or as the interaction or coupling of the two trans- with the same period content. lational directions (longitudinal and transverse (c) Diagram (c) shows story forces whose axis) of response must be considered. For ex- values are obtained by multiplying the story ac- ample, the predominant motion may be skewed celeration by the story mass (or weight). The from the apparent principal axis. This is some- shapes of diagram (c) curves are quite similar what analogous to a Mohr's circle for principal to the shapes of diagram (b) curves because the stresses. Thus, three-dimensional methods of building mass is essentially uniform. analysis are required and each mode shape is (d) Diagram (d) shows story shears, defined in three dimensions by the longitudinal which are a summation of the modal story forces movement, the transverse movement, and the in diagram (c). The higher modes become less angle of rotation. In addition to complicating significant in relation to the first mode because the method of analysis, building irregularities the forces tend to cancel each other due to the complicate the methods used to combine modes. reversal of direction. Except for the top stories, Methods such as SRSS may not be appropriate the SRSS values do not differ substantially from for some three-dimensional methods of dynamic the first mode values. analysis. Procedures for performing three-di-' ) (e) Diagram (e) of figure 2-9 shows the mensional analyses are covered in chapter 5. building overturning moments. Again, the e. Inelastic-nonlinearresponse. In order to 2-12 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A
1ST MODE 2ND MODE 3RD MODE T1 =.880 sec T2 =.288 sec T3 =.164 sec SRSS Roo f 7 2.74 .19 1 .04 2.75 6 2.57 -.11 .00 2.57 5 2.30 -.04 2.30 -. 04 4 1.93 -.18I 1.93 3 1.46 -. I8 .00 1.48 2 .96 .02 I .97 - . -~~-.12~ .04 Ground .52 _ .,, (a) MODAL LATERAL DISPLACEMENTS (inches)
l?~~~~~~? Wt (K) Roo f 1410 .360 .234 71 .121 -7o .446 7 1460 .338 .129 -.007,, 1 .362 6 1460 .303 .324 5 1460 .254 -.148 -.112 .315 -.22 5 4 1460 .193 .1 .004- .297 -. 219 .275 1460 .127 3 -. 146 :.20 [) .201 2 1830 .068 Ground 0 (b) MODAL STORY ACCELERATIONS (g's)
Ht Roo f 65.7' 508 330 170 -- 6759 7 57. 0' 494 18B -10-*' A 529 6 48.3' 443 -1 -166 473 5 6' 371 39. -329 , -6 433 .0 4 30.9' 282 I.3 22.2' 185 -7 -319 156 W4oo D -267 219 - 367 Ground 13.5' 125 0 (c) MODAL STORY FORCES (kips)
^h Roof 8.7 508 330 170 , 629 7 160 -6 1139 6 8.7 1002 51 8 8.7 1445 499 1529 5 1846 4 8.7 1816 283 r I -169 2098 -16 2106 3 0 -175 2312 8.7 2283 -365 2 rd-632 200 2498 trrnun,d 13.5 2408 (d) MODAL STORY SHEARS (kips)
Roo f 0 0 0 0 7 4420 -2871 1479 5474 6 13137 -7378 2871 15338 54 25709 -11719 2819 28394 4 ~~41508 -14181 1349 28394 3 59761 -13781 -174 4368L 2 '79623 ______-10605 -339 61330 Ground 211231 I23. -- 20371 112175 (e) MODAL STORY OVERTURNING MOMENTS (kip-ft)
US Arm)y Corps of Engineers Figure 2-10. Sample modal analysis of a 7-story building.
2-13 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 estimate the behavior of a structure that may of motion corresponding to the building's nat- be subjected to a major, damaging-type earth- ural periods of vibration tend to be magnified. quake, it is necessary to investigate its inelastic In other words, a response spectrum of a build- response characteristics and capacity. The gen- ing's floor motion will have predominant peaks eral procedures discussed in paragraphs a at the participating periods of the building. If through d above are on the basis of elastic-lin- elements are rigid and rigidly attached to the ear distortions of the building's structural ele- structure, the maximum accelerations will be ments. When one major structural element the same as the maximum floor accelerations, begins to yield, changes will begin to occur in such as those shown in the SRSS curve of dia- the structure's behavioral characteristics. For gram (b) in figures 2-9 and 2-10. But, if the example, force distribution, periods of vibra- elements are flexible and have periods of vibra- tion, and mode shapes will be altered as parts tion close to any of the predominant building of various elements yield. Dynamic analysis pro- vibration modes, these elements will experience cedures for nonlinear systems can be very com- accelerations substantially greater than the floor plex, requiring step-by-step, time-history-forcing- accelerations. Generally, a time-history analy- functions, and inelastic force-distortion prop- sis is required to determine the peak response erties of all the structural elements and their of flexible or flexibly attached equipment at up- connections. However, approximate methods per levels of a building. A time-history of the have been developed that give rough approxi- ground motion is used to calculate a time-his- mations as to the inelastic response or capacity tory of the floor motion. The floor motion time- of structures. Post-yield analysis procedures are history is then used to construct a floor response discussed in chapter 5 and illustrative examples spectra. This procedure is illustrated in figure are included in appendix E. 2-11. In chapter 6, an approximate method is shown for constructing design floor response 2-6. Nonstructural elements. spectra. Illustrative examples are included in Elements that are housed in the building, as well appendix F. as portions of the building that are not part of b. Elements attached to adjacent floors. the structural system, must also be investigated Elements extending vertically from floor to floor for their response to earthquake motion. These (e.g., full-height partitions, exterior panels, pip- ) elements are generally categorized as architec- ing) will be subjected to two types of dynamic tural, mechanical, or electrical (refer to Basic motion. One type is the response motion de- Design Manual, chaps 9 and 10). scribed in paragraph a above. The other type is a. Elements attached to floors of build- due to the distortion resulting from the inters- ings. These elements (e.g., mechanical equip- tory displacements between two adjacent story ment, free-standing partitions, storage racks, levels. Interstory displacements for each mode suspended fixtures) respond to floor motion in can be obtained by finding the difference be- much the same manner as a building responds tween adjacent modal lateral story displace- to ground motion. However, the floor motion ments (diagram (a) in figs 2-9 and 2-10). may vary substantially from the ground motion. Interstory displacements for a multi-mode sys- The high-frequency components that make the tem can be approximated by combining the modal ground motion complex tend to be filtered out interstory displacements by the SRSS or other at the higher floor levels, while the components methods.
2-14 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A
Sa Flexibly Mounted FlAor'IRespne Equipment m
Floor Response Floor Response IX Spectrum
__=_~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
WN O.-W WI OIW WIN Al%"% - ~ Ground lotion
US Anny Corps of Engineers Figure 2-11. Response of flexibly-mounted equipment in buildings.
2-15 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A CHAPTER 3 SPECIFICATION OF GROUND MOTION
Section I BASIC STEPS FOR SPECIFICATION OF GROUND MOTION 3-1. Introduction. described in table 3-5 and having locations out- The purpose of this chapter is to describe the side of the limits of paragraph 3-1 b( 1), the ATC methodologies for determining site dependent 3-06 method of section III, paragraph 3-8 of this earthquake ground motions for sites anywhere manual may be used. in the United States. The objective is to develop (3) For sites in the WUS having exceptional design parameters from the available informa- soil conditions conforming to the soil descrip- tion and seismic ground motion. The principal tion of soil profile S3 as described in table 3-5, method of describing these ground motions will the selection of the corresponding site specific be in the form of acceleration response spectra response spectrum shape will consider and em- for input in the dynamic analysis of a given ploy the recommendations of paragraphs 3-6c(3) structure. or 3-6f(3) as directed by the responsible agency. a. Selected method of description. There are If this WUS site location is outside of the limits several methods of arriving at a description of of paragraph 3-1b(1), then the selected spec- future earthquake loading. These are described trum shape may be scaled by the appropriate briefly along with their advantages and disad- site acceleration coefficient A, given in paragraph vantages in appendix C, paragraph C-3. The 3-8. method employing an attenuated site severity (4) For sites in the EUS having the soil con- factor (such as peak ground acceleration, PGA) ditions conforming to soil profile S3 and outside which is used to scale a normalized site spectral of the limits of paragraph 3-lb(1), the method shape (Dynamic Amplification Factor, DAF) is of paragraph 3-8 may be used. judged to be the most appropriate and practical (5) In all cases where methods other than input for the dynamic analysis of building struc- those of paragraph 3-8 are employed, the re- tures and therefore will be the principal method sults will be compared with those from para- for this manual. However, this empirical method graph 3-8, and any significant differences will may be supplemented by available results from be justified and resolved. All final recommen- the other methods; particularly any findings from dations shall be subject to approval by the re- a site soil column response study, as described sponsible agency. in appendix C, paragraph C-3. c. Scope. The scope of this part of the Man- b. Procedures. The following selection pro- ual includes the description of the essential steps cedures will be followed for the evaluation of and related procedures necessary for the spec- site dependent earthquake ground motions, (see ification of site specific ground motion. These fig 3-1). These procedures are dependent upon are listed in paragraph 3-3 for the Western three conditions: the geotectonic regions of the United States (WUS) and the Eastern United Western United States (WUS) and the Eastern States (EUS), and for the deterministic and United States (EUS) as defined in paragraph 3- probabilistic procedures. 4a, the proximity of seismic sources, and the site d. Current state-of-the-art. It is important soil conditions as described in table 3-S5. to recognize that the field of ground motion (1) For sites located within 20 kilometers specification is in a state of evolution. The gen- from a fault or area source in the WUS, or within eral steps and input variables as outlined in this a tectonic province in the EUS, where the source manual are reasonably well accepted by most of or province has a maximum local magnitude of the researchers and users. However, because of 6.0 or greater, the detailed procedures of para- the very active state of development, it is not graphs 3-3 through 3-7 will be considered and possible to outline a step by step procedure which employed as directed by the responsible agency. will remain the same with time as well as from (2) For sites in either the WUS or EUS hav- region to region. Thus, the steps outlined in this ing normal site soil conditions conforming to manual are to be viewed as guidelines rather the description of soil profile types SI or S2 as than as one universally accepted and recom- mended procedure. Nothing in this chapter will prevent substantiated alter- e. Format of results. Various methods for native methods or time history procedures if approved by the evaluation of the level of ground motion and the agency command. its time history or frequency content are de- 3-1 TM 5-809-1O-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986
Western United States (WUS)
Source to Site Surface Distance
g _ _~~ II Soil he 20 Kilometers More than Type or less* 20 Kilometers
S1 or Site Specific Hazard ATC 3-06 Method S Analysis (Para 3-8) 2 ~(para 3-3 to 3-7) (aa38
S3 Same as above Site Specific Spectra Development (para 3-6). Site Specific Hazard Analysis not Required.
* If line fault or area source, then source must have maximum Mmax greater than 6.0, otherwise use Column II.
Eastern United States (EUS)
Soil I II Type Within a province All regions having MmaxŽ 6.0 other than in Column I
S or Site Specific Hazard ATC 3-06 Method 1 Analysis (para 3-8) S2 (para 3-3 to 3-7) (
S3 Same as above Same as above
US Army Corps of Engineers Figure 3-1. Selection procedure.
3-2 27 February 1986 TM 5-809-1O-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A scribed in appendix C, paragraph C-3. Of all these 3-2. Definition of Terms, Glossary, and methods, the empirical method consisting of an Symbols. PGA scaling factor for ground motion severity The methodologies of determining ground mo- at a given risk level, and an effective DAF spec- tion are based on the following disciplines: ge- tral shape, has been selected for the typical con- ology, seismology, dynamics and vibrations, ditions and design objectives of this manual. An probability and statistics. Because of this rather effective response spectrum will be specified for extensive range of subject matter, it is neces- each of the two levels of structural perform- sary to provide both symbols and a glossary of ance. Unless specified by the appropriate agency terms used in this manual along with the related the acceptable risk of exceedance will corre- terminology commonly used in the references spond to: and necessary bibliography. These are given in (1) A fifty percent risk of exceedance in fifty appendix A, Symbols and Notations; and in the years (EQ-I), and Glossary. (2) A ten percent risk of exceedance in one 3-3. General Overview of Seismic Hazard hundred years, (EQ-I1). Analysis and Specification of Ground Table 3-1 shows the relationship between the Motion. exposure time (or economic life of the facility), For engineering design and planning purposes, the probability of exceedance and the return pe- the future earthquake loadings at a site of in- riod. terest must be known. The procedures and steps
Table 3-1. Return period as a function of exposure time and probability of non-exceedance
Exposure Time Years 10 20 30 40 50 100 "Hazard" or Probability of exceeding A
5 195 390 585 780 975 1950 10 95 190 285 390 475 950 20 45 90 135 180 225 449 30 29 57 84 113 140 281 40 20 40 59 79 98 196 50 1 5 29 44 58 72 145 60 1 1 22 33 44 55 110 70 9 17 25 34 42 84 80 7 13 19 25 31 63 90 5 9 14 18 22 44 95 4 7 11 14 18 34 99 3 5 7 9 11 22 99..5 2 4 6 8 10 19
------L_
[IS Ainmy Corps of Liigi rseer±s 3-3 TM 5-809-10-1INAVFAC P-355.1/AFM 88-3, Chapter 13, SectIon A 27 February 1986 for estimating this future loading comes under parameters (such as magnitude, intensity, and the general category of seismic hazard analysis. spectra) to be employed are dependent upon the It should be recognized that there are two dif- type of information available to the analyst and ferent approaches: deterministic and probabi- the needs of the designer. The procedures and listic. Deterministic approaches do not take into the models selected depend on the type, quan- account the uncertainty in the size, the location, tity, and quality of information as well as the and the frequency of seismic events. Probabilis- goal of the analysis. The general procedures for tic approaches incorporate uncertainty in all the evaluating seismic ground motion in the West- above quantities. An overview of the procedures ern United States do not differ greatly from those for deterministic and probabilistic approaches in the Eastern United States. However, since is given in this paragraph. Steps are outlined by the tectonic setting and the available seismic means of flow diagrams and illustrative for- information varies greatly between those two mats. These are shown in figure 3-2 for the two geographic regions, the elements of the proce- main tectonic regions; the Western United States dures are different. A discussion related to se- (WUS) and the Eastern United States (EUS). lection of deterministic or probabilistic a. Algorithm ofBasicStepsofSeismicHazard procedures will be given in paragraph 3-3c. The Analysis. Various earthquake severity param- five basic steps required for the evaluation of a eters at the source and site are described in ap- site specific seismic ground motion are described pendix C, paragraph C-1. The particular below (see fig 3-3). The region-specific flow dia-
Seismic Hazard Analysis Procedure (See Figure 3-3)
I)1
WESTERN UNITED STATES EASTERN UNITED STATES (See Figures 3-4 and 3-5) (See Figures 3-6 and 3-7)
DETERMINISTIC PROBABILISTIC DETERMINISTIC PROBABILISTIC PROCEDURE PROCEDURE PROCEDURE PROCEDURE
US Anmy Corps of Engineers Figure 3-2. General flow diagram selection chart. 3-4 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Sectlon A
IDENTIFICATION AND MODELING OF Step 1 (para 3-4a and 3-4b) SEISMIC SOURCES
DEF1NE SIZE PROBABILISTIC BTAIN PARAMETER SOURCE _ / ~~~~~~~SEISMICITY (para 3-3d) DETERMLNISTIC INFORMATION APPROACH
SELECT LEVEL(S) ESTIMATE LARGEST OF EARTHQUAKE MAGNITUDE MAGNITUDE(s) POSSIBLE FOR Step 11 (para 3-4c ) J ~~~~~~~THESOURCE SELECT MOST CONSERVATIVE RECURRENCE DISTANCE FROM RELATIONSHIP I - SOURCE TO SITE
_. _ _SELECT FORECASTING ATTENUATE ONE OF THE SITE MODEL FOR Step III (para 3-4d ) SEVERITY PARAMETERS THE SOURCE FROM SOURCE TO SITE ,1.~- ATTENUATE THE DEVELOP SPECTRAL SELECTED SITE SHAPE FOR THE SEVERITY PARAMETER Step IV (para 3-5) SITE FROM SOURCE TO SITE. PROBABILISTIC GROUND MOTION INFORMATION F FIXED SEISMIC L DESIGN INPUT CRITERIA FOR THE SITE IDEVELOP EFFECTIVE RESPONSE . . | .~ - - - Step V SPECTRUM FOR GIVEN RISK (para 3-o ) LEVEL AND SITE CONDITION
US Army Corps of Engineers Figure 3-3. General flow chart. 3-5 TM 5-809-1Ol-1NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 grams and illustrations of related procedures formation from experts. The purpose of this step are shown in figures 3-4 and 3-5 for the WUS is to assemble the information required to de- and figures 3-6 and 3-7 for the EUS. Each figure lineate faults and regions within which seismic shows the parallel basic steps as required in the activity can be considered homogeneous. See deterministic and probabilistic procedures. paragraph 3-4b for a detailed discussion and ap- (1) Step I is to identify and model seismic pendix D for examples. sources. The selected type and accuracy of this (2) Step II is to define the size or severity modeling depends on the available geologic, geo- parameter of the seismic event at the source and tectonic, geomorphic, historic, and subjective in- the related recurrence relation. The size will be
DETERMINISTIC PROBABILISTIC
Selection of Earthquake Selection of Earthquake Data Base Data Base
- _F~~~ l Step Irldentify and Model I |Earthquake History and Identify and Mcdel Seismic Sources _ -oGeological Information Seismic Sources I
Determine Largest Earthquake I. Unit Sys Ami I. m- - ¶bL- Incompleteness in Records for Each Source Adjust Data Base Use. Earthquake History and Geological nformation Step II Determine Recurrence Earthquake History and Relationship for each IGeological Informatlon L Select Sround Ile ionio Seismic Source; m 9S max Attenuation Re-ationship J ) T Step 'I:Select Probabilistic Determine FGA at Site Model for Earthquake Occurenc due to Largest event on Each Source and Using I Shortest Distance IStep IV:Select Ground Motion _ Strong Motion Records in Attenuation Relationship Similar Tectonic Setting i
.se Largest PGA
as Design Level Determine Probability of IN IEarthquake Occurence Model I Exceedence of Different and Attenuation Relationship I PCA Levels; Fazard Curve U'se this Largest ?GA to . Scale the Appropriate Site i Response Spectrum Shape Determine PGA Level IType of Facility and Corresponding to Specified I ' AcptableRisk Probability of Exceedence
I
Step V:Select Appropriate - IRegional Attenuation Effects Response Spectrum Shapet *_ and ~Site _Soil C~onditions I Anchor at PCA Level
11SArmy Corps of Engineers Figure 3-4. Flow diagram for the Western United States. 3-6 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1I/AFM 88-3, Chapter 13, Section A
DETERMINISTIC FROBABILISTIC
Step I
Su:-NCESZ SOURCZANDS IvKCIlvE VISIAJICE
DE7ERY!.NE- YAXIPrJM YAG::'1CED A471 C'~n;t7: EC,. T Step II :-r'!"ACE FO' EAC' 11'.'4''
I.ECU"E.DC
UELRCZ PROBABILISTIC PORECASTIC MO0DEL III (PARA 3-Jid) Step
PGA on F,
Step IV VA on A I
I' K I I
ATIEKUAIOI AlTENUATI O4 (PARA 3-5c)
USE AflTIMIATIW4 1-0 LE-EEhKIN Km'AX!Y.UY -11-A AT SITE USE MAXINUM. n SPECThFWIPF,~ SIiTE
SITE HAZARD CURVE l i~ARA )-4d) S Step V I~~~-
SITt ALSPOr., 5PtCliltM SITE RESPONSE SPECTRUM (P'ARA3-6)
UIS Aru y Caorps of Engineers Figure 3-5. Hazard evaluation of WUS.
3-7 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986
DETERMININST IC PROBABILISTIC
Selection of Larthquake Selection of Earthquake |Data Base Data Base T - ~~~___
l Define Tectonic | S-ep I Define Tectonic |Earthquake History anc I Provinces Provinces and Model Tectonic Infcrm2t or ! Seismic Sources _,F _ ___ ~~~~~~~~~~~~~~. I Determine Largest Historical t . Eeent ir each Tectonic Provincej Adjust Data Base Unit System y V.~- -- rl I Incompleteness in Reccrds I |Se'eet Intensity |Azten ua'ion Relationship; Step :I:Determine Recurrence _ Earthq{ ke History anf Relatiorship for each |Tectonic Information I Seismic Source; ,p Imax Determine Site Intensity due to Largest Event in each Psovincee and Using Shortest Distance 'I~ste Ielect ProbabiliSt ic Mollfor Earthuake OccurencJ
Ilse Largest Intetisity as . . r Site Des~gn Intensity |Step IV:Select Intensity Iso Seismal Maps I ) Attenuation Relationship I Corr-late Site Iiesign Intensity i w1 PGth I Cotrelate Attenuated Intensity Intensity and Related I with PCA at Site Ground Motion Data
Use *trt! ILA IC- Scale the Appropriate Site Determine Probability of Earthquake Occurence Model Response Spectrum LEceedence of Different and Attenuation Relations'.Fp M1A Levels; Hlazard Curve _-~ Determine PGA Level IType of Facility andI Ac eris ta l Correspornding to Specified Probability of Exceedence
- Step VaSelect Appropriate Regional Attenuation Effects I Response Spectrum Shape: and Site Soil Conditions Anchor at PGA Level
US Army Corps of Engineers Figure 3-6. Flow diagram for the Eastern United States.
3-8 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A
DETERMI NI ST IC PROBABILISTIC
ECTONI C PROVINCE A,
Step I
SJOUCES SOURCEMODEL (PARA 3-'.b.(2))
DETEMbI!E MAXIM'jW I!r*!wSITY ANI. SITE 1S-AW;E PO.-, Step II
FOR A2
I I RECURIENCERELATIONS (PASIA 3-1.c.(I)(fI
SELEC!? PRC~hABIUS?IC Step III POAECAST'ING MODEL ('ARA 3-4id) (Slot Illustrat...)
ICA
KiA Step IV
AYEII'JAI1 IW ATTEHIIUTI UN4 (YAkp 3-5d)
USE AT'kWL'A'J1ONITL;
USEMAXIMUM TO SCALE APP'PMIAIE
SITE HAZARD CURVE IPARA 3-4d) S a 5 / C .(I4~~A i CIJAY) Step V KA a \ S,-(PCA 1)IDAP)
-- ? CITE RESPCAI:L TI J? SITE RESPONSE SPECTRUM (PARA 3-6)
US Army Corps of Engincers
Figure 3-7. Hazardevaluation of EUS.
3-9 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 one of the magnitude scales (ML, mb, MS) or celeration for a given magnitude event on the epicentral intensity I., or seismic moment, Mo. source could be used. The selection of the pa- The most commonly used size or severity param- rameter used for representing the severity and eter at the source is the Richter magnitude ML- the form of its attenuation relation depends on For the deterministic approach, the frequency the region where the analysis is performed and (or number per unit of time) of occurrences of the type of available data. See paragraph 3-5 various magnitude events need not be deter- for a detailed discussion and appendix D for ex- mined, and the assessment of ground motion at amples. a site will be governed only by the maximum (5) Step V is to represent the effects of dis- level of earthquake magnitude. For the proba- tance, local soil conditions, the magnitude of the bilistic approach,the parameters describing the seismic event, and the structural foundation size source seismicity must be obtained. This infor- and mass on the frequency content of the ground mation usually is in the form of a "recurrence motion. This is represented by the shape (DAF) a relationship," and an upper magnitude or in- of the effective response spectrum for the site tensity cut-off. The recurrence relationship pro- and its formulation is described in paragraphs vides information on magnitude or intensity and 3X6 and 3-7. The final specified spectrum is of the corresponding rate of occurrence or exceed- course scaled down by the forecasted site se- ence of that magnitude anywhere on the source verity. See paragraph 3-8g for examples. under consideration. The upper magnitude cu- b. Use of Results. This available informa- toff consists of the largest (maximum) possible tion on ground motion is utilized for design and/ event that the source can generate. The method or analysis of structures. Chapter 4 shows this of obtaining the above information depends on utilization for prescribed structural perform- the type of region and the data base available ance and selected risk levels. for the region. See appendix C, paragraph C-i c. Selection of Method. The deterministic for background, paragraph 3-4c for a detailed procedures as outlined in the flow diagrams are discussion and appendix D for examples. used exclusively for those important structures (3) Step III is to project the recurrence in- where the consequences of failure are cata- formation from regional information and past strophic; such as nuclear power plants, liquified data into forecasts concerning future occur- natural gas facilities, and dams. These proce- rence. This step is needed in the probabilistic dures tend to compound conservatism (cer- ) approach only. The forecasting model depends tainty of occurrence, largest magnitude and on the type and reliability of the data base. The closest distance from epicenter to the site) and most commonly used forecasting model is the will generally result in extremely large design Homogeneous Poisson probability model. Ho- requirements. For most structures, these highly mogeneous implies a memory-less occurrence of conservative design values cannot be justified events in time and location. When this homo- economically for use. This disadvantage of ex- geneity in time does not appear applicable, Semi- treme conservatism has actually resulted in the Markov and Markov chain models are used (see adoption of probabilistic procedures even for Patwardhen et al. (Biblio 50), Vagliente (Biblio some critical facilities. Deterministic proce- 68), Nishioka and Shah, (Biblio 45). These models dures, therefore, will not be discussed further allow inclusion of memory or time since last event in this manual. and are more involved and require substantially d The STASHA program. The purpose of this more information than the Poisson model. A manual is to provide the user with an over-all simple extension of the Homogeneous Poisson understanding of the procedures, assumptions, model, known as the Non-homogeneous Poisson and computational methods of ground motion model, may be adapted to incorporate time- hazard analysis. However, it is most important dependent information such as the rate of stress to recognize that any actual site hazard evalu- build-up and the time since last event, see Savy ation would require the use of the computer for and Shah (Biblio 52). Another model, usually a development of the various empirical relations uniform probability function, may be employed and the multiple calculations required for prob- to represent the random location of event oc- abilistic accuracy, and prediction uncertainties. currence on the source. See paragraph 3-4d for In order to perform these calculations in an or- a detailed discussion and appendix D for ex- derly manner for each step of the hazard anal- amples. ysis, the STASHA Program has been developed (4) Step IV involves the attenuation of the by the John A. Blume Earthquake Engineering severity parameter from its location on the source Center at Stanford University. Both the user's to the site. Either intensity or peak ground ac- manual and computer program tapes for 3-10 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A STASHA are available at the Corps of Engi- tained in the STASHA user's manual (Stanford neers Office. In the text of this manual, the University, Technical Report No. 36). A descrip- STASHA Program will be referenced whenever tion of STASHA and examples are given in ap- there is a need for extensive computational ef- pendix D. fort or for the representative examples con- Section 11. PROCEDURE FOR SITE SPECIFIC GROUND MOTION 3-4. Determination of Source Seismicity. However, since the future event could occur Each of the probabilistic hazard analysis pro- anywhere over the tectonic province and, there- cedures as presented in paragraph 3-3, and in fore, could be very near the site, the attenuation figures 3-4 to 3-7 is described in this paragraph distances (R) can therefore be short. Also, even and in the following paragraphs 3-5 to 3-8. though there are considerable variations in seis- a. Geotectonic and seismotectonic environ- mic severity patterns in the (EUS), these are ment. In the United States, two general re- not as well defined as in the (WUS). There is a gions are defined which are dependent upon the general smoothing effect over each entire tec- available geologic, geotectonic, geomorphic, his- tonic province and the boundaries between torical, and subjective expert information. It will provinces are often controversial. Also, the rel- be shown that each of the steps for seismic haz- atively low rate of seismic activity in the East ard analysis are region dependent. These re- makes the recurrence estimation over small areas gions are the Western United States (WUS) and very difficult. Further, because most Eastern the Central and Eastern United States (EUS). events have occurred in "pre-instrument" times, The boundary between these regions can be de- their source severity data are in terms of the fined by the eastern boundary of the Rocky more subjective value of intensity rather than Mountains, (Biblio 5). magnitude. Finally, the almost complete lack of (1) Regional Approaches. Due to the in- strong motion recordings makes the direct em- herent difference in the geologic structure in pirical development of attenuation relation- the two regions, two major approaches are used ships in terms of acceleration or velocity in defining seismic sources and assessing future impossible. However, both historical reports and seismic activity. In the Western United States seismological studies indicate significantly lower (WUS) and in many other parts of the world, rates of attenuation in the EUS. A summary of earthquakes occur on faults that extend to the regional differences is given in figure 3-8. surface of the earth. However, in intraplate re- b. Source modelling. Step I in seismic hazard gions, such as the Eastern United States (EUS), analysis is to identify and model seismic sources. this is not necessarily true, and it is difficult to This step depends on the following information recognize and delineate active faults. The two (see figs 3-9 and 3-12): major approaches are (see Biblio 17): -Type and amount of historic seismic oc- -Active Fault Approach currence data base. -Tectonic Province Approach -Geologic, geotectonic, and geomorphic data (2) Procedurefor each approach. The two base. regional approaches require different proce- dures for seismic hazard evaluation. In the ac- -Subjective opinions of experts concerning tive fault approach, seismic sources are relatively the seismicity of the region. well defined along plate boundaries or faults and, The process of source modelling provides two hence, the concentration of seismic events and essential portions of information for site hazard the resulting level of seismicity per unit length analysis: of the source or unit area of the source is rel- atively high. Also, because of the definite loca- -First, the configuration of the source and tion of the source, the source-to-site attenuation its size establishes the number and loca- distances (R) for the seismic severity parame- tion of seismic events for the evaluation ters are reasonably well defined. In the tectonic of source seismicity in paragraph 3-4c. province approach, the seismicity is diffused over a large area because no specific faults are iden- -Second, the configuration and location of tified. Each identified source area is assumed to sources relative to the site determines the have homogeneous (uniform) seismicity, and, attenuation distances (R) for ground mo- therefore, the seismicity per unit area is small. tion severity in paragraph 3-5. 3-11 TM 5-809-10-1INAVFAC P-355.1/AFM 88-3, Chapter 13, Sectlon A 27 February 1986 WESTERN UNITED STATES _WJS9
- Well defined sources
- Significant amounts of data in the fonn of historical repoits.
accelerograms, and geological creep measurements.
- Attenuation data in the form of records at different distance
and soil conditions.
- Relatively high occurence rates.
- High attenuation of ground motion severity
mainly within 100 kilometers.
EASTERN UNITED STATES (EUS) ) - Vague description of source provinces. I' - Some historical reports, and very few strong motion records.
- Relatively low occurence rates.
- Low attenuation of ground motion severity
with significant values at 200 to 300 kilometers.
US Army Corps of Engineers Figure 3-8. Regional differences.
It should be mentioned that currently, the USGS for source identification in each region are de- researchers are attempting to define seismic scribed as follows: source zones for five interior regions of the (1) Source modelling in the Western United United States, preparatory to the construction States. In this region (see fig 3-9), seismic of new national probabilistic ground motion sources are identified and modelled in the fol- seismic hazard maps. The five regions are the lowing ways: Great Basin, the Northern and Southern Rocky (a) Point source. This source charac- Mountains, the Central Interior and the North terizes a small region where repeated past Eastern United States (see Biblio 67). Since this earthquakes have occurred. However, no geo- work is not yet complete, this manual will de- logically identifiable fault exists. Typically, the velop procedures based on the two regions, the size of the region is small compared to the dis- WUS and the EUS. The particular approaches tance from this source to the site. Occasionally, 3-12 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Sectlon A
OBTAIN EPICENTRAL MAP OF REGION FROM DATA BASE, COE UP TO 1980 AND NOAA AFTER DEC. 1980
OBTAIN FAULT MAP OF REGION FROM USCS, STATE AND/OR NRC
INDEPENDENT ASSESSMENT OF THE POSSIBLE PRESENCE OF UNMAPPED FAULTS
MO)DEA, KNUIWN FAULTS BY LINEAR D1PI'lNG PLANE SOURCES TO MArCH TIIE EPICENTRAL AND FAULT MAPS
ASSIGN LTHE SOURCE, AND POINT SOURCE AS APPROPRIATE __Z11______.
ASSIGN AREA SOURCE IF EPICENTERS DO NOT MATCH THE KNOWN FAULTS AND THEN AFTER CONSIDERING GEOLOGY
ASSIGN BACKGROUND SEISMIC SOURCE FOR ALL UNASSIGNED EVENTS (AREA SOURCE)
US Army Corps of Engineers Figure 3-9. Flow chart for step I source identificationand modelling for the WUS.
3-13 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 volcanic sources can be identified as point planes and line sources are used: (Biblio 41), sources. and (Biblio 70). An example to demonstrate as (b) Line sources. Fault traces are taken to how sources are modelled is given in appendix as lines at a certain fixed depth below the ground D. surface. In California, this depth is usually be- (2) Source modellingin the Eastern United tween 5 to 35 kilometers. This "active fault" States. In the Eastern United States, the tec- modelling approach is used wherever the tec- tonic province approach is used (see fig 3-12). tonic structure is more or less evident at the There are various reasons for adopting such an surface. approach; the most important being that the (c) Area sources. This source model is degree of fault and seismic activity in the East- used when the occurrence of earthquakes in a ern United States is low, resulting in very little region cannot be correlated with known faults geologic and historic evidence. Also, in large areas or the geologic structure of the region. There of the Eastern United States, there is a scarcity are also cases where the number of small faults, of geologically recent deposits that would re- or a source of clustered activity, may be consid- cord evidence of recent fault activity. In addi- ered together as an area source. tion, the heavy vegetation covers the faults and prevents their detection. Finally, the recent de- (d) Dippingplane. This source model is velopments of evaluating fault activity in the used when one geologic plate thrusts under an- (WUS) have not been applied in the east due to other plate so as to create a distributed source excessive cost and time involvement; except in of earthquakes. This feature is called a Benioff a few regions such as New Madrid where fault- Zone, and can be modelled by means of dipping ing evidence has been substantiated (Biblio 71). planes upon which earthquakes have variable epicentral depths. Geological conditions such as (a) Area source configuration. One of this occur in Alaska and in Central America. the key features of tectonic province approach is to delineate these provinces as area sources (e) Background area source. In gen- that have a uniform potential to generate earth- eral, events that occur somewhat randomly quakes. Within that area, the future earthquake throughout the region and that cannot be as- activity should be homogeneous. Due to lack of sociated to any fault or source are treated as sufficient historical and geological evidence, background seismicity. They are considered to there is no unique and generally consistent way be part of a large area source with uniformly of delineating these area sources. Two examples low seismicity that extends over the area not on area source configurations for the Eastern covered by the other sources. The earthquake United States are shown in figures 3-13 and 3- location, if not included within one of the pre- 14. viously defined sources, is then in the back- ground zone to account for the possible (b) Using subjective input as furnished occurrence of the random or "floating" earth- by interviews from ten experts, Mortgat (Biblio quake. The effect of the background zone is gen- 63 and 64), has developed homogeneous area erally small since the contribution of the other sources as shown in figure 3-15. With respect to sources are governing the hazard. In some par- this method of using expert opinion, it is well ticular cases however, where the hazard is low, to recognize that experts form their objective the background contribution may be non neg- biases from the particular data and other geo- ligible. logic and seismologic evidence that they may have seen. Since most of the experts work with a sim- (f) Western source conditions. The point, ilar data and information base, the variability line, and area source models are shown in figure in their individual source configuration is due 3-10, and the dipping plane model in figure 3- to their personal biases. Barstow et al (Biblio 11. In source modelling, historical records and 5) have studied statistical techniques to provide the knowledge of geotectonic features of the a methodology for the production of working region play an important role. Due to the high tectonic province and tectonic structure maps seismic activity in the Western United States, for the Eastern and Central United States, iden- and the relatively good geological evidence of tifying areas of uniform seismic hazard. faults, surface rupture, and other tectonic fea- tures, line sources are used most extensively. c. Source seismicity. Step II in seismic haz- Area sources are common in the Pacific North- ard analysis is to evaluate the seismicity of each ' west. In regions such as Alaska, both dipping of the modelled source (see fig 3-16). Evalua-
3-14 27 February 1986 TM 5.-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A
Fault Length L1 IL II_I Line Source
;1
~Site R2 -1 '
I-, Ri N-1 \ ~~~~~~~~~I
\ R3
Line Source \ 6 Point Source
US Army Corps of Engineers
Figure 3-10. Point, line and area sources.
3-15 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 L'w
Ruptured Zone
Latitude distance
Deep Boundary
Dipping A Shallo Planes SiSt te- .0/Boundary
Longitude distance
Ground
US Army Corps of Engineers
Figure 3-11. Dipping plane source.
3-16 27 February 1986 TM 5-809-1O-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A
OBTAIN EPICENTRAL MAP OF REGION FROM DATA BASE
HQ-USACE* UP TO 1980 AND NOAA** AFTER DECEMBER 1980 i IDENTIFY THE TECTONIC PROVINCES AND THE RESULTING AREA SOURCES, INCLUDE FAULT MAP IF AVAILABLE FROM USGS, STATE AND/OR NRC. OBTAIN EXPERT
OPINION ON SOURCE LOCATIONS
ASSIGN BACKGROUND SEISMICITY
FOR EVENTS THAT
CANNOT BE ASSIGNED A SPECIFIC AREA SOURCE
* DAEN-ECE-D Washington, D.C. 20314
** NOAA/NGSDC/TGB 325 Broadway, mail code D-623 Boulder, CO 80303 NOTE: If at a future date, specific faults are identified, then they can be modeled by means of line or dipping plane sources.
US Army Corps of Engineers
Figure 3-12. Flow chart for step I source identification and modelling for the EUS.
3-17 TM 5-809-10-1/NAVFAC P-355.1I/AFM 88-3, Chapter 13, Section A 27 February 1986
I I I
)
RE-printed from "Effects of Uncertainty in Seismicity on Estimates of Seismic Hazard for the East of the United States," McGuire, R. K., Bulletin of the Seismological Society of America, Vol. 67, No. 3, 1977, with permission from the Seismological Society of America.
Figure 3-13. Seismic sources after Algermisson and Perkins (1976).
3-18 13, Section A 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter
~~~ -
Reprinted from "Effects of Uncertainty in Seismicity on Estimates of Seismic Hazard for the East of the United States," McGuire, R. K., Bulletin of the Seismological Society of America, Vol. 67, No. 3, 1977, with permission from the Seismological Society of America. Figure 314. Seismic sources after Hadley and Devne (1974).
Ki'
3-19 . . ,r. ILLINOIS INDIAVIM: _ 1/I t~r~vua /? i I iO
MIS&WAI H~~~~~~~~~~~~~~VIAGINIA ~~UT .
. @X~~~/ it- '" ";W /g n
A7-~~~~~~~~~~~~~~~~~-
L. . , y k: UA "'..# KZ5f~ p I "' L A'\ &WAAS _ _ Io T.__t " ,&&,.,.& .. g \ ( 3
0 J ~ ~ ~~~~ ~~~~~~~~~~~~~~~~ Le
Reprinted from "Seismic Hazard Mnalysis- a Solicitation of Expert Opinion," Nuclear a Regulatory Commission, Tera Corporation, NURFG/CR-1582, Vol. 3, 1980. Figure 3-15. Seismic source after Tera (1981). Z 0. K C~~~~~~~~~~~~~~~~~~~O 27 February 1986 TM 5-809-1 0-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A
FOR A GIVEN SOURCE AND CORRESPONDING HISTORICAL PROVIDE INPUT FROM EVEtITS, PLOT MAGNITUDE EXPERT OPINION INCLUDING GEOLOGICAL OR INTENSITY VS I FREQUENCY OF INPUT OCCURRENCE I
I - I
SELECT AN ANALYTICAL FROM GEOLOGICAL EXPERT FORM, COMMlONLY THE LOG-LINEAR FORM FOR THE RECURRENCE OPINION, ASSIGN RELAT!ONSHIP AND FIT IT I TO THE ABOVE MAXIM4UM MAGNITUDE INFORMAT ION
OBTAIN NORMALIZED
~ECURRENCE RELATIONSHIP
FOR THE SOURCE I REPEAT THE ABOVE
STEPS FOR ALL SOURCES
OF INTEREST
US Army Corps of Engineers Figure 3-16. Flow chart for step If source seismicity and recurrence relationshipfor WUS and EUS.
3-21 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 tion of seismicity involves the following com- The main problem with this approach is that it ponents: reduces the size of the useful sample and mean- ingful statistical averages of large earthquakes -Collection and processing of occurrence cannot be obtained because of their infrequent ) data and formulation of the recurrence occurrences (Biblio 6). An alternative is to cor- relationship. rect for incomplete reporting by a random sim- -Determination of the size of the maximum ulation of missing data (see STASHA). The earthquake a given source is capable of Gutenberg-Richter relationship is given by generating. equation (3-1). (1) Collection of data and formulation of In N(m) = a + f3m (eq 3-1) the recurrencerelationship. The data base for where seismic events on a given source is often incom- In = Natural log to the base e plete, nonhomogeneous in time, and lacking in refinement. The appropriate processing of this N(m) = Average Number of events occurrence information is very important be- greater than or equal to the cause the reliability of results of the hazard magnitude m. analysis are strongly dependent on the consist- a, 1 = constants. ency and the completeness of the input data base. The magnitude-frequency or recurrence rela- Very often, this relationship is used in a slightly tionship is formulated from the number of different format where logarithm to the base 10 earthquakes that a source has generated and is used instead of to the base e. their respective magnitudes. The most common logioN(m) = a + bm (eq 3-2) method of determining this relationship is from historic data. Occasionally, other information One would convert the equation from base e to sources, such as geological evidence and slip rate base 10 by means of the following simple con- of the fault, are used to supplement this histor- version: ical data base. Statistical regression analysis is a = 0.43429a (eq 3-3) commonly used to obtain the best line fit with the "least squared" error. Expert subjective b = 0.43429P (eq 3-4) ) opinion can also be incorporated in order to sup- Such magnitude-frequency relationships are plement the historical data base. The most com- called "recurrence relationships" in the litera- monly used magnitude-frequency relationship ture and a general example is shown in figure is the one suggested by Gutenberg and Richter 3-17. After the recurrence relationship is ob- (Biblio 26). In this relationship, the source se- tained, the following normalization process can verity parameter could be either magnitude or be performed. epicentral intensity. The type of parameter and (a) Normalization to unit length and the constants of the magnitude-frequency re- time. Let T be the time-period over which the lationship vary from one region to the other. recurrence data has been obtained. If the source Data adjustment is usually necessary before us- is a line source, let L be the length of this source. ing the data to determine the parameters of the Then, N(m) = average number of events equal magnitude-frequency relationship. It has been to or greater than magnitude m during the time observed that the completeness of earthquake period T and on source length L for the line records varies with time. In the past, due to low source. population density and lack of interest in earth- quake activity, only large events were recorded. Let With increased instrumental coverage, inter- N0(mW = N(m) mediate and lesser earthquakes have been re- N'() - LT corded with more frequency, producing an apparent increase in seismic activity with time then which biases the statistics from uncorrected ln(N'(m)) = In N(m) - InN(m) -ln(LT) catalogs of data. In recognition of this time bias, LT the evaluation of the degree of completeness of ln(N'(m)) = a + 1 m - In(LT) the available earthquake record is an important =a - In(LT) + 1 m step in the analysis of data. One possibility is to confine analysis to sections of the record that or are complete for the earthquakes of interest. ln(N'(m)) = a' + 1 m (eq 3-5) 3-22 27 February 1986 TM 5-809-1-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A where or N'(m) average number of events equal ln(N'(m)) = a' + ( m, to or greater than magnitude m with a' = a - ln(AT) per unit time and unit of source (eq 3-6) length Where N'(m) and a' are now normalized with a' = a - ln(LT) respect to the source area A. Note that the value of a does not change when (c) Sampling uncertainty. For a given the recurrence relationship is normalized. This magnitude m, the fitted line gives the average step of normalizing the recurrence relationship value of N(m) or N'(m), and this average or is usually done by the seismic hazard analysis expected rate value is required for the proba- computer program. The purpose of presenting bilistic forecasting model in paragraph 3-4d. this step here is to indicate that in the normal- However, there is considerable scatter of the ization, it is assumed that for a given source, the actual recorded number of events. To take this number of events equal to or greater than a scatter into account, a probability distribution given magnitude is homogeneous in time and function is generally assumed for the number space. Thus, the mean rate of occurrence does of events equal to or greater than a given m. not change with time or along the given source. Further, the fitted recurrence line, because of More will be discussed on this topic when the limited data base and the largely subjective eval- probabilistic-forecasting models are presented. uation of the maximum magnitude, has a sam- (b) Normalization to unit area and pling error. This sampling error is an indicator time. If the area source with area A was con- of the difference between the sample fitted line sidered instead of the line source, the relation- from the limited data source and the true line ship would have a simlar format: that would be obtained from a very large data source, figure 3-17. The STASHA, (Stanford N'(m) N(m)T University, Technical Report No. 36) program gives a probabilistic representation for this ln(N'(m)) = a - ln(AT) + p m sampling uncertainty in the N(m) value.
lnN (m)
L \ Fitted 'line -lnN (m) = a + Sm
.. Ic bound for sampling \-. X. error of fitted line
Data IN
m
US Army Corps of Engineers Figure 3-17. Linear Recurrence Relationship. 3-23 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 (d) Non-linear relations. Other forms of the recurrence constant for most of the WUS recurrence relationships have been used by re- sources lies between about 1.1 to approximately searchers. Dalal (Biblio 21) has used Gaussian 2.5. Figure 3-19 shows the recurrence relation- and log-Normal probability distribution models. ship for the northern section of the San Andreas Mortgat et al (Biblio 42) have used a bilinear Fault in California. It should be mentioned here relationship as shown in figure 3-18. Here, two that one large fault such as the San Andreas lines are fitted to the data. The point where the may be broken down into two or more homo- two lines meet is usually determined subjec- geneous segmental sources and the recurrence tively from the geologic considerations concern- relationship may be determined for each of these ing capabilities or rates of large magnitudes on segmental sources. This use of homogeneous the source. Cornell and Merz (Biblio 39) have segments is quite common in California where used a quadratic form for their recurrence re- there is evidence of varying degrees of seismic- lationship. Recently, Dong et al. (Biblio 23) have ity on the large sources. The source severity pa- applied the maximum entropy concept to obtain rameter employed in developing these recurrence minimally biased recurrence relationships. (See relationships in the WUS is usually the Richter app D for some examples). magnitude (which can be considered to be the (e) Recurrence relationship for sources same as the local magnitude ML)- In appendix in the Western United States. The "active fault" C, paragraph C-1, these variations magnitude approach is usually employed in this region. scales are defined. Therefore, based on the fault locations and the (f) Recurrence relationshipsfor sources modelling of these faults as line sources, past in Eastern United States. The tectonic prov- seismic events are assigned according to their ince approach is used for modelling sources in relative proximity to the different sources. This the eastern United States. Therefore, all the process of event assignment is usually per- sources are area sources, and these usually cover formed by expert judgement with recognition rather large regions. With respect to the source that epicentral locations are subject to error and severity parameter, most of the historical data that events are more likely to occur on the known in the East is compiled in the form of the Mod- fault rather than on the adjacent area. The ified Mercalli Intensity Scale. However, there are STASHA program has a procedure for event as- cases where the most recent data is in local Mag- ) signment. It has been found that the value of nitude (ML) or body wave magnitude (mb).
lnN(m)
nN(m) = a1 + 8
reak Point
7- 1 nN(m) = a2 + BI
Maximum Possible Magnitude for the Source (l m maM Imax
US Army Corps of Engineers
Figure 3-18. Bilinear Recurrence Relationship. 3-24 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A
1000.0
'100.0010~~~~~~E k \
lnN(H) 8.59 1.18M
10.0 _ N
3 < s s~~~~~~~~~~~~~~~~~~~~~~~1 z N.,~~~~~~~~~
¢1.0 _ \ \ s N? _
3.00 4.00 S.00 6.00 7.00 8.00
Richter Magnitude (M)
US Army Corps of Engineers
Figure 3-19. Recurrence relation for North San Andreas TM 5-809-10-1/NAYFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 Therefore, in order to make all the data con- the term "maximum credible" event. Such a term sistent, one approach consists in converting the should be discouraged from use. Instead, use of magnitude into an epicentral Modified Mercalli a term such as the "maximum earthquake size' Intensity (MMI). should be encouraged. The size of the maximum'n" ) Let 1o be the epicentral MMI earthquake is used in source seismicity deter- mination in two ways: M be the Richter magnitude. -Deterministic use of the maximum earth- In paragraph C-1, appendix C, the relationships quake in the design process (see figures between these two parameters for the Eastern 3-4 and 3-6). United States are given. -Probabilistic use of the maximum earth- Also, Nuttli (Biblio 46) has developed a rela- quake, in the recurrence relationship. Here, tionship between the body wave magnitude mb the value of this earthquake size provides and the epicentral intensity I., the upper cut off magnitude in linear re- currence relationship, or it could be an mb = 1.75 + 0.50Io (eq 3-7) asymptote in the non-linear recurrence Using relationships such as these, the occur- relationship, see figure 3-18. rence data in the form of the source intensity 1. can be obtai- ±d. The Gutenberg-Richter re- The estimate of the size of the maximum earth- quake for a given source is based on the follow- currence formal ',r intensity is then written in the following form: ing factors: 1-Geologic evaluation of the regional tec- ln[N(I )] = aj + 03iI0 0 (eq 3-8) tonic framework. Where atcand pi are regression constants. Ye- gian (Biblio 71) and TERA (Biblio 63,64) have 2-Historical seismicity of the source and given values of pi, for the EUS. A shortcoming the surrounding region. of using epicentral intensity (I.) as a parameter 3-Geologic history of displacement (from is that I, unlike magnitude, is not a direct meas- trenching investigations). ure of a source severity. By definition, intensity ) is a number corresponding to particular ob- 4-Relationship between earthquake mag-__ served effects and these are often influenced by nitude and fault rupture length. both the site condition and the prevailing local 5-Relationship between earthquake mag- types of construction. In order to overcome this nitude and amount of fault displace- shortcoming, an alternative approach involves ment. the estimation of magnitude of the historical events (before instrument records) in terms of Out of the five factors mentioned above, the tec- their estimated epicentral intensity, felt area, tonic province approach in the EUS would per- and fall-off intensity. This requires a large mit the use of only the first three. When the amount of background research effort. How- active fault approach is employed in the WUS, ever, most large events in the EUS have been then all of the five factors will be used for such assigned a magnitude based on this method by an evaluation. Whether one decides to use a spe- different researchers (Nuttli, et al. (Biblio 47) ). cific maximum earthquake value or a probabi- Smaller events of less importance in the anal- listic distribution representation of the maximum ysis can be converted to magnitude using one of earthquake value, the STASHA program can the relationships in appendix C-1, or equation handle both forms of this input information. 3-7. In the formulation of the recurrence rela- (a) Determination of the Size of the Max- tion in the EUS, it is usually assumed (because imum Earthquake-Western United States. In this of lack of data) that the same pi value applies region, seismic sources are usually line sources throughout very large regions and that local (active fault approach). For such sources the variations apply only to the level of seismicity maximum earthquake size is usually based on (parameter al). The range of values for the pa- the fault rupture length or the maximum amount rameter PI' is from 0.80 to 0.92. of displacement that may be associated with the (2) Determination of the maximum earth- causative fault. Not only the historical data base quake. One of the most controversial and im- is used, but also geological data from trenching, portant variables of interest in representing or other geomorphological studies; Sieh, (Biblio "- ) source seismicity is that of the size of the max- 60) can be employed. Recently, (Aki, (Biblio 1); imum earthquake. Past literature has employed Kanamori and Geller, (Biblio 24); Molnar, (Bib- 3-26 27 February 1986 TM 5-809-10-1/tNAVFAC P-355.1/AFM 88-3, Chapter 13, Section A lio 40)) seismic moment has been related to the sponding event can be considered as an upper fault rupture area, along with the fault shear bound. This last approach should include all modulus and average slip. Relating the maxi- available information such as local or regional mum seismic moment Mo.max to moment mag- strain release or stress field data (See para nitude Mm gives the value of the largest moment 3-4c(3)). magnitude. It is useful to note that Mm is equal (3) Use of Seismic Moment to Represent to ML for ML values between 5 and 7. Empirical Source Seismicity. One of the more recent de- relationships between M, fault rupture length velopments in seismic hazard analysis is to use L and fault displacement D are developed from seismic moment (M,) to describe source seis- world wide data (Bonilla and Buchanan, (Biblio micity. Seismologists have introduced a "phys- 11); Slemmons, (Biblio 61)). Paragraph C-4, ap- ical" parameter called seismic moment M, to pendix C gives these relationships. The tables describe size of an earthquake. This develop- and relationships presented in paragraph C-4 ment is relatively new and its practical imple- should not be used exclusively but together with mentation for seismic hazard analysis has not historical and other geologic evidence. The his- been achieved. Paragraph C-4, appendix C, in- torical record of earthquakes in a given region troduces the users of this manual to this new may be one of the few indicators of the potential concept. for future earthquakes. However, extreme cau- d. ProbabilisticForecastingModels. Step III tion must be exercised when extrapolated fore- is to forecast source severity of future earth- casts are made. The time period of records in quakes on each of the identified sources (see fig the United States is relatively short and there- 3-20), once the sources of seismic activity have fore statistical prediction should always be com- been identified (para 3-4b) and the seismicity pared or modified by expert judgment concerning of the identified sources has been determined seismicity. In paragraph C-4, table C-11 shows (para 3-4c). These forecasting models are not the slip rate activity of some of the faults of the based on extrapolation of past data, but are based Western United States, and figure C-10 shows on stochastic models. These models from 'the fault slip versus time. This type of information probability theory field of stochastic processes can also be incorporated probabilistically in as- may however employ data for the evaluation of sessing fault activity and in estimating the size their parameters. The type of stochastic fore- of maximum earthquake events. This will be dis- casting model selected depends on the accept- cussed further in the forecasting paragraph, able type and level of assumptions about the 3-4d. seismic occurrence on each of the sources. The (b) Determination of the Size of Maxi- most widely used model is called the homoge- mum Earthquake-Eastern United States. In this neous Poisson Model. Typical examples of this region, seismic sources are modelled by the tec- approach are given in the following references: tonic province approach. The most commonly Cornell (Biblio 18), Cornell and Van Marcke used method of determining the size of the max- (Biblio 19), Stepp (Biblio 62), Algermissen (Bib- imum earthquake is through historical records. ho 3), McGuire (Biblio 37), Shah et al. (Biblio Very little information (if any) is available on 58), Wiggins (Biblio 69), Der Kiureghian and the fault rupture or fault displacement and hence Ang (Biblio 22), Liu and Fagel (Biblio 34), Kir- these two parameters cannot be related to the emidjian and Shah (Biblio 32). This is normally size. To overcome the problem of limited histor- called a memoryless process because of the as- ical data in estimating the maximum earth- sumption that the probability of occurrence or quake size, the opinions of experts should be nonoccurrence of an earthquake in any given obtained. Two principal methods are used to de- year and for a given source does not depend on termine the maximum earthquake size. The first the time interval since the last occurrence. For one consists of using the size of the largest his- most practical cases where the future time ho- torical event subjectively incremented by a safety rizon is of the order of fifty to one hundred years, factor such as half a magnitude or one intensity this is a reasonable assumption and is suitable unit. The other consists in using the earthquake for the purposes of this manual. A non-homo- size corresponding to a 1000 to 5000 year return geneus Poisson model has also been used to ac- period from the recurrence relationship. Al- count for the dependence of the mean rate of though this last method is somewhat ad hoc it occurrence on time. Savy and Shah 1981, (Biblio is felt that, in the present geologic framework, 52) have shown the use of this model. In order the near future will be similar to the past and to account for the lack of sufficient historic oc- that the 1000 to 5000 year choice represents a currence data and also to take into account geo- low enough probability such that the corre- logical data (such as slip rate, size of past rupture 3-27 TM 809-10IM"NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986
FROM NORMALIZED RECURRENCE RELATIONSIP,
OBTAIN MEAN RATE OF OCCURRENCE FOR
MAGNITUDE OR INTENSITY OF INTEREST
SELECT STOCHASTIC FORECASTING MODEL COMPATIBLE WITH THE GEOLOGICAL AND SEISMOLOGICAL INFORMATION
o Homogeneous Poisson Model (widely used) o Non Homogeneous Poisson Model o Bayesian models
MODIFY THE MEAN RATE OF OCCURRENCE IF. GEOLOGICAL INFORMATION AND/OR EXPERT) SUBJECTIVE OPINION-CAN SIGNIFICANTLY CHANCE THE STATISTICAL ESTIMATE OF THE RATE OF OCCURENCE e i
OBTAIN PROBABILITY OF OCCURRENCE OF DIFFERENT MAGNITUDES OR INTENSITIES FOR FUTURE TIME PERIOD T AND FOR TOTAL SOURCE DIMENSIONS
I X -T~~~~~~~ REPEAT PREVIOUS STEPS FOR ALL THE SOURCES
US Army Corps of Engineers Figure3-20. Flow chart for step III seismic forecasting model.
' _ /)
3-28 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A length and the amount of fault displacement per and event), Bayesian models have been developed. These models assume a Poisson occurrence model x is the mean rate of occurrence. along with a Bernoulli model for the size of each If one defines the return period (TR) as the time occurrence. The STASHA program describes this interval during which the expected number of type of model; see appendix D for an example. occurrences is one, then this much used engi- (1) When the occurrence of a future event neering parameter in risk analysis is obtained is independent of the past occurrences, then the as follows: the expected number of events for homogeneous Poisson model is a reasonable the Poisson process of equation 3-9 is given by model. The Poisson model of occurrence can be written as E(N(t)lk) = Xt (eq 3-11) where E(N(t)l) = Expected number of events t ( t)n (eq PN(nt) = e 3-9) for future time t given X. where PN(n,t) =Probability of having n events If equation 3-11 is equated to one, we get the in a future time period t definition of return period. n =number of events KTR= 1 1 and hence TR-= (eq 3-12) X = mean rate of events per unit of time (years) TB is therefore the average time interval be- (2) If Xis independent of time, then the pro- tween events, and is also the reciprocal of the cess is called homogeneous. If Xvaries with time, annual risk of occurrence. The value of Xis usu- the process is called non-homogeneous. ally obtained from the recurrence relationship (3) For earthquake events to follow the ho- developed in paragraph 3-4c. Let N'(m) = a' + mogeneous poisson model, the following as- p m be the average number or rate of events sumptions must be valid: equal to or greater than magnitude m per unit of time and per unit of source dimension. Then, -Earthquakes are spatially independent; using the Poisson occurrence model, the prob- -Earthquakes are temporarily indepen- ability of n events equal to or greater than mag- dent; nitude m in future time t for source of length L (or area A) is given by -The probability that two seismic events exp(-N'(m)Lt )n (N'(zn)Lt)" will take place at the same place and at P(n,m,t) = the same instant of time approaches zero. (eq 3-13) The first assumption implies that occurrence or Thus, nonoccurrence of a seismic event at one site or location or source does not affect the occur- P(O,m,t) = exp(-N'(m)Lt) (eq 3-14) rence or nonoccurrence of another seismic event or probability of at least one event above mag- at some other location or source. The second nitude m for a source of length L in future time assumption implies that the seismic events do t is given by not have memory in time. The third assumption implies that for a small time interval dt, no more 1 - P(O,m,t) = 1 - exp(-N'(m)Lt) than one seismic event can occur. This assump- (eq 3-15) tion is considered to be realistic and fits the Equation 3-15 provides the most elementary physical phenomenon reasonably well. hazard statement for the occurrence of a given (4) It can be shown that if the arrival of magnitude (or greater) on a given source. The earthquake events follow the Poisson process, probability of exceeding a given level of site in- then the random description of the time interval tensity (such as PGA) needs consideration of between two events follows exponential distri- the location of the event (epicenter or rupture bution. Thus, length) on the source and also the consideration f(t) = ke"' t ¢ 0 (eq 3-10) of all sources affecting the site. This is treated in the next paragraph 3-5. = 0, Otherwise f(t) is the probability distribution func- tion for the interarrival time t between events, 3-29 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 3-5 Selection of the attenuation relation ponent of interest and the distance from the for the determination of seismic source. These types of developments are avail- severity at a site. able for body waves (See Savy, (Biblio 53)). ) Step IV of the seismic hazard analysis deals with However, the most commonly used methods for the methods of evaluating the severity of ground ground motion estimation in engineering and motion at the site where the structure is located, for seismic hazard and risk analysis are the ones given the information developed in the previous based on empirical relationships. In this man- three steps. ual, a short description of these empirical tech- a. Attenuation of ground motion. When a niques will be presented. For a detailed study rupture along a fault plane occurs, vibratory see Idriss (Biblio 30) or the OASES study by ground motions are generated. These motions Woodward-Clyde Consultants (Biblio 70). It is travel out from the source as body and surface commonly accepted by seismologists and geo- waves (See fig C-2). As these waves travel far- physicists that the type, the amount, and the ther out from the source, they are attenuated. geometry of the rupture surface influences the The type and amount of attenuation depends on amplitude and frequency of motion near the many factors, the most important of which are source. Other factors influencing the near-source listed below: motion characteristics are the velocity of rup- ture, the stress drop, the physical properties of -Size or source severity of the event on the the fault plane material, and the pattern of non- source uniformity of rupture on the rupture surface. -Type of fault mechanism The larger the rupture surface, the greater the ground motion. However, there are definite up- -Transmission path of the seismic waves per limits for both the rupture size and the re- from source to the site sulting motion. The wave patterns generated at -Vibration or wave frequency of interest the source travel out in all directions in the form of the seismic ground motion of complex wave forms. The regions through which these wave forms travel from source to -Distance from the source to the site site constitute the "transmission path." It has -Local site soil response effect been observed that the transmission path influ- ences the attenuation of wave forms in both the frequency and amplitude domains. The decaying In estimating the type and severity of ground of amplitude with distance is usually referred motion that would exist at a site due to some to as the "attenuation." In the frequency do- future seismic event, the analyst should incor- main, higher frequency components in the wave porate the above parameters in his model. The form get filtered out as the distance from the current state-of-the-art methods for estimating source to site increases. In this paragraph, only the ground motion can be classified into two the amplitude attenuation will be discussed. groups. Paragraph 3-6 considers the aspects of fre- -Methods based on wave propagation the- quency attenuation and its influence on the re- ories through elastic and non-elastic sponse spectrum shape. media with appropriate damping charac- b. Empirical attenuation relations. Various teristics. empirical relationships are available in the lit- erature to describe the relationship between the -Empirical methods based on past data. size of the event, the distance from the source In the first method, various researchers in re- and the site ground motion parameter of inter- cent years have developed models to study dis- est (see fig 3-21). In working with these rela- placement (or some other ground motion tionships, the question of distance from the parameter) wave forms as a function of the type source to the site arises. The most "realistic" of event and the distance from the source. In distance to be selected could be either the epi- particular, the models for estimating the sur- central distance, hypocentral distance, distance face wave patterns have been quite good and fit from the site to the energy release center, or the the data well (See Boore, (Biblio 12); Frazier, distance from site to the closest rupture loca- (Biblio 4); McCann, (Biblio 35)). There are some tion on the fault. Earlier relationships have used other models which look at the attenuation of epicentral distance; however, with the availa- Fourier spectra with distance. Such models take bility of more data in recent years, it has becomeK ' / into account the damping characteristics of the evident that this distance is not the most rele- transmission media, the wave frequency com- vant. Some studies have used hypocentral dis- 3-30 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A
WLI S EUS REGION OF INTEREST
-
USE PGA AS SITE USUALLY INTENSITY SEVERITY PARAMETER DATA IS USED FOR ATTENUATION EQUATION - ___ __- I CONVERT ATTENUATED SELECT TRANS. PATH INTENSITY FOR THE SITE TO DESIRED GROUND 1- MOTION PARAMETER SELECT THE APPROPRIATE I ATTENUATION RELATIONSHIP DISTANCE WEIGHTING IF DESIRED I _F MAGNITUDE WEIGHTING IF DESIRED INCORPORATE UNCERTAINTY INFORMATION -~~~~I- I SOIL EFFECT, IF INFOR- MATION IS AVAILABLE
I V USE 50 PERCENTILE VALUE INCORPORATE OF PGA FOR SITE L UNCERTAINTY INFORMATION SEVERITY SCALING
USE 50 PERCENTILE VALUE
OF PGA FOR SITE SEVERITY SCALING
US Army Corps of Engineers Figure 3-21. Step IV, attenuation of ground motion from source to site.
tance. The recent relationships use the concept some of the distance definitions used in the of significant distance. This is the shortest dis- literature. tance to the ruptured source. Figure 3-22 shows 3-31 TM 5-809-10--1/NAVFAC P-355.1/AFM 88-3, Chapter 13, SectIon A 27 February 1986 Rf
RC - Distance to energy center Be - Epicentral distance R- Distance to causative fault Rat - Hypocentral distance m- Map distance to energy center as- Significant distance
Reprinted from "Offshore Alaska Seismic Exposure Study (OASES)." 1978, with permission from Woodward-Clyde Con- sultants. Figure 3-22. Attenuation distances.
(1) Recent studies have indicated (OASES, the Western United States. For most earth- Biblio (70) that the transmission path B is very quakes in California and Hawaii, transmission important. Thus, for shallow earthquakes path A should be assumed. Also, there are im- (transmission path A in fig 3-23) there is one portant differences in rates of attenuation for attenuation relationship; whereas for deeper the WUS and EUS regions. These will be dis- earthquakes, (transmission path B in fig 3-23) cussed in the paragraphs for these regions. there is a separate attenuation relationship. This (2) Many empirical attenuation relation-' transmission path dependence has been ob- ships are available in the literature. They all served in data collected in Alaska, Japan and in have their shortcomings in both accuracy and 3-32 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A
Z
0
W-J W
U 'C
I
3 10 30 100 300 1000 DISTANCE EKMI
Reprinted from "Offshore Alaska Seismic Exposure Study (OASES)," 1978, with permission from Woodward-Clyde Con- sultants.
Figure 3-23. OASES attenuation. applicability for a given site. The scatter of data -OASES Model (Biblio 70) with respect to the estimated relationships is Figure 3-24 shows the first two of these rela- considerable. Hence, this scatter should be prop- tionships. The third relationship is given in erly accounted for in the use of the attenuation fig- ure 3-23. relationships. See appendix D for an example. c. Attenuation of ground motion in the West- (1) The mathematical relationship used for ern United States. The abundance of strong modeling the attenuation of peak acceleration motion records in the WUS makes empirical with distance is expressed by Campbell (Biblio 14) by the regression analysis the ideal tool to predict equation: ground motion. A number of assumptions can PGA = aexp(bM)(R + C(M))-dexp(-rR) have a significant impact on the results of such (eq 3-16) regression analyses. The most important ones where PGA are the attenuation mathematical forms, the is the mean of the peak acceleration scaled from the regression techniques (linear, non-linear, horizontal component of the ac- weighted vs. non-weighted), the data base se- celerogram in g units. lection criteria, the definition of magnitude, at- M is the magnitude (M = ML for mag- tenuation, and site soil condition. Three of the nitude less than 6.0) most recent attenuation models developed for (M = Ms for mag- the WUS are given below: nitude greater than -Campbell Model (Biblio 14) 6.0) -Joyner and Boore Model (Biblio 31) 3-33 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 I
zCP 2 0 ti W-J V C) I 0 N -J 0I- w 0B Ld
.01 I 10 100 DISTANCE (kn)
Reprinted from "Near-Source Attenuation of Peak Horizontal Acceleration," Campbell, K. W., and "Peak Horizontal Acceleration and Velocity from Strong Motion Records," Joyn.er, W. B. and Boore, D. M., Bulletin of the Seismological Society of America, Vol. 7I, No. 6, 1981, with permission from the Seismological Society of America. Figure 3-24. Attenuation Relations.
r is the absorption coefficient which af- PGA = 0.22 exp(0.734M)(R + 0.567exp 089 fects the rate of attenuation. (0.345M)Y) 1 exp ( - rR) (eq 3-17) R is the closest distance in kilometers to where the value of r for the WUS is given by the surface projection of the rupture r = 0.0423 - 0.00911M + 0.000573m2 (eq 3-18) zone. The 84th percentile value is obtained by multi- a, b, and d are regression constants. C(M) is a plying equation 3-17 by 1.49. This step assumes function which models possible nonlinear mag- that the natural logarithm of PGA has a stand- nitude and distance scaling effects in the near ard error of 0.40. field that may be supported by the data. Ac- cording to Campbell, The Joyner and Boore relationships (1981) are as follows: C(M) 0.567 exp(0.345M) logA = - 1.02 + 0.249M, - logR1 Substituting this into equation 3-16 along with - 0.00255R, + 0.26P - (eq 3-19) the values for a, b, and d gives the following equation for the median value of peak acceler- where R, = (d2 + 7.3)2 (5.0 < Mo S 7.7) ation: (eq 3-20) 3-34 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A A is the peak horizontal acceleration in g tion model for the EUS is a difficult task for units several reasons. First, there is not much strong motion data available from EUS earthquakes. M. is the moment magnitude. Second, it is generally agreed that one cannot d is the closest distance to the surface pro- directly use a ground motion model developed jection of the fault rupture in kilome- for the Western United States (WUS) because ters. data from a number of sources, e.g., Nuttli (Bib- Ho 48), Chung and Bernreuter (Biblio 15) that P is zero for 50 percentile value and one the attenuation of seismic energy in the EUS is for the 84 percentile value. much different (more gradual) than in the WUS. (2) The OASES (Biblio 70) relationship has Four approaches appear applicable to develop the following mathematical format: an EUS ground motion model. Given the limited amount of intensity data available for the EUS, PGA = biexp(b M)(R + C)b (eq 3-21) 2 3 three of the approaches use intensity as an in- where PGA is the peak horizontal acceleration termediary variable to compare the ground mo- in cm/sec2 . tion between WUS and EUS: bl, b2, and b3 are regression constants Let Is = site intensity R is the closest distance to fault rupture IC = epicentral intensity in kilometers. C is a constant dependent on magnitude R =distance from source to the site M, but independent of transmission M = magnitude path. C = 0.864exp(0.463M.) (eq 3-22) F() and g() functional forms GM =ground motion parameter, such For different transmission paths and soil con- as peak acceleration or peak ditions, values of regression constants b1, b2 and velocity b3 along with the standard deviation of In ( PGA) are given in table 3-2. Use of any one of the Distance Weighting three attenuation relationships should give rea- Is = f (I., R) (EUS Data) sonable results. d. Attenuation of ground motion in the East- Log GM = g(I,,R) and in some cases ern United States. Developing a ground mo- G(ISMLR) (WUS D
Table 3-2. OASES attenuation constants for median PGA values.
Standard Range Deviation of Mag- bI b3 of log nitudes (PGA) Mc_
(Shallow siti 191 0.823 -1.56 0.568 4 to 7.5
Focus ______Events) rock Eve )site 157 1.04 -1.90 0.579 4 to 7.5
Path B (Deep Focus stiff 284 0.587 -1.05 0.70 5 to 8.5 or Subduc- site tion Zone rock 276 0.68 -1.20 0.70 4 to 8.5 Events) site 27 068 1 120 0704t 8. I
Reprinted from "Offshore Alaska Seismic Exposure Study (OASES)," 197S, with permis- sion from Woodward-Clvdc Consultants. 3-35 TM 5-809-10-1/NAVFAC P-355.1IAFM 88-3, Chapter 13, Section A 27 February 1986 Magnitude weighting the future. (EUS Data) (1) Empirical models using an intensity Is = f(I.,R) attenuation data base. The first three ap Log GM = g(l,,M) (WUS Data) proaches require a relation giving the atten-._./ uation of intensity as a function of distance. It No weighting would be ideal to have a number of earthquakes I. = f(I.,R) (EUS Data) with a range of epicentral intensity (IJ) and many (WUS Data) reports of site intensity (Is) for each earth- Log GM = g(15) quake. Then it would be possible to obtain the The fourth method uses a theoretical approach required relation of the form through a simple such as Nuttli's (Biblio 48) model. It combines regression analysis: theoretical modeling with measured regional Q values (damping value of the transmission me- (Is - IW)= C1 + C2R + C3InR (eq 3-23) dium), assumes the near-source ground motion However, no such data set exists in a usable in the EUS is the same as in the WUS, and scales form. Considerable data does exist, but it is in only by magnitude. If it is kept in mind that the the form of isoseismals for given earthquakes. elements of ground motion models are a com- Isoseismals have a number of drawbacks, in- bination of source travel path and local site ef- cluding the fact that they are generally subjec- fects, it can be seen that all four approaches tively determined. Of even greater significance make a common assumption. This is that the set is the fact that isoseismals represent the aver- of WUS earthquakes, making up the strong age distance at which a given intensity was felt, ground motion data set, adequately represents rather than average intensity at a given dis- future earthquakes in the EUS in terms of such tance. Six earthquakes, that have been studied parameters as dynamic stress drop, static stress in enough detail to develop sufficient data for drop, seismic moment, and focal mechanism. determining the required coefficients in equa- Validity of this common assumption can be ver- tion 3-23 by regression analysis, are listed be- ified only as more information is generated in low. Maximum Analysis Name Date Intensity Source Southern Illinois 11-9-1968 Vll G.A. Bollinger ~ (Biblio 9) Cornwall-Massena 9-4-1944 V1l R.J. Holt (Biblio 9) Ossippee 12-20-1940 VlI R.J. Holt (Biblio 9) Giles County 5-31-1897 VII-VIll G.A. Bollinger (Biblio 9) Charleston 8-31-1886 X G.A. Bollinger (Biblio 9) New Madrid 1811-1812 XI-XII 0.Nuttli (Biblio 9) (2) Stronggroundmotiondata base. This to be correlated with spectral amplitude as well data base allows correlation of site intensity with as PGA, the data sets are more limited. The most such information as peak ground acceleration common one consists of the California Institute (PGA), velocity (PGV), distance from recording of Technology (CIT) data tapes, such as those site to the epicenter and/or nearest approach of of Trifunac and Brady or McGuire and Barn- the fault rupture plane, earthquake magnitude, hard. These sets are then used to obtain rela- and information about site geology (See para- tions of the form: graph C-1, appendix C-1). A number of such data bases have been developed, e.g., Murphy In GM = Cl + C21, + C3 In R (eq 3-24a' and O'Brien (Biblio 43), Trifunac and Brady (Biblio 66), McGuire and Barnhard (Biblio 38), or Boore et al., (Biblio 13). If the site intensity is In GM = C1 + C21. + C3M (eq 3-24b) 3-36 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A or Sa of interest. Several models are plotted for PGA in figure 3-25. Based on the methods sug- In GM = Cl + C I1 (eq 3-24c) 2 gested in this section, any one of the following and four attenuation relationships can be used.
In GM = C, + C2M + C3 In R + C4S 1. Gupta and Nuttli model (1976). (Biblio (eq 3-24d) 25). where 2. Bollinger model (1977). (Biblio 8). 3. Ossippee model (1977). (Biblio 64) GM = Ground motion parameter (PGA, PGV, 4. Model developed by Tera Corporation. or spectrum S, at a given period) (b) The Tera Model is based on the first Is = Site intensity three models mentioned above. This model has R = distance measure (epicentral, closest the following format: approach, etc.) log PGA = 0.74 + 1.12 mb - 0.733 In R - M = generally local magnitude 0.0007R (for R > 20 kilometers.) S = Site type parameter (for soil S = 0; for rock S = 1) = -1.47 + 1.12 mb The parameters C. are determined by regression (for R s 20 kilometers.) analysis using an appropriate data set. The val- PGA is in cm/sec2 ues of I., R, and site type for some records differ Mb is the body wave magnitude = significantly between data sets. Thus some (0.98ML - 0.29) choices are involved. R is the epicentral distance in Kms. (3) Site Correction Factor. The ideal way to include a generic correction factor for rock e. Uncertainty associated with ground mo- sites is to perform the required regression anal- tion model applied in the east. One weakness ysis using only the rock subset of the data in of the approach applied in the EUS has to do place of equation (3-23) one could use: with apportioning an attenuation model into submodels. The uncertainty contained in each Is - 1o = C + C R + C In R + C4S 1 2 3 of the submodels increases the uncertainty in (eq 3-25) the final prediction (Cornell, et al., (Biblio 20). where S = site type (S = 0 for soil and S = 1 Although at the present time, there does not for rock), and in place of equations 3-24a, b, c, appear to be any rational alternative to this. one could include a site type in the relation be- This added uncertainty significantly influences tween ground motion, site intensity, and dis- the seismic hazard results. Improved estimates tance or magnitude. Unfortunately, the intensity could be obtained through additional work on attenuation data does not include the site type this topic. When an attenuation model is derived and the intensity assigned is not generally at a directly from recorded ground motion, the sta- site where an accelerograph would be located, tistical uncertainty usually corresponds to a one but rather it is determined from isoseismals or standard deviation confidence level of 1.6-2.0 nearby reports of intensity. This reduces the times the mean. When the uncertainty in mean applicability of the above approach. predictions of intermediate parameters (such (a) Another method consists of intro- as intensity) is rigorously included, this multi- ducing the variation between soil and rock sites plicative factor becomes 2.0-2.9 (Cornell, et al., at the level of equations (3-24) and the general (Biblio 20). A hazard analysis, which results in ground motion model for the EUS is the com- a one standard deviation confidence level equal bination of equation (3-23), the appropriate form to 2 or 3 times the mean predicted value of site of equations (3-24) and the inclusion of the term, severity is being dominated by this multiplica- C4S (S = 0 for soil sites arid S = 1 for rock sites) tive factor. It should be recognized that a large where C4 is obtained from WUS data, (In(GM = part of the uncertainty is due to the use of data C4S + C1 + C2ML + C31nR). The resulting ground representing all possible earthquake types and motion is of the form: all possible travel paths. The necessity for this is to acquire a sufficient statistical sample size ln GM = Cl + C2I0 + C3R + C4ln R for averages and empirical prediction equa- (eq 3-26) tions. However, in most cases the seismic hazard where GM is PGA, PGV or any spectral ordinate at a particular site is largely determined by a
3-37 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986
la
I
)
lu.^ 1X0 * 1000 Epicentral disunce ikm)
"Seismic Hazard Analysis-Solicitation of Expert Opinion," Nuclear Regulatory Commission, reprinted from NUREG/CR-1582, Vol. 3, 1980.
Figure 3-25. Comparison of ground motion models for Mb = 5.5.
particular type of earthquake (e.g., magnitude the random uncertainty associated with predic- range, depth, focal mechanism, etc.), with a par- tion of median PGA levels in the EUS should be ticular path. It is believed that a detailed con- substantially different than in the WUS for given sideration of this specific local knowledge would parameters. Therefore uncertainty measures significantly reduce the attenuation model un- similar to those values obtained in the WUS from certainty. Also, as stated in the next paragraph, direct regression on strong motion data are rec- the median forecasted value of PGA is used for ommended for use in the EUS. scaling the response spectrum shape. The high f. Site severity for scaling the response spec- uncertainty in actual PGA values does not enter trum shape. For the purpose of scaling the ap- into this scaling procedure; only the statistical propriate site response spectrum shape (DAF) sampling uncertainty of predicted median PGA as described in the next paragraph 3-6, it is rec- as it estimates the true (infinite sample size me- ommended that the median or 50 percentile value dian value) median is of concern. Aside from the of PGA be used in the attenuation equation. The use of sub-models (such as conversion of I to mean value shall be used if the median is not M), there is no a priori reason to believe that given by the attenuation equation. For a given 3-38 27 February 1986 TM 5-809-10-1 /NAVFAC P-355.1 /AFM 88-3, Chapter 13, Section A convoluted seismic hazard or return period of 'probability of exceeding a given level of site in- severity at the site, it is judged that the median tensity (such as PGA) involves the convolution value is sufficiently conservative for spectral of the probabilities of all the possible combi- scaling purposes. Note that PGA data used for nations of source intensities (M or I) and at- empirical attenuation relations is the PGA from tenuation distances R that can produce or exceed the principal component of the recorded time the given level of PGA. Figure 3-26 provides a history. Further conservatism due to the spec- simplified illustration of the typical condition tral enveloping property of the specified DAF for a line source and an area source. shape is discussed in paragraph 3-7. (1) On the line source the set of all possible g. Computation of total hazard at the combinations of rupture length location, its cor- site. The process of computing the hazard or responding attenuation distance RHand mag-
iRupture length due to M having random location Line Source on source
jRI
Site having given level of PGA
Area element having random location on source M occurs in this element.
Area Source US Army Corps of Engineers Figure 3-26. Descriptionof sets of M and R required for a given PGA. 3-39 TM 5-809-10-1 /NAVFAC P-355. I/AFM 88-3, Chapter 13, Section A 27 February 1986 nitude Mi are able to produce or exceed the given -Analytical Soil-Column Response PGA at the site. Similarly the set of the M; and area element location Rj produces the PGA at (see para C-3, app C for an overview of all meth- the site from the area source. ods). The results of any or all of these methods (2) The total probability of exceeding the may be combined to define the appropriate spec- PGA is the probability of the union of the oc- tra for structural design and analysis; this is currences of all the sets of Mi and Ri combina- usually done in a rather subjective manner to best represent the quality of information from tions on the line source, and Mj and Rj combinations on the area source. The convolu- each method, (See SEAOC Pamphlet, (Biblio tion operation required for this total probability 55)). However, before proceeding to the descrip- for a selected range of given PGA values can be tion of these methods and the formation of the very lengthy and is best performed by a com- site specific spectra, it is useful to review the puter program such as STASHA (Stanford Uni- major factors that govern the shape and size of versity Technical Report, No. 36). A simple the response spectrum. example of this type of calculation is given in a. Spectral shape factors. It is generally rec- paragraph 3-7c. ognized that the frequency content and corre- sponding response spectrum shape is governed (3) Finally, a sensitivity analysis involving the probable upper and lower bound values of by the following source and site factors. the parameters of the hazard analysis may be -Characteristics of Soil Deposits Underly- performed. For example, when large uncertain- ing the Site ties exist due to sparse data and (or) judge- -Magnitude of Seismic Event producing the mentally assigned values in source locations R, Site Ground Motion recurrence parameters a, I, Mmax, and different but applicable attenuation relations, then sep- -The Source Fault Rupture Characteristics arate runs of PGA evaluations may be per- formed using probable upper and lower bounds -The Source-to-Site Travel Path Charac- for each individual parameter. The results of teristics of Distance and Wave Attenua- this analysis are useful to identify the impor- tion Properties tant factors that significantly effect the calcu- The second and third factors are recognized sub- lated PGA, such that perhaps more information jects of research, but are not generally incor- ) can be obtained to better evaluate these factors porated in site spectra with the exception that of parameters. Also, the resulting probable records for spectral averaging purposes may be' bounds on a PGA for a given return period pro- grouped according to magnitude levels. The first vide a numerical description of the quality or "soil type" factor is well established and used in. stability of the hazard analysis and can assist in most site specific ground motion studies. The the final assignment of the design spectral scal- fourth "travel path" factor is also an estab- ing value for the PGA. lished procedure for both distant sites in all re- gions, and for the representation of the low 3-6. Site specific response spectra, step V. attenuation rates in the Eastern United States. The exact prediction of future ground motions Detailed discussions and procedures for deter- (such as the accelerogram x(t)) at a site is not mination of spectra are given in appendix C, par- possible. Therefore, forecasted response spec- agraphs C-2 and C-3. tra representative of this motion offer the most b. Statistical averages of normalized re- effective method of specifying the future. Hav- sponse spectra. In this first empirical method, ing the value of site severity from step IV of the the shape of the spectrum is determined by a seismic hazard analysis, this value provides the statistical analysis (evaluation of averages and basis for scaling the response spectrum shape standard deviations) of past earthquake strong resulting from step V, treated in this paragraph, motion accelerograms; as classified according to and summarized in figure 3-27. site conditions, distance from the source and In practice, the response spectrum shape may size of the event. All the response spectra for a be obtained by three rather common tech- common set of conditions are normalized by the niques; two of which are empirical, and one ana- recorded PGA, see figure 3-28. lytical method: The mean and standard deviations of the nor- malized spectra (referred to as the Dynamic -Averages of Normalized Spectra Amplification Factor or DAF) are then calcu- -Attenuation of Spectral Ordinates lated. This statistical summary is used to fore-
3-40 27 February 1986 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A
GIVEN FORECASTED MEDIAN PGA FOR EQ-I AND EQ-II (para 3-le) DETERMINE FORECASTED RESPONSE SPECTRA 1 4'I AVERAGES OF ATTENUATION OF ANALYTICAL SOIL- NORMALIZED SPECTRAL ORDINATES COLUMN RESPONSE SPECTRA HAVING FOR APPROPRIATE FOR APPROPRIATE COMMON SITE SOIL TRAVEL PATH, TIME HISTORIES CONDITION MAGNITUDE, AND ON BED ROCK (para 3-6b) SITE SOIL CONDITION AND SOIL COLUMN (para 3-6c) MODEL (para 3-6d)
i
ATC 3-06 SHAPE AS SPECIFIED IN PARA 3-8 ] -I ____
COMPARE AND JUSTIFY FINAL SPECTRAL SHAPE USING INFORMATION FROM A.LL AVAILABLE METHODS (appendix C, para C-4)
SCALE FINAL SPECTRAL SHAPE ACCORDING TO EQ-I AND EQ-II LEVELS OF PGA AND DAMPING VALUES (para 3-7 and 3-8)
US Anny Corp~s of Eiigineers
Figure 3-27. Step V, site specific response spectra.
3-41 TM 5-809-1 0-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986
,Sa DAF 1 a(t)1
B DAMPING
tV,06 v , t T T
Sa 2 DAF2
t T T DAF Sa n a(t) n I n )
'04-00 -I!, A A AA . . t - V V \1 V , T T n iMDAF X DAF. MDAF = 1 1 (For all values of T) n
n 2 E (DAF - I MDAF) Variance (DAY) I (For al 1 (n-i) values of T) T
US Army Corps of Engineers Figure 3-28. Statisticalaveraging of normalized spectra.
3-42 27 February 1986 TM 5-809-1O-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A cast the spectral shape of future events according information is not complete. See appendix C, to the particular site conditions. The method, paragraph C-3. even though widely used for practical applica- e. Site specific earthquakespectra. The pro- tions, has some shortcomings. The procedure of cedures of paragraphs 3-6b, c, and d, have all or normalizing according to PGA creates a large in part lead to generalized versions of earth- coefficient of variation (standard deviation di- quake spectra. Some of the important recom- vided by average), particularly in the long pe- mendations resulting from these procedures are riod region. However, since no better means of given here and in the next paragraphs on shape normalization is yet available, this technique has effects. These include the methods of: provided the primary source of design earth- quake spectral shapes. See Seed et al. (Biblio -Newmark-Hall, (Biblio 44) 56), Kiremedjian and Shah (Biblio 33) and ATC -Seed et al, (Biblio 56) 3-06 (National Bureau of Standards, Special Publication 510). -Kiremidjian and Shah, (Biblio 33) c. Attenuation of spectral ordinates. The -ATC 3-06 second empirical approach of forming a site (1) Newmark-Hall Method of Constructing spectrum is by the use of attenuation equations Elastic Response Spectrum. This is an empirical for spectral ordinates at specific period values method of constructing an elastic spectrum. It for a set of records and then statistically ana- employs the following normalized values for lyzing these attenuated ordinates. This again ground motion: provides a mean and standard deviation descrip- tion of the site spectrum such that an upper Acceleration 1g confidence limit can be given in terms of one or Velocity 48 in/sec. more standard deviations. This method has the advantage of avoiding a normalization method Displacement 36" with its inherent creation of large spectral var- Thus, for a peak ground acceleration of interest, iability. This advantage is offset, however, by as forecasted for the site, construct the ground the need for the use of spectral attenuation re- motion parameters on the tripartite plot. As an lations that have large prediction error. Also, example, let the PGA value be 0.35g. For this the development of these relations requires a case, ground motion values are: sufficient set of records applicable for a common seismic region; the method is therefore limited Acceleration A = 0.35g to these regions (see app C, para C-3d). This (Og x 0.35) method, however, may find increased applica- Velocity V = 16.8 in/sec. bility in the Eastern United States (see Nuttli: (48 in/sec x 0.35) Biblio 47), not because of the availability of data for that region, but because the method can in- Displacement D = 12.6" corporate expert opinion and theories for wave (36" x 0.35) transmission peculiar to the region and its pos- Draw this ground motion spectrum on the tri- tulated sources of seismicity. The most current partite paper. (fig 3-29). application of this technique is given by NUREG/ (a) The second step is to construct an CR-1582, Vol. 3 and 4, (Biblio 63, 64). "elastic" response spectrum. To construct this d. Analytical soil column response. The third spectrum, a table of amplification factors, based or analytical method of obtaining a spectral shape on the study of past spectra, is available. See is based on a site specific study of the strong table 3-3 from (Biblio 44). motion accelerogram. If the acceleration time These amplification factors are functions of history at the bedrock level for a given site can damping ratios, and the described confidence be formulated, then using the overlying soil lay- level. As an example, consider a 5 percent damp- ers as a filter, the response on the surface can ing ratio, and the median level. be determined. Thus, the transfer function of (b) The lines of constant acceleration, ve- the soil layer and the motion at the bedrock level locity, and displacement representing the elastic determines the time history and corresponding response spectrum are given by the correspond- spectral shape at the surface. The problem with ing ground motion values times the appropriate this method is that a time history at the bedrock factors from the table. level has to be formulated. This may not be an easy task for a region where the seismotectonic S.= (.35g)(2.12) = 0.74g
3-43 TM 5-809-10-1/NAVFAC P-355.1I/AFM 88-3, Chapter 13, Section A 27 February 1986
)
FREQUENCY - Hz 300.0 30.0 IC0 0I litIlEI I I .,, I I I r L , ...... - 100.0 -______
U
10.0
4-
L- U 0
z 0 1.0 a. a W )
a- 4 a 0S a
P in 0.1
0.01 0.ol 0.1 I 0 FuRIOD-SEC Reprinted from "Earthquake Spectra and Persig,." Newmark, N. M. and Ilal), W. J., EERI Monograph Series, 1982, with permission from the niarthqtaakr Engineering Research Institute. Figure 3-29. Newmark-Hall Spectrum.
3-44 27 February 1986 TM 5-.809-1O0-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A Table 3-3. Spectrum amplification factors for horizontal elastic response.
Damping, One Sigma (S4.1o) Median (50%) 9eCritical A V D A V D
0.5 5.10 3.84 3.04 3.68 2.59 2.01 1 4.38 3.38 2.73 3.21 2.31 1.82 2 3.66 2.92 2.42 2.74 2.03 1.63 3 3.24 2.64 2.24 2.46 1.86 1.52 5 2.71 2.30 2.01 2.12 1.65 1.39 7 2.36 2.08 1.85 1.89 - 1.51 1.29 10 1.99 1.84 1.69 1.64 1.37 1.20 20 1.26 1.37 1.38 1.17 1.08 1.01
Reprinted from "Earthquake Spectra and Design," Newmark, N. M. and Hall, W. J., EERI Honograph Series, 1982, with permis- sion from the Earthquake Engineering Research Institute.
S, = (16.8 in/sec) (1.65) = 27.7 in/sec response spectrum. This total risk must involve Sd= (12.6')(1.39) = 17.5 in. the convolution of probability functions for both the forecasted PGA scaling factor and the DAF (c) These constant levels are plotted on spectral shape. See Kiremidjian and Shah (Bib- the tri-partite paper, and along with recom- lio 33) for examples. A more simplified reliabil- mended connecting lines as given in (Biblio 44), ity calculation is given in paragraph 3-7. the complete spectrum is defined. This New- (4) The ATC 3-06 method uses much of mark/Hall method provides a direct procedure methods (1) and (2) as background justifica- of forming a spectrum, and also has the advan- tion. It, however, goes further to provide sim- tage of constructing inelastic yield force and de- plified DAF shapes for not only the soil types formation spectra in terms of structural ductility but also the tectonic region. Because of this sim- factors (see Biblio 44). Also the site soil con- ple, yet representative quality, it is recom- ditions can be represented by either the known mended that these ATC 3-06 shapes be used for forecasted peak ground velocity or prescribed the appropriate site conditions and tectonic re- relations between peak ground acceleration, ve- gion. Therefore, unless there are special site locity, and displacement. However, since the conditions, close active sources, or high risk f a- representation and description of site soil con- cilities, these shapes as scaled by the forecasted ditions are not as detailed as in the following site severity values can provide the input spec- methods, the use of this Newmark/Hall method tra for design and analysis. The complete ATC is not recommended except for general compar- 3-06 method for site severity and response spec- ison with other methods. tra is given in paragraph 3-8. In order to rep- (2) SeedetaL. This method provides mean resent the particular regional attenuation effects DAF, and mean plus one standard deviation that are indicated when the A, value exceeds the shapes for different categories of site condi- Aa value on the contour maps* given in para- tions, see figures 3-30 and 3-31. These DAF graph 3-8, the spectral shape should be found shapes may be scaled to the forecasted PGA value using the respective contour map values of A, having a given risk value at the site. and A., then this shape should be scaled by the (3) Kiremidjian and Shah. This method is ratio of the forecasted PGA to the contour map similar to method (2), and a definite listing of value of A.. The PGA value corresponds to the the data base and the site soil conditions is pro- hazard level or return period of EQ-I or EQ-Il. vided. Also, in addition to mean and mean plus f. Factors affecting response spectral one standard deviation shapes, probability func- shapes. As mentioned in paragraph 3-6a there tions are given for the random DAF values as are several important conditions or factors that they are scattered about the mean value. This can alter the shape or frequency content of the probability information is most useful for cal- response spectrum. culating the total risk of exceeding a specified 3-45 TM 5-809-1O0-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986
0
0
*I
U
0 E a E a
0 0.5 ID L5 2.0 25 3.0 Period - seconds
Reprinted from "Site-Dependent Spectra for Earthquake Resistant Design," Seed, It.B. et al, Report No. EERC 74-12, University of California at Berkeley, 1974.
Figure3-30. Average acceleration spectra for different site conditions.
(1) Type and duration of fault rup- source to the site increases. In other words, the ture. Generally the type and duration of the predominant period of motion increases with fault rupture affects the frequency content of distance and size of the seismic event. The en- the seismic wave. Various seismological papers gineering implication of this observation is ob- are available which describe the theoretical for- vious. Taller structures are affected more by large mulation of the above mentioned dependence. distant earthquakes than are the shorter (or (Haskell, (Biblio 28, 29); Savage, (Biblio 51)). stiffer) structures at the same location. According to these models, the seismic wave (3) Local site soil conditions. The effects characteristic in the time and frequency domain of local site soil conditions on the frequency con- is a function of the radiation pattern (source tent can be very significant. The response of a and propagating geometry), seismic moment given layered soil media to a seismic bedrock (size of the event or energy release level) and motion depends heavily on the transfer function the source mechanism. of the soil. Thus, stiffer soils transfer higher (2) Size of event in terms of magnitude or frequency components whereas softer soils seismic moment and distance from source to transfer lower frequency components. Exten- site. Based on the recorded ground motion sive studies of the available strong motion ac- characteristics, many empirical relationships are celerograms by many researchers have shown available to show the dependence of the re- that the shape of the Response Spectrum changes sponse spectrum shape on the size of an event, with the site condition. There are usually three the distance from the source to the site and the classifications of soils: soft alluvium deposits predominant period (or frequency) of the mo- (soil class 0), intermediate stiff soils (soil class tion. Figures 3-32 and 3-33 show such empirical 1) and firm soils or rocks (soil class 2). These results. It can be seen from these figures that classifications could be made on the basis of shear the higher frequency components are filtered wave velocities. As a guide to such a possible out from seismic waves as the distance from the 3-46 27 February 1986 TM 5-809-10-1 /NAVFAC P-355.1 /AFM 88-3, Chapter 13, Section A
C 0 C
0I'l 4.1
U C
E X0 EI WI aU I
Period - seconds
Reprinted from "Site-Dependent Spectra for Earthquake Resistant Design." Seed, H. B. et al, Report No. EERC 74-12, University of California at Berkeley, 1974. Figure 3-31. 64 Percentile acceleration spectra for different site conditions.
3-47 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 )
1.4
1.2
1.0 A~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
i0.8
J0.6 AA
04 LEGEND ) _;_ 0 0 ~ ~~~~~~~~~~.0SO(SmWUtmY PD 0 Auaf.tu Pod 02 A COOMbnwat
0 150 200 300 350 400 450 500 Uo DO 00 Epicantr Difv4c-km
Reprinted from "Characteristics of Rock Motions During Earthquakes," Seed, R. B. et al, Report EERC 68-5, University of California at Berkeley, 1968. Figure 3-32. Predominantperiods for motions in rock-earthquakemagnitude = 7.
3-48 27 February 1986 TM 5-809-10-1/lNAVFAC P-355. IAFM 88-3, Chapter 13, Section A
5.2
M.?
0.6 M.6~~~~~~~~~~~~~~~~~..
10 40 8M10- 20 20 .5326 ICA_.AA._.
N~ow fromn CousWir F&A -kem
0 25 50 75 BOO 125 ISO 175 200 2wz Distonceforom Cavwivev Podt-fmies
Reprinted from "Characteristics of Rock Motions During Earthquakes," Seed, H. B. et al. Report EERC 68-5, University of California at Berkeley, 1968. Figure 3-33. Predominantperiods for maximum accelerationsin rock.
3-49 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986 classification, the following procedure is rec- soils of class 2. Under very special conditions, ommended (Vs is the shear wave velocity): (such as in Mexico City, where the city is on an old lake bed), the spectral peak could occur at Firm Site: Vs - 450 meters/sec. a period as long as 1.5 to 2.5 seconds. Intermediate Stiff: 250 - Vs < 450 meters/ (4) Regional geology. This is a most im- sec. portant effect, not only for the Western United States where there is a reasonable amount of Soft Alluvium Deposits: Vs < 250 metersl strong motion records, but for the Eastern sec. United States where data is sparse and predic- Also, for the purpose of this manual, these soil tions of future ground motion must be based classes 0,1,2 may be considered to correspond to upon geological features. The future develop- the soil types SI, S2, S3 respectively, as described ments in ground motion prediction will depend in table 3-5. Figures 3-34, 3-35 and 3-36 taken strongly upon inferred behavior of possible from Kiremidjian and Shah, (Biblio 33), show earthquake source mechanisms, and the corre- the effect of the soil conditions on the frequency sponding propagation of effects in the general content of ground motion. It can be seen from geological structure. One of the most prominent these figures that for soil class 0, the spectral characteristics of Eastern United States seis- peak occurs at higher period than for stiffer micity is the exceptional transmission of peak
mean + one standard deviation 3.' I.' 0 u *~ -- Kiremidj ian and Shah ) 0
-4 'a. '44
.-4
U 1.0 C,
0. 0.0 0.5 1.0 1.5 2.0 2.5 3.0 Period T in sec.
Reprinted from "Probabilistic Site- Dependent Response Spectra," Kiremidjian, A. i. and Shah, H. C., Journal of the Structural Division, Proceedings of the ASCE, %ol. 10D, No. STI, January 198o, with permission from the American Society of C:ivil Engineers.
Figure3-34. Comparison of DAF from Kiremidjian and Shah to Seed et al, soil class - 0, damping = 5%. 3-50
Fu 27 FebruarY 1986 TM 5"09-1 -1NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A
V q4 I I 1.' 11 Kjiremidjian and Shah 03 . I"~ I ……- -- Seed et al C. ~I, to
mean + one standard deviation
U -4
(0 sUI~
mean
0. 0.5 0.0 1.0 1.5 2.0 2.5 3.0 Period T in sec.
Reprinted from "Probabilistic Site- Dpeendent Response Spectra, K;rnfiifiaIt A. S. and Shah. it. C., Journal of the Structural Division, Proceedings of the ASCE. Vol. 106. No. ST1. January 1980, wit permiSsion from the American SocietY of
Civil Engineers. Co present study to DAF from Seed e Figure 3-3. Comnparisonof DAF from rsn td oDFfo ede l olcas 1 apn %
3-51 TM 5-809-10-1/NAVFAC P-355.1/AFM 88-3, Chapter 13, Section A 27 February 1986