- Maps - General Information
- Map Types
- Global Earthquake ShakeMaps
- ShakeMap Atlas
- Scenario Earthquakes
For additional background information, download the ShakeMap Manual
A ShakeMap is a representation of ground shaking produced by an earthquake. The information it presents is different from the earthquake magnitude and epicenter that are released after an earthquake because ShakeMap focuses on the ground shaking produced by the earthquake, rather than the parameters describing the earthquake source. So, while an earthquake has one magnitude and one epicenter, it produces a range of ground shaking levels at sites throughout the region depending on distance from the earthquake, the rock and soil conditions at sites, and variations in the propagation of seismic waves from the earthquake due to complexities in the structure of the Earth's crust.
Part of the strategy for generating rapid-response ground motion maps is to determine the best format for reliable presentation of the maps given the diverse audience, which includes scientists, businesses, emergency response agencies, media, and the general public. In an effort to simplify and maximize the flow of information to the public, we have developed a means of generating not only peak ground acceleration and velocity maps, but also an instrumentally-derived, estimated Modified Mercalli Intensity map. This map makes it easier to relate the recorded ground motions to the expected felt and damage distribution. The Instrumental Intensity map is based on a combined regression of recorded peak acceleration and velocity amplitudes. (see Intensity Maps)
Even with the current seismic station distribution in Calfornia, data gaps are common, particularly for events outside the densely-instrumented metropolitan regions surrounding Los Angeles and San Francisco. In order to stabilize contouring and minimize the misrepresentation of the ground motion pattern due to data gaps, we augment the data with predicted values in areas without data. Given the epicenter and magnitude (and for larger earthquakes, fault geometry if avaliable), peak motion amplitudes in spare regions are estimated from the ground motion prediction equations.
The instrumental intensity map shows station symbols as open triangles in order to see the underlying intensity value. The legend bar at the bottom explains the colors (and see Intensity Maps below). For the intensity map as with other maps, station locations are the best indicator of where the map is most accurate: Near seismic stations the shaking is well constrained by data; far from such stations, the shaking is estimated using standard seismological inferences and interpolation.
Note: ShakeMaps are generated automatically following moderate and large earthquakes. These are preliminary ground shaking maps, normally posted within several minutes of the earthquake origin time. The acceleration and velocity values are raw and are at least initially, NOT checked by humans. Further, since ground motions and intensities typically can vary significantly over small distances, these maps are only APPROXIMATE. At small scales, they should be considered unreliable. Finally, the input data is raw and unchecked, and may contain errors. (See Disclaimer)
Maps - General Information
- The Red star usually located near the center of the map is the epicenter.
- Small unfilled circles, shown on some maps, are points where strong motion values were estimated and used to fill gaps in the station distribution.
- The colored triangles indicate reporting stations. In California,
- Red triangles are stations from the Caltech/USGS digital telemetered network.
- Blue triangles represent California Geological Survey (CGS) and California Strong Motion Instrumentation Program (CSMIP) dial-up stations.
- Yellow triangles represent stations from the ANZA Regional Network.
- Green triangles represent National Strong Motion Project (NSMP) dial-up stations or, on historical maps, non-digital stations from which strong motion records were digitized.
- On the Intensity maps, station symbols are open, allowing the underlying intensity color to show thru. On the other maps station symbols (circles, triangles) are color coded according to the data type mentioned above.
The station information includes the station code and name, the agency that manages the station, the station location coordinates in degrees latitude and longitude, and the peak acceleration, velocity, and spectral acceleration for each component of ground motion (when available). Spectral acceleration maps are only made for larger earthquakes (normally, magntidue greater than 5.5). When the peak ground motion maps are made, the value from the peak horizontal component of ground motion is used as the value for the station. This value is highlighted in bold in the station information.
Components from many stations are defined by three letter codes. The last letter indicates the orientation (Z = vertical, N = horizontal north, E = horizontal east). The first two letters indicate the instrument class:
|VL||low gain channels on the analog network|
|VH||high gain channels on the analog network|
|AS||FBA's on the analog network|
|HL||FBA's on the digital network|
|HN||FBA's on the digital network|
|BH||broadband data streams|
|HH||broadband data streams|
FBA's (force balance accelerometers) are designed to record extremely large ground motions and can accurately record waves from very large earthquakes. However, ground motions from small and moderate earthquakes are often too small to trigger these instruments or rise above instrument noise. On the other hand, Broadband seismic sensors can record extremely small ground motions and accurately record waves from earthquakes that range from very small up to moderately large. A number of stations have both FBA and broadband sensors. For ShakeMap, the network tends to emphasize FBA recordings for large ground motions and broadband recordings for small ground motions.
Occassionally, station channels will be flagged due to problems with the station or possibly anomalous peak values. In this case, the popup window of station information will indicate the flagging with the following codes:
|G||Glitch (clipped or below noise)|
|I||Incomplete time series|
|N||Not in list of known stations|
Peak Acceleration Maps
Peak horizontal acceleration at each station is contoured in units of percent-g (where g = acceleration due to the force of gravity = 981 cm/s/s). The peak values of the vertical components are not used in the construction of the maps because they are, on average, lower than the horizontal amplitudes and ground motion prediciton equations used to fill in data gaps between stations are based on peak horizontal amplitudes. The contour interval varies greatly and is based on the maximum recorded value over the network for each event.
For moderate to large events, the pattern of peak ground acceleration is typically quite complicated, with extreme variability over distances of a few km. This is attributed to the small scale geological differences near the sites that can significantly change the high-frequency acceleration amplitude and waveform character. Although distance to the causative fault clearly dominates the pattern, there are often exceptions, due to local focussing and amplification. This makes interpolation of ground motions at one site to a nearby neighbor somewhat risky. Peak acceleration pattern usually reflects what is felt from low levels of shaking up to to moderate levels of damage.
Peak Velocity Maps
Peak velocity values are contoured for the maximum horizontal velocity (in cm/sec) at each station. As with the acceleration maps, the vertical component amplitudes are disregarded for consistency with the regression relationships used to estimate values in gaps in the station distribution. Typically, for moderate to large events, the pattern of peak ground velocity reflects the pattern of the earthquake faulting geometry, with largest amplitudes in the near-source region, and in the direction of rupture (directivity). Differences between rock and soil sites are apparent, but the overall pattern is normally simpler than the peak acceleration pattern. Severe damage, and damage to flexible structures is best related to ground velocity. For reference, the largest recorded ground velocity from the 1994 Northridge (Magnitude 6.7) earthquake made at the Rinaldi Receiving station, reached 183 cm/sec.
Spectral Response Maps
Following earthquakes larger than magnitude 5.5, spectral response maps are made. Response spectra portray the response of a damped, single-degree-of-freedom oscillator to the recorded ground motions. This data representation is useful for engineers determining how a structure will react to ground motions. The response is calculated for a range of periods. Within that range, the Uniform Building Code (UBC) refers to particular reference periods that help define the shape of the "design spectra" that reflects the building code.
ShakeMap spectral response maps are made for the response at three UBC reference periods: 0.3, 1.0, and 3.0 seconds. For each station, the value used is the peak horizontal value of 5% critically damped pseudo-acceleration.
Rapid Instrumental Intensity Maps
As an effort to simplify and maximize the flow of information to the public, we have developed a means of generating estimated Modified Mercalli Intensity maps based on instrumental ground motion recordings (Wald et al., 1999). These "Instrumental Intensities" are based on a combined regression of peak acceleration and velocity amplitudes vs. observed intensity for eight significant California earthquakes (1971 San Fernando, 1979 Imperial Valley, 1986 North Palm Springs, 1987 Whittier, 1989 Loma Preita, 1991 Sierra Madre, 1992 Landers, and 1994 Northridge).
From the comparison with observed intensity maps, we find that a regression based on peak velocity for intensity > VII and on peak acceleration for intensity < VII is most suitable. This is consistent with the notion that low intensities are determined by felt accounts (sensitive to acceleration). Moderate damage, at intensity VI-VII, typically occurs in rigid structures (masonry walls, chimneys, etc.) which also are sensitive to high-frequency (acceleration) ground motions. As damage levels increase, damage also occurs in flexible structures, for which damage is proportional to the ground velocity, not acceleration. By relating recorded ground motions to Modified Mercalli intensities, we can now estimate shaking intensities within a few minutes of the event based on the recorded peak motions made at seismic stations.
A descriptive table of Modified Mercalli Intensity is available in the education section of the website. A table of intensity descriptions with the corresponding peak ground acceleration (PGA) and peak ground velocity (PGV) values used in the ShakeMaps is given below. ShakeMap uses PGA to estimate intensities lower than V, it linearly combines PGA & PGV for intensities greater than V and less than VII, and it uses PGV for intensities greater than VII (See Wald et al., 1999b, for more details).
ShakeMaps are computed as the uncertainty-weighted combination of ground motion amplitudes from a Ground Motion Prediction Equation (GMPE), seismic data, and (optionally) reports of macro seismic intensity. This weighted-averaging process allows us to compute an uncertainty at each grid point in a ShakeMap. Since the GMPE also provides an estimate of ground motion uncertainty at each point, we can compute the ratio of the final ShakeMap uncertainty to the GMPE uncertainty. This ratio lets us know at each grid point if the ShakeMap is more or less uncertain than a purely predictive map generated by the GMPE.
We utilize the uncertainty ratio to produce a graded map of uncertainty. Where the ratio is 1.0 (meaning the ShakeMap is purely predictive), the map is colored white. Where the ratio is greater than 1.0 (meaning that the ShakeMap uncertainty is high because of unknown fault geometry) the map shades toward dark red, and where the uncertainty is less than 1.0 (because the presence of data decreases the uncertainty) the map shades toward dark blue. These maps provide a quick visual summary of quality of the ground motion estimates over the area of interest. ShakeMaps are also given a letter grade, based on the mean uncertainty ratio within the area of the MMI 6 contour (on the theory that this is the area most important to accurately represent). A ratio of 1.0 is given a grade of "C". Maps with mean ratios greater than 1.0 get grades of "D" or "F". Ratios less than 1.0 earn grades of "B" or "A". If the map does not contain areas of MMI >= 6, no grade is assigned.
Global Earthquake ShakeMaps
For regions around the world were there are insufficient near-real time strong motion seismic stations to generate an adequate, strong-ground-motion data-controlled ShakeMap, we can still provide a very useful estimate of the shaking distribution using the ShakeMap software.
Initially, a point source approximation (hypocenter and magnitude) is used to constrain region-specific empirical ground motion estimation. Site amplification is approximated from a relationship developed between topographic gradient and shear-wave velocity. Additional constraints for these predictive maps come primarily from three important sources, the availability of which varies depending on the region in which the earthquake occurred, as well as a function of time after the earthquake occurrence. These constraints include: (1) additional earthquake source information, particularly fault rupture dimensions, (2) observed macroseismic intensities (including via the USGS "Did You Feel It?" system, and (3) observed strong ground motions, where and when available.
When intensity data are used directly, the peak ground motion parameters are inferred from the macroseismic observations using the equations of Wald et al, (1999a). This is the opposite approach normally used in ShakeMaps where numerous seismic recording are available, and from them the intensities are then inferred.
Input data for global ShakeMaps are depicted with circles for intensity reports and triangles for strong motion data. Intensities are further separated by color: Blue circles indicate data collected via the USGS "Did You Feel It?" Community Internet Intensity Map system and Yellow circles indicate traditional Modifed Mercalli Intensity assignments. For the intensity maps, circles are open, allowing the underlying intensity color to show thru. On the other maps observation location symbols (circles, triangles) are color coded according to the data type mentioned above.
An atlas of maps of peak ground motions and intensity “ShakeMaps” has been developed for approximately 1,000 recent and historical global earthquakes (Allen and others, 2008). These maps are produced using established ShakeMap methodology (Wald and others, 1999c, 2005) and constraints from macroseismic intensity data, instrumental ground motions, regional topographically-based site amplifications (Wald and Allen, 2007), and published earthquake rupture models. The Atlas of ShakeMaps provides a consistent and quantitative description of the distribution of shaking intensity for recent global earthquakes (January 1973 – September 2007 for Version 1.0). We anticipate that the Atlas will be regularly updated with more data constraints for historical events and the addition of future significant events as time progresses.
The Atlas was developed specifically for calibrating global earthquake loss estimation methodologies to be used in the USGS Prompt Assessment of Global Earthquakes for Response (PAGER) Project. PAGER will employ these loss models to rapidly estimate the impact of global earthquakes as part of the USGS National Earthquake Information Center’s earthquake response protocol. Though developed primarily for PAGER, we anticipate many other uses for the historical ShakeMap Atlas, including disaster response planning, capacity building and outreach programs, in addition to calibration of other global loss methodology approaches.
The primary sources for instrumental data are:
- Pacific Earthquake Engineering Research Center’s Next Generation Attenuation (NGA) database
- Consortium of Organizations for Strong-Motion Observation Systems (COSMOS)
- European Strong-Motion Database (ESD)
- Kyoshin Network (K-NET), Japan
Macroseismic intensity observations are gathered from several online sources, including:
- National Geophysical Data Center, National Oceanic and Atmospheric Administration (NOAA)
- Centro Regional de Sismología para América del Sur (CERESIS)
- Istituto Nazionale di Geofisica e Vulcanologia (INGV), Italy
Community Internet Intensity Maps data obtained from the USGS Did You Feel It? (DYFI?) system (Wald and others, 1999a; Atkinson and Wald, 2007) are a valuable ground-shaking constraint and are used for U.S. events since 1999 and for global events since 2003. Additional intensity data were gathered or digitized from numerous earthquake reconnaissance reports and peer reviewed publications, or through the generous contribution from colleagues around the world.
An important, consolidated source for finite fault models is provided by Martin Mai of the Swiss Seismological Service, Zurich. Other fault rupture dimensions were digitized from observations of surface displacement, finite-fault source inversions and slip distribution determined from teleseismic observations, or more recently, from InSAR observations. Less well-constrained faults have been estimated from earthquake aftershock distributions.
The ShakeMap Atlas aims to provide the best estimate of the shaking distribution for historical earthquakes. One important feature of the Atlas is that data constraints are not limited to those data produced by the National Earthquake Information Center, and we will openly accept data contributions that help to improve the representation of the shaking distribution for any of the events. References to all data sources in the Atlas of ShakeMaps (as of August 2008) are provided in Allen and others (2008). Please contact Trevor Allen or David Wald with any comments or data contributions.
Below we indicate the default values and predictive equations used in constructing the Atlas of ShakeMaps. Variations in ground-motion prediction equations (GMPEs) and other configurations may change on an event-by-event basis. Please see the event info.xml on individual ShakeMap download pages or Allen and others (2008) for more information.
Date Range: January 1973 – September 2007
Active crust GMPE: Boore and others (1997)
Subduction zone GMPE (intraslab and interface): Youngs and others (1997)
Stable continent GMPE: Atkinson and Boore (2006)
Seismic site-conditions: Wald and Allen (2007) <http://earthquake.usgs.gov/vs30/>
Peak ground motion to macroseismic intensity conversions: Wald and others (1999b)
- Allen, T.I., Wald, D.J., Hotovec, A.J., Lin, K., Earle, P.S., and Marano, K.D., 2008, An Atlas of ShakeMaps for selected global earthquakes: U.S. Geological Survey Open-File Report, 2008-1236, 34 p.
- Atkinson, G.M., and Boore, D.M., 2006, Earthquake ground-motion predictions for eastern North America: Bull. Seism. Soc. Am., v. 96, p. 2181-2205.
- Atkinson, G.M., and Wald, D.J., 2007, "Did You Feel It?" intensity data: A surprisingly good measure of earthquake ground motion: Seism. Res. Lett., v. 78, no. 3, p. 362-368.
- Boore, D.M., Joyner, W.B., and Fumal, T.E., 1997, Equations for estimating horizontal response spectra and peak acceleration from Western North American earthquakes: A summary of recent work: Seismological Research Letters, v. 68, no. 1, p. 128-153.
- Wald, D.J., and Allen, T.I., 2007, Topographic slope as a proxy for seismic site conditions and amplification: Bull. Seism. Soc. Am., v. 97, no. 5, p. 1379-1395.
- Wald, D.J., Quitoriano, V., Dengler, L., and Dewey, J.W., 1999a, Utilization of the Internet for rapid Community Intensity Maps: Seism. Res. Lett., v. 70, p. 680-697.
- Wald, D.J., Quitoriano, V., Heaton, T.H., and Kanamori, H., 1999b, Relationship between Peak Ground Acceleration, Peak Ground Velocity, and Modified Mercalli Intensity in California: Earthquake Spectra, v. 15, no. 3, p. 557-564.
- Wald, D.J., Quitoriano, V., Heaton, T.H., Kanamori, H., Scrivner, C.W., and Worden, B.C., 1999c, TriNet "ShakeMaps": Rapid generation of peak ground-motion and intensity maps for earthquakes in southern California: Earthquake Spectra, v. 15, no. 3, p. 537-556.
- Wald, D.J., Worden, B.C., Quitoriano, V., and Pankow, K.L., 2005, ShakeMap manual: technical manual, user's guide, and software guide: U.S. Geological Survey, 132 p.
- Youngs, R.R., Chiou, S.-J., Silva, W.J., and Humphrey, J.R., 1997, Strong ground motion attenuation relationships for subduction zone earthquakes: Seism. Res. Lett., v. 68, no. 1, p. 58-73.
Earthquake Scenarios describe the expected ground motions and effects of specific hypothetical large earthquakes. In planning and coordinating emergency response, utilities, emergency responders, and other agencies are best served by conducting training exercises based on realistic earthquake situations, ones that they are most likely to face. Scenario earthquakes can fill this role; they can be generated for any potential hypothetical future or past historic earthquake by the following steps.
First, assume a particular fault or fault segment will rupture over a certain length relying on consensus-based information about the potential behavior of the fault. For historic events, the actual rupture dimensions may be constrained based on existing observations or models. Second, estimate ground motions at all locations in a chosen region surrounding the causative fault.
These earthquake scenarios are not earthquake predictions. That is, no one knows in advance when or how large a future earthquake will be. However, if we make assumptions about the size and location of a hypothetical future earthquake, we can make a reasonable prediction of the effects of the assumed earthquake, particularly the way in which the ground will shake. This knowledge of the potential shaking effects is the main benefit of the earthquake scenario for planning and preparedness purposes.
Choosing An Appropriate Earthquake Scenario
The scenario earthquakes compiled on the Northern and Southern California ShakeMap web pages represent 236 different earthquakes anticipated for California: this set of scenario earthquakes was determined by the Unified California Earthquake Rupture Forecast (UCERF3). In 2008, UCERF2 concluded that the likelihood of one or more large (M≥6.7) earthquakes in over the next 30 years is 99%. The likelihoods of one or more large (M≥6.7) earthquakes over the next 30 years in Northern or Southern California are 93% and 97% respectively. This regional probabilities consider both large earthquakes occurring on the major (mapped) fault systems as well as large earthquakes occurring in the background, that is, on unmapped faults or on faults considered unlikely to rupture. As part of the UCERF process, the major California fault systems were subdivided into individual fault segments. These fault segments are the shortest fault sections deemed capable of rupturing in large earthquakes. Each fault was divided into segments by evaluating paleoseismic data as well as structural changes along the fault systems such as bends, offsets, and different rock types. UCERF allowed for the possibility of either single segment or multi-segment earthquakes on these fault systems. Each possible combination of segments was called a rupture source and each rupture source represents a possible future scenario earthquake. UCERF3 estimated the likelihood of occurrence of all combinations of the possible rupture sources for each fault system. The resulting rupture sources, 236 in all, along with their mean magnitude, probability of occurring in the next 30 years, and the associated uncertainties in the probabilities are listed in Table 1, divided into Northern and Southern California faults. The scenario ShakeMaps graphically illustrate the strength and regional extent of shaking that can be expected from future earthquakes in Northern and Southern California. It is important to note that the predicted shaking is a median estimate: when a large earthquake actually occurs, the ground shaking will exceed these estimates in some places and will be lower in others. Users interested in specific scenarios for planning purposes are encouraged to make such a request by filling out a Comment form
Scenario Earthquakes for California as determined by UCERF3.
|Northern California Faults||Lat||Long||Magnitude|
|Bartlett Springs M7.3 Scenario||39.31||-122.83||7.3|
|Battle Creek M6.7 Scenario||40.40||-122.11||6.7|
|Big Lagoon-Bald Mtn M7.5 Scenario||41.11||-123.79||7.5|
|CalaverasCC M6.4 Scenario||37.44||-121.80||6.4|
|CalaverasCC+CS M6.5 Scenario||37.43||-121.79||6.5|
|CalaverasCN M6.9 Scenario||37.78||-121.98||6.9|
|CalaverasCN+CC M7.0 Scenario||37.76||-121.97||7.0|
|CalaverasCN+CC+CS M7.0 Scenario||37.74||-121.95||7.0|
|CalaverasCS M5.8 Scenario||36.86||-121.41||5.8|
|Cedar Mtn-Mahogany Mtn M7.1 Scenario||41.71||-121.88||7.1|
|Collayomi M6.7 Scenario||38.82||-122.72||6.7|
|Fickle Hill M7.1 Scenario||40.85||-123.83||7.1|
|Fish Slough M6.8 Scenario||37.37||-118.29||6.8|
|Gillem-Big Crack M6.8 Scenario||41.72||-121.51||6.8|
|Great Valley 1 M6.8 Scenario||39.60||-122.39||6.8|
|Great Valley 2 M6.5 Scenario||39.25||-122.37||6.5|
|Great Valley 3 Mysterious Ridge M7.1 Scenario||38.99||-122.35||7.1|
|Great Valley 4a Trout Creek M6.6 Scenario||38.61||-122.19||6.6|
|Great Valley 4b Gordon Valley M6.8 Scenario||38.46||-122.16||6.8|
|Great Valley 5 Pittsburg Kirby Hills M6.7 Scenario||38.23||-121.95||6.7|
|Great Valley 7 M6.9 Scenario||37.63||-121.51||6.9|
|Great Valley 8 M6.8 Scenario||37.34||-121.21||6.8|
|Great Valley 9 M6.8 Scenario||37.01||-121.01||6.8|
|Great Valley 10 M6.5 Scenario||36.73||-120.81||6.5|
|Great Valley 11 M6.6 Scenario||36.56||-120.67||6.6|
|Great Valley 12 M6.4 Scenario||36.42||-120.50||6.4|
|Great Valley 13 Coalinga M7.1 Scenario||36.20||-120.47||7.1|
|Great Valley 14 Kettleman Hills M7.2 Scenario||35.90||-120.28||7.2|
|Concord-Green Valley M6.8 Scenario||38.31||-122.16||6.8|
|Greenville M7.0 Scenario||37.51||-121.55||7.0|
|Hartley Springs M6.8 Scenario||37.70||-118.90||6.8|
|Hat Creek-McArthur-Mayfield M7.2 Scenario||41.24||-121.56||7.2|
|Hayward-Rodgers CreekHN M6.6 Scenario||37.88||-122.26||6.6|
|Hayward-Rodgers CreekHN+HS M7.0 Scenario||37.88||-122.26||7.0|
|Hayward-Rodgers CreekHS M6.8 Scenario||37.59||-121.99||6.8|
|Hayward-Rodgers CreekRC M7.1 Scenario||38.19||-122.51||7.1|
|Hayward-Rodgers CreekRC+HN M7.2 Scenario||38.43||-122.68||7.2|
|Hayward-Rodgers CreekRC+HN+HS M7.3 Scenario||38.53||-122.75||7.3|
|Hilton Creek M6.9 Scenario||37.52||-118.66||6.9|
|Honey Lake M7.0 Scenario||40.01||-120.06||7.0|
|Hosgri M7.3 Scenario||34.89||-120.75||7.3|
|Hunting Creek-Berryessa M7.1 Scenario||38.55||-122.26||7.1|
|Likely M7.0 Scenario||40.92||-120.29||7.0|
|Little Salmon On-Offshore M7.5 Scenario||40.66||-123.87||7.5|
|Little Salmon Offshore M7.3 Scenario||40.91||-124.17||7.3|
|Little Salmon Onshore M7.1 Scenario||40.66||-123.89||7.1|
|Los Osos M7.0 Scenario||35.24||-120.87||7.0|
|Maacama-Garberville M7.4 Scenario||38.75||-122.88||7.4|
|Mad River M7.2 Scenario||40.85||-123.79||7.2|
|McKinleyville M7.2 Scenario||40.86||-123.75||7.2|
|Mono Lake M6.8 Scenario||38.03||-119.06||6.8|
|Monte Vista-Shannon M6.5 Scenario||37.31||-122.12||6.5|
|Monterey Bay-Tularcitos M7.3 Scenario||36.37||-121.54||7.3|
|N. San AndreasSAN M7.5 Scenario||37.82||-122.60||7.5|
|N. San AndreasSAN+SAP+SAS M7.8 Scenario||39.16||-123.83||7.8|
|N. San AndreasSAO M7.4 Scenario||39.22||-123.86||7.4|
|N. San AndreasSAO+SAN M7.7 Scenario||40.13||-124.19||7.7|
|N. San AndreasSAO+SAN+SAP M7.9 Scenario||40.16||-124.24||7.9|
|N. San AndreasSAO+SAN+SAP+SAS M7.9 Scenario||40.13||-124.19||7.9|
|N. San AndreasSAP M7.2 Scenario||37.27||-122.11||7.2|
|N. San AndreasSAP+SAS M7.5 Scenario||37.79||-122.57||7.5||N. San AndreasSAS M7.1 Scenario||36.82||-121.50||7.1|
|North Tahoe M6.7 Scenario||39.09||-119.95||6.7|
|Ortigalita M7.1 Scenario||37.25||-121.27||7.1|
|Point Reyes M6.9 Scenario||38.01||-122.89||6.9|
|Quien Sabe M6.6 Scenario||36.79||-121.24||6.6|
|Rinconada M7.5 Scenario||36.63||-121.70||7.5|
|Robinson Creek M6.7 Scenario||38.19||-119.22||6.7|
|Round Valley M7.1 Scenario||37.27||-118.52||7.1|
|S. San AndreasCH M7.1 Scenario||35.37||-119.93||7.1|
|S. San AndreasCH+CC M7.4 Scenario||35.66||-120.22||7.4|
|S. San AndreasPK M6.1 Scenario||35.80||-120.35||6.1|
|S. San AndreasPK+CH M7.1 Scenario||35.80||-120.35||7.1|
|S. San AndreasPK+CH+CC M7.4 Scenario||35.85||-120.40||7.4|
|S. San AndreasPK+CH+CC+BB M7.6 Scenario||35.90||-120.46||7.6|
|S. San AndreasPK+CH+CC+BB+NM M7.7 Scenario||35.75||-120.30||7.7|
|San Gregorio M7.5 Scenario||36.36||-121.90||7.5|
|San Luis Range So Margin M7.2 Scenario||34.91||-120.20||7.2|
|Surprise Valley M7.2 Scenario||41.17||-119.94||7.2|
|Table Bluff M7.2 Scenario||40.76||-124.19||7.2|
|Trinidad M7.5 Scenario||40.90||-123.74||7.5|
|West Napa M6.7 Scenario||38.19||-122.26||6.7|
|West Tahoe M7.1 Scenario||38.87||-119.94||7.1|
|White Mountains M7.4 Scenario||37.04||-118.17||7.4|
|Zayante-Vergeles M7.0 Scenario||36.80||-121.55||7.0|
|Southern California Faults||Lat||Long||Magnitude|
|Anacapa-Dume Alt 1 M7.2 Scenario||34.09||-118.72||7.2|
|Anacapa-Dume Alt 2 M7.2 Scenario||34.09||-118.61||7.2|
|Birch Creek M6.6 Scenario||37.03||-118.25||6.6|
|Blackwater M7.1 Scenario||35.09||-117.11||7.1|
|Burnt Mtn M6.8 Scenario||34.11||-116.47||6.8|
|Calico-Hidalgo M7.4 Scenario||34.22||-116.18||7.4|
|Casmalia Orcutt Frontal M6.7 Scenario||34.90||-120.58||6.7|
|Channel Islands Thrust M7.3 Scenario||34.15||-119.34||7.3|
|Chino Alt 1 M6.7 Scenario||33.96||-117.78||6.7|
|Chino Alt 2 M6.8 Scenario||33.97||-117.77||6.8|
|Clamshell-Sawpit M6.7 Scenario||34.30||-117.91||6.7|
|Cleghorn M6.8 Scenario||34.29||-117.39||6.8|
|Coronado Bank M7.4 Scenario||33.28||-117.92||7.4|
|Cucamonga M6.7 Scenario||34.23||-117.46||6.7|
|Death Valley Black Mtns Frontal M7.3 Scenario||36.48||-116.93||7.3|
|Death Valley Connected M7.8 Scenario||37.77||-118.19||7.8|
|Death Valley No M7.3 Scenario||36.57||-116.90||7.3|
|Death Valley No Of Cucamongo M7.2 Scenario||37.31||-117.67||7.2|
|Death Valley So M6.9 Scenario||35.63||-116.44||6.9|
|Deep Springs M6.8 Scenario||37.44||-118.04||6.8|
|Earthquake Valley M6.8 Scenario||33.17||-116.57||6.8|
|Elmore Ranch M6.7 Scenario||33.19||-115.69||6.7|
|ElsinoreCM M6.9 Scenario||32.83||-116.09||6.9|
|ElsinoreGI M6.9 Scenario||33.82||-117.57||6.9|
|ElsinoreGI+T M7.3 Scenario||33.82||-117.56||7.3|
|ElsinoreGI+T+J+CM M7.7 Scenario||32.82||-116.09||7.78|
|ElsinoreT M7.1 Scenario||33.64||-117.33||7.1|
|ElsinoreT+J+CM M7.6 Scenario||32.79||-116.00||7.6|
|ElsinoreW Combined Wmerge2.f M7.0 Scenario||33.86||-117.59||7.0|
|Elysian Park Upper M6.7 Scenario||34.14||-118.09||6.7|
|Eureka Peak M6.7 Scenario||33.99||-116.34||6.7|
|GarlockGC M7.3 Scenario||35.42||-117.75||7.3|
|GarlockGC+GW M7.6 Scenario||34.82||-118.91||7.6|
|GarlockGE M6.9 Scenario||35.60||-116.83||6.9|
|GarlockGE+GC M7.5 Scenario||35.59||-116.45||7.5|
|GarlockGE+GC+GW M7.7 Scenario||35.59||-116.47||7.7|
|GarlockGW M7.3 Scenario||34.82||-118.92||7.3|
|Gravel Hills-Harper Lk M7.1 Scenario||34.88||-116.95||7.1|
|Helendale-So Lockhart M7.4 Scenario||34.36||-116.84||7.4|
|Hollywood M6.7 Scenario||34.16||-118.30||6.7|
|Holser Alt 1 M6.8 Scenario||34.36||-118.75||6.8|
|Hunter Mountain Connected M7.6 Scenario||35.63||-116.93||7.6|
|Hunter Mountain-Saline Valley M7.2 Scenario||36.49||-117.48||7.2|
|Imperial M7.0 Scenario||32.74||-115.40||7.0|
|Independence M7.2 Scenario||36.55||-118.01||7.2|
|Johnson Valley No M6.9 Scenario||34.33||-116.47||6.9|
|Laguna Salada M7.3 Scenario||32.19||-115.18||7.3|
|Landers M7.4 Scenario||34.17||-116.42||7.4|
|Lenwood-Lockhart-Old Woman Springs M7.5 Scenario||34.42||-116.67||7.5|
|Lions Head M6.8 Scenario||34.73||-120.32||6.8|
|Little Lake M6.9 Scenario||35.69||-117.71||6.9|
|Los Alamos-West Baseline M6.9 Scenario||34.64||-120.33||6.9|
|Malibu Coast Alt 1 M6.7 Scenario||34.05||-118.62||6.7|
|Malibu Coast Alt 2 M7.0 Scenario||34.07||-118.62||7.0|
|Mission Ridge-Arroyo Parida-Santa Ana M6.9 Scenario||34.41||-119.87||6.9|
|Newport Inglewood Connected Alt 1 M7.5 Scenario||32.59||-117.15||7.5|
|Newport Inglewood Connected Alt 2 M7.5 Scenario||32.59||-117.15||7.5|
|Newport-Inglewood Alt 1 M7.2 Scenario||33.65||-117.97||7.2|
|Newport-Inglewood Alt 2 M7.2 Scenario||33.66||-117.98||7.2|
|Newport-Inglewood Offshore M7.0 Scenario||33.51||-117.80||7.0|
|North Channel M6.8 Scenario||34.36||-119.45||6.8|
|North Frontal East M7.0 Scenario||34.20||-116.78||7.0|
|North Frontal West M7.2 Scenario||34.23||-117.23||7.2|
|Northridge M6.9 Scenario||34.33||-118.74||6.9|
|Oak Ridge Connected M7.4 Scenario||34.18||-119.60||7.4|
|Oak Ridge Offshore M7.0 Scenario||34.17||-119.60||7.0|
|Oak Ridge Onshore M7.2 Scenario||34.20||-119.15||7.2|
|Owens Valley M7.3 Scenario||36.52||-118.03||7.3|
|Owl Lake M6.7 Scenario||35.63||-116.83||6.7|
|Palos Verdes Connected M7.7 Scenario||33.93||-118.52||7.7|
|Palos Verdes M7.3 Scenario||33.30||-117.93||7.3|
|Panamint Valley M7.4 Scenario||35.61||-116.91||7.4|
|Pinto Mtn M7.3 Scenario||34.06||-116.70||7.3|
|Pisgah-Bullion Mtn-Mesquite Lk M7.3 Scenario||34.14||-116.01||7.3|
|Pitas Point Connected D2.1 M7.3 Scenario||34.39||-119.16||7.3|
|Pitas Point Lower West M7.3 Scenario||34.49||-119.82||7.3|
|Pitas Point Lower-Montalvo M7.3 Scenario||34.47||-119.56||7.3|
|Pitas Point Upper M6.9 Scenario||34.36||-119.59||6.9|
|Pleito M7.1 Scenario||34.86||-119.25||7.1|
|Puente Hills Coyote Hills M6.9 Scenario||34.06||-117.91||6.9|
|Puente Hills LA M7.0 Scenario||34.13||-118.08||7.0|
|Puente Hills M7.1 Scenario||34.05||-117.92||7.1|
|Puente Hills Santa Fe Springs M6.7 Scenario||34.08||-118.09||6.7|
|Raymond M6.8 Scenario||34.18||-118.01||6.8|
|Red Mountain M7.4 Scenario||34.41||-119.36||7.4|
|Rose Canyon M6.9 Scenario||32.59||-117.15||6.9|
|S. San AndreasBB M7.1 Scenario||34.92||-119.37||7.1|
|S. San AndreasBB+NM+SM M7.6 Scenario||34.92||-119.36||7.6|
|S. San AndreasBB+NM+SM+NSB+SSB M7.8 Scenario||34.92||-119.37||7.8|
|S. San AndreasBB+NM+SM+NSB+SSB+BG+CO M7.9 Scenario||33.71||-116.14||7.9|
|S. San AndreasBG M7.1 Scenario||33.88||-116.31||7.1|
|S. San AndreasBG+CO M7.4 Scenario||33.72||-116.12||7.4|
|S. San AndreasCC M7.2 Scenario||35.19||-119.74||7.2|
|S. San AndreasCC+BB+NM+SM M7.7 Scenario||35.19||-119.74||7.7|
|S. San AndreasCC+BB+NM+SM+NSB M7.8 Scenario||35.22||-119.77||7.8|
|S. San AndreasCC+BB+NM+SM+NSB+SSB M7.9 Scenario||33.99||-116.85||7.9|
|S. San AndreasCC+BB+NM+SM+NSB+SSB+BG M7.9 Scenario||33.82||-116.85||7.9|
|S. San AndreasCC+BB+NM+SM+NSB+SSB+BG+CO M8.0 Scenario||33.62||-116.03||8.0|
|S. San AndreasCH+CC+BB+NM+SM M7.8 Scenario||35.66||-120.22||7.8|
|S. San AndreasCH+CC+BB+NM+SM+NSB+SSB M7.9 Scenario||33.97||-116.83||7.9|
|S. San AndreasCO M7.0 Scenario||33.70||-116.14||7.0|
|S. San AndreasNM M6.9 Scenario||34.81||-118.89||6.9|
|S. San AndreasNM+SM M7.5 Scenario||34.33||-117.59||7.5|
|S. San AndreasNM+SM+NSB M7.6 Scenario||34.15||-117.22||7.6|
|S. San AndreasNM+SM+NSB+SSB M7.7 Scenario||34.81||-118.89||7.7|
|S. San AndreasNM+SM+NSB+SSB+BG M7.8 Scenario||33.85||-116.35||7.8|
|S. San AndreasNM+SM+NSB+SSB+BG+CO M7.8 Scenario||33.36||-115.71||7.8|
|S. San AndreasNSB M6.9 Scenario||34.28||-117.47||6.9|
|S. San AndreasNSB+SSB M7.2 Scenario||34.28||-117.47||7.2|
|S. San AndreasNSB+SSB+BG M7.5 Scenario||33.85||-116.31||7.5|
|S. San AndreasNSB+SSB+BG+CO M7.6 Scenario||33.71||-116.13||7.6|
|S. San AndreasPK+CH+CC+BB+NM+SM M7.8 Scenario||35.85||-120.40||7.8|
|S. San AndreasPK+CH+CC+BB+NM+SM+NSB M7.9 Scenario||35.90||-120.46||7.9|
|S. San AndreasPK+CH+CC+BB+NM+SM+NSB+SSB M7.9 Scenario||33.98||-116.84||7.9|
|S. San AndreasPK+CH+CC+BB+NM+SM+NSB+SSB+BG M8.0 Scenario||33.81||-116.27||8.0|
|S. San AndreasPK+CH+CC+BB+NM+SM+NSB+SSB+BG+CO M8.0 Scenario||33.53||-115.92||8.0|
|S. San AndreasSM M7.3 Scenario||34.61||-118.27||7.3|
|S. San AndreasSM+NSB M7.4 Scenario||34.64||-118.35||7.4|
|S. San AndreasSM+NSB+SSB M7.6 Scenario||34.64||-118.35||7.6|
|S. San AndreasSM+NSB+SSB+BG M7.7 Scenario||33.81||-116.27||7.7|
|S. San AndreasSM+NSB+SSB+BG+CO M7.8 Scenario||33.54||-115.92||7.8|
|S. San AndreasSSB M6.9 Scenario||34.15||-117.22||6.9|
|S. San AndreasSSB+BG M7.3 Scenario||33.87||-116.34||7.3|
|S. San AndreasSSB+BG+CO M7.5 Scenario||33.37||-115.70||7.5|
|San Cayetano M7.2 Scenario||34.55||-118.78||7.2|
|San Gabriel M7.3 Scenario||34.37||-118.25||7.3|
|San JacintoA+C M7.5 Scenario||33.85||-117.06||7.5|
|San JacintoB M6.8 Scenario||33.18||-116.16||6.8|
|San JacintoB+SM M7.1 Scenario||33.20||-116.19||7.1|
|San JacintoCC M7.0 Scenario||33.44||-116.50||7.0|
|San JacintoCC+B M7.2 Scenario||33.44||-116.50||7.2|
|San JacintoCC+B+SM M7.3 Scenario||33.47||-116.54||7.3|
|San JacintoSBV M7.1 Scenario||34.18||-117.40||7.1|
|San JacintoSBV+SJV M7.3 Scenario||34.18||-117.40||7.3|
|San JacintoSBV+SJV+A+C M7.7 Scenario||34.25||-117.50||7.7|
|San JacintoSJV M7.0 Scenario||34.01||-117.22||7.0|
|San JacintoSM M6.7 Scenario||32.98||-115.89||6.7|
|San Joaquin Hills M7.1 Scenario||33.53||-117.95||7.1|
|San Jose M6.7 Scenario||34.13||-117.73||6.7|
|San Juan M7.1 Scenario||35.54||-120.24||7.1|
|Santa Cruz Island M7.2 Scenario||34.06||-119.95||7.2|
|Santa Monica Alt 1 M6.6 Scenario||34.10||-118.42||6.6|
|Santa Monica Alt 2 M6.8 Scenario||34.16||-118.36||6.8|
|Santa Monica Connected Alt 1 M7.3 Scenario||34.15||-118.45||7.3|
|Santa Monica Connected Alt 2 M7.4 Scenario||34.16||-118.39||7.4|
|Santa Rosa Island M6.9 Scenario||34.00||-119.94||6.9|
|Santa Susana Alt 1 M6.9 Scenario||34.40||-118.49||6.9|
|Santa Ynez Connected M7.4 Scenario||34.48||-120.28||7.4|
|Santa Ynez East M7.2 Scenario||34.46||-119.58||7.2|
|Santa Ynez West M7.0 Scenario||34.49||-120.26||7.0|
|Sierra Madre Connected M7.3 Scenario||34.20||-117.72||7.3|
|Sierra Madre M7.2 Scenario||34.20||-117.72||7.2|
|Sierra Madre San Fernando M6.7 Scenario||34.37||-118.28||6.7|
|Simi-Santa Rosa M6.9 Scenario||34.34||-118.73||6.9|
|So Emerson-Copper Mtn M7.1 Scenario||34.17||-116.20||7.1|
|So Sierra Nevada M7.5 Scenario||35.31||-117.93||7.5|
|Superstition Hills M6.8 Scenario||33.02||-115.84||6.8|
|Tank Canyon M6.4 Scenario||35.76||-117.30||6.4|
|Ventura-Pitas Point M7.0 Scenario||34.34||-119.23||7.0|
|Verdugo M6.9 Scenario||34.20||-118.11||6.9|
|White Wolf M7.2 Scenario||35.08||-118.87||7.2|
Estimating Ground Motions for Scenario Earthquake ShakeMaps
At present, ground motions are estimated using an empirical attenuation relationship, which is a predictive relationship that allows the estimation of the peak ground motions at a given distance and for an assumed magnitude. Thus, ground motions are estimated for a given magnitude earthquake, and at a particular distance from the assumed fault, in a manner consistent with recordings of past earthquakes under similar conditions. For ShakeMap, we use the relationship of Boore et al. (1997) for peak and spectral acceleration, and we use Joyner and Boore's (1988) relationship for peak velocity. We use these predictive relationships to estimate peak ground motions on rock sites, and then correct the amplitude at that location based on the site soil conditions as we do in the general ShakeMap interpolation scheme. Site conditions come directly from the Statewide Site Conditions Map for California (Wills et al., 2000) and we correct for site amplification with the amplitude and frequency-dependent factors determined by Borcherdt (1994).
Attributes and Limitations of Current Maps
Our approach is simple and approximate. We account for fault finiteness by measuring the distance to the surface projection of the fault location (Joyner and Boore's distance definition), but we do not consider the direction of rupture nor do we modify the peak motions by a directivity term. With this approach, the location of the earthquake epicenter does not have any effect on the resulting ground motions; only the location and dimensions of the fault matter. If we were to add directivity to the calculations, than different choices of epicentral location would result in significantly different motions for the same magnitude earthquake and fault segment. Rather, our approach here is to show the average effect since it is difficult to show results for every possible epicentral location.
Our empirical predictive approach also only gives average peak ground motions values so it does not account for all the expected variability in motions, other than the aforementioned site amplification variations. Actual ground motions show significant variability for a given distance, magnitude, and site condition and, hence, the scenario ground motions are more uniform than would be expected for an actual earthquake. The true variations are partially attributable to 2D and 3D wave propagation, path effects (such as basin edge amplification and focusing), differences in motions among earthquakes of the same magnitude, and complex site effects not accounted for by our method.
Earthquake scenarios are used heavily in emergency response planning. Primary users for response planning include city, county, state and federal government agencies (e.g., the California Office of Emergency Services, FEMA, the Army Corp of Engineers), emergency response planners and managers for utilities, businesses, and other large organizations. Scenarios are also used for loss-estimation by utilities, governments, and industry.
Scenarios are of fundamental interest to the community and scientific audiences interested in the nature of the ground shaking likely experienced in past earthquakes as well as the possible effects due to rupture on known faults in the future.
In addition, more detailed and careful analysis of the ground motion time histories (seismograms) produced by such scenario earthquakes is highly beneficial for earthquake engineering considerations. Engineers require site-specific ground motions for detailed structural response analysis of existing structures and future structures designed around specified performance levels. In the future, with these scenarios we will also provide synthetic time histories of strong ground motions that include rupture directivity effects.
Future Scientific Advances
While current earthquake modeling techniques are sufficient for providing useful motion time histories and scenario ShakeMaps based on empirical means (e.g., ground motion attenuation relations), substantial improvement will require developing cost-effective numerical tools for proper theoretical inclusion of known complex ground motion effects. These efforts are underway and must continue in order to obtain site, basin and deeper crustal structure, to characterize and test 3D earth models (including attenuation and nonlinearity), and to improve numerical wave propagation methods to obtain useful, site-specific, ground motion time histories.