Search USGS⌕ Search

USGS · 70211190

Preliminary report on the 22 December 2003, M 6.5 San Simeon, California earthquake

Abstract

The M w 6.5 San Simeon earthquake struck the central California coast on 22 December 2003 at 19:15:56 UTC (11:15:56 am local time.) The epicenter was located 11 km northeast of the town of San Simeon, and 39 km west-northwest of Paso Robles ( Figure 1 ), as reported by the California Integrated Seismic Network (CISN, the California region of the Advanced National Seismic System [ANSS]). The mainshock nucleated at 35.702°N, 121.108°W and a depth of 7.1 km, and the rupture propagated unilaterally to the southeast. The strong directivity of the rupture resulted in a concentration of damage and aftershock activity to the southeast of the hypocenter. The worst earthquake damage occurred in Paso Robles, where two people died in the collapse of an unreinforced masonry building. The accurate and rapid earthquake information provided in near real-time by CISN/ANSS to the Governor's Office of Emergency Services made it possible to focus emergency response in the source area, although the earthquake was felt from San Francisco to Los Angeles. The San Simeon earthquake occurred on a reverse fault striking northwest and most likely dipping to the northeast. Although motion along the Pacific-North America plate boundary in California is dominantly strike-slip, there is a small compressional component through central California. Repeated thrust earthquakes such as the San Simeon event accommodate this compression and build the Coast Ranges. Other recent thrust earthquakes in central California include the 1983 Coalinga (M 6.4) and the 1985 Kettleman Hills (M 6.0) earthquakes. Prior earthquakes in the vicinity of the San Simeon event include a M 5-6 earthquake in 1853, a M 5.7 earthquake in 1906, and the M L 6.2 Bryson earthquake of 1952 ( Figure 1 ) ( McLaren and Savage, 2001 .) The San Simeon earthquake occurred on a previously unknown blind thrust fault. No surface rupture associated with the earthquake has been identified. A number of roads, including State Highway 46, buckled due to the earthquake, but this deformation appears mainly to be failure of road fill due to ground shaking and not the result of tectonic surface rupture. Extrapolation of the fault plane to the surface would roughly align with the surface trace of the Oceanic Fault, but this is thought to be a vertical strike-slip fault. Two models for the kinematics of the region have previously been proposed. The first is a fault-propagation fold model developed by Namson and Davis ( 1990 ) for the Santa Lucia mountains ∼30 km to the southeast of the San Simeon sequence. The mainshock geometry is similar to, although more steeply dipping than, the main blind thrust of this model, implying that this model may be applicable to the San Simeon region as well. The second is the model of McLaren and Savage ( 2001 ), in which the region is dominated by strike-slip faulting with shortening on high-angle reverse faults. This model also may be applicable, although the dip of the San Simeon mainshock is shallower than predicted. The San Simeon earthquake was followed by a vigorous aftershock sequence, with 165 events above M 3 reported by CISN within the first week of the mainshock. Although the event triggered many aftershocks, it did not significantly impact the seismicity rates of other nearby faults such as the San Andreas Fault and the San Simeon-Hosgri fault zone. The only triggered seismicity seems to be a few small events within the mainshock coda at the Geysers geothermal area, north of San Francisco. The San Simeon earthquake did, however, trigger shallow creep on the San Andreas Fault at Parkfield and hydrologic changes in hot springs in Paso Robles.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 35.22767235493586° to 35.782170703266075° latitude; -121.26434326171875° to -120.63812255859375° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Jeanne L. Hardebeck, John Boatwright, D. Dreger, Rakesh Goel, V. Graizer, Kenneth W. Hudnut, Ji Chen, Lucile M. Jones, John O. Langbein, Jian Lin, Evelyn A. Roeloffs, Robert W. Simpson, K. Stark, Ross S. Stein, John Tinsley. 2004. Preliminary report on the 22 December 2003, M 6.5 San Simeon, California earthquake. https://doi.org/10.1785/gssrl.75.2.155

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Time for an updated macroseismic scale

Macroseismology plays a crucial role in earthquake hazard and risk analyses, linking past earthquake occurrences and impacts with those of the present and future. The popularity of internet-based community contributions and the usefulness of intensities for communicating seismic hazards, early warning, and impacts to the public have all bolstered the use of macroseismology beyond its traditional utility in characterizing sources and shaking from historical and modern earthquakes (Wald et al., 2023). Nevertheless, many modern approaches for collecting building damage data and assigning macroseismic intensities suffer limitations. Many countries do not or no longer routinely collect such data, and even after significant and even catastrophic earthquakes, building damage data are often not collected or shared systematically, nor are field-based macroseismic intensities professionally assigned. Additionally, the United States (U.S.), Canada, and other nations still use variations of the Modified Mercalli Intensity (MMI) scale (Wood and Neumann, 1931) for assigning macroseismic intensities. At nearly 100 years old, the MMI scale predates most modern building types. As a result, the scale is silent about and blind to building damage that now matters to most domestic occupants. It says nothing about single- or multi-unit wood-frame buildings, the predominant type of residential structure in the US and Canada, whose damage is a particularly important indicator of shaking intensity, nor about the collapse or safety of modern building types. Importantly, the dependence on the U.S. Geological Survey’s “Did You Feel It?” (DYFI) for domestic intensity data collection limits their accuracy type at higher, damaging intensity levels, where attention to building type is requisite. The MMI and some of its widely used variants are often poorly documented; there is no official version, nor is it even clear which version of the scale was employed in many studies (Musson, 2002). The MMI scale is consistent with—but is more limited than—the 1998 European Macroseismic Scale (EMS-98, Grünthal, 1998). The new International Macroseismic Scale of 2025 (IMS-25; Grünthal et al., 2025) builds on EMS-98—a well-documented, modern scale already widely used in Europe and well beyond. IMS-25 further expands EMS-98 internationally by enhancing the EMS-98 vulnerability table and damage grade collections to include more worldwide building types and their evolving design codes. Because the IMS requires well-defined vulnerability classes and damage grades specific to various structural types, its worldwide adoption will help unify post-event data collection on building damage. Similarly, the IMS requires statistical analyses of damage levels for specific structural types, ensuring that building damage sampling is sufficient for other important engineering analyses, such as the derivation of fragility functions, calibration of loss models, and the assessment of losses and damages after significant events

Seismological Research Letters↗

Intensity prediction equations for intraslab earthquakes in the Cascadia Subduction Zone

In the U.S. Pacific Northwest (PNW), intraslab earthquakes located within the subducting Juan de Fuca slab are responsible for most of the region’s moderate‐magnitude seismicity. Historically, these earthquakes have caused significant damage, such as the 2001 M 6.8 Nisqually, 1965 M 6.5 Seattle, and 1949 M 6.8 Olympia events. Intensity prediction equations (IPEs) aim to relate source‐to‐site distance and earthquake magnitude to the resultant felt shaking intensity at a particular location, and help to estimate hazard from earthquakes. We find that a commonly used IPE developed by Atkinson et al. (2014) does not accurately describe shaking intensity from most intraslab events in the PNW, likely because it was developed using primarily crustal fault earthquakes. In this study, we develop two new IPEs that provide an improved fit to observed intensities from Cascadia subduction zone intraslab earthquakes. To develop our IPEs, we compute linear mixed‐effects regressions using community decimal intensity data from the U.S. Geological Survey’s ‘Did You Feel It?’ program from 36 intraslab events that occurred in British Columbia, Canada, Washington, Oregon, and northern California from 2000 to 2025. We also investigate the relationship between the new intraslab IPEs and select parameters related to near‐surface conditions and deeper geologic structure that may be related to ground‐motion intensities, but find no correlations. These new IPEs provide a better fit to observations of intraslab earthquakes in the PNW than Atkinson et al. (2014) , particularly within 300 km of the epicenter. They also provide a better fit to the 2018 M 7.1 Anchorage, Alaska, intraslab earthquake, suggesting that they may have broader applicability to global subduction zones.

California, Oregon, Washington↗

Amploc: A complete Python workflow for amplitude-based locations

Amplitude‐based locations are well suited for characterizing long‐duration seismic sources that may migrate with time, such as debris flows, avalanches, rockfall, and volcanic tremor. Although amplitude‐based location methods have existed for decades, their adoption in real‐time monitoring systems remains limited. We present a complete Python‐based workflow that includes station corrections, amplitude‐based location algorithms, and interactive review tools that are appropriate for real‐time use. The algorithms themselves are not new; however, Amploc integrates several previously published methods into a unified operational framework. The workflow also introduces a browser‐based graphical interface that allows users to review real‐time location results using formats familiar to many observatories. By integrating calibration, location algorithms, and visualization tools, Amploc lowers the barrier to operational use of amplitude‐based locations and improves characterization of phenomena that inform hazard assessments such as long‐period earthquakes, volcanic tremor, and debris flows.

British Columbia, Oregon, Washington↗