Site-Response Maps for the Los Angeles Region
Based on Earthquake Ground Motions


Open-File Report 96-723

Version 1.1

This report is preliminary and has not been reviewed for conformity with U.S. Geological Survey editorial standards. Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the U.S. Government.

Denver, Colorado 80225-00046

Abstract.--Ground-motion records from aftershocks of the 1994 Northridge earthquake and main-shock records from the 1971 San Fernando, 1987 Whittier Narrows, 1991 Sierra Madre, and 1994 Northridge earthquakes are used to estimate site response in the urban Los Angeles, California, area. Two frequency bands are considered, 0.5-1.5 Hz and 2.0-6.0 Hz. Instrument characteristics prevented going to lower frequencies, and frequencies above 6.0 Hz are less important to the building inventory. Site response determined at the instrumented locations is associated with the surficial geology and contoured to produce a continuous spatial estimation of site response. The maps in this report are preliminary and will evolve as more data become available and more analysis is done.


Site response, or the site specific amplification of ground motion as a function of frequency, is a fundamental component of earthquake hazard assessment. Site response is also a desired component of any regional seismic hazard map. Previous site response studies in Los Angeles include Campbell (1976), who showed that damage during the 1933 Long Beach earthquake was greater for sites on unconsolidated soils than at sites on consolidated soils, for a given distance from the Newport- Inglewood fault zone. Rogers and others (1979) calculated site response in Long Beach using spectral ratios of Nevada Test Site (NTS) nuclear events. They found amplification factors as high as 11 in the frequency band from 0.16 to 5.0 Hz for alluvium-to-rock ratios. This work was continued by Rogers and others (1980 ,1984, 1985) to include additional sites in the Los Angeles area. They confirmed that low strain records from NTS nuclear events gave similar site response estimates as strong-motion records from the 1971 San Fernando earthquake. Their results show elevated site response at frequencies greater than 2.0 Hz with increased mean void ratio in the near-surface layers, unconsolidated sediment thickness, and depth to basement rock. At frequencies less than 2.0 Hz, the most significant factors were reported to be thickness of Quaternary sediments and the depth to basement. The 1994 Northridge earthquake greatly increased the available ground-motion data for the urban Los Angeles area from which site-response estimates can be made. Recordings of this earthquake and its aftershocks have also greatly increased our understanding of the variabiliy in site response over widely different length scales; from one city block to across entire sedimentary basins. In addition, instrument deployment for the Northridge earthquake targeted several specific areas of severe damage and compared them with low or nondamaged areas. This expanded data set has led to the feasibility of creating site-response maps for the Los Angeles area.


The data for this study are recorded earthquake ground motions. We use both aftershock and main-shock records to increase the size of our database. Deployment of portable digital recorders after the 1994 Northridge earthquake yielded an extensive collection of aftershock records for a wide variety of surficial geologic units (Meremonte, 1996; Glassmoyer, 1996). Figure 1 shows the distribution of portable recorders and table 1 gives their locations. Over 1,300 records from 61 sources and 87 sites are utilized. The main-shock records are SMA-1 strong motion accelerograms for the four events: 1971 San Fernando (M6.6), 1987 Whittier Narrows (M5.8), 1991 Sierra Madre (M5.8), and 1994 Northridge (M6.7). These records are available for 223 sites and give coverage to the south and west not represented in the Northridge aftershock data set. These records come from four primary sources: U.S. Geological Survey, California Strong Motion Instrumentation Program, University of Southern California, and Los Angeles Department of Water and Power (Porcella and others, 1994; Shakal and others, 1994, 1995; Trifunac and others, 1994). The main-shock station locations are shown on figure 2 and listed in table 2.

Figure 1
Click for larger image.

FIGURE 1. Aftershock station locations, indicated by triangles and listed in table 1, used in the determination of site response. Surficial geology is from Tinsley and Furnal (1985).


The method used to calculate site response is described in detail in two papers: Hartzell and others (1996) for the aftershock records, and Harmsen (1996) for the main-shock records. These two data sets present their own individual processing problems that are handled separately and described in the above papers. The general process, however, is relatively simple and is summarized here. The S waves from the above data records are first transformed to the frequency domain using a fast Fourier transform. Then, following Andrews (1986), the S-wave spectrum, U(f), for source i and site j is decomposed as

Si(f).Rj(f).Pij(f,r) = Uij(f),

were S(f) is the source spectrum, R(f) is the site spectrum, P(f) is the path spectrum, f is frequency, and r is the source to station distance. General assumptions are made about the nature of P(f,r), including 1/r attenuation with distance and an elastic attenuation, Q(f). Incorporating these propagation path effects and taking the logarithm of both sides gives the linear expression

log Si(fk ) + log Rj(fk ) = log U'ij(fk),

where U'(f) is the path-corrected S-wave spectrum. Linear least squares analysis can then be used to solve for the best-fitting source and site terms. Averages are then taken over frequency bands of R(f) for each site to obtain the site amplification factors. Two frequency bands are selected, 0.5-1.5 Hz and 2.0-6.0 Hz, on the basis of the bandwidth of the available ground-motion data. Main-shock records with suspected nonlinear effects have been omitted from the analysis.

The aftershock and strong-motion data results must be normalized with respect to each other. In the above analysis, a reference site is selected to remove a fundamental degree of freedom in the problem. The site response at this reference site must be specified. To accomplish the normalization, we select a common reference site for both the aftershock and main-shock record analysis. We choose one of the most competent rock sites in our database and assume its amplification is approximately 1 at all frequencies. This site is the Encino Reservoir Dam, ER1 in the aftershock data set (table 1) and ERD in the main-shock data set (table 2). These two stations are separated by approximately 20 ft.

Click for larger image.

FIGURE 2. Strong-motion station locations, indicated by triangles and listed in table 2, used in the determination of site response.


Given the point determinations of site response at each of the recording sites, a method is needed to spatially continue these values to produce a map. Surficial geology and shallow S-wave velocity are logical means of interpolating and extrapolating our point values of site response. The correlation of larger ground motion with less competent surficial material and lower S-wave velocity has long been observed and used in studies of site response (Borcherdt, 1970, 1994; Boore and others, 1993, 1994; Joyner and Fumal, 1984; and Harmsen, 1996). In the work of Hartzell and others (1996) and Harmsen (1996), however, significant variations in site response are observed for the same mapped surficial geologic unit, particularly within the young basin sediments (factor of 2 in amplitude within a few hundred meters). These results reflect the fact that shallow S-wave velocities vary significantly within these sediments (Williams and others, 1996; Harmsen, 1996) and buried structures such as folds and mini-basins focus and defocus energy at the surface. Because of these variations in S-wave velocity and the sparseness of the S-wave data set, S-wave velocity is presently not a feasible means of spatially continuing our site amplification values. For this report we have chosen a compromise by calculating average site amplification values for Mesozoic and Tertiary rocks and Pleistocene and Holocene alluvium, as well as contouring the amplification values within the more spatially variable sedimentary basins. With minor corrections to the extent of the hard-rock units, we use the surficial geology of Tinsley and Fumal (1985).

Table 3 lists the average amplification values for the four ages (Mesozoic, Tertiary, Pleistocene, and Holocene) and for the two frequency bands 0.5-1.5 Hz and 2.0-6.0 Hz. Two different averages are given: (1) the geometric average of all sites within the indicated surficial geologic unit, and (2) the value predicted by a regression relationship between shallow S-wave velocity and site amplification obtained by Harmsen (1996). For the second case, we use the arithmetic average of the S-wave velocity in each geologic unit. Because of variability within a rock type and possible topographic amplification at some sites, the average Mesozoic values are not all 1. We therefore normalize the amplifications contoured on the maps and lists in tables 1, 2, and 3 so that the average Mesozoic values are 1. (The pre-normalization values for Mesozoic units are: 1.4 (2.0-6.0Hz, average #1), 1.9 (2.0-6.0Hz, average #2), 1.0 (0.5-1.5Hz, average #1), 1.3 (0.5-1.5Hz, average #2).)

Contouring of amplification values within the sedimentary basins is done using the inverse distance-weighted interpolation algorithm in the Arc/Info GIS software package. An interpolated value is a linearly weighted combination of a set of sample points. The weight is a function of inverse distance. This weighted distance average cannot be greater than the highest or less than the lowest input. Therefore, it cannot create highs or lows if these extremes have not already been sampled. In this process all the Pleistocene and Holocene data points are contoured together, without distinction. After contouring, the Mesozoic and Tertiary units are layed down as masks over the contours. The result is discontinuous contours at the contacts of the sedimentary basins with the Mesozoic

and Tertiary outcrops, reflecting the discontinuous change in material properties and amplification factors. In addition, some of the contours are given hachures to clarify the down-gradient direction.



Graphics files of the site-response maps are available to those who have the capability to plot large-format, color encapsulated-PostScript or HPGL2 files. The files can be accessed from

For information on how to download files from the internet, click here.

The site response map files are in subdirectory \downloads\.


The site-response maps have been approved for release and publication by the Director of the U.S. Geological Survey (USGS). However, the USGS reserves the right to revise these maps pursuant to further analysis and review. Furthermore, these maps are released on condition that neither the USGS nor the United States Government may be held liable for any damages resulting from their authorized or unauthorized use.


Many organizations and individuals contributed to the collection of the data using in this report. Strong motion records came from The California Strong Motion Instrumentation Program (CSMIP), the Civil Engineering Department of the University of Southern California (USC), the Los Angeles Department of Water and Power (LADWP), and the USGS. Aftershock data records were collected by Tom Bice, Paul Bodin, David Carver, Edward Cranswick, Arthur Frankel, Mark Meremonte, Jack Odum, Dee Overturf, Rob Williams, and David Worley of the USGS, Golden, Colo.; Chris Dietel, Gene Sembera, and Leif Wennerberg of the USGS, Menlo Park, Calif.; and Sue Hough of the USGS, Pasadena, Calif.

Table 1. Aftershock station locations and amplification factors.

Table 2. Main-shock station locations and amplification factors

Table 3.--Average amplification


Andrews, D.J., 1986, Objective determination of source parameters and similarity of earthquakes of different size, in Earthquake Source Mechanics: Geophysical Monograph 37, v. 6, p. 259-267.

Boore, D.M., Joyner, W.B., and Fumal, T.E., 1993, Estimation of response spectra and peak accelerations from western North American earthquakes--An interim report: U.S. Geological Survey Open-File Report 93-509, 72 p.

Boore, D.M., Joyner, W.B., and Fumal, T.E., 1994, Estimation of response spectra and peak accelerations from western North American earthquakes--An interim report, Part 2: U.S. Geological Survey Open-File Report 94-127, 40 p.

Borcherdt, R.D., 1970, Effects of local geology on ground motion near San Francisco Bay: Bulletin of the Seismological Society of America, v. 60, p. 29-61.

Borcherdt, R.D., 1994, Estimates of site-dependent response spectra for design (methodology and justification): Earthquake Spectra, v. 10, p. 617-654.

Campbell, K.W., 1976, A note on the distribution of earthquake damage in Long Beach, 1933: Bulletin of the Seismological Society of America, v. 66, p. 1001-1006.

Glassmoyer, G., 1996, Digital recordings of aftershocks of the 17 January 1994 Northridge, California, earth quake: U.S. Geological Survey Open-File Report 96-0000 (in preparation).

Harmsen, S., 1996, Determination of site amplification in the Los Angeles urban area from inversion of strong-motion records: Bulletin of the Seismological Society of America (submitted).

Hartzell, S., Leeds, A., Frankel, A., and Michael, J., 1996, Site response for urban Los Angeles using aftershocks of the Northridge earthquake: Bulletin of the Seismological Society of America, v. 86, p. S168-S192.

Joyner, W., and Fumal, T.E., 1984, Use of measured S-wave velocity for predicting geologic site effects on strong ground motion: Proceedings of the Eight World Conference on Earthquake Engineering, San Francisco, Calif., v. 2, p. 777-783.

Meremonte, M., Frankel, A., Cranswick, E., Carver, D., and Worley, D., 1996, Urban seismology-Northridge aftershocks recorded by multi-scale arrays of portable digital seismographs: Bulletin of the Seismological Society of America, v. 86, p. 1350-1363.

Porcella, R.L., Etheredge, E.C., Maley, R.P., and Acosta, A.V., 1994, Accelerograms recorded at USGS national strong-motion network stations during the Ms=6.6 Northridge, California, eathquake of January 17, 1994: U.S. Geological Survey Open-File Report 94-141,

100 p.

Rogers, A.M., Tinsley, J.C., Hays, W.H., and King, K.W., 1979, Evaluation of the relation between near-surface geological units and ground response in the vicinity of Long Beach, California: Bulletin of the Seismological Society of America, v. 69, p.1603-1622.

Rogers, A.M., Covington, P.A., Park, R.B., Borcherdt, R.D., and Perkins, D.M., 1980, Nuclear event time histories and computed site transfer functions for locations in the Los Angeles region: U.S. Geological Survey Open-File Report 80-1173, 207 p.

Rogers, A.M., Borcherdt, R.D., Covington, P.A., and Perkins, D.M., 1984, A comparative ground response study near Los Angeles using recordings of Nevada nuclear tests and the 1971 San Fernando earthquake: Bulletin of the Seismological Society of America, v. 74, p. 1925-1949.

Rogers, A.M., Tinsley, J.C., and Borcherdt, R.D., 1985, Predicting relative ground response, in Evaluating earthquake hazards in the Los Angeles region-an earth-science perspective: U.S. Geological Survey Professional Paper 1360, p. 221-247.

Shakal, A., Huang, M., Darragh, R., Cao, T., Sherburne, R., Malhotra, P., Cramer, C., Sydnor, R., Graizer, V., Maldonado, G., Petersen, C., and Wampole, J., 1994, CSMIP strong-motion records from the Northridge, California earthquake of 17 January 1994: California Strong Motion Instrumentation Program Report 94-07, 308 p.

Shakal, A., Huang, M., Darragh, R., Brady, A., Trifunac, M., Lindvall, C., Wald, D., Heaton, T., and Mori, J., 1995, Recorded ground and structure motions, in Northridge earthquake of January 17, 1994, reconnaissance report: Earthquake Spectra, v. 11, Supplement C, p. 13-96.

Tinsley, J.C., and Fumal, T.E., 1985, Mapping Quaternary sedimentary deposits for areal variations in shaking response, in Evaluating earthquake hazards in the Los Angeles region-An earth-science perspective: U.S. Geological Survey Professional Paper 1360, p. 101-125.

Trifunac, M.D., Todorovska, M.I., and Ivanovic, S.S., 1994, A note on distribution of uncorrected peak ground accelerations during the Northridge, California, earthquake of 17 January 1994: Soil Dynamics and Earthquake Engineering, v. 13, p. 187-196.

Williams, R.A., Stephenson, W.J., Odum, J.K., and Worley, D.M., 1996, Shallow P- and S-wave velocities at eleven aftershock recording stations of the Northridge earthquake, San Fernando valley, California: U.S. Geological Survey Open-File Report 96-261, 10 p.

Return to Engineering Seismology Page logo