See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/28784819

Response envelope spectrum and interpretation of strong

earthquake ground motion

Article in Bulletin of the Seismological Society of America · April 1971

DOI: 10.1785/BSSA0610020343 · Source: OAI

CITATIONS READS

110 249

1 author:

Mihailo Trifunac
University of Southern California
284 PUBLICATIONS 9,647 CITATIONS

SEE PROFILE

Some of the authors of this publication are also working on these related projects:

3D Wave Diffraction View project

Nonlinear response of soft soils during strong earthquakes: full-scale observations and numerical simulations using finite differences
View project

All content following this page was uploaded by Mihailo Trifunac on 22 August 2014.

The user has requested enhancement of the downloaded file.

 Bulletin of the Seismological Society of America. Vol. 61, No. 2, pp. 343-356. April, 1971

RESPONSE ENVELOPE SPECTRUM AND I N T E R P R E T A T I O N OF STRONG

EARTHQUAKE Gt~OUND ~/[OTION

BY M. D. TRIFUNAC*

ABSTRACT

A multiple filter technique is developed using a single-degree-of-freedom

lightly-damped oscillator as a narrow band-pass filter. It provides for a physi-
cally simple approach to the analysis of accelerograph records from the point
of view of structural engineering.
The analysis of several typical accelerograph records indicates that a signifi-
cant portion of strong earthquake ground motion consists of surface waves. It is
concluded that the duration of intense shaking will be determined predominantly by
the dispersion properties of the ground.

INTRODUCTION
An understanding of the phenomena associated with strong earthquake motion is
important for many areas of research and application. This motion of destructive in-
tensity endangers the safety of all man-made structures and has, on several occasions,
led to complete or partial devastation of whole cities. The recording and interpreting
of strong-motion data is also essential for various studies of earthquake mechanism,
for high resolution of the patterns of earthquake energy release and, in general, for
better understanding of the close field of earthquake ground motion.
The basic information related to this phenomenon comes from strong-motion ac-
celerographs. The first instruments of this kind were installed in the field some 40 years
ago, and the first strong-motion record was obtained during the 1933 Long Beach,
California, earthquake. Since that time, several hundred earthquakes have been re-
corded, predominantly in California and Japan.
Although the history of modern recording of teleseismic earthquake waves is not
much older than the measurements of strong, near-field, ground motion, the nature of
teleseismic waves and the great number of recorded shocks resulted in early interpre-
tations of the different wave forms displayed in those records. This was not the case
with the strong-motion records.
The inhomogeneities and discontinuities in the Earth's crust lead to a wide variety
of wave velocities. For distant recordings different frequency-wave components are
dispersed and cause natural separation of various wave types. This separation, which
helps in the record interpretation, is further enhanced by the great distances and hence
the long travel times which allow dispersion effects to take place.
Strong earthquake ground motion close to the source of energy release is less ob-
vious to interpret. In principle, it is composed of essentially the same types of waves,
but often they are not clearly separated in time because of the source proximity. In
addition, high frequency motions, whose amplitudes decay rapidly with distance and
so cannot be observed on most teleseismic records, are also present. Also source size,
relative to the distance to the recording station, and the spatial and temporal distribu-
tion of energy release become important. All of these complexities are probably the

* Present address: Lamont-Doherty Geological Observatory, of Columbia University,

Palisades, New York.
343
3Axz~ BULLETIN OF THE SEISMOLOGICAL SOCIETY OF AMERICA

reason why, so far, no notable attempts have been made to determine the detailed
character of strong earthquake ground motion.
Because the collection of recorded strong-motion accelerograms for engineering ap-
plications is still far from complete, many research workers have made valuable at-
tempts to generate artificial accelerograms. Some of these are based on white noise
(Housner, 1947; Rosenblueth, 1956; Rosenbleuth et al., 1962; Bycroft, 1960), Sta-
tionary Gaussian process (Tajimi, 1960; Housner and Jennings, 1964; Barstein, 1960)
and various nonstationary processes (Bolotin, 1960; Bogdanoff et al., 1961; Cornell,
1964; Amin and Ang, 1966; Shinozuka and Sato, 1967; Jennings et al., 1968). Although
these approaches represent a significant step forward in the engineering use of ac-
celerograms, they are often based on simplified interpretations of the properties of
strong motion. Usually only the simplest features of the recorded accelerograms are
employed to model the "rapid build-up" of the motion, the "central portion of strong
shaking" and the "decaying tail." The nonstationary time dependence of the intensity
and frequency content of strong motion are usually neglected in such models. The
reasons for these simplifications come, partly, from the mathematical complexity in
generating these processes but are mainly the result of the lack of our knowledge of
the actual character of strong-motion aceelerograms.
The area of earthquake engineering research concerned with the effects of local
geology on the amplitudes of strong-motion waves is also directly affected by the lack
of detailed interpretation of strong-motion accelerograms. Most of these studies are
based on models in which earthquake waves are assumed to propagate vertically
toward the set of horizontal surface layers, irrespective of the question of whether such
motions are predominant or not. However, the simple qualitative analysis of the param-
eters which govern the properties of strong-motion waves leads to the conclusion that
a predominant part of the near-field ground shaking may be composed of surface
waves associated with energy propagating horizontally through soft surface layers as a
wave guide (Trifunae, 1969). The following analysis of some typical strong-motion
records supports such an interpretation.

RESPONSE ENVELOPE SPECTRUM (RES)

The need for improved methods of assessment of information from seismograms
stimulated the development of the moving window Fourier analysis and multiple filter
techniques. Alexander (1963) used frequency and velocity windows to extract informa-
tion on wave modes. Archambeau, Flinn and Lambert (1966) studied the dispersion of
body-wave arrivals by using digital filters. Landisman, Dziewonski, and Sato (1969)
developed a process called "moving window analysis" to display the amplitudes and
phases of transient seismic signals as functions of period and group velocity. Dzie-
wonski, Bloch and Landisman (1969) reported on the "multiple filter technique" as an
efficient method for analyzing multiple dispersed signals. These approaches based on
different numerical methods are aimed at the same objective, which is to extract a
maximum of information present in the seismograms. This information typically con-
sists of the dispersion as indicated by arrival times of different frequency components
in the seismic surface waves. It further consists of the phase data, ellipticity of the
surface particle motion, and other features from which the character of energy release
at the source can be determined. Trifunac and Brune (1970), for example, used a simple
form of moving window Fourier analysis to study the extended time behavior of seis-
mic energy release.
In engineering seismology and in particular in the analysis of strong-motion ac-
 INTERPRETATION OF STRONG EARTHQUAKE GROUND MOTION 345

celerograms, similar information is needed. The knowledge of the time dependent

amplitudes and the frequency content of strong-motion records is required for pre-
dicting the vibration behavior of engineering structures.
Small motion of a single-degree-of-freedom, viscously-damped, linear oscillator is
described by the well-known differential equation
2
2 + 2con~2 ÷ co, x = - a ( t ) (1)

where
x = relative motion of the mass,

= fraction of critical viscous damping,

co~ = natural frequency of vibration,

a(t) = absolute acceleration of the ground.

The vibration of several types of structures such as elevated water tanks and small
buildings can be described by the single-degree-of-freedom oscillator. For tall build-
ings, chimneys, and towers, when contribution of the higher modes of vibration can-
not be neglected, each mode of vibration may be represented by the equation of the
same form as (1). Then, by combining the response from all higher modes, the result-
ing response of the multi-degree-of-freedom system can be determined (Merchant and
Hudson, 1962). Hence, from the engineering point of view, the most valuable in-
formation about the structural vibrations during earthquake motion will be ob-
tained by studying the response of the system described by equation (1). Similar
ideas led to the introduction of the concept of the response spectrum by Benioff
(1934) and Blot (1941) which was further developed and applied by Housner, Martel
and Alford (1953), Hudson (1956), and others.
In order to find a physically simple method for analyzing aeeelerograph records
from the point ;of view of structural engineering, we consider again equation (1).
It is well known that if for a(t) is substituted -Aco 2 sin cot and if it is assumed that the
vibration will go on for a long time, so that transient vibrations generated by the ini-
tim conditions have decayed, the steady state part of the relative response x(t) is
determined by the following expression

(co)2
x(t) = ~ A-sin (cot - ,1~) (2)

\con/_1 COn

where

(2co_. r_
-- tan-1 k,co2 _ co2/" (3)

The factor 1 / % / - in equation (2) represents a dynamic amplieation factor and gives
the ratio between the dynamic and static deflections. It is easily seen that for small
/', this factor is maximum when w = con and equals ½f. For most engineering structures,
346 BULLETIN OF THE SEISMOLOGICAL SOCIETY OF AMERICA

¢ the equivalent viscous damping is in the range between 1 and 5 per cent of critical.
For structures oscillating in the nonlinear and plastic range, the equivalent /" in-
creases to about 10 per cent and more. Also when co = co~ from equation (3), it follows
that the phase shift • is equal to rr/2. Thus we see that a typical lightly-damped
single-degree-of-freedom system acts like a narrow band filter which amplifies the
input frequencies, centered around co., by 1/2f times and follows the input wave
form with 7r/2 delay. When A = A(t) is a slowly changing function of time, the re-
sponse x(t) will also have a slowly-varying amplitude. The envelope of the response of
the single-degree-of-freedom oscillator is thus approximately A (t)/2~-.
The physical appropriateness of the single-degree-of-freedom system as a narrow
band filter suggests that the response of a series of such systems can be used to infer
the character of the input acceleration. This is, of course, the basic idea behind the

{D ~r--

LJ
©

~ O ~
_m
C3

> - -

SECONDS
L I I I I I I I I I 210 I
0 tO

FIG. I. The relative displacement and its envelope for an oscillator with natural frequency
o~ = i0 rad/sec and fraction of critical damping ~" = 0.10, for the NS component of the El Centro,
1940, strong motion accelerogram (Figure 5).

Response Spectrum, and the following analysis is simply a means of introducing more
detail into the use of such spectra. Computing the response envelopes for many filters
for equally-spaced frequencies will enable one to construct a three-dimensional dia-
gram which indicates the time and frequency dependence of the instantaneous en-
velope A(t, co) of the harmonic component of ground motion with frequency co at
time t. We will call A(t, co) the Response Envelope Spectrum (RES) of the ground
displacement. Likewise we will call coA(t, co) the RES of the ground velocity and c02A(t, co)
the RES of the ground acceleration.
The relationship between the RES and the response spectrum is considered next.
A point on the relative displacement spectrum for the frequency co corresponds to
the maximum of A (t, co) during the time of excitation. Thus RES contains all of the
information required to construct the relative response spectrum (i.e., relative dis-
placement, and approximately velocity or acceleration spectrum), and, in addition,
retains information about the time at which various responses occur.
The RES technique is closely related to the "multiple filter technique" or the "mov-
ing window- analysis." It is, in fact, simply a multiple filter technique with filter prop-
erties specified by the single-degree-of-freedom, viscously-damped oscillator. There-
fore the RES can be used as a tool to study dispersion properties of seismic waves and
complicated multiple arrivals. Like multiple filter techniques or moving window
 INTERPRETATION OF STRONG EARTHQUAKE GROUND MOTION 347

analysis, it enables one to determine group-velocity curves for those cases for which
the signal-to-noise ratio is so small that the peak and the trough method cannot be
applied. Another important merit of the RES is that it conveys the information about
the length of large amplitude response of structures.
The RES in Figures 3 to 6 were computed in the following way. Recorded ground
acceleration, shown on the side of each figure, is substituted in equation (1) as a(t).
The a(t) for computation was in the form of equally-spaced data points with At =
0.02 see, interpolated from the original unequally-spaced data. Equation (1) was then
integrated b y the Adams-Moulton method, for each frequency co and damping f =
0.10 of critical. The envelope of the relative response x(t) was approximated by con-
neeting the successive peaks in x(t) by a straight line. This is illustrated in Figure 1,
where the response and the envelope are computed for the NS component of the E1
Centro 1940 aeeelerogram for an oscillator with co~ = 10 rad/sec and f = 0.10. The
computations were repeated for 100 oscillators with frequencies equally spaced at
aco~ = 0.166 rad/see and ranging from 0.166 rad/sec to 16.66 rad/see. All of these

TABLE 1
INFORMATION RELATED TO FIGURES 3, 4, 5 AND 6

Figure No. Event* Component Peak Amplitude w2A EpicentraI Distances

(crn/sec2) A (kin)

3 9 NS 13.58 19
4 9 EW 19.80 19
1A 10.0
5 1B NS 143.29 12.5
1C 15.0
1A 10.0
6 1B EW 89.78 12.5
1C 15.0

* These events were recorded during the Imperial Valley California earthquake, 1940. The
names: Event IA, IB, IC and Event 9, are derived from our previous study (Trifunac and Brune,
1970).

envelopes together define J A (t, co). The RES for co2A(t, co) was chosen for this work
because, as will be seen later, our main interest in this paper is concerned with studying
the properties of accelerograms in which higher frequencies are emphasized. To form a
rectangular mesh for Figures 3 to 6, it is necessary to interpolate equally-spaced
envelope amplitudes along the time coordinates. This was done by interpolating
amplitudes every 0.25 sec. For simplicity of visual presentation, all of the amplitudes
were normalized so that the peak equals 5. Intermediate amplitudes were determined
in the following way. When the normalized amplitude An was k < A,, < k + ½, for
k = 0, 1, 2, 3, 4, 5, A~ was rounded off to k and that value was printed in the cor-
responding slot. When lc ÷ ½ < A n ~ k ~- 1, An was substituted by a blank and
nothing was printed on the output paper. In this way elevation contours of J A (t, co)
became clear and easy to interpret. Because of the photographic reduction of the
original scale, all of the contours of J A (t, co) in Figures 3 to 6 are also traced by hand.
The actual amplitudes of the RES may be obtained by multiplying normalized am-
plitude values by the "peak value" given in Table 1 for each figure. Arrivals of the
surface-wave groups were computed using the theoretical dispersion curves (Figure 2)
and the epicentral distances shown in Table 1.
 348 BULLETIN OF THE SEISMOLOGICAL SOCIETY OF AMERICA

ANALYSIS OF ACCELEROGRAMS FOR THE IMPERIAL VALLEY IV~AY 18 1940

CALIFORNIA, EARTHQUAKE

Several typical accelerograms from the E1 Centro accelerograph record of May 18

1940 were chosen as the optimal samples for this analysis for the following reasons.
The dispersion of the surface waves can be theoretically estimated, provided the layer
thicknesses and velocities of wave propagation are well known in the area between
the station and the source. Detailed studies of the local travel times and layer thick-
ness in the Imperial Valley were made by Biehler (1964). Trifunac and Brune (1970)
used the data from the Westmoreland seismic profile (Biehler, 1964) to compute the
group-velocity curves for several Love and Rayleigh modes (Figure 2). These dis-

4.0 I

5 2

>:
9.0
4/ j /

> j~4 ~ 2

o I
LO ~ ~

LOVEWAVES
----- RAYLEIGHWAVES

0 I
I 2 B 4 5
PERIOD, seconds
FIG. 2. Love- and Rayleigh-wave dispersion curves for the structure corresponding to the West-
moreland profile (Biehler, 1964) in the Imperial Valley.

persion curves are based on the five-layer model. Layer properties were taken as
shown in Table 2 (Trifunac and Brune, 1970).
The flat character of the fundamental-mode group velocity indicates that most of
the energy in the short-period range (shorter than about 5 see) will travel with a
velocity of about 1 km/sec. In the same paper, approximate distances and tentative
locations for all shocks recorded on the E1 Centro strong-motion accelerograph were
determined. Trifunac and Brune (1970) also interpreted all obvious (P and S)
phases for all recorded events. Therefore, these accelerograms permit a comparison
of an interpretation based on RES with that carried out in the previous study.
We first analyze an aftershock referred to as Event 9 (Trifunac and Brune, 1970).
Figure 3 shows the trace and RES of the NS acceleration component and Figure 4
shows the same for the EW acceleration component. Epicentral distance for this
event was estimated (Trifunac and Brune, 1970) to be about 13 to 19 km from El
Centro. If it is assumed that the actual distance was about 19 km, it is possible to
 INTERPRETATION OF STRONG EARTHQUAKE GROUND MOTION 349

calculate the expected arrival times for all surface phases for which theoretical dis-
persion curves were calculated (Figure 2). It should be pointed out that our interpre-
tation here is only qualitative since we have estimated only the theoretical group
velocities in the absence of actual measurements. Nevertheless, when the expected
arrivals for the surface phases are superimposed on RES, one can easily associate
various bursts of arriving energy in the aecelerogram with the estimated arrivals.
In this way, we interpret the RES high amplitudes between 17 and 22 sec (Figure 3)
to be associated with the fundamental Love and Rayleigh waves. Likewise, waves
arriving during the time interval from about 6 to 16 see are interpreted to be the higher
surface wave modes, as may be seen in the same Figures 3 and 4.
I t may be difficult and perhaps useless to draw a clear line in the RES separating
surface waves from body waves. This is clear from the Figures 3 and 4 where it may
be seen that some surface-wave phases (in RES from 6 to 10 sec and for frequency
co < 8 rad/see, Figures 3 and 4) do arrive simultaneously with the higher frequency
energy (in RES from 6 to about 10 see and for co > 8 tad/see, Figures 3 and 4) which

TABLE 2
LAYER PROPERTIES
Layer Thickness (kin) a (km/sec) fl (km/sec) ~" (gr/cm~)

1 0.18 1.70 0.98 1.28

2 0.55 1.96 1.13 1.36
3 0.98 2.71 1.57 1.59
4 1.19 3.76 2.17 1.91
5 2.68 4.69 2.71 2.19
6 ~ 6.40 3.70 2.71

we may call S wave. The arrivals after 22 see (Figures 3 and 4) cannot be interpreted
in a simple way and probably represent reflections which traveled along the indirect
paths and hence arrived later at the recording station.
In the light of these interpretations, based on the calculated arrivals of surface
waves, we can conclude that, except for the high frequency arrivals from about 6 to
8 see (Figures 3 and 4), all other motions recorded during E v e n t 9 are represented b y
surface waves. The first burst of higher frequencies which we may call the S-wave
arrival is the prominent feature of an accelerogram in almost every case. However,
from the RES, one can find that, from the point of view of the vibration amplitudes
of the wide class of simple oscillators, surface waves may be more significant. During
surface-wave accelerations, structures may experience the longest and biggest oscilla-
tions. In Figure 3, for example, the response spectrum peak would occur for a fre-
quency 15.3 rad/sec at a time t = 13.25 sec which is definitely in the region of surface-
wave excitation.
The duration of strong earthquake ground motion is one of the most important
parameters governing destruction of man-made structures. In terms of duration, the
body, P and S, waves are recorded on most accelerograms in the form of relatively
higher frequency bursts of energy, arriving at the station for several seconds. The
surface waves, on the other hand, usually composed of somewhat longer period waves,
arrive during considerably longer time intervals. This extended arrival of the energy
in the form of the surface waves is the consequence of the dispersion through the hori-
zontal waves guides which are usually characterized by the low velocity of wave
 eL C~M~G EVE~I S ~S PE~ ~LbE is i~,~
~ELAII~E ~ESFCNS~ ~CCEL~R. =NVbLLPE Sp~LI,~M. RbSP~AS~ ApPLJI&CE I ~ C~. IS L ~ I ; I ~ C ~Y ~bLIIPLYIAC
SPELI,~L V~LbE bY I p E PE~K ~ALUEo ~AMPII~G IS O.ZOO CF I h E CRIIICAL

F R ~ ~ U £ ~ C ~t - R A L i A ~ b I S E ~ L i~ ~
o.I ,]

CCC~CC~C~G4~CC~C~Cc~C~c~c~G~C~c~u~c~c~C~GC~¢~CGCC~C(~C~cG~c -- O

cc~coc~cc~ccc~cccccccccccco~ccc~cccccc~ccc~acc~ccccccc¢c~ccccc~cccc~cc~accccc~c~ec~cccc
~ c ccceccccc ccc~ccccccccccc~ee~ac~cce~c~oc~e~c~e~c~coe~¢o~c~cc~cc~cccc~c~Gcc~o~c~c~cccoccc
c c c c c c c c c c c c c c c o c c CCCCCCc~e~eOCOCCCuCCOOaOoOOOCeOCCOa~GocceecCOooCccOCCGCCCCeCCOOCCaOCa¢ CCCC~CCCCCCC
c ccccacacccccccec cc~ccccccc~a~c~cc~cc~cc~ca~occ~eco~ccc~z~cGc~a~cc~occ~c~acc~cc~cccu~c
c c c c o o ¢ CCCCCCC CCC¢ ~ c c ~ c ~ c c c c ~ c G c c c ~ c c ~ c c ~ c ~ c c c o ~ a ~ a ~ c ~ c c ~ a c c ~ e ~ c c ~ c c c ~ a ~ e ~ c c a ~ c c ~ c ~ e c c c ~ e c c
~ CCCOOOC~ccccCCCCCCCCcc CcccceooCUCCCe~eOaaCOOeCCOOCOoua¢oaOCCCCaoCC~COCCOnCCOCOOaOCCCeCeCCCCeOcoeoCCC
c~ccc~cc~ccc~cc~c~cccc~c~ca~cec~cc~cc~cc~cc~ec~cc~ooc~ecc~c~e~c~ccc~ce~cea~cccce~cc~eccc
c c c c e e c c o c c o c u c c c e CCCC C~cCO~o~OaaeCaOCOOCCCOC~COOCOOOOO~eCCCCCCOOCOaaCCOOCCCCOOCZCCOCOOCCOCOoCCCO~a~
CCCCOCC~c¢ ¢C¢CCCC¢C~CC CCCCCC oeCOCCCeCCeOCCGOOeCCeC~OOOCOe e c c c c c c e c o c e c e o c c c c c e e e o c c e c ~ e c c o c c c c c c e c c c
~ CCCCOCCCCeCCCCCCCCCCCC CCcCOoeOCaCCCOCCr~OCCOaa~OOCOOOCCOOCOCCOOOCCOC~CCOCCC ¢ C c c c o c c ~ a c o c c c c c c c c c a o o c
CeeCCOCCuCCCCC COcCOCGCCCCOUO~aeeOCCCCCeCCCCecacoecoeoeuoC cCcCCCOeaOG~CCCCCCCCaCeeccOOCOC CCCCOC~CCC~C
CCCCOOCCCCCCCgCCCCCCeCcCCOaO a u c a c c c c o c c c G ~ e ~ c o o o c c a o a o c ~ u a a o o c c o c c o c o c c c c c c c c c c c c c c e c c o c c c c c u c c c o c c c

--5

--I0

--15

--20

E~Ec~C~ccccc~ccc ~cco / ooocGoooooo cccccGc~ccccccc~ccc

ccccoococccc~cccccc cc¢¢ccoc -- ooooccooocoocooooo ccoccc
CCCCCCCCCCCOCCCOGCC~C£CCO00 / OCGOOCOOOCCaOOOOOCOJOG COCCCCCCCC
22 OCCCOOCCOOCCCCCCCCCC (CCCOC CCO00CCO0000GO000COCOC OCO000COGCOCCOOOOCCCOCCCCGCCCCCCOCOC
OCCCOOCGCCCCCCCCCCCCC~.G~O OCO00eCOOCCCOCCOOOOOOOCCCOOOOOCOCOGCOCCCCCOOOOCCUOCO¢COCOOCCCOCCC
CCCCO0000CCCCCCC(CCCCCCCCOCO0000000GOCO00COOC~OOCO0000000000 OOOCOCOCCCCCCOCOCOOCCOCCOCCCCC~COCO00C
C¢CCCCCCCCCCCCCOCCCCCC(CCCOODCOCOC CGCCCCOOCCOGCGGGGC¢ OCCOCOCCCCCCCCCGGU~OOOO0~CCCCCCCC~COC
CCCCOCCCOCCCCCCOCCCCOCCCCCOO00QOCC 000 COOCO OCCO
CCCCOCCCCCCCC(CCCCCCCCCCCCCGCCOCC CC
CCCCCCCOOCOCOCCCCCCCCCCCCOCO0000CC CCCCOOCOCO00C
24 OCCCCOCCOCCCCCCCCCCCCCC60OO 0 0 0 0 GOGCOOCCO000COCCCCCOC~OOCOO000COCOOCCCOOCC
cccccccooocoocoocccccccccoc ooo cooccocco¢ccoco
CC(CCCCC~CCCCCCCC£CCCCCOCG C GCCCCCCCC¢OCCCG
occccoccccccccocoocccccoG ooocoooo ooo
25 CCCCCGCCCCCCCCCCCCCCCCCC Oooco 11111111 . . . . . . ,, • --2_5
cccccocccccccccoaccccGco ~ ~
CCCCCOCCCCCCCCCCCCCCCCOC 0 L) l ~ l l l l l ] ~ OCCCOOGOCCOCCCCGCC
cccccoccccccccccccccccccc o cccc c
z~ cccccoccccccccccococuccoo ococoooo~o o~ocoo
cccococccccoocccccoccccoo ooooccccoccooo CCCCGcocccoccocccccccccoooc
CCCCOOCCCC(CCCCCCCCCOCCGO OOCOOCGCOOCC~flOUOOOCOOOOGOOCO000 GCOCOUCCOCGCCCCCOOGCCOOCC
cccccoccccccccccccccccoco ococccoccoocccoooc c c c o c c c c c c ~ o o c ccccoooccoooocccccocccococ
27 cccccoccccccc~(cccolccoco ~COOGCCOOOJ G
lf OGOCC CCCOCCOOCGCCCC
CCCCCCCCCCCC~CCCCCCCCCCO0 CCCCCC ~
OCCC~OCCCCCCCCCCCCCCCCCCCC CGOOC ~ CCOCOC
OCCCCOCOCCCCOOCOOCCCCCCOOC OttO00 OCGCOOOC"
~ cccccoccoccccccGcccccccoccc cccccccc cocccccccccco~c
occcooccccccccccccccccccccccoocococooco o¢oococcc
ccccooccccc~ccccccccccccooooooooccoocco oucoccccccco
GCCCCCCCOCCOCCCCCCCCCCCCCCCOGOCOCCCOCC CGOCCCOCCCCOC GCCCCCGCCCCCCCCCCCCC
~ C(CCC~CCCCCCCCCOCOCCOCCCCCOOOO~OOGO00C OO00OOOOCGCCCCCCOOOCOCOOCOOCCCCCCOOCCOCOOCCCCOOCCCOCOC
occcoococccccccocccocccoccccooco&cooooo COCGOOOOCOOOCCCCOCCCCCOOCOOCOOCDOOO COCCCCOOGOCUCGC
CCCCOOCCCCCOCCCOCCCCCCCCCCCCU~CCCCCOCGO0 OCCCCOCCOC~ OCCCCCCCC CCCOOCC
CCCCO~CCCCEOCCCOCCC~CCCCCCCO0000CCCCCCOOOC OOCCO00COO oocc
~c ccccc~cccccccccccccccccccccoooccccooocuoooooo CCGGGGC --30
c(ccoo,~coccocccccccccgccoocooococcgocco~flcoo ~ OOCCC OCGCGOOC
C£CCCCCCO(CCCC~CCCCCCCCCCCCOG~CCCCCGCCGC
C
lfO0 CCCCOOCCCCCCCCCCCCC
CCCCOCCCCCCGCCCCCCCCCCCCCCCCOC CCCGCCO00CO0 oooooooo~o
~z ccccooccccccccccccccccccccogoc ccocc~c~cco cccoccoccc

R.E,S, NORMALIZED AMPLITUDES

0 I 2 3 ,4 5
F~e. 8. The c~2A(t, c~) Response Enve]ope Spectrum (P~ES) for E v e n t 9, NS component,.
350
 EL C~D,TRCJ 6V~hT ~ eN PE~. '~LbL IS I.~,eO

)ELAIIV~ I ~ $ ~ [ h S ~ ACCELEH. EAVI.LCF6 5P~CTHLM. RESPGhS6 AMPLLFbL~ IN C~. IS 661~1{N61~ uY #UIIIPLYIh(~

SPECIR~L V@~.~£ EY l h ~ p E ~ VALUE. ~A@~P~G IS C*lOO U~r THE CHIIICAL

h R E ~w b i: h C y - R A C I A h S i S E C CN C o.I
x ~ ) ~ 5 6 ) ~ ~ 10 11 12 13 i~ I~ le

ccccocccccccc¢ C c ~ c ~ c ¢ c c ~ c ~ c ~ ¢ ~ ¢ ~ G ~ C ~ c ~ ¢ ~ c G c ~ C ~ C c ~ c ~ c ~ c c G ~ c
-0
G C C ~ c ~ c c ~ c ~ C c c ¢ ~ ~ ~ u ~ ¢ ~ ¢ ~ C c ~ ¢ ~ G ¢ ~
cccccocc¢cc¢cc C~c~cc¢¢~¢~c~cCccc~¢c~c~cc~G~c~c¢~c~Ccccccc¢~c¢cc~c~ccccccc~cccc
cccc~ccccc~c¢cc~ccccC¢cc~c~c~c~ccc~cc~cc~c~c~c~c~c~c¢~c~¢~cc~cc~c~
~ c ~ G ~ c ~ ¢ ~ ¢ ~ c c ~ c ~ C ~ o ~ c ~ G ~ u c ~ c u C c ~ c ~ c c ~ ¢ ~ ¢ ~ G ~ c ~ c c ~
Cccc~c~cccc~cCcGCcc¢¢¢cccc~c~cc~c~c~cc~¢~cc~c~¢c~c~c~c¢~ccc~¢~c(l~ccc~ccc¢~
2 c cccocccccc¢c¢ccoccc¢ cccoooo~ooocoooc~ooocoooooooooooo¢ooooooccoocooocco¢ccecoccuocooooocooooovGoooo
ccc~¢cc~cccccc~Cc¢cCc¢ccc~c~¢c~c~cc~c~c~cc~c~cc~c~c¢cc~a¢~e~cc~ccc~
c c c C ~ c ~ ¢ c c ~ c ~ c ~ c ~ c C c c c ~ c ~ c ~ c ~ c ~ c ~ C ~ c ~ c c ~ c c ~ c ~ c ~ c ~ c c ~
~ccc~cc~cc~c~c~cc~cG~u~c~¢~c~c~c~c~¢~¢~c~c~c~c~cc~c~c~ccc~cc
~c~c~Cc~cC~cc~cCcc~c~C~c~c~c¢~G~G~cc~`~¢~cc~c~c~c~c¢~c~c~c
c c c ~ ¢ ~ c ~ c c c c C c c C ~ c G c c c c C ~ c ~ c ~ c ~ ¢ ~ c ~ G ~ c ~ c ~ c c ~ c ~ c c ~ c ~ c ~ c ~ c
~c¢c~C~G~c~c~c~cccc~cc~c~c~c~c~~c~c~cCc~c¢¢¢~G~cc¢c~c~¢~
¢¢c¢~cc~cc~cCc~ccccccc¢G~ac~ccc~caacc~c~cc~a~ccc~c¢~cc~ca¢c~cc¢cG~c~c~c~cc~c~cc~c~e
~, ccccccG ccoccc c c c c c C c c c c ~ c ~ c ~ c ~ ~ c ~ c ~ c ~ c ~ c ~ c G ~

-5

--I0

--15

OCO~OOOCGGCQC (CCCCCCCCCCO00 0000 ~ l l l Z ~ 1 1 1 1

~ccccooooocooo,~,ccoocccooooo ooooooo o o o o ~ o i l i i ~ ii i i i 11 i oo
-20
cccccocoocccc¢c~cc ~ c c c c c coo~cccoccooc
c c o o c o cc'ococo OGG
ocoooccooccococc~cccccoooeooococoGoooo~o~-'i~oooooo coooooooo
~1 occoooccccccccccc ccoccocococococcoooGemoo
ocooooccocco¢ccoc~ ooccoccccocc occccoooooco~c
cccc~ocoooocoocc~,ocooo oooooo¢ ooo~cco~c
Q¢ccooccccccocccocc~ccccccoocococ.c~occoooc occccoooo ooooo¢o
O£CO00COOCCOCCCOCGC CCCCCCO00OC~4~GOOOOOOOOOQOQOOGQOUO OCOCO0000CO000 OOOCOOOOCCOOOOOOOO
a2 ccccccccoococcccocco~ococco~occo~ooooocooooooooco oooooooooccccocccooooo
occcooccocoooccccccc cccooo ooooococoooccocccoooo
occcooccocccccccccccoc~oococococcooocoooooooooooo cooooo
ococooccccccoccoccccocccccooooooocooocoo~ocoooouooooo coooccoocococcccccocoo
~) occoocccoccooococooooccGoooooo¢oocooocoooooooooooooooo oooooooooocccoccooocccocoooo
~c~cc~c~cc~c~cccc~cGc~cc~c~cc~cc~C~cc~cc~cc~c~cc~c~ccc~d
O¢OCOOCOOCCOCOCCCCCOCCCCCCQOOCOO¢COOOOOOOCOOOOCOOOCOOOCOOCOOGOOOOOCOOCCOOCoOCOO oocoo
~¢¢c~c~c~c~c~cccc~¢~c~¢~c~c~cc~c~Cc~c~c~c~cccccc¢c~c~cc~c~
2~ c¢c¢c¢~c~c¢~¢ccc¢cc~c~c~¢~cc~¢cc~cccc~¢~ccc¢c~¢~ccccc~¢~
~ c c ~ c c c c ~ c c ~ c c ~ c c ~ c ~ c ~ c ~ c ~ c ~ c ~ c C c ~ c c ~ c c c c ~ c ~ OOC
CCOeOOC¢CCCCCCCCCCCOC¢CCCOCOOOCOOCOOOOOOOOOOOOOOOGO000000 O0000C 0
CC¢COOCCCCCCC¢(COCCC¢CCCCCO000COCOOOOGOOOCOOGOOOOCOOOOCOOGO00CO0
~ ~ c c ~ c ~ c ~ c ~ c ~ c ~ c ~ ~ c ¢ ~ c ~ c ~ c c c ~ c ~ c ~ c ~ c c ~
OCCCCOCCC(CCCC¢CCCCCCCOCCO0000000COOCCOOCCO000000~O000CO00COCOCOCCOCOOCOOCOC --25
occccccccccocccooccc~cccocoooococccoc¢ooocooooo ~coocc o
~CCOOCCOCCOO¢¢COOOCOCCCCOOOOOCOOCCOOCGOOCOO COOOCOOOCC
2e ¢coc¢ocGcccocccocccccccc¢¢oooocoooooooooooo Q o
ccccooccoccocccoccccoccc¢ooooococccocoooo oocc~ccoc
cccccoc(cc¢¢¢ccooccoocccoooouoooooo~ocooo oooooocccccocoocoooooo
occccoccoccoccccoccc¢coooooooococooooooooooooooocoooooooocoocooooco¢occoooo
~ ccoccocooccoooccoc¢¢occcccoooocooco¢ocooocoooocoooooo cooc.~ooooocccc~oc¢o o~ooooo
occooocoo¢~ooococcococoooooooooocoooocooccccoooooo coocooco cooocooooc
CCCCCOCCOCCCCCCCOCCGCCCCCCOCOOCOOCOOOOOOOOOOOOOOOOOOOCC
ccoooocooccccc(occccccccccooooccccoocoooocoooc¢oococoooooo ~
o
~a ococcocooocooc¢ooccoooccoooooooo~oooooooccoooocoocooo¢oooo
ocOC~OCOcCCOCGCOCcocccccccooo~cocCCCOCOOOCOOOOCOOC~OOCCO~O
occcoococcccccccccccoccccoooooco~cocccooooooo ooooooooooo
oocoooccoocooococccooccooooooooooooooooooooo ooocooooccccooooooooo
2~ ccccccccocccoccoc(ccccccccoocococcc~ccooo cccccccooocoooooccccooccco©oc
cccccooccccococ~cccooccccoooooooooo oooooocoooooooooooo o ococoooooooocooooooooo
cccccoccocccccccccccccccccocooooo ooocoooccooooo cccccccoccc coccc
ccccoocooccoccccccccccccccccooco ooooocoocooooo oooooo~ooocc~cooco¢oooc
~o ococooccocccccccoc¢cccoccooooooooo
occcoococcccccccocccccccoocGooco~cooo
ooooocoooccooooo
ocoooo¢oooooooo
oco~ooooooooooooo~gg~ooo
ooococccocoooooooooocc¢c~o~oooccc --:30
cccccoccoccoocccocooocccccoooocooooco o¢oooococcooooooccoooccccccccoGooooooooococooccooooo
OCCCCOCOGCCOCC¢CCCCOOC¢CGGO000CO000CO OOCO00CCOOCO00COOCCOCOOCCCCCCCCCOOCC¢O00C¢OCOCOCO00¢O
~ OCCO00CCCCCCCCCO00CCOCCOCCO000CO00000 O000000COOOCOCOCOOCOCOC¢C¢CCOCOOOOCO00COCCCO00CCOC~CC

R. E. S. NORMALIZED AMPLITUDES

0 I 2 3 4 5
FIG. 4. The co2A(t, ~o) Response Envelope Spectrum (RES) for E v e n t 9, E W component.
351
 352 BULLETIN OF THE SEISMOLOGICAL SOCIETY OF AMERICA

N~LAIIVE RESFChSE ~CCEL~, E~V~LLP~ ~P~CIHbN, RE~PChS~ AbPLIIUD6 I N CV* [S [EI$1kED BY PbLIIPLYXKG
SB[CI~AL ~AL~[ BY 3H~ ~ E ~ WlLUh. D~gi~G i~ ~ * k ~ OF I~E C~IIICAL
0.3
F R E Gb E h C y - R A C i A N ~ / S E C 0 N C r
Z Z . ~ 5 6 , . ~ io - zz , z* m5 .
-0

4
J

-5

j 1

ZO -I0

IX

12

1.3

~. - 15

~o --20

21

occcooccoccccc cccccccoG~cooocooooo cooccooocoooooooooccococoooo

cccccocccococcc(k~ccooc~cocoooccoooooo oooooo~oooooooooccoo¢oooccoccooeooo
23 ¢OCOCOCGOCCOCOCOC~C~CCOCCOOGOOOOCOOO0 CO O00000GOOOCCO000000000CO00¢OCC~OOC
cc(ccccoccooccocr'~oocrcccoooocococ~oooocoooooooooooocoocooccoooo occo
cccoocooccoocoocoo c¢occoooo occooooooooooocoooooooooocoeocooocococccoccoco
ccccocooccooccoc6 ~co~ooooooo~ooooeooocoooooooooooooooooooocoo¢ooococoooo oooo
ccco~cccccco~cocc~cocoo~ooooocooooooooooooooooooooooooooocoo~ooooo
cccoococccoccccccc~cccc~ccooo~occcooooooooooocoocoooo~oo~ooooooccocccoocc
z~ cccooc~o~coo¢coccccc~o~oooooocoooooocoooo oooooooocoococoooooooo ~ ~*,~,,t~,izzzzzl -25
(CCCCCCOCCOOCCOCCCOCCCO0000DOQOOCCO000000 oo Illll lllllllll)Xki)
CCCOOCCCCCCOCCGCCCCOC~CCCO00000CCCO~O0 ~OOCOOO~CGCO00 Lkkzl
cooooccoccocccccccooococoooooooo oooooooooo zz
26 CCOOOCCCCCOCCCOC£CCCC£CGOGO00C o OOOOOO000000COOCCO0000000006 ~ ~r~Zl
CCOOOCOGCOGOOOOOCOOOCQOCO0000OO O000000COOOOOO llllll

Z£(t£~Ctt£tCCCt~£CCCOCC6 o0 T I'
ccccoccccccoe¢ooccoocooooooo ooGooooooooooo
ccccccccccccccccccccocccccoo ocooccoooooo ocooaoocoeooccc~oc
28 ccoccooocoooocooo¢ooc~ooooo oooooooooooooooooooooooo oooooooo ooooooococooooooooo
~cccoccccccoccccccccccccooo oooccoooooo~ooooooo~ooooooocooocooocoocoo oo~occoocoo
occcoccocccooccooccccccooooo ooooocooooooooooooooooooooooooooo¢oooooocooooooooooocooooooooooo
oo~¢o~o~cc~ccccoccooc~cooo ooooooooooooooo~oooooooooocooocococcoccoocooooccooooccoooo~oooooo
2S ccccocccccccc(ccccccccoooooooo ooooooooooooooooooocoooooooooocoooccoccoocooooooooooccc~oocoooooo
cccoocooccooocooooocccoooooooooooooocooocooooooooooooooooooooooo¢ooo~cooc°°c° ooooo
cccococcocccccccocccoccc~ooooocoooooocoooccoooooccoooocoooccooooo ooooo
CCCO0000CCOOCCOCCOCOCCOOOO000COGGO00000000000000CO000OOOOO000000 OOO000000
30 occccocooccooccoococccooooooooooooooooooocoooooooooooooooooooooooo
OCCO00CO00COOOCCOCCCOCCOCO0000COOCO000000CO00000000000000000OO0000
C~OCC O00000000CO0~CO00OO
O0000000GCOCO0~OOO000
--30
~coooo~occcccc~occcccc~cocooooooocooooooooooocooooooooooooooooo O00000COCOOOO00000
ccccoocoocoooccoccoocc~ooooooocoooooooooooooooooooooooooooooocoo
~z ccccooccoccccccoccccccccooooooooo~ooGcooocoooocoocoooccoocoooooooo

R. E.S. NORMALIZED AMPLITUDES

0 I 2 3 4 5

FIG. 5. The co2A (t, ~) Response Envelope Spectrum (RES) for Event~ 1A, 1B, and 1C,
NS component.
 INTERPRETATION OF STRONG EARTHQUAKE GROUND ~[OTION 353

~L CENTRD MAIN EVENT EW pEAK VALUE IS 89.78

RELATIVE RESPONSE aCCELER. ENVELOPE SPECTRUM, R~$PDNSE ~*#4PLITUDE IN CM. IS OBTAINED By MULTIPLYING
5PECTR~t VALUE BY THE PEAK VALUE, D~MPING IS 0,).00 OF THE CRITICC~L
03
F R E Q U E N C Y -- R A D I A N S ] $ E C 0 N O

--0

--5

--I0
L

--15

--20

o¢oooooo~coooc~oo cooooooooooooooooooo ooooooooo ooooo o

oooooocoooocoo~cooooocooocoooooocoo oooo oo ~ ¢I
z3 oooooooooooooooo~ooooooooooooocooooo~o ooooooooooooc oocooo "~.
oooooooooooooooooo~p~oooooooooooooo o~'oo ocooooooooooooo oooooooooocooonoooo ~o
ooooooocoooooeoooo~Koooooooooooc~oooooooooooooooocoo ooooooooeooooooocooooo
O00000GOOOOGO0000 CO000000000 O00000000000GODO00 O00000000000OO0000000CO0000
2¢ coooooeoooocoooooooo),4~oog~%oooooooococoo ocooooooo ocoonoo ooo0oo
o...................oo~ooo . . . . . . . 22Z~22Z2Z~Z2~2
oooooovo~ocooooooooooooooooooooooooo ~
oooooooooooo°oooooeooooooooooooooooooo 111~I ~ P
zs o~oocoooooocoeoooooooooooooooooooooooooooooooco oo~o0oo --~
oooococooooooooooocoooooooooooooo oooooooooo ~ ~ ~,7"rrr~,'~,~ A
ocooooocooocooooooooc~ ~ - ~ ~ ~-#-rrm"rr1~ ~....... Ii
ooooooooooo~oooooooo
ooo~ocooocoooo ooo ocoo ooooo
occoooooocoooo ooo oooocoooo oooo
cO
ocoocoooooooooo ooo ooooo oo IH, ......... ,,,,,,,

ccooo~cooocooooooooo oo ~ oooooo O
o~oooooooccoooooooooo ooooooooo~ ~ oooooooo oooooooo
zB ocoocoocoocooooooooo ooooooooooooocoo ooooooooooo ooooooooooeooo o
OOCGO00~O000CO000000 O00000000000000OO0000QO~OnOOOOO00000000 0000000000000000000000 u)
oooooooooooooooooooo ooooooooooooooooooooocooooooooooooooooooo~ooooooooooooooocoooooooooooo
ooooo°ooooocooooooco o~oooooooooooooooooooooooooooooooooooooooooooooooooooooooo~eoooooooooooo
z~ oooooooocooooooooooo oooooooooo~oooooooooooocoo~oooooooooooooooooooooooooooocoooooooooooooeoo
onoooooooocooo~oooooo oooooooooooooooooooooooooc~ooooooooooooooooooooooooooooooooooooooooooooooo
oooooooocooooooooooooo oooooooooooooooooooooooooooooooooooooooooooooooooocooooooooooooooooooooooooo
~ ~ ~ ~ ~ ~ ~
3O O00000000000000000CO000000 O000000noOnO00000000000000000000000000000000000000000000000000oO000000000 - - 3 0
O0000000CO00OO0OOCO000000 O0000000000OO0000000000000000 o
ocooooooccooooooooooooooo ocooooooooooooooooooooooooooooo
ooooooooooooococooooeoooooooooooooooooooooooooocooooooooo
~ ¢ooooooooocoooooooooooooooooooooooooooooooooooooooooooooooo ooooo

R. E .S. NORMALIZED AMPLITUDES

0 I 2 3 4 5

FIG. 6. The ~,2A(t, ~) Response Envelope Spectrum (RES) for Events] 1A, 1B, and 1C,
EW component.
354 BULLETIN OF THE SEISMOLOGICAL SOCIETY OF AMERICA

propagation. Dispersion curves in Figure 2 can serve as one example of this. The sur-
face-wave motions characterized by such a dispersion will thus last for

(1
A V rain
1)
V max

time units, where A is the epicentral distance and V is the group velocity.
The main part of the energy release during the Imperial Valley, 1940 earthquake
was recorded on the accelerograph at E1 Centro for about 25 sec (Figures 5 and 6).
The approximate analysis of this record (Trifunac and Brune, 1970) indicates that
at least four distinct shocks (Events 1A, 1B, 1C, and 2) may be associated with indi-
vidual bursts of S-wave arrivals on the NS and EW accelerograph records at about 2,
8, 13 and 23.5 sec. Corresponding P-wave arrivals for each of these events cannot be
unequivocally identified because of the high level of short-period vibrations produced
by the previous events. For this reason, the estimated distances based on the S-P
times may be in error. In this study, we interpret RES for these events in an attempt
to find the best fit of the observed and expected arrivals of surface-wave energy. The
results of this fit are shown in Figures 5 and 6 which indicate that our previous esti-
mates (Trifunac and Brune, 1970) for the Events 1A and 1C were probably within
the correct range. The tentative interpretatiou based on the surface-wave arrivals
indicates that Event 1B could have occurred within a smaller epicentral distance
than estimated previously. Surface-wave arrivals for Event 2 are not interpreted
here separately. The epicentral distance for this event is probably 35 to 44 km which
places it at the southernmost end of the observed surface fault. For this distance,
fundamental Love and Rayleigh modes would arrive more than 50 sec after the strong
motion was triggered and would be mixed with the arrivals from Event 3. If we as-
sume the epicentral distances for Events 1A, 1B and 1C to be about 10, 12.5 and 15
kin, then the expected surface-wave arrivals are as indicated in Figures 5 and 6.
Although the accuracy in the above interpretations depends on a knowledge of the
correct epieentral distances, the applicability of our theoretical dispersion curves,
and on the properties of the RES, the most important tentative conclusion of this
work seems well established. This conclusion is, as already mentioned, that the pre-
dominant part of strong earthquake ground motion consists of surface waves. The
form of these waves is considerably more complicated than for teleseismic surface
waves because of the proximity and size of the source of the seismic energy release
and the frequency domain considered. In our interpretation these waves are surface
waves in the sense that their arrivals are governed by the travel-time-distance rela-
tionships based on a simplified model for surface-wave dispersion.

CONCLUSIONS

The main results of this study can be summarized as follows:

The Response Envelope Spectrum (RES), a multiple filter technique, represents a
physically simple and informative tool for the analysis of strong earthquake ground
motion accelerograph records. It contains all of the information necessary to deter-
mine the relative response spectra and is a logical extension of response spectrum tech-
niques. The most important property of this method is that it displays the time de-
pendence of the amplitude and frequency changes in the accelerograph record. It is
particularly useful in detecting arrivals of different wave forms and may be used as a
 INTERPRETATION OF STRONG EARTHQUAKE GROUND MOTION 355

tool in dispersion analysis. Another important property of RES is that it conveys the
information about the length of large amplitude response of structures.
The analysis of several strong-motion accelerograph records using the RES indi-
cates that the predominant part of the strong earthquake ground motion, generated
by a near-shallow earthquake energy release is represented by surface waves. This is
important from the point of view of destructiveness and also from the point of view of
the establishment of the correct simplified wave-propagation models for studying the
effects of local geology. The analysis indicates that the duration of the intense shaking
will be mainly governed by the dispersion properties of the ground.
ACKNOWLEDG3/IENTS
I am indebted to Professor D. E. Hudson and Mr. F. E. Udwadia for a critical reading of the
manuscript.
This study was supported in part by the National Science Foundation grant (Engineering
Mechanics Program) and the Earthquake Research Affiliates Program at the California Institute
of Technology.
REFERENCES
Alexander, S. S. (1963). Surface wave propagation in the western United States, Ph.D. Thesis,
California Institute of Technology.
Amin, M. and A. H. Ang (1966). A nonstationary stochastic model for strong-motion earthquakes,
Structural Research Series No. 806, University of Illinois, Department of Civil Engnieering.
Archambeau, C. B., E. A. Flinn, and D. G. Lambert (1966). Detection, analysis and interpretation
of teleseismic signals, 1. Compressional phases from the Salmon event, J. Geophys. Res. 71,
3483-3501.
Barstein, M. F. (1960). Application of probability methods for designing the effects of seismic
forces on engineering structures, Proceedings of the Second World Conference on Earthquake
Engineering, Japan, 1467-1482.
Benioff, H. (1934). The physical evaluation of seismic destructiveness. Bull. Seism. Soc. Am.
24,398 403.
BiehIer, S. (1964). Geophysical study of the Salton Trough of Southern California (Ph.D. Thesis),
California Institute of Technology.
Biot, M. A. (1941). A mechanical analyzer for the prediction of earthquake stresses, Bull. Seism.
Soc. Am. 31,151-171.
Bogdanoff, J. L., J. E. Goldberg, and M. C. Bernard (1961). Response of a simple structure to a
random earthquake-like disturbance, Bull. Seism. Soc. Am. 51,293-310.
Bolotin, V. V. (1960). Statistical theory of aseismic design of structures, Proceedings of the Second
World Conference on Earthquake Engineering, Japan, 1365-1374.
Bycroft, G. N. (1960). White noise representation of earthquake, J. Eng. Mech. Div., Am. Soc.
Civil Engrs. 86, 1-16.
Cornell, C. A. (1964). Stochastic process models in structural engineering. Technical Report No.
84, Stanford University, Department of Civil Engineering.
Dziewonski, S., S. Bloeh, and M. Landisman (1969). A technique for the analysis of transient
seismic signals, Bull. Seism. Soe. Am. 59,427-444.
Housner, G. W. (1947). Characteristics of strong-motion earthquakes, Bull. Seism. Soc. Am. 37,
19-31.
Housner, G. W., R. R. Martel, and J. L. Alford (1953). Spectrum analysis of strong-motion earth-
quakes, Bull. Seism. Soc. Am. 43, 9~119.
Housner, G. W. and P. C. Jennings (1964). Generation of artificial earthquakes, J. Eng. Mech.
Div., ASCE, 90, 113-150.
Hudson, D. E. (1956). Response spectrum techniques in engineering seismology, Proceedings of
the World Conference on Earthquake Engineering, Earthquake Research Institute and the
University of California, Berkeley, 4, 1.
Jennings, P. C., G. W. Housner, and N. C. Tsai (1968). Simulated Earthquake Motions, Earth-
quake Engineering Research Laboratory, California Institute of Technology.
Landisman, M., A. Dziewonski and Y. Sato (1969). Recent improvements in the analysis of surface
wave observations, Geophys. or. 17, 369-403.
 356 BULLETIN OF THE SEIS~[OLOGICAL SOCIETY OF A~IERICA

Merchant, H. C. and D. E. Hudson (1962). Mode superposition in multi-degree of freedom systems

using earthquake Response Spectrum data, Bull. Seism. Soc. Am. 52,405-416.
Rosenblueth, E. (1956). Some applications of probability theory in aseismic design, Proceedings
of the World Conference on Earthquake Engineering, Berkeley, 8, 1-18.
Rosenblueth, E. and Bustamante, J. E. (1962). Distribution of structural response to earthquakes,
J. Eng. Mech. Div., ASCE, 88, 75-106.
Shinozuka, M. and Y. Sato (1967). Simulation of nonstationary random processes, J. Eng. Mech.
Div., ASCE, 93, 11-40.
Tajimi, H. (1960). A statistical method of determining the maximum response of a building struc-
ture during an earthquake, Proceedings of the Second World Conference on Earthquake Engineer-
ing, Japan, 781-797.
Trifunac, M. D. (1969). Investigation of Strong Earthquake Ground Motion, Earthquake Engineer-
ing Research Laboratory, California Institute of Technology.
Trifunae, M. D. and J. N. Brune (1970). Complexity of energy release during the imperial valley,
California, earthquake of 1940, Bull. Seism. Soc. Am. 60,137-160.

CALIFORNIA INSTITUTE OF TECHNOLOGY

I~ASADENA, CALIFORNIA

Manuscript received August 31 1970

View publication stats